Development and parametrization of continuum solvent models. I. Models based on the boundary element method

A series of different simplifications of the boundary element method (BEM) for solving the Poisson–Boltzmann equation is investigated in an effort to obtain an accurate and fast enough treatment of electrostatic effects to be incorporated in Monte‐Carlo and molecular dynamics simulation methods. The tested simplifications include increasing the size of Boundary Elements, decreasing the surface dot density, and ignoring the interactions between the polarization charges. Combined with terms describing the nonelectrostatic solvation effects, the simplified BEM polarization terms were built into expressions for the solvation potential. The solvation potential is treated as empirical consistent force field equations. The intervening parameters, including atomic and probe radii, are derived by different fitting strategies of calculated vs experimental vacuum to water transfer energies of 173 charged, polar, and nonpolar small molecules. These fits are shown to yield very good correlations (rms ∼1.4 kcal/mol), e...

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

[2]  Thomas Simonson,et al.  Solvation Free Energies Estimated from Macroscopic Continuum Theory: An Accuracy Assessment , 1994 .

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

[4]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[5]  M Karplus,et al.  Simulation analysis of the stability mutant R96H of T4 lysozyme. , 1991, Biochemistry.

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

[7]  B. Montgomery Pettitt,et al.  Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids , 1985 .

[8]  Shoshana J. Wodak,et al.  Computer simulations of liquid water: treatment of long-range interactions , 1990 .

[9]  Barry Honig,et al.  A local dielectric constant model for solvation free energies which accounts for solute polarizability , 1992 .

[10]  Thomas E. Creighton,et al.  Stability of folded conformations , 1991 .

[11]  Alexander A. Rashin,et al.  Hydration phenomena, classical electrostatics, and the boundary element method , 1990 .

[12]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[13]  John D. Roberts,et al.  Basic principles of organic chemistry , 1965 .

[14]  Alexander A. Rashin,et al.  A simple method for the calculation of hydration enthalpies of polar molecules with arbitrary shapes , 1987 .

[15]  J. Rullmann,et al.  “Ensemble” iterative relaxation matrix approach: A new NMR refinement protocol applied to the solution structure of crambin , 1993, Proteins.

[16]  Shoshana J. Wodak,et al.  Molecular dynamics simulation of polarizable water by an extended Lagrangian method , 1992 .

[17]  J. C. Martin,et al.  Ground-state analogs of transition states for attack at sulfonyl, sulfinyl, and sulfenyl sulfur: a sulfuranide dioxide (10-S-5) salt, a sulfuranide oxide (10-S-4) salt, and a sulfuranide (10-S-3) salt , 1985 .

[18]  A. T. Hagler,et al.  Consistent force field studies of intermolecular forces in hydrogen-bonded crystals. 2. A benchmark for the objective comparison of alternative force fields , 1979 .

[19]  B. Lee,et al.  Solvent reorganization contribution to the transfer thermodynamics of small nonpolar molecules , 1991, Biopolymers.

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

[21]  C. Tanford,et al.  Theory of Protein Titration Curves. I. General Equations for Impenetrable Spheres , 1957 .

[22]  S J Wodak,et al.  Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Donald G. Truhlar,et al.  AM1-SM2 and PM3-SM3 parameterized SCF solvation models for free energies in aqueous solution , 1992, J. Comput. Aided Mol. Des..

[24]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

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

[26]  H. Berendsen,et al.  The electric potential of a macromolecule in a solvent: A fundamental approach , 1991 .

[27]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[28]  E. Clementi,et al.  Analytical potentials from ab initio computations for the interaction between biomolecules. VII. Polar amino acids and conclusions , 1981 .

[29]  M. Karplus,et al.  Electrostatic effects of charge perturbations introduced by metal oxidation in proteins. A theoretical analysis. , 1988, Journal of molecular biology.

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

[31]  F E Cohen,et al.  Protein model structure evaluation using the solvation free energy of folding. , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[32]  B. Honig,et al.  Accurate First Principles Calculation of Molecular Charge Distributions and Solvation Energies from Ab Initio Quantum Mechanics and Continuum Dielectric Theory , 1994 .

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

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

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

[36]  D. Chandler,et al.  Theory of the hydrophobic effect , 1977 .

[37]  M. Gilson,et al.  Prediction of pH-dependent properties of proteins. , 1994, Journal of molecular biology.

[38]  J L Sussman,et al.  Acetylcholinesterase: electrostatic steering increases the rate of ligand binding. , 1993, Biochemistry.

[39]  Julia M. Goodfellow Computer Modelling in Molecular Biology , 1992 .

[40]  H. Scheraga,et al.  Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[41]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[42]  Cyrus Chothia,et al.  Molecular structure of a new family of ribonucleases , 1982, Nature.

[43]  D. Bashford Electrostatic effects in biological molecules , 1991 .

[44]  Richard Wolfenden,et al.  Comparing the polarities of the amino acids: side-chain distribution coefficients between the vapor phase, cyclohexane, 1-octanol, and neutral aqueous solution , 1988 .

[45]  S J Wodak,et al.  Calculations of electrostatic properties in proteins. Analysis of contributions from induced protein dipoles. , 1987, Journal of molecular biology.

[46]  B Honig,et al.  Salt effects on ligand-DNA binding. Minor groove binding antibiotics. , 1994, Journal of molecular biology.

[47]  K Wüthrich,et al.  Hydration of proteins. A comparison of experimental residence times of water molecules solvating the bovine pancreatic trypsin inhibitor with theoretical model calculations. , 1993, Journal of molecular biology.

[48]  Anna Tempczyk,et al.  Electrostatic contributions to solvation energies: comparison of free energy perturbation and continuum calculations , 1991 .

[49]  M. Karplus,et al.  Hidden thermodynamics of mutant proteins: a molecular dynamics analysis. , 1989, Science.

[50]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[51]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[52]  Norbert Muller,et al.  Search for a realistic view of hydrophobic effects , 1990 .