Browse
Submit Data
Databases
API
Help

Dataset Information

17 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Accelerated Molecular Dynamics Applied to the Peptaibol Folding Problem.

ABSTRACT: The use of enhanced sampling molecular dynamics simulations to facilitate the folding of proteins is a relatively new approach which has quickly gained momentum in recent years. Accelerated molecular dynamics (aMD) can elucidate the dynamic path from the unfolded state to the near-native state, "flattened" by introducing a non-negative boost to the potential. Alamethicin F30/3 (Alm F30/3), chosen in this study, belongs to the class of peptaibols that are 7-20 residue long, non-ribosomally synthesized, amphipathic molecules that show interesting membrane perturbing activity. The recent studies undertaken on the Alm molecules and their transmembrane channels have been reviewed. Three consecutive simulations of ~900 ns each were carried out where N-terminal folding could be observed within the first 100 ns, while C-terminal folding could only be achieved almost after 800 ns. It took ~1 μs to attain the near-native conformation with stronger potential boost which may take several μs worth of classical MD to produce the same results. The Alm F30/3 hexamer channel was also simulated in an E. coli mimicking membrane under an external electric field that correlates with previous experiments. It can be concluded that aMD simulation techniques are suited to elucidate peptaibol structures and to understand their folding dynamics.

SUBMITTER: Tyagi C

PROVIDER: S-EPMC6747184 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Accelerated molecular dynamics simulations of protein folding.

Project description:Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.

| S-EPMC4487363 | biostudies-literature

Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration.

Project description:Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential. In conventional AMD the reweighting procedure has an inherent statistical problem in systems with large acceleration, where the points with the largest biases will dominate the reweighted result and reduce the effective number of data points. We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling. The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

| S-EPMC2651661 | biostudies-literature

Using Selectively Applied Accelerated Molecular Dynamics to Enhance Free Energy Calculations.

Project description:Accelerated molecular dynamics (aMD) has been shown to enhance conformational space sampling relative to classical molecular dynamics; however, the exponential reweighting of aMD trajectories, which is necessary for the calculation of free energies relating to the classical system, is oftentimes problematic, especially for systems larger than small poly peptides. Here, we propose a method of accelerating only the degrees of freedom most pertinent to sampling, thereby reducing the total acceleration added to the system and improving the convergence of calculated ensemble averages, which we term selective aMD. Its application is highlighted in two biomolecular cases. First, the model system alanine dipeptide is simulated with classical MD, all-dihedral aMD, and selective aMD, and these results are compared to the infinite sampling limit as calculated with metadynamics. We show that both forms of aMD enhance the convergence of the underlying free energy landscape by 5-fold relative to classical MD; however, selective aMD can produce improved statistics over all-dihedral aMD due to the improved reweighting. Then we focus on the pharmaceutically relevant case of computing the free energy of the decoupling of oseltamivir in the active site of neuraminidase. Results show that selective aMD greatly reduces the cost of this alchemical free energy transformation, whereas all-dihedral aMD produces unreliable free energy estimates.

| S-EPMC2976571 | biostudies-literature

Accelerated Molecular Dynamics Simulation for Helical Proteins Folding in Explicit Water.

Project description:In this study, we examined the folding processes of eight helical proteins (2I9M, TC5B, 1WN8, 1V4Z, 1HO2, 1HLL, 2KFE, and 1YYB) at room temperature using the explicit solvent model under the AMBER14SB force field with the accelerated molecular dynamics (AMD) and traditional molecular dynamics (MD), respectively. We analyzed and compared the simulation results obtained by these two methods based on several aspects, such as root mean square deviation (RMSD), native contacts, cluster analysis, folding snapshots, free energy landscape, and the evolution of the radius of gyration, which showed that these eight proteins were successfully and consistently folded into the corresponding native structures by AMD simulations carried out at room temperature. In addition, the folding occurred in the range of 40~180 ns after starting from the linear structures of the eight proteins at 300 K. By contrast, these stable folding structures were not found when the traditional molecular dynamics (MD) simulation was used. At the same time, the influence of high temperatures (350, 400, and 450 K) is also further investigated. Study found that the simulation efficiency of AMD is higher than that of MD simulations, regardless of the temperature. Of these temperatures, 300 K is the most suitable temperature for protein folding for all systems. To further investigate the efficiency of AMD, another trajectory was simulated for eight proteins with the same linear structure but different random seeds at 300 K. Both AMD trajectories reached the correct folded structures. Our result clearly shows that AMD simulation are a highly efficient and reliable method for the study of protein folding.

| S-EPMC6691143 | biostudies-literature

Tripleurin XIIc: Peptide Folding Dynamics in Aqueous and Hydrophobic Environment Mimic Using Accelerated Molecular Dynamics.

