Abstract:We describe a method for efficient sampling of the energy landscape of a protein in atomic molecular dynamics simulations. A simulation is divided into alternatively occurring relaxation phases and excitation phases. In the relaxation phase (conventional simulation), we use a frequently updated reference structure and deviations from this reference structure to mark whether the system has been trapped in a local minimum. In that case, the simulation enters the excitation phase, during which a few slow collective modes of the system are coupled to a higher temperature bath. After the system has escaped from the minimum (also judged by deviations from the reference structure) the simulation reenters the relaxation phase. The collective modes are obtained from a coarse-grained Gaussian elastic network model. The scheme, which we call ACM-AME (amplified collective motion-assisted minimum escaping), is compared with conventional simulations as well as an alternative scheme that elevates the temperature of all degrees of freedom during the excitation phase (amplified overall motion-assisted minimum escaping, or AOM-AME). Comparison is made using simulations on four peptides starting from non-native extended or all helical structures. In terms of sampling low energy conformations and continuously sampling new conformations throughout a simulation, the ACM-AME scheme demonstrates very good performance while the AOM-AME scheme shows little improvement upon conventional simulations. Limited success is achieved in producing structures close to the native structures of the peptides: for an S-peptide analog, the ACM-AME approach is able to reproduce its native helical structure, and starting from an all-helical structure of the villin headpiece subdomain (HP-36) in implicit solvent, two out of three 150 ns ACM-AME runs are able to sample structures with 3-4 Angstrom backbone root-mean-square deviations from the nuclear magnetic resonance structure of the protein. (C) 2003 American Institute of Physics.

Efficient and Reproducible Folding Simulations of the Trp-cage Protein with Multiscale Molecular Dynamics

Trp-cage: Folding Free Energy Landscape in Explicit Water

Folding Trp-cage to NMR resolution native structure using a coarse-grained model

Folding Simulation of Trp-cage Utilizing a New AMBER Compatible Force Field with Coupled Main Chain Torsions

Molecular Dynamics Characterisations of the Trp-Cage Folding Mechanisms: in the Absence and Presence of Water Solvents

Simulation of the thermodynamics of folding and unfolding of the Trp-cage mini-protein TC5b using different combinations of force fields and solvation models

The Folding Thermodynamics and Kinetics of Crambin Using an All-Atom Monte Carlo Simulation

Folding Thermodynamics and Mechanism of Five Trp-Cage Variants from Replica-Exchange MD Simulations with RSFF2 Force Field.

Computational Model for Protein Unfolding Simulation.

Replica exchange molecular dynamics method for protein folding simulation.

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Folding Simulations for Proteins with Diverse Topologies Are Accessible in Days with A Physics-Based Force Field and Implicit Solvent

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Application of the Accelerated Molecular Dynamics Simulations to the Folding of A Small Protein

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

The "Beacon" Structural Model of Protein Folding: Application for Trp-Cage in Water

Efficiently Explore The Energy Landscape Of Proteins In Molecular Dynamics Simulations By Amplifying Collective Motions

Ratcheted molecular-dynamics simulations identify efficiently the transition state of protein folding

Atomic-detailed milestones along the folding trajectory of protein G

Unraveling Low-Resolution Structural Data Of Large Biomolecules By Constructing Atomic Models With Experiment-Targeted Parallel Cascade Selection Simulations

Thermodynamics and folding pathways of trpzip2: an accelerated molecular dynamics simulation study.