Free-Energy Profiles for Membrane Permeation of Compounds Calculated Using Rare-Event Sampling Methods

The free-energy profile of a compound is an essential measurement in evaluating the membrane permeation process by means of theoretical methods. Computationally, molecular dynamics (MD) simulation allows the free-energy profile calculation. However, MD simulations frequently fail to sample membrane permeation because they are rare events induced in longer timescales than the accessible timescale of MD, leading to an insufficient conformational search to calculate an incorrect free-energy profile. To achieve a sufficient conformational search, several enhanced sampling methods have been developed and elucidated the membrane permeation process. In addition to these enhanced sampling methods, we proposed a simple yet powerful free-energy calculation of a compound for the membrane permeation process based on originally rare-event sampling methods developed by us. Our methods have a weak dependency on external biases and their optimizations to promote the membrane permeation process. Based on distributed computing, our methods only require the selection of initial structures and their conformational resampling, whereas the enhanced sampling methods may be required to adjust external biases. Furthermore, our methods efficiently search membrane permeation processes with simple scripts without modifying any MD program. As demonstrations, we calculated the free-energy profiles of seven linear compounds for their membrane permeation based on a hybrid conformational search using two rare-event sampling methods, that is, (1) parallel cascade selection MD (PaCS-MD) and (2) outlier flooding method (OFLOOD), combined with a Markov state model (MSM) construction. In the first step, PaCS-MD generated initial membrane permeation paths of a compound. In the second step, OFLOOD expanded the unsearched conformational area around the initial paths, allowing for a broad conformational search. Finally, the trajectories were employed to construct reliable MSMs, enabling correct free-energy profile calculations. Furthermore, we estimated the membrane permeability coefficients of all compounds by constructing the reliable MSMs for their membrane permeation. In conclusion, the calculated coefficients were qualitatively correlated with the experimental measurements (correlation coefficient (R²) = 0.8689), indicating that the hybrid conformational search successfully calculated the free-energy profiles and membrane permeability coefficients of the seven compounds.

Protein Structure Validation Derives a Smart Conformational Search in a Physically Relevant Configurational Subspace

Since proteins perform biological functions through their dynamic properties, molecular dynamics (MD) simulation is a sophisticated strategy for investigating their functions. Analyses of trajectories provide statistical information about a specific protein as a free-energy landscape (FEL). However, the timescale of normal MD is shorter than that of biological functions, resulting in statistically insufficient conformational sampling, finally leading to unreliable FEL calculation. To search for a broad configurational subspace, an external bias is imposed on a target protein as biased sampling. However, its regulation is challenging because the optimal strength of the perturbation is unknown. Furthermore, a physically irrelevant configurational subspace was searched when imposing an inappropriate external bias. To address this issue, we newly proposed an external biased regulation scheme known as the G-factor external bias limiter (GERBIL). In GERBIL, protein configurations generated by external bias are structurally validated by an indicator (G-factor), enabling the search for a physically relevant subspace. In addition to biased sampling, nonbiased sampling might search for a physically irrelevant configurational subspace because repeating multiple MD simulations from several initial structures tends to search for an overly broad configurational subspace. For this issue, the structural qualities of configurations generated by nonbiased sampling have not been investigated. Therefore, we confirmed whether the G-factor screened the collapsed (low-quality) configurations generated by nonbiased sampling. To address this issue, the outlier flooding method (OFLOOD) was adopted in GERBIL as a nonbiased sampling method, which is referred to as OFLOOD-GERBIL. OFLOOD rapidly expands a configurational subspace by resampling the rarely occurring states of a given protein and tends to search an overly broad subspace. Thus, we considered that GERBIL might improve the excessive conformational search of OFLOOD for a physically irrelevant configurational subspace. As a demonstration, OFLOOD and OFLOOD-GERBIL were applied to a globular protein (T4 lysozyme) and their conformational search qualities were assessed. Based on our assessment, normal OFLOOD without the outlier validation frequently sampled low-quality configurations, whereas OFLOOD-GERBIL with the outlier validation intensively sampled high-quality configurations. In conclusion, OFLOOD-GERBIL derives a smart conformational search in a physically relevant configurational subspace, indicating that protein structure validation works in both nonbiased and biased sampling methods.

The role of ATP in solubilizing RNA-binding protein fused in sarcoma

