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.

ELViM: Exploring Biomolecular Energy Landscapes through Multidimensional Visualization

ELViM: Exploring Biomolecular Energy Landscapes through Multidimensional Visualization

Understanding the Energy Landscape of Intrinsically Disordered Protein Ensembles

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

A Mode Evolution Metric to Extract Converged Reaction Coordinates for Biomolecular Conformational Transitions

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

Volumetric Visualization And Analysis Of Molecular Interaction Field

MD DaVis: interactive data visualization of protein molecular dynamics

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

PROTERAN: animated terrain evolution for visual analysis of patterns in protein folding trajectory.

MolSieve: A Progressive Visual Analytics System for Molecular Dynamics Simulations

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

HTMoL: full-stack solution for remote access, visualization, and analysis of Molecular Dynamics trajectory data

Free energy along drug-protein binding pathways interactively sampled in virtual reality

Exploring Molecular Energy Landscapes by Coupling the DFTB Potential with a Tree-Based Stochastic Algorithm: Investigation of the Conformational Diversity of Phthalates

Exploring the Folding Pathways of Proteins through Energy Landscape Sampling: Application to Alzheimer's E-Amyloid Peptide #

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

New Features in Visual Dynamics 3.0

visGReMLIN: graph mining-based detection and visualization of conserved motifs at 3D protein-ligand interface at the atomic level

Efficient Exploration of Membrane-Associated Phenomena at Atomic Resolution

Efficient sampling of high-dimensional free energy landscapes using adaptive reinforced dynamics