Project description:Peptaibols are a special class of fungal peptides with an acetylated N-terminus and a C-terminal 1,2-amino alcohol along with non-standard amino acid residues. New peptaibols named tripleurins were recently identified from a strain of the filamentous fungal species Trichodermapleuroti, which is known to cause green mould disease on cultivated oyster mushrooms. To understand the mode of action of these peptaibols, the three-dimensional structure of tripleurin (TPN) XIIc, an 18-mer peptide, was elucidated using an enhanced sampling method, accelerated MD, in water and chloroform solvents. Non-standard residues were parameterized by the Restrained Electrostatic Potential (RESP) charge fitting method. The dihedral distribution indicated towards a right-handed helical formation for TPN XIIc in both solvents. Dihedral angle based principal component analysis revealed a propensity for a slightly bent, helical folded conformation in water solvent, while two distinct conformations were revealed in chloroform: One that folds into highly bent helical structure that resembles a beta-hairpin and another with an almost straight peptide backbone appearing as a rare energy barrier crossing event. The hinge-like movement of the terminals was also observed and is speculated to be functionally relevant. The convergence and efficient sampling is addressed using Cartesian PCA and Kullback-Leibler divergence methods.

| S-EPMC6359335 | biostudies-literature

Accelerated Molecular Dynamics for Peptide Folding: Benchmarking Different Combinations of Force Fields and Explicit Solvent Models.

Project description:Accelerated molecular dynamics (aMD) protocols were assessed on predicting the secondary structure of eight peptides, of which two are helical, three are β-hairpins, and three are disordered. Protocols consisted of combinations of three force fields (ff99SB, ff14SB, ff19SB) and two explicit solvation models (TIP3P and OPC), and were evaluated in two independent aMD simulations, one starting from an extended conformation, the other starting from a misfolded conformation. The results of these analyses indicate that all three combinations performed well on helical peptides. As for β-hairpins, ff19SB performed well with both solvation methods, with a slight preference for the TIP3P solvation model, even though performance was dependent on both peptide sequence and initial conformation. The ff19SB/OPC combination had the best performance on intrinsically disordered peptides. In general, ff14SB/TIP3P suffered the strongest helical bias.

| S-EPMC10207344 | biostudies-literature

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics.

Project description:An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038-2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pK(a) values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pK(a) values in closer agreement to those obtained experimentally for the more slowly converging residues.

| S-EPMC2817915 | biostudies-literature

Gaussian Accelerated Molecular Dynamics in NAMD.

Project description:Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that provides efficient free energy calculations of biomolecules. Like the previous accelerated molecular dynamics (aMD), GaMD allows for "unconstrained" enhanced sampling without the need to set predefined collective variables and so is useful for studying complex biomolecular conformational changes such as protein folding and ligand binding. Furthermore, because the boost potential is constructed using a harmonic function that follows Gaussian distribution in GaMD, cumulant expansion to the second order can be applied to recover the original free energy profiles of proteins and other large biomolecules, which solves a long-standing energetic reweighting problem of the previous aMD method. Taken together, GaMD offers major advantages for both unconstrained enhanced sampling and free energy calculations of large biomolecules. Here, we have implemented GaMD in the NAMD package on top of the existing aMD feature and validated it on three model systems: alanine dipeptide, the chignolin fast-folding protein, and the M3 muscarinic G protein-coupled receptor (GPCR). For alanine dipeptide, while conventional molecular dynamics (cMD) simulations performed for 30 ns are poorly converged, GaMD simulations of the same length yield free energy profiles that agree quantitatively with those of 1000 ns cMD simulation. Further GaMD simulations have captured folding of the chignolin and binding of the acetylcholine (ACh) endogenous agonist to the M3 muscarinic receptor. The reweighted free energy profiles are used to characterize the protein folding and ligand binding pathways quantitatively. GaMD implemented in the scalable NAMD is widely applicable to enhanced sampling and free energy calculations of large biomolecules.

| S-EPMC5743237 | biostudies-literature

Implementation of Accelerated Molecular Dynamics in NAMD.

Project description:Accelerated molecular dynamics (aMD) is an enhanced-sampling method that improves the conformational space sampling by reducing energy barriers separating different states of a system. Here we present the implementation of aMD in the parallel simulation program NAMD. We show that aMD simulations performed with NAMD have only a small overhead compared with classical MD simulations. Through example applications to the alanine dipeptide, we discuss the choice of acceleration parameters, the interpretation of aMD results, as well as the advantages and limitations of the aMD method.

| S-EPMC3115733 | biostudies-literature

Gaussian Accelerated Molecular Dynamics: Theory, Implementation, and Applications.

Project description:A novel Gaussian Accelerated Molecular Dynamics (GaMD) method has been developed for simultaneous unconstrained enhanced sampling and free energy calculation of biomolecules. Without the need to set predefined reaction coordinates, GaMD enables unconstrained enhanced sampling of the biomolecules. Furthermore, by constructing a boost potential that follows a Gaussian distribution, accurate reweighting of GaMD simulations is achieved via cumulant expansion to the second order. The free energy profiles obtained from GaMD simulations allow us to identify distinct low energy states of the biomolecules and characterize biomolecular structural dynamics quantitatively. In this chapter, we present the theory of GaMD, its implementation in the widely used molecular dynamics software packages (AMBER and NAMD), and applications to the alanine dipeptide biomolecular model system, protein folding, biomolecular large-scale conformational transitions and biomolecular recognition.

| S-EPMC5927394 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data