Incorporating variable dielectric environments into the generalized Born model.

A generalized Born (GB) model is proposed that approximates the electrostatic part of macromolecular solvation free energy over the entire range of the solvent and solute dielectric constants. The model contains no fitting parameters, and is derived by matching a general form of the GB Green function with the exact Green's function of the Poisson equation for a random charge distribution inside a perfect sphere. The sphere is assumed to be filled uniformly with dielectric medium epsilon(in), and is surrounded by infinite solvent of constant dielectric epsilon(out). This model is as computationally efficient as the conventional GB model based on the widely used functional form due to Still et al. [J. Am. Chem. Soc. 112, 6127 (1990)], but captures the essential physics of the dielectric response for all values of epsilon(in) and epsilon(out). This model is tested against the exact solution on a perfect sphere, and against the numerical Poisson-Boltzmann (PB) treatment on a set of macromolecules representing various structural classes. It shows reasonable agreement with both the exact and the numerical solutions of the PB equation (where available) considered as reference, and is more accurate than the conventional GB model over the entire range of dielectric values.

[1]  Yong Duan,et al.  Distinguish protein decoys by Using a scoring function based on a new AMBER force field, short molecular dynamics simulations, and the generalized born solvent model , 2004, Proteins.

[2]  Ronald M. Levy,et al.  AGBNP: An analytic implicit solvent model suitable for molecular dynamics simulations and high‐resolution modeling , 2004, J. Comput. Chem..

[3]  Anthony K. Felts,et al.  Free Energy Surfaces of -Hairpin and -Helical Peptides Generated by Replica Exchange Molecular Dynamics with the AGBNP Implicit Solvent Model , 2004 .

[4]  Velin Z. Spassov,et al.  Introducing an Implicit Membrane in Generalized Born/Solvent Accessibility Continuum Solvent Models , 2002 .

[5]  F. Young Biochemistry , 1955, The Indian Medical Gazette.

[6]  J. Kirkwood,et al.  Theory of Solutions of Molecules Containing Widely Separated Charges with Special Application to Zwitterions , 1934 .

[7]  C. Brooks,et al.  An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins. , 2003, Biophysical journal.

[8]  J. Apostolakis,et al.  Continuum Electrostatic Energies of Macromolecules in Aqueous Solutions , 1997 .

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

[10]  W. C. Still,et al.  The GB/SA Continuum Model for Solvation. A Fast Analytical Method for the Calculation of Approximate Born Radii , 1997 .

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

[12]  Gregory D. Hawkins,et al.  Pairwise solute descreening of solute charges from a dielectric medium , 1995 .

[13]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

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

[15]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[16]  A. Bondi van der Waals Volumes and Radii , 1964 .

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

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

[19]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[20]  T. Simonson,et al.  Protein molecular dynamics with the generalized born/ACE solvent model , 2001, Proteins.

[21]  Tomasz Grycuk,et al.  Deficiency of the Coulomb-field approximation in the generalized Born model: An improved formula for Born radii evaluation , 2003 .

[22]  C. Cramer,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. , 1999, Chemical reviews.

[23]  Donald Bashford,et al.  An Object-Oriented Programming Suite for Electrostatic Effects in Biological Molecules , 1997, ISCOPE.

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

[25]  Yutaka Ishikawa,et al.  Scientific Computing in Object-Oriented Parallel Environments , 1997, Lecture Notes in Computer Science.

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

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

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

[29]  M. Karplus,et al.  A Comprehensive Analytical Treatment of Continuum Electrostatics , 1996 .

[30]  D. Case,et al.  Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model , 2000 .

[31]  P Beroza,et al.  Calculations of proton-binding thermodynamics in proteins. , 1998, Methods in enzymology.

[32]  Ray Luo,et al.  Comparison of generalized born and poisson models: Energetics and dynamics of HIV protease , 2000 .

[33]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[34]  Richard A. Friesner,et al.  Solvation Free Energies of Peptides: Comparison of Approximate Continuum Solvation Models with Accurate Solution of the Poisson−Boltzmann Equation , 1997 .

[35]  H. Nymeyer,et al.  Simulation of the folding equilibrium of α-helical peptides: A comparison of the generalized Born approximation with explicit solvent , 2003, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

[39]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[40]  Ting Wang,et al.  Implicit solvent models for flexible protein–protein docking by molecular dynamics simulation , 2002, Proteins.

[41]  David A. Case,et al.  Effective Born radii in the generalized Born approximation: The importance of being perfect , 2002, J. Comput. Chem..