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.

Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Free energy simulations with the AMOEBA polarizable force field and metadynamics on GPU platform

Editorial: Advanced Sampling and Modeling in Molecular Simulations for Slow and Large-Scale Biomolecular Dynamics

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Advances in enhanced sampling molecular dynamics simulations for biomolecules

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

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

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

Development of an anisotropic polarizable model for the all-atom AMOEBA force field

Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Enhanced Sampling of Biomolecular Slow Conformational Transitions Using Adaptive Sampling and Machine Learning

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis.

Evaluating polarizable biomembrane simulations against experiments

Assessing Ion-Water Interactions in the AMOEBA Force Field Using Energy Decomposition Analysis of Electronic Structure Calculations.

Incorporating Neural Networks into the AMOEBA Polarizable Force Field