Influence of the solvent representation on vibrational entropy calculations: Generalized born versus distance‐dependent dielectric model

The harmonic model is the most popular approximation for estimating the “configurational” entropy of a solute in molecular mechanics/Poisson‐Boltzmann solvent accessible surface area (MM/PBSA)‐type binding free energy calculations. Here, we investigate the influence of the solvent representation in the harmonic model by comparing estimates of changes in the vibrational entropies for 30 trypsin/ligand complexes on ligand binding. Second derivatives of Amber generalized Born (GB) solvation models are available in the nucleic acid builder code. They allow one to use these models for the calculation of vibrational entropies instead of using a simpler solvation model based on a distance‐dependent dielectric (DDD) constant. Estimates of changes in the vibrational entropies obtained with a DDD model are systematically and significantly larger, by on average, 6 kcal mol−1 (at T = 300 K), than estimates obtained with a GB model and so are more favorable for complex formation. The difference becomes larger the more the vibrational entropy contribution disfavors complex formation, that is, the larger the ligand is (for the complexes considered here). A structural decomposition of the estimates into per‐residue contributions reveals polar interactions between the ligand and the surrounding protein, in particular involving charged nitrogens, as a main source of the differences. Snapshots minimized with the DDD model showed a structural deviation from snapshots minimized in explicit water that is larger by, on average, 0.5 Å RMSD compared to snapshots that were minimized with GBHCT. As experimental vibrational entropies of biomacromolecules are elusive, there is no direct way to establish a solvent model's superiority. Thus, we can only recommend using the GB harmonic model for vibrational entropy calculations based on the reasoning that smaller structural deviations should point to the implicit solvent model that closer approximates the energy landscape of the solute in explicit solvent. © 2012 Wiley Periodicals, Inc.

[1]  A Kitao,et al.  Energy landscape of a native protein: Jumping‐among‐minima model , 1998, Proteins.

[2]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[3]  T. Šolmajer,et al.  Electrostatic screening in molecular dynamics simulations. , 1991, Protein engineering.

[4]  Paul A. Bates,et al.  Can MM‐PBSA calculations predict the specificities of protein kinase inhibitors? , 2006, J. Comput. Chem..

[5]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[6]  E. Mehler,et al.  Electrostatic effects in proteins: comparison of dielectric and charge models. , 1991, Protein engineering.

[7]  Holger Gohlke,et al.  Converging free energy estimates: MM‐PB(GB)SA studies on the protein–protein complex Ras–Raf , 2004, J. Comput. Chem..

[8]  D. Case,et al.  Rescoring docking hit lists for model cavity sites: predictions and experimental testing. , 2008, Journal of molecular biology.

[9]  Jacob Kongsted,et al.  An improved method to predict the entropy term with the MM/PBSA approach , 2009, J. Comput. Aided Mol. Des..

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

[11]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[12]  A. Roitberg,et al.  All-atom structure prediction and folding simulations of a stable protein. , 2002, Journal of the American Chemical Society.

[13]  David A. Case,et al.  Second derivatives in generalized Born theory , 2006, J. Comput. Chem..

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

[15]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[16]  Olivier Michielin,et al.  Comparison between computational alanine scanning and per‐residue binding free energy decomposition for protein–protein association using MM‐GBSA: Application to the TCR‐p‐MHC complex , 2007, Proteins.

[17]  D. Case.  Normal mode analysis of protein dynamics , 1994 .

[18]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[19]  N. Go,et al.  New method for calculating the conformational entropy of a regular helix. , 1974, Macromolecules.

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

[21]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[22]  R W Pickersgill,et al.  A rapid method of calculating charge-charge interaction energies in proteins. , 1988, Protein engineering.

[23]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[24]  Irwin D Kuntz,et al.  Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods. , 2003, Biopolymers.

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

[26]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[27]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[28]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[29]  C. Tanford.  Macromolecules , 1994, Nature.

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

[31]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[32]  Chandra Verma,et al.  Dissecting the Vibrational Entropy Change on Protein/Ligand Binding: Burial of a Water Molecule in Bovine Pancreatic Trypsin Inhibitor , 2001 .

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

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

[35]  J. Richardson,et al.  Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation. , 1999, Journal of molecular biology.

[36]  R. Zhou.  Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

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

[38]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[39]  M. Whitlow,et al.  Conformations of trypsin-bound amidine inhibitors of blood coagulant factor Xa by double REDOR NMR and MD simulations. , 1999, Journal of medicinal chemistry.

[40]  R. Levy,et al.  Free energy surfaces of β‐hairpin and α‐helical peptides generated by replica exchange molecular dynamics with the AGBNP implicit solvent model , 2004 .

[41]  V. Hornak,et al.  Investigation of Salt Bridge Stability in a Generalized Born Solvent Model. , 2006, Journal of chemical theory and computation.

[42]  Harold A. Scheraga,et al.  Analysis of the Contribution of Internal Vibrations to the Statistical Weights of Equilibrium Conformations of Macromolecules , 1969 .

[43]  A. W. Schüttelkopf,et al.  PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. , 2004, Acta crystallographica. Section D, Biological crystallography.

[44]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.