Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis

In this report we describe an accurate numerical method for calculating the total electrostatic energy of molecules of arbitrary shape and charge distribution, accounting for both Coulombic and solvent polarization terms. In addition to the solvation energies of individual molecules, the method can be used to calculate the electrostatic energy associated with conformational changes in proteins as well as changes in solvation energy that accompany the binding of charged substrates. The validity of the method is examined by calculating the hydration energies of acetate, methyl ammonium, ammonium, and methanol. The method is then used to study the relationship between the depth of a charge within a protein and its interaction with the solvent. Calculations of the relative electrostatic energies of crystal and misfolded conformations of Themiste dyscritum hemerythrin and the VL domain of an antibody are also presented. The results indicate that electrostatic charge‐solvent interactions strongly favor the crystal structures. More generally, it is found that charge‐solvent interactions, which are frequently neglected in protein structure analysis, can make large contributions to the total energy of a macromolecular system.

[1]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[2]  J. Randles.  The real hydration energies of ions , 1956 .

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

[4]  D. R. Rosseinsky,et al.  Electrode Potentials and Hydration Energies. Theories and Correlations , 1965 .

[5]  J. Hine,et al.  Structural effects on rates and equilibriums. XIX. Intrinsic hydrophilic character of organic compounds. Correlations in terms of structural contributions , 1975 .

[6]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

[7]  W. Hol,et al.  Structure of bovine liver rhodanese. I. Structure determination at 2.5 A resolution and a comparison of the conformation and sequence of its two domains. , 1978, Journal of molecular biology.

[8]  W G Hol,et al.  Structure of bovine liver rhodanese. I. Structure determination at 2.5 A resolution and a comparison of the conformation and sequence of its two domains. , 1978, Journal of Molecular Biology.

[9]  M Karplus,et al.  Picosecond dynamics of tyrosine side chains in proteins. , 1979, Biochemistry.

[10]  Wayne A. Hendrickson,et al.  Structure of the hydrophobic protein crambin determined directly from the anomalous scattering of sulphur , 1981, Nature.

[11]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

[12]  C. Paul,et al.  Building models of globular protein molecules from their amino acid sequences. I. Theory. , 1982, Journal of molecular biology.

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

[14]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[15]  M. L. Connolly.  Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

[16]  B. Honig,et al.  On the environment of ionizable groups in globular proteins. , 1984, Journal of molecular biology.

[17]  L. Sieker,et al.  BINUCLEAR IRON COMPLEXES IN METHEMERYTHRIN AND AZIDOMETHEMERYTHRIN AT 2. 0-A RESOLUTION. , 1984 .

[18]  M. Karplus,et al.  An analysis of incorrectly folded protein models. Implications for structure predictions. , 1984, Journal of molecular biology.

[19]  M. Sternberg,et al.  Electrostatic interactions in globular proteins: calculation of the pH dependence of the redox potential of cytochrome c551. , 1985, Journal of molecular biology.

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

[21]  F. Quiocho,et al.  Sulphate sequestered in the sulphate-binding protein of Salmonella typhimurium is bound solely by hydrogen bonds , 1985, Nature.

[22]  Kurt Wüthrich,et al.  Secondary structure of the α-amylase polypeptide inhibitor Tendamistat from Streptomyces tendae determined in solution by 1H nuclear magnetic resonance , 1985 .

[23]  Barry Honig,et al.  Reevaluation of the Born model of ion hydration , 1985 .

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

[25]  J. Moult,et al.  An algorithm for determining the conformation of polypeptide segments in proteins by systematic search , 1986, Proteins.

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

[27]  Y. Satow,et al.  Phosphocholine binding immunoglobulin Fab McPC603. An X-ray diffraction study at 2.7 A. , 1985, Journal of molecular biology.

[28]  J. Warwicker,et al.  Continuum dielectric modelling of the protein-solvent system, and calculation of the long-range electrostatic field of the enzyme phosphoglycerate mutase. , 1986, Journal of theoretical biology.

[29]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[30]  B. Honig,et al.  Calculation of electrostatic potentials in an enzyme active site , 1987, Nature.

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

[32]  Harold A. Scheraga,et al.  Free energies of hydration of solute molecules. 3. Application of the hydration shell model to charged organic molecules , 1987 .

[33]  M. Karplus,et al.  Prediction of the folding of short polypeptide segments by uniform conformational sampling , 1987, Biopolymers.

[34]  R. Bruccoleri,et al.  Criteria that discriminate between native proteins and incorrectly folded models , 1988, Proteins.

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

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