Efficient simulation method for polarizable protein force fields: Application to the simulation of BPTI in liquid

被引:58
作者
Harder, E
Kim, BC
Friesner, RA
Berne, BJ
机构
[1] Columbia Univ, Dept Chem, New York, NY 10027 USA
[2] Columbia Univ, Ctr Bimol Simulat, New York, NY 10027 USA
关键词
D O I
10.1021/ct049914s
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
A methodology for large scale molecular dynamics simulation of a solvated polarizable protein, using a combination of permanent and inducible point dipoles with fluctuating and fixed charges, is discussed and applied to the simulation of water solvated bovine pancreatic trypsin inhibitor (BPTI). The electrostatic forces are evaluated using a generalized form of the P3M Ewald method which includes point dipoles in addition to point charge sites. The electrostatic configuration is propagated along with the nuclei during the course of the simulation using an extended Lagrangian formalism. For the system size studied, 20000 atoms, this method gives only a marginal computational overhead relative to nonpolarizable potential models (1.23-1.45) per time step of simulation. The models employ a newly developed polarizable dipole force field for the protein(1) with two commonly used water models TlP4P-FQ and RPOL. Performed at constant energy and constant volume (NVE) using the velocity Verlet algorithm, the simulations show excellent energy conservation and run stably for their 2 ns duration. To characterize the accuracy of the solvation models the protein structure is analyzed. The simulated structures remain within 1 A of the experimental crystal structure for the duration of the simulation in line with the nonpolarizable OPLS-AA model.
引用
收藏
页码:169 / 180
页数:12
相关论文
共 66 条
[61]   EFFECTIVE POTENTIALS FOR LIQUID WATER USING POLARIZABLE AND NONPOLARIZABLE MODELS [J].
WALLQVIST, A ;
BERNE, BJ .
JOURNAL OF PHYSICAL CHEMISTRY, 1993, 97 (51) :13841-13851
[62]   DYNAMICS OF MOLECULES WITH INTERNAL DEGREES OF FREEDOM BY MULTIPLE TIME-STEP METHODS [J].
WATANABE, M ;
KARPLUS, M .
JOURNAL OF CHEMICAL PHYSICS, 1993, 99 (10) :8063-8074
[63]   THE EFFECT OF LONG-RANGE ELECTROSTATIC INTERACTIONS IN SIMULATIONS OF MACROMOLECULAR CRYSTALS - A COMPARISON OF THE EWALD AND TRUNCATED LIST METHODS [J].
YORK, DM ;
DARDEN, TA ;
PEDERSEN, LG .
JOURNAL OF CHEMICAL PHYSICS, 1993, 99 (10) :8345-8348
[64]   ATOMIC-LEVEL ACCURACY IN SIMULATIONS OF LARGE PROTEIN CRYSTALS [J].
YORK, DM ;
WLODAWER, A ;
PEDERSEN, LG ;
DARDEN, TA .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 1994, 91 (18) :8715-8718
[65]  
ZHOU R, 1995, J CHEM PHYS, V103
[66]   Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems [J].
Zhou, RH ;
Harder, E ;
Xu, HF ;
Berne, BJ .
JOURNAL OF CHEMICAL PHYSICS, 2001, 115 (05) :2348-2358