Discrimination between native and intentionally misfolded conformations of proteins: ES/IS, a new method for calculating conformational free energy that uses both dynamics simulations with an explicit solvent and an implicit solvent continuum model

A new method for calculating the total conformational free energy of proteins in water solvent is presented. The method consists of a relatively brief simulation by molecular dynamics with explicit solvent (ES) molecules to produce a set of microstates of the macroscopic conformation. Conformational energy and entropy are obtained from the simulation, the latter in the quasi‐harmonic approximation by analysis of the covariance matrix. The implicit solvent (IS) dielectric continuum model is used to calculate the average solvation free energy as the sum of the free energies of creating the solute‐size hydrophobic cavity, of the van der Waals solute‐solvent interactions, and of the polarization of water solvent by the solute's charges. The reliability of the solvation free energy depends on a number of factors: the details of arrangement of the protein's charges, especially those near the surface; the definition of the molecular surface; and the method chosen for solving the Poisson equation. Molecular dynamics simulation in explicit solvent relaxes the protein's conformation and allows polar surface groups to assume conformations compatible with interaction with solvent, while averaging of internal energy and solvation free energy tend to enhance the precision. Two recently developed methods—SIMS, for calculation of a smooth invariant molecular surface, and FAMBE, for solution of the Poisson equation via a fast adaptive multigrid boundary element—have been employed. The SIMS and FAMBE programs scale linearly with the number of atoms. SIMS is superior to Connolly's MS (molecular surface) program: it is faster, more accurate, and more stable, and it smooths singularities of the molecular surface. Solvation free energies calculated with these two programs do not depend on molecular position or orientation and are stable along a molecular dynamics trajectory. We have applied this method to calculate the conformational free energy of native and intentionally misfolded globular conformations of proteins (the EMBL set of deliberately misfolded proteins) and have obtained good discrimination in favor of the native conformations in all instances. Proteins 32:399–413, 1998. © 1998 Wiley‐Liss, Inc.

[1]  Craig C. Douglas,et al.  Multigrid methods in science and engineering , 1996 .

[2]  C. Tanford,et al.  Empirical correlation between hydrophobic free energy and aqueous cavity surface area. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

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

[4]  B. Berne,et al.  Molecular Dynamics Study of the Dependence of Water Solvation Free Energy on Solute Curvature and Surface Area , 1995 .

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

[6]  Peter A. Kollman,et al.  Free energy calculation methods: A theoretical and empirical comparison of numerical errors and a new method qualitative estimates of free energy changes , 1997, J. Comput. Chem..

[7]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[8]  L. Greengard,et al.  A Fast Poisson Solver for Complex Geometries , 1995 .

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

[10]  P. Claverie,et al.  Improvements of the continuum model. 1. Application to the calculation of the vaporization thermodynamic quantities of nonassociated liquids , 1988 .

[11]  Volkhard Helms,et al.  Free energies of hydration from thermodynamic integration: Comparison of molecular mechanics force fields and evaluation of calculation accuracy , 1997, J. Comput. Chem..

[12]  Dusanka Janezic,et al.  Harmonic analysis of large systems. II. Comparison of different protein models , 1995, J. Comput. Chem..

[13]  Frederick P. Brooks,et al.  Computing smooth molecular surfaces , 1994, IEEE Computer Graphics and Applications.

[14]  M. Sippl.  Calculation of conformational ensembles from potentials of mena force , 1990 .

[15]  Y. Guissani,et al.  A computer‐simulation study of hydrophobic hydration of rare gases and of methane. I. Thermodynamic and structural properties , 1991 .

[16]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

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

[18]  S Vajda,et al.  Empirical potentials and functions for protein folding and binding. , 1997, Current opinion in structural biology.

[19]  W. L. Jorgensen,et al.  AN EXTENDED LINEAR RESPONSE METHOD FOR DETERMINING FREE ENERGIES OF HYDRATION , 1995 .

[20]  K. Sharp,et al.  Finite difference Poisson‐Boltzmann electrostatic calculations: Increased accuracy achieved by harmonic dielectric smoothing and charge antialiasing , 1997 .

[21]  Harold A. Scheraga,et al.  A combined iterative and boundary-element approach for solution of the nonlinear Poisson-Boltzmann equation , 1992 .

[22]  Robert B. Hermann,et al.  Theory of hydrophobic bonding. II. Correlation of hydrocarbon solubility in water with solvent cavity surface area , 1972 .

