Protein molecular dynamics with electrostatic force entirely determined by a single Poisson‐Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson‐Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500‐ps molecular dynamics (MD) simulation with explicit water and a 500‐ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant. Proteins 2002;48:497–504. © 2002 Wiley‐Liss, Inc.

[1]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[2]  Wei Zu Chen,et al.  A Stochastic Dynamics Simulation Study Associated with Hydration Force and Friction Memory Effect , 2000 .

[3]  Yunyu Shi,et al.  Investigating the dielectric effects of channel pore water on the electrostatic barriers of the permeation ion by the finite difference Poisson-Boltzmann method , 1998, European Biophysics Journal.

[4]  B Honig,et al.  Reconciling the magnitude of the microscopic and macroscopic hydrophobic effects. , 1991, Science.

[5]  J. Andrew McCammon,et al.  Computation of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation , 1993 .

[6]  T L Blundell,et al.  The structure of 2Zn pig insulin crystals at 1.5 A resolution. , 1988, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[7]  J. Andrew McCammon,et al.  Comparison of Continuum and Explicit Models of Solvation: Potentials of Mean Force for Alanine Dipeptide , 1996 .

[8]  A. Li,et al.  Identification and characterization of the unfolding transition state of chymotrypsin inhibitor 2 by molecular dynamics simulations. , 1996, Journal of molecular biology.

[9]  Ray Luo,et al.  Comparison of generalized born and poisson models: Energetics and dynamics of HIV protease , 2000 .

[10]  J. Andrew McCammon,et al.  The Low Dielectric Interior of Proteins is Sufficient To Cause Major Structural Changes in DNA on Association , 1996 .

[11]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[12]  Michael K. Gilson,et al.  The inclusion of electrostatic hydration energies in molecular mechanics calculations , 1991, J. Comput. Aided Mol. Des..

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

[14]  R. Zauhar,et al.  The incorporation of hydration forces determined by continuum electrostatics into molecular mechanics simulations , 1991 .

[15]  S. Cannistraro,et al.  Molecular dynamics of copper plastocyanin: simulations of structure and dynamics as a function of hydration , 1994 .

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

[17]  J. Andrew McCammon,et al.  Conformational sampling with Poisson–Boltzmann forces and a stochastic dynamics/Monte Carlo method: Application to alanine dipeptide , 1997 .

[18]  James Andrew McCammon,et al.  Molecular dynamics simulation with a continuum electrostatic model of the solvent , 1995, J. Comput. Chem..

[19]  T. Simonson,et al.  Protein molecular dynamics with the generalized born/ACE solvent model , 2001, Proteins.

[20]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[21]  J. Andrew McCammon,et al.  Application of Poisson—Boltzmann solvation forces to macromolecular simulations , 1997 .

[22]  Kim A. Sharp,et al.  Incorporating solvent and ion screening into molecular dynamics using the finite‐difference Poisson–Boltzmann method , 1991 .

[23]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[24]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[25]  W. F. Gunsteren,et al.  β-Hairpin stability and folding: Molecular dynamics studies of the first β-hairpin of tendamistat , 2000 .

[26]  Tony J. You,et al.  Conformation and hydrogen ion titration of proteins: a continuum electrostatic model with conformational flexibility. , 1995, Biophysical journal.

[27]  Martin Karplus,et al.  pH-Dependence of Protein Stability: Absolute Electrostatic Free Energy Differences between Conformations† , 1997 .