Protein molecular dynamics with the generalized born/ACE solvent model
- 22 August 2001
- journal article
- research article
- Published by Wiley in Proteins-Structure Function and Bioinformatics
- Vol. 45 (2) , 144-158
- https://doi.org/10.1002/prot.1134
Abstract
Implicit solvent models are increasingly important for the study of proteins in aqueous solution. Here, the generalized Born (GB) solvent polarization model as implemented in the analytical ACE potential [Schaefer and Karplus (1996) J Phys Chem 100:1578] is used to perform molecular dynamics simulations of two small, homologous proteins: the immunoglobulin‐binding domain of streptococcal protein G and the Ras binding domain of Raf. Several model parameterizations are compared through more than 60 ns of simulation. Results are compared with two simpler solvent models—an accessible surface area model and a distant‐dependent dielectric model, with finite‐difference Poisson calculations, with existing explicit solvent simulations, and with experimental data. The simpler models yield stable but distorted structures. The best GB/ACE implementation uses a set of atomic Voronoi volumes reported recently, obtained by averaging over a large database of crystallographic protein structures. A 20% reduction is applied to the volumes, compensating in an average sense for an excessive de‐screening of individual charges inherent in the ACE self‐energy and for an undersolvation of dipolar groups inherent in the GB screening function. This GB/ACE parameterization yields stable trajectories on the 0.5–1‐ns time scale that deviate moderately (∼1.5–2.5 Å) from the X‐ray structure, reproduce approximately the surface distribution of charged, polar, and hydrophobic groups, and reproduce accurately backbone flexibility as measured by amide NMR‐order parameters. Over longer time scales (1.5–3 ns), some of the protein G runs escape from the native energy basin and deviate strongly (3 Å) from the native structure. The conformations sampled during the transition out of the native energy basin are overstabilized by the GB/ACE solvation model, as compared with a numerical treatment of the full dielectric continuum model. Proteins 2001;45:144–158.Keywords
This publication has 68 references indexed in Scilit:
- Environment and exposure to solvent of protein atoms. Lysozyme and insulinPublished by Elsevier ,2004
- Generalized Born Model Based on a Surface Integral FormulationThe Journal of Physical Chemistry B, 1998
- Solution conformations and thermodynamics of structured peptides: molecular dynamics simulation with an implicit solvation modelJournal of Molecular Biology, 1998
- A fast cluster counting algorithm for percolation on and off latticesComputer Physics Communications, 1998
- Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian Dynamics programComputer Physics Communications, 1995
- A 500 ps molecular dynamics simulation study of interleukin-1β in waterJournal of Molecular Biology, 1992
- A precise analytical method for calculating the electrostatic energy of macromolecules in aqueous solutionJournal of Molecular Biology, 1990
- Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparisonProteins-Structure Function and Bioinformatics, 1988
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Volumen und Hydratationswärme der IonenThe European Physical Journal A, 1920