A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems

An efficient molecular dynamics (MD) algorithm is presented in this paper for biomolecular systems, which incorporates a novel variation on the fast multipole method (FMM) coupled to the reversible reference system propagator algorithm (r‐RESPA). A top‐down FMM is proposed which calculates multipoles recursively from the top of the box tree instead of from the bottom in Greengard’s original FMM, in an effort to be more efficient for noncubic or nonuniform systems. In addition, the use of noncubic box subdivisions of biomolecular systems is used and discussed. Reversible RESPA based on a Trotter factorization of the Liouville propagator in generating numerical integration schemes is coupled to the top‐down FMM and applied to a MD study of proteins in vacuo, and is shown to be able to use a much larger time‐step than the standard velocity Verlet method for a comparable level of accuracy. Furthermore, by using the FMM it becomes possible to perform MD simulations for very large biomolecules, since memory and...

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

[2]  C. Tanford.  Macromolecules , 1994, Nature.

[3]  H. Trotter.  On the product of semi-groups of operators , 1959 .

[4]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[5]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

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

[7]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[8]  W. F. van Gunsteren,et al.  Effect of constraints on the dynamics of macromolecules , 1982 .

[9]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[10]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[11]  J. Haile,et al.  A multiple time-step method for molecular dynamics simulations of fluids of chain molecules , 1984 .

[12]  J. Skilling,et al.  Maximum entropy image reconstruction: general algorithm , 1984 .

[13]  V. Rokhlin.  Rapid solution of integral equations of classical potential theory , 1985 .

[14]  Andrew W. Appel,et al.  An Efficient Program for Many-Body Simulation , 1983 .

[15]  B. Jönsson,et al.  Vectorizing a general purpose molecular dynamics simulation program , 1986 .

[16]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[17]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

[18]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[19]  S. Toxvaerd.  Comment on constrained molecular dynamics of macromolecules , 1987 .

[20]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[21]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[22]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[23]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with long range forces , 1991 .

[24]  Mark E. Tuckerman,et al.  Molecular dynamics in systems with multiple time scales: Systems with stiff and soft degrees of freedom and with short and long range forces , 1991 .

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

[26]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[27]  J. W. Causey,et al.  Accelerated molecular dynamics simulation with the parallel fast multipole algorithm , 1992 .

[28]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

[29]  Masakatsu Watanabe,et al.  Dynamics of molecules with internal degrees of freedom by multiple time-step methods , 1993 .

[30]  Mark E. Tuckerman,et al.  Vibrational relaxation in simple fluids: Comparison of theory and simulation , 1993 .

[31]  D. Chandler,et al.  Diabatic surfaces and the pathway for primary electron transfer in a photosynthetic reaction center , 1993 .

[32]  Martin Head-Gordon,et al.  Derivation and efficient implementation of the fast multipole method , 1994 .

[33]  B. Berne,et al.  A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules , 1994 .

[34]  Piero Procacci,et al.  Computer simulation of solid C60 using multiple time‐step algorithms , 1994 .

[35]  Jiro Shimada,et al.  Performance of fast multipole methods for calculating electrostatic interactions in biomacromolecular simulations , 1994, J. Comput. Chem..

[36]  Scott Wilson.  The Protein Folding Problem and Tertiary Structure Prediction , 1994, Birkhäuser Boston.

[37]  M. Karplus,et al.  SIMULATIONS OF MACROMOLECULES BY MULTIPLE TIME-STEP METHODS , 1995 .

[38]  R. Jaenicke.  The protein folding problem and tertiary structure prediction , 1995 .

[39]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.