Solvent models for protein–ligand binding: Comparison of implicit solvent poisson and surface generalized born models with explicit solvent simulations

Solvent effects play a crucial role in mediating the interactions between proteins and their ligands. Implicit solvent models offer some advantages for modeling these interactions, but they have not been parameterized on such complex problems, and therefore, it is not clear how reliable they are. We have studied the binding of an octapeptide ligand to the murine MHC class I protein using both explicit solvent and implicit solvent models. The solvation free energy calculations are more than 103 faster using the Surface Generalized Born implicit solvent model compared to FEP simulations with explicit solvent. For some of the electrostatic calculations needed to estimate the binding free energy, there is near quantitative agreement between the explicit and implicit solvent model results; overall, the qualitative trends in the binding predicted by the explicit solvent FEP simulations are reproduced by the implicit solvent model. With an appropriate choice of reference system based on the binding of the discharged ligand, electrostatic interactions are found to enhance the binding affinity because the favorable Coulomb interaction energy between the ligand and protein more than compensates for the unfavorable free energy cost of partially desolvating the ligand upon binding. Some of the effects of protein flexibility and thermal motions on charging the peptide in the solvated complex are also considered. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 591–607, 2001

[1]  Ronald M. Levy,et al.  Prediction of pKa Shifts without Truncation of Electrostatic Interactions: An Explicit Solvent Calculation for Succinic Acid , 1996 .

[2]  D. Beveridge,et al.  A MODIFICATION OF THE GENERALIZED BORN THEORY FOR IMPROVED ESTIMATES OF SOLVATION ENERGIES AND PK SHIFTS , 1998 .

[3]  Gerhard Hummer,et al.  Molecular Theories and Simulation of Ions and Polar Molecules in Water , 1998 .

[4]  Thomas Simonson,et al.  Microscopic Dielectric Properties of Cytochrome c from Molecular Dynamics Simulations in Aqueous Solution , 1995 .

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

[6]  R. Friesner,et al.  Generalized Born Model Based on a Surface Integral Formulation , 1998 .

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

[8]  William L. Jorgensen,et al.  Temperature and size dependence for Monte Carlo simulations of TIP4P water , 1985 .

[9]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

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

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

[12]  Donald G. Truhlar,et al.  PARAMETRIZED MODEL FOR AQUEOUS FREE ENERGIES OF SOLVATION USING GEOMETRY-DEPENDENT ATOMIC SURFACE TENSIONS WITH IMPLICIT ELECTROSTATICS , 1997 .

[13]  Bruce Tidor,et al.  Optimizing electrostatic affinity in ligand-receptor binding: Theory, computation, and ligand properties , 1998 .

[14]  Terry P. Lybrand,et al.  Hydration of chloride and bromide anions: determination of relative free energy by computer simulation , 1985 .

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

[16]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[17]  William L. Jorgensen,et al.  Monte Carlo simulations of the hydration of substituted benzenes with OPLS potential functions , 1993, J. Comput. Chem..

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

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

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

[21]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

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

[23]  I. Tamm,et al.  Fundamentals of the theory of electricity , 1979 .

[24]  B. Berne,et al.  Large scale simulation of macromolecules in solution: Combining the periodic fast multipole method with multiple time step integrators , 1997 .

[25]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[26]  B Tidor,et al.  Computation of electrostatic complements to proteins: A case of charge stabilized binding , 1998, Protein science : a publication of the Protein Society.

[27]  J. W. Causey,et al.  Accelerated molecular dynamics simulation with the parallel fast multipole algorithm , 1992 .

[28]  R. Levy,et al.  Computer simulations with explicit solvent: recent progress in the thermodynamic decomposition of free energies and in modeling electrostatic effects. , 1998, Annual review of physical chemistry.

[29]  P. A. Peterson,et al.  Crystal structure of an H-2Kb-ovalbumin peptide complex reveals the interplay of primary and secondary anchor positions in the major histocompatibility complex binding groove. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

[33]  R. Zauhar,et al.  A new method for computing the macromolecular electric potential. , 1985, Journal of molecular biology.

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

[35]  J A McCammon,et al.  Poisson-Boltzmann analysis of the lambda repressor-operator interaction. , 1992, Biophysical journal.

[36]  M. Karplus,et al.  pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model. , 1990, Biochemistry.

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

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

[39]  Alan E. Mark,et al.  Dielectric properties of trypsin inhibitor and lysozyme calculated from molecular dynamics simulations , 1993 .

[40]  William L. Jorgensen,et al.  Cis-trans energy difference for the peptide bond in the gas phase and in aqueous solution , 1988 .

[41]  Fumio Hirata,et al.  Conserving energy during molecular dynamics simulations of water, proteins, and proteins in water , 1990 .

[42]  M. B. Pinto,et al.  Optimized δ expansion for relativistic nuclear models , 1997, nucl-th/9709049.

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

[44]  B. Honig,et al.  New Model for Calculation of Solvation Free Energies: Correction of Self-Consistent Reaction Field Continuum Dielectric Theory for Short-Range Hydrogen-Bonding Effects , 1996 .

[45]  Bhyravabhotla Jayaram,et al.  Solvation Free Energy of Biomacromolecules: Parameters for a Modified Generalized Born Model Consistent with the AMBER Force Field , 1998 .

[46]  J. Mccammon,et al.  Effect of artificial periodicity in simulations of biomolecules under Ewald boundary conditions: a continuum electrostatics study. , 1999, Biophysical chemistry.

[47]  Heather A. Carlson,et al.  Accuracy of free energies of hydration for organic molecules from 6‐31g*‐derived partial charges , 1993, J. Comput. Chem..

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

[49]  Martin Karplus,et al.  SOLVATION. A MOLECULAR DYNAMICS STUDY OF A DIPEPTIDE IN WATER. , 1979 .

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

[51]  Richard A. Friesner,et al.  An automatic three‐dimensional finite element mesh generation system for the Poisson–Boltzmann equation , 1997 .

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

[53]  J. G. Goodwin,et al.  Surface concentrations and residence times of intermediates on samarium oxide during the oxidative coupling of methane , 1990 .

[54]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[55]  O. Steinhauser,et al.  Calculation of the dielectric properties of a protein and its solvent: theory and a case study. , 1997, Journal of molecular biology.

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

[57]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.