Biophysical Journal 65: 48-55 (1993)
© 1993 the Biophysical Society

This Article Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Oberoi, H Articles by Allewell, N M Search for Related Content PubMed PubMed Citation Articles by Oberoi, H Articles by Allewell, N M

Multigrid solution of the nonlinear Poisson-Boltzmann equation and calculation of titration curves.

Department of Biochemistry, University of Minnesota, St. Paul 55108.

ABSTRACT

Although knowledge of the pKa values and charge states of individual residues is critical to understanding the role of electrostatic effects in protein structure and function, calculating these quantities is challenging because of the sensitivity of these parameters to the position and distribution of charges. Values for many different proteins which agree well with experimental results have been obtained with modified Tanford-Kirkwood theory in which the protein is modeled as a sphere (reviewed in Ref. 1); however, convergence is more difficult to achieve with finite difference methods, in which the protein is mapped onto a grid and derivatives of the potential function are calculated as differences between the values of the function at grid points (reviewed in Ref. 6). Multigrid methods, in which the size of the grid is varied from fine to coarse in several cycles, decrease computational time, increase rates of convergence, and improve agreement with experiment. Both the accuracy and computational advantage of the multigrid approach increase with grid size, because the time required to achieve a solution increases slowly with grid size. We have implemented a multigrid procedure for solving the nonlinear Poisson-Boltzmann equation, and, using lysozyme as a test case, compared calculations for several crystal forms, different refinement procedures, and different charge assignment schemes. The root mean square difference between calculated and experimental pKa values for the crystal structure which yields best agreement with experiment (1LZT) is 1.1 pH units, with the differences in calculated and experimental pK values being less than 0.6 pH units for 16 out of 21 residues. The calculated titration curves of several residues are biphasic.

This article has been cited by other articles:

				X. Liu, Z.-R. Zhang, M. D. Fuller, J. Billingsley, N. A. McCarty, and D. C. Dawson CFTR: A Cysteine at Position 338 in TM6 Senses a Positive Electrostatic Potential in the Pore Biophys. J., December 1, 2004; 87(6): 3826 - 3841. [Abstract] [Full Text] [PDF]

				S. Varma and E. Jakobsson Ionization States of Residues in OmpF and Mutants: Effects of Dielectric Constant and Interactions between Residues Biophys. J., February 1, 2004; 86(2): 690 - 704. [Abstract] [Full Text] [PDF]

				S. Bransburg-Zabary, E. Nachliel, and M. Gutman Gauging of the PhoE Channel by a Single Freely Diffusing Proton Biophys. J., December 1, 2002; 83(6): 2987 - 3000. [Abstract] [Full Text] [PDF]

				X. Yuan, V. J. LiCata, and N. M. Allewell Effects of Assembly and Mutations Outside the Active Site on the Functional pH Dependence of Escherichia coli Aspartate Transcarbamylase J. Biol. Chem., January 19, 1996; 271(3): 1285 - 1294. [Abstract] [Full Text] [PDF]

				B Honig and A Nicholls Classical electrostatics in biology and chemistry Science, May 26, 1995; 268(5214): 1144 - 1149. [Abstract] [PDF]

				T. Fukamizo, A. H. Juffer, H. J. Vogel, Y. Honda, H. Tremblay, I. Boucher, W. A. Neugebauer, and R. Brzezinski Theoretical Calculation of pKaReveals an Important Role of Arg205 in the Activity and Stability of Streptomyces sp. N174 Chitosanase J. Biol. Chem., August 11, 2000; 275(33): 25633 - 25640. [Abstract] [Full Text] [PDF]