A Poisson–Boltzmann dynamics method with nonperiodic boundary condition

We have developed a well-behaved and efficient finite difference Poisson–Boltzmann dynamics method with a nonperiodic boundary condition. This is made possible, in part, by a rather fine grid spacing used for the finite difference treatment of the reaction field interaction. The stability is also made possible by a new dielectric model that is smooth both over time and over space, an important issue in the application of implicit solvents. In addition, the electrostatic focusing technique facilitates the use of an accurate yet efficient nonperiodic boundary condition: boundary grid potentials computed by the sum of potentials from individual grid charges. Finally, the particle–particle particle–mesh technique is adopted in the computation of the Coulombic interaction to balance accuracy and efficiency in simulations of large biomolecules. Preliminary testing shows that the nonperiodic Poisson–Boltzmann dynamics method is numerically stable in trajectories at least 4 ns long. The new model is also fairly e...

[1]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[2]  M. Marchi,et al.  A dielectric continuum molecular dynamics method , 2001 .

[3]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[4]  D. D. Yue,et al.  Theory of Electric Polarization , 1974 .

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

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

[7]  J. Andrew McCammon,et al.  Electrostatic energy calculations by a Finite‐difference method: Rapid calculation of charge–solvent interaction energies , 1992 .

[8]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[9]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[10]  M. Gilson,et al.  Nucleic acid base-pairing and N-methylacetamide self-association in chloroform: affinity and conformation. , 1999, Biophysical chemistry.

[11]  Kim A. Sharp,et al.  Electrostatic interactions in macromolecules , 1994 .

[12]  Tony You,et al.  An analytical algorithm for the rapid determination of the solvent accessibility of points in a three‐dimensional lattice around a solute molecule , 1995, J. Comput. Chem..

[13]  M. Karplus,et al.  A Comprehensive Analytical Treatment of Continuum Electrostatics , 1996 .

[14]  K. Sharp,et al.  On the calculation of pKas in proteins , 1993, Proteins.

[15]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[16]  Michael K. Gilson,et al.  Synthetic Adenine Receptors: Direct Calculation of Binding Affinity and Entropy , 2000 .

[17]  Wei Zu Chen,et al.  Protein molecular dynamics with electrostatic force entirely determined by a single Poisson‐Boltzmann calculation , 2002, Proteins.

[18]  J. A. Grant,et al.  A Gaussian Description of Molecular Shape , 1995 .

[19]  Barry Honig,et al.  Focusing of electric fields in the active site of Cu‐Zn superoxide dismutase: Effects of ionic strength and amino‐acid modification , 1986, Proteins.

[20]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

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

[22]  S. Sridharan,et al.  A rapid method for calculating derivatives of solvent accessible surface areas of molecules , 1995, J. Comput. Chem..

[23]  B. Dominy,et al.  Development of a generalized Born model parameterization for proteins and nucleic acids , 1999 .

[24]  H. Scheraga,et al.  A fast adaptive multigrid boundary element method for macromolecular electrostatic computations in a solvent , 1997 .

[25]  D. Case,et al.  Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model , 2000 .

[26]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[27]  Richard A. Friesner,et al.  Numerical solution of the Poisson–Boltzmann equation using tetrahedral finite‐element meshes , 1997 .

[28]  W. Im,et al.  Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation , 1998 .

[29]  Randy J. Zauhar,et al.  Computing the electric potential of biomolecules: Application of a new method of molecular surface triangulation , 1990 .

[30]  P A Kollman,et al.  Free energy calculations on dimer stability of the HIV protease using molecular dynamics and a continuum solvent model. , 2000, Journal of molecular biology.

[31]  J. Andrew McCammon,et al.  Dielectric boundary smoothing in finite difference solutions of the poisson equation: An approach to improve accuracy and convergence , 1991 .

[32]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[33]  J. A. McCammon,et al.  Calculating electrostatic forces from grid‐calculated potentials , 1990 .

[34]  K. Sharp,et al.  Macroscopic models of aqueous solutions : biological and chemical applications , 1993 .

[35]  Yunyu Shi,et al.  Stochastic dynamics simulation of alanine dipeptide: Including solvation interaction determined by boundary element method , 1997 .

[36]  Barry Honig,et al.  Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation , 1990 .

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

[38]  J. Andrew Grant,et al.  A smooth permittivity function for Poisson–Boltzmann solvation methods , 2001, J. Comput. Chem..

[39]  Donald Bashford,et al.  An Object-Oriented Programming Suite for Electrostatic Effects in Biological Molecules , 1997, ISCOPE.

[40]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[41]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

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

[43]  Marcia O. Fenley,et al.  Fast Boundary Element Method for the Linear Poisson-Boltzmann Equation , 2002 .

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

[45]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .

[46]  J. Warwicker,et al.  Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. , 1982, Journal of molecular biology.

[47]  S. Subramaniam,et al.  Treatment of electrostatic effects in proteins: Multigrid‐based newton iterative method for solution of the full nonlinear poisson–boltzmann equation , 1994, Proteins.

[48]  Michael J. E. Sternberg,et al.  Regular representation of irregular charge distributions , 1984 .

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

[50]  K. Sharp,et al.  Protein folding and association: Insights from the interfacial and thermodynamic properties of hydrocarbons , 1991, Proteins.

[51]  M. Gilson,et al.  The determinants of pKas in proteins. , 1996, Biochemistry.

[52]  Michael J. Holst,et al.  Multigrid solution of the Poisson—Boltzmann equation , 1992, J. Comput. Chem..

[53]  Huan-Xiang Zhou,et al.  Macromolecular electrostatic energy within the nonlinear Poisson–Boltzmann equation , 1994 .

[54]  Ruhong Zhou,et al.  Poisson−Boltzmann Analytical Gradients for Molecular Modeling Calculations , 1999 .

[55]  Klaus Schulten,et al.  Molecular Dynamics Simulations in Heterogeneous Dielectrica and Debye-Hückel Media - Application to the Protein Bovine Pancreatic Trypsin Inhibitor , 1992 .

[56]  Nathan A. Baker,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation I. Algorithms and examples , 2000 .

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

[58]  M. Gilson,et al.  Strength of Solvent-Exposed Salt-Bridges† , 1999 .

[59]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[60]  Nathan A. Baker,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation II. Refinement at solvent‐accessible surfaces in biomolecular systems , 2000 .

[61]  M K Gilson,et al.  Modeling Molecular Recognition: Theory and Application , 2000, Journal of biomolecular structure & dynamics.

[62]  M K Gilson,et al.  Theory of electrostatic interactions in macromolecules. , 1995, Current opinion in structural biology.

[63]  J. A. McCammon,et al.  Solving the finite difference linearized Poisson‐Boltzmann equation: A comparison of relaxation and conjugate gradient methods , 1989 .

[64]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

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

[66]  M. Gilson.  Modeling protonation equilibria in biomolecules , 1997 .

[67]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[68]  B Honig,et al.  On the pH dependence of protein stability. , 1993, Journal of molecular biology.

[69]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[70]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[71]  B. Honig,et al.  On the calculation of electrostatic interactions in proteins. , 1985, Journal of molecular biology.

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

[73]  M K Gilson,et al.  Energetics of charge–charge interactions in proteins , 1988, Proteins.

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

[75]  E. Alexov,et al.  Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties. , 1997, Biophysical journal.

[76]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.