Intrinsically disordered protein (IDP) plays an important role in liquid-liquid phase separation (LLPS). RNA-binding protein fused in sarcoma (FUS) is a well-studied IDP that induces LLPS since its low-complexity core region (FUS-LC-core) is essential for droplet formation through contacts between FUS-LC-cores. Several experimental studies have reported that adenosine triphosphate (ATP) concentrations modulate LLPS-driven droplet formation through the dissolution of FUS. To elucidate the role of ATP in this dissolution, microsecond-order all-atom molecular dynamics (MD) simulations were performed for a crowded system of FUS-LC-cores in the presence of multiple ATP molecules. Our analysis revealed that the adenine group of ATP frequently contacted the FUS-LC-core, and the phosphoric acid group of ATP was exposed to the external solvent, which promoted both hydration and solubilization of FUS.

Structural Validation by the G-Factor Properly Regulates Boost Potentials Imposed in Conformational Sampling of Proteins

Free energy landscapes (FELs) of proteins are indispensable for evaluating thermodynamic properties. Molecular dynamics (MD) simulation is a computational method for calculating FELs; however, conventional MD simulation frequently fails to search a broad conformational subspace due to its accessible timescale, which results in the calculation of an unreliable FEL. To search a broad subspace, an external bias can be imposed on a protein system, and biased sampling tends to cause a strong perturbation that might collapse the protein structures, indicating that the strength of the external bias should be properly regulated. This regulation can be challenging, and empirical parameters are frequently employed to impose an optimal bias. To address this issue, several methods regulate the external bias by referring to system energies. Herein, we focused on protein structural information for this regulation. In this study, a well-established structural indicator (the G-factor) was used to obtain structural information. Based on the G-factor, we proposed a scheme for regulating biased sampling, which is referred to as a G-factor-based external bias limiter (GERBIL). With GERBIL, the configurations were structurally validated by the G-factor during biased sampling. As an example of biased sampling, an accelerated MD (aMD) simulation was adopted in GERBIL (aMD-GERBIL), whereby the aMD simulation was repeatedly performed by increasing the strength of the boost potential. Furthermore, the configurations sampled by the aMD simulation were structurally validated by their G-factor values, and aMD-GERBIL stopped increasing the strength of the boost potential when the sampled configurations were regarded as low-quality (collapsed) structures. This structural validation is regarded as a "Brake" of the boost potential. For demonstrations, aMD-GERBIL was applied to globular proteins (ribose binding and maltose-binding proteins) to promote their large-amplitude open-closed transitions and successfully identify their domain motions.

Ligand Binding Path Sampling Based on Parallel Cascade Selection Molecular Dynamics: LB-PaCS-MD

Parallel cascade selection molecular dynamics (PaCS-MD) is a rare-event sampling method that generates transition pathways between a reactant and product. To sample the transition pathways, PaCS-MD repeats short-time MD simulations from important configurations as conformational resampling cycles. In this study, PaCS-MD was extended to sample ligand binding pathways toward a target protein, which is referred to as LB-PaCS-MD. In a ligand-concentrated environment, where multiple ligand copies are randomly arranged around the target protein, LB-PaCS-MD allows for the frequent sampling of ligand binding pathways. To select the important configurations, we specified the center of mass (COM) distance between each ligand and the relevant binding site of the target protein, where snapshots generated by the short-time MD simulations were ranked by their COM distance values. From each cycle, snapshots with smaller COM distance values were selected as the important configurations to be resampled using the short-time MD simulations. By repeating conformational resampling cycles, the COM distance values gradually decreased and converged to constants, meaning that a set of ligand binding pathways toward the target protein was sampled by LB-PaCS-MD. To demonstrate relative efficiency, LB-PaCS-MD was applied to several proteins, and their ligand binding pathways were sampled more frequently than conventional MD simulations.

A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing

The Protein Folding Problem (PFP) is a big challenge that has remained unsolved for more than fifty years. This problem consists of obtaining the tertiary structure or Native Structure (NS) of a protein knowing its amino acid sequence. The computational methodologies applied to this problem are classified into two groups, known as Template-Based Modeling (TBM) and ab initio models. In the latter methodology, only information from the primary structure of the target protein is used. In the literature, Hybrid Simulated Annealing (HSA) algorithms are among the best ab initio algorithms for PFP; Golden Ratio Simulated Annealing (GRSA) is a PFP family of these algorithms designed for peptides. Moreover, for the algorithms designed with TBM, they use information from a target protein’s primary structure and information from similar or analog proteins. This paper presents GRSA-SSP methodology that implements a secondary structure prediction to build an initial model and refine it with HSA algorithms. Additionally, we compare the performance of the GRSAX-SSP algorithms versus its corresponding GRSAX. Finally, our best algorithm GRSAX-SSP is compared with PEP-FOLD3, I-TASSER, QUARK, and Rosetta, showing that it competes in small peptides except when predicting the largest peptides.

