Performance of fast multipole methods for calculating electrostatic interactions in biomacromolecular simulations

The fast multipole method proposed by Greengard and Rokhlin (GR) is applied to large biomacromolecular systems. In this method, the system is divided into a hierarchy of cells, and electric field exerted on a particle is decomposed into two parts. The first part is a rapidly varying field due to nearby cells, so that it needs rigorous pairwise calculations. The second part is a slowly varying local field due to distant cells; hence, it allows rapid calculations through a multipole expansion technique. In this work, two additional possibilities for improving the performance are numerically examined. The first is an improvement of the convergence of the expansion by increasing the number of nearby cells, without including higher‐order multipole moments. The second is an acceleration of the calculations by the particle–particle and particle–mesh/multipole expansion (PPPM/MPE) method, which uses fast Fourier transform instead of the hierarchy. For this purpose, the PPPM/MPE method originally developed by the authors for a periodic system is extended to a nonperiodic isolated system. The advantages and disadvantages of the GR and PPPM/MPE methods are discussed for both periodic and isolated systems. It is numerically shown that these methods with reasonable costs can reduce the error in potential felt by each particle to 0.1–1 kcal/mol, much smaller than the 30‐kcal/mol error involved in conventional simple truncations. © 1994 by John Wiley & Sons, Inc.

[1]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[2]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[3]  A. Ladd.  Long-range dipolar interactions in computer simulations of polar liquids , 1978 .

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

[5]  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.

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

[7]  Anthony J. C. Ladd,et al.  Monte-Carlo simulation of water , 1977 .

[8]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[9]  K. Schmidt,et al.  Implementing the fast multipole method in three dimensions , 1991 .

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

[11]  Arieh Warshel,et al.  The extended Ewald method: A general treatment of long‐range electrostatic interactions in microscopic simulations , 1988 .

[12]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

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

[14]  Antonio Rizzo,et al.  Selected topics in ab initio computational chemistry in both very small and very large chemical systems , 1991 .

[15]  O. Steinhauser,et al.  Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. , 1992, Biochemistry.

[16]  W. F. V. Gunsteren,et al.  Moleküldynamik‐Computersimulationen; Methodik, Anwendungen und Perspektiven in der Chemie , 1990 .

[17]  Haruki Nakamura,et al.  Numerical Calculations of Reaction Fields of Protein-Solvent Systems , 1988 .

[18]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions , 1992 .

[19]  Jiro Shimada,et al.  Efficient calculations of coulombic interactions in biomolecular simulations with periodic boundary conditions , 1993, J. Comput. Chem..

[20]  Toshikazu Ebisuzaki,et al.  A special-purpose computer for gravitational many-body problems , 1990, Nature.

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

[22]  Leslie Greengard,et al.  The fast mutipole method for gridless particle simulation. , 1987 .

[23]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[24]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

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

[26]  Bernard Pettitt,et al.  Peptides in ionic solutions: A comparison of the Ewald and switching function techniques , 1991 .

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

[28]  H. Berendsen,et al.  COMPUTER-SIMULATION OF MOLECULAR-DYNAMICS - METHODOLOGY, APPLICATIONS, AND PERSPECTIVES IN CHEMISTRY , 1990 .

[29]  Peter A. Kollman,et al.  Conformational and energetic effects of truncating nonbonded interactions in an aqueous protein dynamics simulation , 1993, J. Comput. Chem..