[23]  Bhyravabhotla Jayaram,et al.  Free energy calculations of ion hydration: an analysis of the Born model in terms of microscopic simulations , 1989 .

[24]  T. Hansson,et al.  On the Validity of Electrostatic Linear Response in Polar Solvents , 1996 .

[25]  Dusanka Janezic,et al.  Harmonic analysis of large systems. I. Methodology , 1995, J. Comput. Chem..

[26]  B Honig,et al.  On the calculation of binding free energies using continuum methods: Application to MHC class I protein‐peptide interactions , 1997, Protein science : a publication of the Protein Society.

[27]  Steven W. Rick,et al.  The Aqueous Solvation of Water: A Comparison of Continuum Methods with Molecular Dynamics , 1994 .

[28]  S. Wodak,et al.  Factors influencing the ability of knowledge-based potentials to identify native sequence-structure matches. , 1994, Journal of molecular biology.

[29]  Robert E. Bruccoleri,et al.  Finite difference Poisson-Boltzmann electrostatic calculations: Increased accuracy achieved by harmonic dielectric smoothing and charge antialiasing , 1997, J. Comput. Chem..

[30]  Harold A. Scheraga,et al.  A fast adaptive multigrid boundary element method for macromolecular electrostatic computations in a solvent , 1997, J. Comput. Chem..

[31]  R. Jernigan,et al.  Structure-derived potentials and protein simulations. , 1996, Current opinion in structural biology.

[32]  E. Lattman,et al.  Experimental measurement of the effective dielectric in the hydrophobic core of a protein. , 1997, Biophysical chemistry.

[33]  K. Dill,et al.  Statistical potentials extracted from protein structures: how accurate are they? , 1996, Journal of molecular biology.

[34]  S. Sriharan,et al.  The fast multipole boundary element method for molecular electrostatics: An optimal approach for large systems , 1995, J. Comput. Chem..

[35]  Akbar Nayeem,et al.  MSEED: A program for the rapid analytical determination of accessible surface areas and their derivatives , 1992 .

[36]  H. Berendsen,et al.  A comparison of techniques for calculating protein essential dynamics , 1997 .

[38]  G. Phillips,et al.  A global model of the protein-solvent interface. , 1994, Biophysical journal.

[39]  I. Tuñón,et al.  Molecular surface area and hydrophobic effect. , 1992, Protein engineering.

[40]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[41]  Ronald M. Levy,et al.  Gaussian fluctuation formula for electrostatic free‐energy changes in solution , 1991 .

[42]  K. B. Ward,et al.  Occluded molecular surface: Analysis of protein packing , 1995, Journal of molecular recognition : JMR.

[43]  C. Chothia,et al.  Hydrophobic bonding and accessible surface area in proteins , 1974, Nature.

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

[45]  H A Scheraga,et al.  Coupling between folding and ionization equilibria: effects of pH on the conformational preferences of polypeptides. , 1996, Journal of molecular biology.

[46]  Lev Davidovich Landau,et al.  Mechanics, Vol. I of Course on Theoretical Physics , 1972 .

[47]  Narendra S. Goel,et al.  Electrostatic field in inhomogeneous dielectric media I.: indirect boundary element method , 1995 .

[48]  J Moult,et al.  Comparison of database potentials and molecular mechanics force fields. , 1997, Current opinion in structural biology.

[49]  Alexandre Varnek,et al.  A fast and Space-Efficient boundary element method for computing electrostatic and hydration effects in large molecules , 1996, J. Comput. Chem..

[50]  Dusanka Janezic,et al.  Harmonic analysis of large systems. III. Comparison with molecular dynamics , 1995, J. Comput. Chem..

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

[52]  J. Moult,et al.  Determination of the conformation of folding initiation sites in proteins by computer simulation , 1995, Proteins.

[53]  R. Linford,et al.  Comparison of the Gibbs energy and entropy of interfaces water-n-hexane and water-perfluorotributylamine , 1970 .

[54]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

[55]  R. A. Scott,et al.  Discriminating compact nonnative structures from the native structure of globular proteins. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Adding backbone to protein folding: why proteins are polypeptides. , 1995, Folding & design.

[57]  William L. Jorgensen,et al.  Free energy of TIP4P water and the free energies of hydration of CH4 and Cl- from statistical perturbation theory , 1989 .

[58]  Gerhard Hummer,et al.  The hydration free energy of water , 1995 .

