Is Poisson-Boltzmann theory insufficient for protein folding simulations?

The Poisson-Boltzmann theory has been widely used in the studies of energetics and conformations of biological macromolecules. Recently, introduction of the efficient generalized Born approximation has greatly extended its applicability to areas such as protein folding simulations where highly efficient computation is crucial. However, limitations have been found in the folding simulations of a well-studied beta hairpin with several generalized Born implementations and different force fields. These studies have raised the question whether the underlining Poisson-Boltzmann theory, on which the generalized Born model is calibrated, is adequate in the treatment of polar interactions for the challenging protein folding simulations. To address the question whether the Poisson-Boltzmann theory in the current formalism might be insufficient, we directly tested our efficient numerical Poisson-Boltzmann implementation in the beta-hairpin folding simulation. Good agreement between simulation and experiment was found for the beta-hairpin equilibrium structures when the numerical Poisson-Boltzmann solvent and a recently improved generalized Born solvent were used. In addition simulated thermodynamic properties also agree well with experiment in both solvents. Finally, an overall agreement on the beta-hairpin folding mechanism was found between the current and previous studies. Thus, our simulations indicate that previously observed limitations are most likely due to imperfect calibration in previous generalized Born models but not due to the limitation of the Poisson-Boltzmann theory.

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

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

[3]  Alexander D. MacKerell,et al.  Force field influence on the observation of π-helical protein structures in molecular dynamics simulations , 2003 .

[4]  V. Muñoz,et al.  Folding dynamics and mechanism of β-hairpin formation , 1997, Nature.

[5]  V S Pande,et al.  Molecular dynamics simulations of unfolding and refolding of a beta-hairpin fragment of protein G. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[6]  W. Im,et al.  Optimized atomic radii for protein continuum electrostatics solvation forces. , 1999, Biophysical chemistry.

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

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

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

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

[11]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[12]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

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

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

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

[16]  J. Onuchic,et al.  Funnels, pathways, and the energy landscape of protein folding: A synthesis , 1994, Proteins.

[17]  Wilfred F. van Gunsteren,et al.  The importance of solute-solvent van der Waals interactions with interior atoms of biopolymers. , 2001 .

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

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

[20]  Ray Luo,et al.  A Poisson–Boltzmann dynamics method with nonperiodic boundary condition , 2003 .

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

[22]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[23]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[24]  Nathan A. Baker,et al.  Solvation forces on biomolecular structures: A comparison of explicit solvent and Poisson–Boltzmann models , 2004, J. Comput. Chem..

[25]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[26]  Emilio Gallicchio,et al.  On the nonpolar hydration free energy of proteins: surface area and continuum solvent models for the solute-solvent interaction energy. , 2003, Journal of the American Chemical Society.

[27]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[28]  L. Serrano,et al.  A short linear peptide that folds into a native stable β-hairpin in aqueous solution , 1994, Nature Structural Biology.

[29]  Ruhong Zhou,et al.  Comment on: Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? The Poisson-Boltzmann equation by Zhou, R. et al. , 2004 .

[30]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.