Selection Rules for Outliers in Outlier Flooding Method Regulate Its Conformational Sampling Efficiency

The outlier flooding method (OFLOOD) has been proposed as an enhanced conformational sampling method of proteins. In OFLOOD, rarely occurring states of proteins are detected as sparse conformational distributions (outliers) with a clustering algorithm. The detected outliers are intensively resampled with short-time molecular dynamics (MD) simulations. As a set of cycles, OFLOOD repeats selections of outliers and their conformational resampling. Herein, as an essential issue to be tackled to perform OFLOOD efficiently, a selection rule for outliers should be carefully specified. Generally, many outliers are detected from distributions on conformational subspaces with the clustering. Judging from its computational costs, it is unreasonable to select all the detected outliers upon the conformational resampling. Therefore, it is important to consider which outliers should be selected from the sparse distributions when restarting their short-time MD simulations with limited computational costs. In this sense, we investigated the conformational sampling efficiency of OFLOOD by changing the selection rules for outliers. To address the conformational sampling efficiency of OFLOOD depending on its selection rules, outliers to be resampled were selected by focusing their probability occurrences (populations of outliers). As a comparison, a random selection rule for outliers was also considered. Through the present assessment, the random selection of outliers showed the most efficient conformational sampling efficiency compared to the other OFLOOD trials using the biased selection rules, indicating that a variety of outliers should be selected and resampled during the OFLOOD cycles. In conclusion, the random outlier selection rule is the best strategy to perform OFLOOD efficiently.

Selection rules on initial structures in parallel cascade selection molecular dynamics affect conformational sampling efficiency

Parallel cascade selection molecular dynamics (PaCS-MD) is a conformational sampling method for generating transition pathways from a given reactant to a product. In PaCS-MD, initial structures relevant to conformational transitions of proteins are selected and resampled by short-time MD simulations. As a general reaction coordinate, a root-mean-square deviation measured from the product (RMSD) is employed to rank the resampled configurations. Quantitatively, n initial structures are randomly selected from among the top X % of highly ranked configurations and resampled again. In PaCS-MD, the selection of initial structures and their conformational resampling are repeated as a cycle to promote the essential conformational transitions. Therefore, rules for selecting the initial structures might affect the conformational sampling efficiency. In the present study, to address the conformational sampling efficiency depending on the selection rule, the open-closed transition of di-ubiquitin was reproduced by PaCS-MD based on the resampling from the top X = 0.1, 1.0, 2.0, 5.0, 10.0, 25.0, 30.0, 40.0, and 50.0% of highly ranked configurations. Judging from broadness of sampled conformational area and required cycles, we conclude that the resampling from the top ∼2.0% of highly ranked configurations might be the most efficient for generating a set of transition pathways in PaCS-MD.

Efficient Conformational Search Based on Structural Dissimilarity Sampling: Applications for Reproducing Structural Transitions of Proteins

Structural Dissimilarity Sampling (SDS) is proposed as an efficient conformational search method to promote structural transitions essential for the biological functions of proteins. In SDS, initial structures are selected based on structural dissimilarity, and conformational resampling is repeated. Conformational resampling is performed as follows: (I) arrangement of initial structures for a diverse distribution at the edge of a conformational subspace and (II) promotion of the structural transitions with multiple short-time molecular dynamics (MD) simulations restarting from the diversely distributed initial structures. Cycles of (I) and (II) are repeated to intensively promote structural transitions because conformational resampling from the initial structures would quickly expand conformational distributions toward unvisited conformational subspaces. As a demonstration, SDS was first applied to maltodextrin binding protein (MBP) in explicit water to reproduce structural transitions between the open and closed states of MBP. Structural transitions of MBP were successfully reproduced with SDS in nanosecond-order simulation times. Starting from both the open and closed forms, SDS successfully reproduced the structural transitions within 25 cycles (a total of 250 ns of simulation time). For reference, a conventional long-time (500 ns) MD simulation under NPT (300 K and 1 bar) starting from the open form failed to reproduce the structural transition. In addition to the open-closed motions of MBP, SDS was applied to folding processes of the fast-folding proteins (chignolin, Trp-cage, and villin) and successfully sampled their native states. To confirm how the selections of initial structures affected conformational sampling efficiency, numbers of base sets for characterizing structural dissimilarity of initial structures were addressed in distinct trials of SDS. The parameter searches showed that the conformational sampling efficiency was relatively insensitive with respect to the numbers of base sets, indicating the robustness of SDS for actual applications.