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.

Automated Event Detection and Activity Monitoring in Long Molecular Dynamics Simulations

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

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics.

Running and analyzing massively parallel molecular simulations

MD-TASK: a software suite for analyzing molecular dynamics trajectories

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

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

drMD: Molecular Dynamics for Experimentalists

MDAnalysis: A toolkit for the analysis of molecular dynamics simulations

Discovering Protein Conformational Flexibility Through Artificial Intelligence Aided Molecular Dynamics.

Studying functional dynamics in bio-molecules using accelerated molecular dynamics

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

Data Driven Classification of Ligand Unbinding Pathways

Systematically Constructing Kinetic Transition Network in Polypeptide from Top to Down: Trajectory Mapping

Molecular Dynamics Simulation for All

EasyAmber: A comprehensive toolbox to automate the molecular dynamics simulation of proteins

Enhancing MD simulations: ASGARD's automated analysis for GROMACS

A Workflow for Exploring Ligand Dissociation from a Macromolecule: Efficient Random Acceleration Molecular Dynamics Simulation and Interaction Fingerprints Analysis of Ligand Trajectories

Ligand Binding and Release Investigated by Contact-Guided Iterative Multiple Independent Molecular Dynamics Simulations

StreaMD: the toolkit for high-throughput molecular dynamics simulations

Elucidation of Protein-Ligand Interactions by Multiple Trajectory Analysis Methods