Normal mode partitioning of Langevin dynamics for biomolecules
- 14 April 2008
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 128 (14), 145101
- https://doi.org/10.1063/1.2883966
Abstract
We propose a novel normal mode multiple time stepping Langevin dynamics integrator called NML. The aim is to approximate the kinetics or thermodynamics of a biomolecule by a reduced model based on a normal mode decomposition of the dynamical space. Our basis set uses the eigenvectors of a mass reweighted Hessian matrix calculated with a biomolecular force field. This particular choice has the advantage of an ordering according to the eigenvalues, which have a physical meaning of being the square of the mode frequency. Low frequency eigenvalues correspond to more collective motions, whereas the highest frequency eigenvalues are the limiting factor for the stability of the integrator. In NML, the higher frequency modes are overdamped and relaxed near their energy minimum while respecting the subspace of low frequency dynamical modes. Our numerical results confirm that both sampling and rates are conserved for an implicitly solvated alanine dipeptide model, with only 30% of the modes propagated, when compared to the full model. For implicitly solvated systems, NML gives a twofold improvement in efficiency over plain Langevin dynamics for sampling a small 22 atom (alanine dipeptide) model and in excess of an order of magnitude for sampling an 882 atom (bovine pancreatic trypsin inhibitor) model, with good scaling with system size subject to the number of modes propagated. NML has been implemented in the open source software PROTOMOL.Keywords
This publication has 69 references indexed in Scilit:
- Normal mode based flexible fitting of high-resolution structure into low-resolution experimental data from cryo-EMJournal of Structural Biology, 2004
- Dynamics of large proteins through hierarchical levels of coarse‐grained structuresJournal of Computational Chemistry, 2001
- All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of ProteinsThe Journal of Physical Chemistry B, 1998
- Essential dynamics of proteinsProteins-Structure Function and Bioinformatics, 1993
- Multiple Time Steps: Limits on the Speedup of Molecular Dynamics Simulations of Aqueous SystemsMolecular Simulation, 1993
- Reversible multiple time scale molecular dynamicsThe Journal of Chemical Physics, 1992
- On the convergence of the SHAKE algorithmComputer Physics Communications, 1991
- Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range InteractionsMolecular Simulation, 1991
- Protein normal-mode dynamics: Trypsin inhibitor, crambin, ribonuclease and lysozymeJournal of Molecular Biology, 1985
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983