A fast and accurate computational approach to protein ionization
- 1 November 2008
- journal article
- Published by Wiley in Protein Science
- Vol. 17 (11), 1955-1970
- https://doi.org/10.1110/ps.036335.108
Abstract
We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.Keywords
This publication has 75 references indexed in Scilit:
- Accurate, conformation-dependent predictions of solvent effects on protein ionization constantsProceedings of the National Academy of Sciences of the United States of America, 2007
- The dominant role of side‐chain backbone interactions in structural realization of amino acid code. ChiRotor: A side‐chain prediction algorithm based on side‐chain backbone interactionsProtein Science, 2007
- pK values of the ionizable groups of proteinsProtein Science, 2006
- Molecular Dynamics Simulations of Large Integral Membrane Proteins with an Implicit Membrane ModelThe Journal of Physical Chemistry B, 2005
- Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientationJournal of Molecular Biology, 1999
- Structural Origins of pH and Ionic Strength Effects on Protein Stability: Acid Denaturation of Sperm Whale ApomyoglobinJournal of Molecular Biology, 1994
- Validation of the general purpose QUANTA ®3.2/CHARMm® force fieldJournal of Computational Chemistry, 1992
- Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsinJournal of Molecular Biology, 1992
- 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