PERFORMANCE OF FAST MULTIPOLE METHODS FOR CALCULATING ELECTROSTATIC INTERACTIONS IN BIOMACROMOLECULAR SIMULATIONS

被引:59
作者
SHIMADA, J
KANEKO, H
TAKADA, T
机构
[1] Fundamental Research Laboratories, NEC Corporation, Tsukuba, Ibaraki, 305, 34, Miyukigaoka
关键词
D O I
10.1002/jcc.540150105
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
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. (C) 1994 by John Wiley and Sons, Inc.
引用
收藏
页码:28 / 43
页数:16
相关论文
共 37 条
[11]   SELECTED TOPICS IN ABINITIO COMPUTATIONAL CHEMISTRY IN BOTH VERY SMALL AND VERY LARGE CHEMICAL-SYSTEMS [J].
CLEMENTI, E ;
CORONGIU, G ;
BAHATTACHARYA, D ;
FEUSTON, B ;
FRYE, D ;
PREISKORN, A ;
RIZZO, A ;
XUE, W .
CHEMICAL REVIEWS, 1991, 91 (05) :679-699
[12]   SIMULATION OF ELECTROSTATIC SYSTEMS IN PERIODIC BOUNDARY-CONDITIONS .1. LATTICE SUMS AND DIELECTRIC-CONSTANTS [J].
DELEEUW, SW ;
PERRAM, JW ;
SMITH, ER .
PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1980, 373 (1752) :27-56
[13]   ATOMIC LEVEL SIMULATIONS ON A MILLION PARTICLES - THE CELL MULTIPOLE METHOD FOR COULOMB AND LONDON NONBOND INTERACTIONS [J].
DING, HQ ;
KARASAWA, N ;
GODDARD, WA .
JOURNAL OF CHEMICAL PHYSICS, 1992, 97 (06) :4309-4315
[14]   A FAST ALGORITHM FOR PARTICLE SIMULATIONS [J].
GREENGARD, L ;
ROKHLIN, V .
JOURNAL OF COMPUTATIONAL PHYSICS, 1987, 73 (02) :325-348
[15]  
GREENGARD L, 1989, CHEM SCRIPTA, V29A, P139
[16]  
GREENGARD L, 1988, 602 YAL U DEP COMP S
[17]   CONFORMATIONAL AND ENERGETIC EFFECTS OF TRUNCATING NONBONDED INTERACTIONS IN AN AQUEOUS PROTEIN DYNAMICS SIMULATION [J].
GUENOT, J ;
KOLLMAN, PA .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1993, 14 (03) :295-311
[18]  
HANSEN JP, 1986, MOL DYNAMICS SIMULAT, P89
[19]   THE OPLS POTENTIAL FUNCTIONS FOR PROTEINS - ENERGY MINIMIZATIONS FOR CRYSTALS OF CYCLIC-PEPTIDES AND CRAMBIN [J].
JORGENSEN, WL ;
TIRADORIVES, J .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1988, 110 (06) :1657-1666
[20]   THE EXTENDED EWALD METHOD - A GENERAL TREATMENT OF LONG-RANGE ELECTROSTATIC INTERACTIONS IN MICROSCOPIC SIMULATIONS [J].
KUWAJIMA, S ;
WARSHEL, A .
JOURNAL OF CHEMICAL PHYSICS, 1988, 89 (06) :3751-3759