[59]  David Chandler,et al.  Statistical mechanics of chemical equilibria and intramolecular structures of nonrigid molecules in condensed phases , 1976 .

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

[61]  Benoît Roux,et al.  Molecular basis for the Born model of ion solvation , 1990 .

[62]  Michael L. Connolly,et al.  Computation of molecular volume , 1985 .

[63]  K. Sharp,et al.  Entropy in protein folding and in protein-protein interactions. , 1997, Current opinion in structural biology.

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

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

[66]  M J Sippl,et al.  Structure-derived hydrophobic potential. Hydrophobic potential derived from X-ray structures of globular proteins is able to identify native folds. , 1992, Journal of molecular biology.

[67]  Iñaki Tuñón,et al.  GEPOL: An improved description of molecular surfaces. III. A new algorithm for the computation of a solvent‐excluding surface , 1994, J. Comput. Chem..

[68]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[69]  Y N Vorobjev,et al.  SIMS: computation of a smooth invariant molecular surface. , 1997, Biophysical journal.

[70]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[71]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[72]  J M Thornton,et al.  SIRIUS. An automated method for the analysis of the preferred packing arrangements between protein groups. , 1990, Journal of molecular biology.

[73]  D. Case,et al.  Thermodynamics of a reverse turn motif. Solvent effects and side-chain packing. , 1997, Journal of molecular biology.

[74]  Randy J. Zauhar,et al.  SMART: A solvent-accessible triangulated surface generator for molecular graphics and boundary element applications , 1995, J. Comput. Aided Mol. Des..

[75]  R. Jernigan,et al.  Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. , 1996, Journal of molecular biology.

[76]  S. Bryant,et al.  Critical assessment of methods of protein structure prediction (CASP): Round II , 1997, Proteins.

[77]  R. Pierotti,et al.  A scaled particle theory of aqueous and nonaqueous solutions , 1976 .

[78]  A Tropsha,et al.  Energetic decomposition of the alpha-helix-coil equilibrium of a dynamic model system. , 1996, Biopolymers.

[79]  B. J. Yoon,et al.  A boundary element method for molecular electrostatics with electrolyte effects , 1990 .

[80]  Michael Levitt,et al.  Finite‐difference solution of the Poisson–Boltzmann equation: Complete elimination of self‐energy , 1996, J. Comput. Chem..

[81]  Michael J. Holst,et al.  Numerical solution of the nonlinear Poisson–Boltzmann equation: Developing more robust and efficient methods , 1995, J. Comput. Chem..

[82]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[83]  K. Merz,et al.  Analysis of a large data base of electrostatic potential derived atomic charges , 1992 .

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

[85]  A. H. Juffer,et al.  Comparison of atomic solvation parametric sets: Applicability and limitations in protein folding and binding , 1995, Protein science : a publication of the Protein Society.

[86]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[87]  A. Varnek,et al.  A fast and Space‐efficient boundary element method for computing electrostatic and hydration effects in large molecules , 1996 .

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

[89]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[90]  Anders Wallqvist,et al.  Computer Simulation of Hydrophobic Hydration Forces on Stacked Plates at Short Range , 1995 .

[91]  R. Zauhar,et al.  The rigorous computation of the molecular electric potential , 1988 .

[92]  A. Ben-Naim.  Solvent effects on protein association and protein folding , 1990, Biopolymers.

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

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

[95]  Howard Reiss,et al.  Aspects of the Statistical Thermodynamics of Real Fluids , 1960 .

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

[97]  R. Levy,et al.  Thermodynamics of the Hydration Shell. 1. Excess Energy of a Hydrophobic Solute , 1994 .

[98]  M J Sternberg,et al.  Application of scaled particle theory to model the hydrophobic effect: implications for molecular association and protein stability. , 1994, Protein engineering.

[99]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[100]  H. Scheraga,et al.  Entropy changes accompanying association reactions of proteins. , 1963, The Journal of biological chemistry.

[101]  Robert B. Hermann,et al.  Modeling hydrophobic solvation of nonspherical systems: Comparison of use of molecular surface area with accessible surface area , 1997, J. Comput. Chem..

[102]  A. Beyer,et al.  An improved pair potential to recognize native protein folds , 1994, Proteins.

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

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

[105]  A. Y. Meyer.  The size of molecules , 1987 .

[106]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[107]  Bruce J. Berne,et al.  A Monte Carlo simulation of the hydrophobic interaction , 1979 .