New ways to boost molecular dynamics simulations

We describe a set of algorithms that allow to simulate dihydrofolate reductase (DHFR, a common benchmark) with the AMBER all‐atom force field at 160 nanoseconds/day on a single Intel Core i7 5960X CPU (no graphics processing unit (GPU), 23,786 atoms, particle mesh Ewald (PME), 8.0 Å cutoff, correct atom masses, reproducible trajectory, CPU with 3.6 GHz, no turbo boost, 8 AVX registers). The new features include a mixed multiple time‐step algorithm (reaching 5 fs), a tuned version of LINCS to constrain bond angles, the fusion of pair list creation and force calculation, pressure coupling with a “densostat,” and exploitation of new CPU instruction sets like AVX2. The impact of Intel's new transactional memory, atomic instructions, and sloppy pair lists is also analyzed. The algorithms map well to GPUs and can automatically handle most Protein Data Bank (PDB) files including ligands. An implementation is available as part of the YASARA molecular modeling and simulation program from www.YASARA.org. © 2015 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

[1]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[2]  R. Dror,et al.  Systematic Validation of Protein Force Fields against Experimental Data , 2012, PloS one.

[3]  Gert Vriend,et al.  Making optimal use of empirical energy functions: Force‐field parameterization in crystal space , 2004, Proteins.

[4]  Diwakar Shukla,et al.  OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. , 2013, Journal of chemical theory and computation.

[5]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[6]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[7]  Axel T. Brunger,et al.  X-PLOR Version 3.1: A System for X-ray Crystallography and NMR , 1992 .

[8]  Gert Vriend,et al.  YASARA View—molecular graphics for all devices—from smartphones to workstations , 2014, Bioinform..

[9]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[10]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[11]  D. J. Price,et al.  A modified TIP3P water potential for simulation with Ewald summation. , 2004, The Journal of chemical physics.

[12]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[13]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[14]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[15]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[16]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[17]  Keehyoung Joo,et al.  Improving physical realism, stereochemistry, and side‐chain accuracy in homology modeling: Four approaches that performed well in CASP8 , 2009, Proteins.

[18]  Helmut Grubmüller,et al.  Multiple time step algorithms for molecular dynamics simulations of proteins: How good are they? , 1998 .

[19]  A M Gronenborn,et al.  Improvements and extensions in the conformational database potential for the refinement of NMR and X-ray structures of proteins and nucleic acids. , 1997, Journal of magnetic resonance.

[20]  Bernard R. Brooks,et al.  Pressure-Based Long-Range Correction for Lennard-Jones Interactions in Molecular Dynamics Simulations: Application to Alkanes and Interfaces , 2004 .

[21]  D. Tildesley,et al.  Multiple time-step methods in molecular dynamics , 1978 .

[22]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[23]  Kresten Lindorff-Larsen,et al.  Evaluating the Effects of Cutoffs and Treatment of Long-range Electrostatics in Protein Folding Simulations , 2012, PloS one.

[24]  K. Lindorff-Larsen,et al.  How robust are protein folding simulations with respect to force field parameterization? , 2011, Biophysical journal.

[25]  Klaus Schulten,et al.  Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions , 1991 .

[26]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[27]  J Hermans,et al.  Microfolding: Conformational probability map for the alanine dipeptide in water from molecular dynamics simulations , 1988, Proteins.

[28]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[29]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[30]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[31]  R. Unger,et al.  Chaos in protein dynamics , 1997, Proteins.

[32]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[33]  Tetsuya Morishita,et al.  Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath , 2000 .