热门搜索词：
Structure‐Based pKa Calculations Using Continuum Electrostatics Methods
关键词： structure based pka calculations来源： 互联网
- Abstract
- Table of Contents
- Figures
- Literature Cited
Abstract
Electrostatic free energy is useful for correlating structure with function in proteins in which ionizable groups play essential functional roles. To this end, the pK a values of ionizable groups must be known and their molecular determinants must be understood. Structure?based calculations of electrostatic energies and pK a values are necessary for this purpose. This unit describes protocols for pK a calculations with continuum electrostatics methods based on the numerical solution of the linearized Poisson?Boltzmann equation by the method of finite differences. Critical discussion of key parameters, approximations, and shortcomings of these methods is included. Two protocols are described for calculations with methods modified empirically to maximize agreement between measured and calculated pK a values. Applied judiciously, these methods can contribute useful and novel insight into properties of surface ionizable groups in proteins.
Keywords: pKa calculations; continuum electrostatics; finite difference; Poisson?Boltzmann; UBHD
GO TO THE FULL PROTOCOL: PDF or HTML at Wiley Online Library Table of Contents
- Basic Protocol 1: Calculating pKa Values Using the FDPB Method and the Single‐Site Charge Model (FDPB/SS)
- Alternate Protocol 1: Calculating pKa Values Using the FDPB Method and the Full Charge Model (FDPB/F)
- Guidelines for Understanding Results
- Commentary
- Literature Cited
- Figures
- Tables
GO TO THE FULL PROTOCOL: PDF or HTML at Wiley Online Library Materials
GO TO THE FULL PROTOCOL: PDF or HTML at Wiley Online Library Figures
-
View ImageFigure 8.11.1 Thermodynamic cycle for p K a calculations. Thermodynamic cycle used in the FDPB/SS method for p K a calculations. p K a model represents the p K a of an ionizable group in a model compound. p K a prot is the p K a of the group in the protein. The transfer free energies, ΔG i tr , are the calculated electrostatic free energy changes for transferring the ionizable group from water to the protein environment in the neutral ( q = 0) and ionized ( q = 1) states. Larger circles denote ionizable groups in the protein; smaller circles denote the polar atoms of the protein, which are treated in these methods in terms of partial charges. -
View ImageFigure 8.11.2 Model of the protein‐water system used for calculation of electrostatic potentials with FDPB methods. The solid line represents the van der Waal's envelope of the protein. The dashed line describes the water‐accessible surface that constitutes the boundary between the water phase with high dielectric constant and the protein phase with low dielectric constant ɛin . The dotted line represents the ion exclusion surface. A single Asp side chain is represented, with partial charges given for the atoms of the group. The grid is necessary for the solution of the Poisson‐Boltzmann equation by the method of finite differences. -
View ImageFigure 8.11.3 H+ titrations of three acidic groups calculated with the FDPB/SS method. The curves were calculated with the 1stn.pdb structure with FDPB/SS method using ɛin = 20 and ionic strength = 100 mM. The p K a value [listed under p K a (app) in Table ] represents the pH where the extent of protonation is 0.5. -
View ImageFigure 8.11.4 Overall H+ titration calculated by FDPB methods. Plot of the average charge ( Q ) of 1stn.pdb calculated with FDPB in 100 mM ionic strength: (solid line) FDPB/SS with ɛin = 20; (dashed‐dot) FDPB/F with ɛin = 4; (dashed) FDPB/SS‐HH with ɛin = 20; (dotted) calculated with the Henderson‐Hasselbalch equation using the p K a values of model compounds. -
View ImageFigure 8.11.5 Comparison of p K a (app) values of an acidic residue calculated with different FDPB methods. The set of p K a values for a representative group in staphylococcal nuclease were calculated with nine different implementations of FDPB methods to illustrate the range of values and their sensitivity to different parameters. The effects of different values of ɛin (4 versus 20), different ionic strengths (100 mM versus 1 M), different charge distribution methods (FDPB/SS versus FDPB/F), different atomic charge sets (PARSE versus CHARMm), different tautomeric states, different structures (static versus MD relaxed), and different definition of the dielectric boundary (water accessible versus van der Waal's), are compared. -
View ImageFigure 8.11.6 Energetic contributions to p K a (app) values calculated with different FDPB methods. These data illustrate how the calculated p K a values are parsed into Born (solid), background (gray), and Coulomb (white) energies by different implementations of FDPB methods.
Videos
Literature Cited
| Alexov, E. 2003. Role of the protein side‐chain fluctuations on the strength of pair‐wise electrostatic interactions: Comparing experimental with computed pK(a)s. Proteins 50:94‐103. | |
| Alexov, E.G., and Gunner, M.R. 1997. Incorporating protein conformational flexibility into the calculation of pH‐dependent protein properties. Biophys. J. 74:2075‐2093. | |
| Alexov, E.G. and Gunner, M.R. 1999. Calculated protein and proton motions coupled to electron transfer: Electron transfer from QA‐ to QB in bacterial photosynthetic reaction centers. Biochemistry 38:8253‐8270. | |
| Antosiewicz, J., McCammon, J.A., and Gilson, M.K. 1994. Prediction of pH‐dependent properties of proteins. J. Mol. Biol. 238:415‐436. | |
| Antosiewicz, J., Briggs, J.M., Elcock, A.H., Gilson, M.K., and McCammon, A. 1996a. Computing ionization states of proteins with a detailed charge model. J. Comput. Chem. 17:1633‐1644. | |
| Antosiewicz, J., McCammon, A.J., and Gilson, M.K. 1996b. The determinants of pKas in proteins. Biochemistry 35:7819‐7833. | |
| Archontis, G. and Simonson, T. 2005. Proton binding to proteins: A free‐energy component analysis using a dielectric continuum model. Biophys. J. 88:3888‐3904. | |
| Baker, N., Holst, M., and Wang, F. 2000. Adaptive multilevel finite element solution of the Poisson‐Boltzmann equation II. Refinement at solvent‐accessible surfaces in biomolecular systems. J. Comput. Chem. 21:1343‐1352. | |
| Baptista, A.M. and Soares, C.M. 2001. Some theoretical and computational aspects of the inclusion of proton isomerism in the protonation equilibrium of proteins. J. Phys. Chem. B 105:293‐309. | |
| Bashford, D. and Gerwert, K. 1992. Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsin. J. Mol. Biol. 224:473‐486. | |
| Bashford, D. and Karplus, M. 1990. pKa's of ionizable groups in proteins: Atomic detail from a continuum electrostatic model. Biochemistry 29:10219‐10225. | |
| Bashford, D. and Karplus, M. 1991. Multiple‐site titration curves of proteins: An analysis of exact and approximate methods for their calculation. J. Phys. Chem. 95:9556‐9561. | |
| Bashford, D., Case, D., Dalvit, C., Tennant, L., and Wright, P. 1993. Electrostatic calculations of side‐chain pKa values in myoglobin and comparison with NMR data for histidines. Biochemistry 32:8045‐8056. | |
| Beroza, P. and Case, D.A. 1996. Including side chain flexibility in continuum electrostatic calculations of protein titration. J. Phys. Chem. 100:20156‐20163. | |
| Beroza, P. and Case, D.A. 1998. Methods to address the change in conformation resulting from ionization process OR fluctuations inherent at a particular pH or in a particular structure. Methods Enzymol. 295:170‐189. | |
| Beroza, D., Fredkin, D.R., Okamura, M.Y., and Feher, G. 1991. Protonation of interacting residues in a protein by a Monte Carlo method: Application to lysozyme and the photosynthetic reaction center of Rhodobacter sphaeroides. Proc. Natl. Acad. Sci. U.S.A. 88:5804‐5808. | |
| Brooks, B.R., Bruccoleri, R.E., Olafson, B.D., States, D.J., Swaminathan, S., and Karplus, M. 1983. CHARMM: A Program for macromolecular energy, minimization, and dynamics calculations. J. Comput. Chem. 4:187‐217. | |
| Davis, M.E., Madura, J.D., Luty, B.A., and McCammon, J.A. 1991. Electrostatics and diffusion of molecules in solution‐Simulations with the University of Houston Brownian Dynamics program. Comp. Phys. Commun. 62:187‐197. | |
| Demchuk, E. and Wade, R.C. 1996. Improving the continuum dielectric approach to calculating pKa s of ionizable groups in proteins. J. Phys. Chem. 100:17373‐17387. | |
| Demchuk, E., Genick, U.K., Woo, T.T., Getzoff, E.D., and Bashford, D. 2000. Protonation states and pH titration in the photocycle of photoactive yellow protein. Biochemistry 39:1100‐1113. | |
| Dillet, V., Dyson, H.J., and Bashford, D. 1998. Calculations of electrostatic interactions and pK(a)s in the active site of Escherichia coli thioredoxin. Biochemistry 37:10298‐10306. | |
| Dimitrov, R.A. and Crichton, R.R. 1997. Self‐consistent field approach to protein structure and stability .1. pH dependence of electrostatic contribution. Proteins 27:576‐596. | |
| Dwyer, J.J., Gittis, A.G., Karp, D.A., Lattman, E.E., Spencer, D.S., Stites, W.E., and García‐Moreno E., B. 2000. High apparent dielectric constants in the interior of a protein reflect water penetration. Biophys. J. 79:1610‐1620. | |
| Elcock, A.H. 1999. Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability. J. Mol. Biol. 294:1051‐1062. | |
| Fitch, C.A., Karp, D.A., Lee, K.K., Stites, W.E., Lattman, E.E., and García‐Moreno E., B. 2002. Experimental pKa values of buried residues: analysis with continuum methods and role of water penetration. Biophys. J. 82:3289‐3304. | |
| Fitch, C.A., Whitten, S.T., Hilser, V.J., and García‐Moreno E., B. 2005. Molecular mechanism of pH‐driven conformational transitions of proteins: Insights from continuum electrostatics calculations of acid unfolding. Proteins 63:113‐126. | |
| Garcia‐Moreno E., B. and Fitch, C.A. 2004. Structural interpretation of pH and salt‐dependent processes in proteins with computational methods. In Energetics Of Biological Macromolecules, Pt E. (M. J.M. Holt, M.L. Johnson, and G.K. Ackers, eds.) pp. 20‐51. Academic Press Inc., San Diego. | |
| Georgescu, R.E., Alexov, E.G., and Gunner, M.R. 2002. Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. Biophys. J. 83:1731‐1748. | |
| Gibas, C.J. and Subramaniam, S. 1996. Explicit solvent models in protein pKa calculations. Biophys. J. 71:138‐147. | |
| Gilson, M.K. 1993. Multiple‐site titration and molecular modeling: Two rapid methods for computing energies and forces for ionizable groups in proteins. Proteins 15:266‐282. | |
| Gilson, M.K. 1995. Theory of electrostatic interactions in macromolecules. Curr. Biol. 5:216‐223. | |
| Gilson, M.K., Sharp, K.A., and Honig, B.H. 1988. Calculating the electrostatic potential of molecules in solution‐Method and error assessment. J. Comput. Chem. 9:327‐335. | |
| Gorfe, A.A., Ferrara, P., Caflisch, A., Marti, D.N., Bosshard, H.R., and Jelesarov, I. 2002. Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase. Proteins 46:41‐60. | |
| Harvey, S.C. and Hoekstra, P. 1972. Dielectric relaxation spectra of water adsorbed on lysozyme. J. Phys. Chem. 76:2987‐2994. | |
| Havranek, J.J. and Harbury, P.B. 1999. Tanford‐Kirkwood electrostatics for protein modeling. Proc. Natl. Acad. Sci. U.S.A. 96:11145‐11150. | |
| Holst, M., Baker, N., and Wang, F. 2000. Adaptive multilevel finite element solution of the Poisson‐Boltzmann equation I. Algorithms and examples. J. Comput. Chem. 21:1319‐1342. | |
| Jorgensen, W.L. and Tirado‐Rives, J. 1988. The OPLS potential functions for proteins. Energy minimizations for crystals of cyclic peptides and crambin. J. Am. Chem. Soc. 110:1657‐1666. | |
| Karshikoff, A. 1995. A simple algorithm for the calculation of multiple‐site titration curves. Protein Eng. 8:243‐248. | |
| Khare, D., Alexander P., Antosiewicz J., Bryan P., Gilson M., and Orban J. 1997. pKa measurements from nuclear magnetic resonance for the B1 and B2 immunoglobin G‐binding domains of protein G: Comparison with calculated values for nuclear magnetic resonance and x‐ray structures. Biochemistry 36:3580‐3589. | |
| Klapper, I., Hagstrom, R., Fine, R., Sharp, K., and Honig, B. 1986. Focussing of electric fields in the active site of Cu‐Zn superoxide dismutase: Effects of ionic strength and amino‐acid modification. Proteins 1:47‐59. | |
| Koumanov, A., Spitzner, N., Rüterjans, H., and Karshikoff, A.D. 2001. Ionization properties of titratable groups in ribonuclease T1 II. Electrostatic analysis. Eur. Biophy. J. 30:198‐206. | |
| Koumanov, A., Ruterjans, H., and Karshikoff, A. 2002. Continuum electrostatic analysis of irregular ionization and proton allocation in proteins. Proteins 46:85‐96. | |
| Krishtalik, L.I., Kuznetsov, A.M., and Mertz, E.L. 1997. Electrostatics of proteins: Description in terms of two dielectric constants simultaneously. Proteins 28:174‐182. | |
| Kuhlman, B., Luisi, D., Young, P., and Raleigh, D. 1999. pKa values and the pH dependent stability of the N‐terminal domain of L9 as probes of electrostatic interactions in the denatured state: Differentiation between local and nonlocal interactions. Biochemistry 38:4896‐4903. | |
| Kundrotas, P.J. and Karshikoff, A. 2002. Modeling of denatured state for calculation of the electrostatic contribution to protein stability. Prot. Sci. 11:1681‐1686. | |
| Lee, K.K., Fitch, C.A., Lecomte, J.T.J., and Garcia‐Moreno E., B. 2002. Electrostatic effects in highly charged proteins: Salt sensitivity of pKa values of histidines in staphylococcal nuclease. Biochemistry 41:5656‐5667. | |
| Li, H., Robertson, A.D., and Jensen, J.H. 2005. Very fast empirical prediction and rationalization of protein pKa values. Proteins 61:704‐721. | |
| Linderstrøm‐Lang, K. 1924. On the ionization of proteins. C R Trav. Lab. Carlsberg 15:1‐29. | |
| MacKerell, A.D., Bashford, D., Bellott, M., Dunbrack, R.L., Evanseck, J.D., Field, M.J., Fischer, S., Gao, J., Guo, H., Ha, S., Joseph‐McCarthy, D., Kuchnir, L., Kuczera, K., Lau, F.T.K., Mattos, C., Michnick, S., Ngo, T., Nguyen, D.T., Prodhom, B., Reiher, W.E., Roux, B., Schlenkrich, M., Smith, J.C., Stote, R., Straub, J., Watanabe, M., Wiorkiewicz‐Kuczera, J., Yin, D., and Karplus, M. 1998. All‐atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 102:3586‐3616. | |
| Madura, J.D., Briggs, J.M., Wade, R.C., Davis, M.E., Luty, B.A., Ilin, A., Antosiewicz, J., Gilson, M.K., Bagheri, B., Scott, L.R., and McCammon, J.A. 1995. Electrostatics and diffusion of molecules in solution‐Simulations with the University of Houston Brownian Dynamics program. Comp. Phys. Commun. 91:57‐95. | |
| Matthew, J.B., Gurd, F.R.N., García‐Moreno E., B., Flanagan, M.A., March, K.L., and Shire, S.J. 1985. pH‐dependent properties in proteins. CRC Crit. Rev. Biochem. 18:91‐197. | |
| Mehler, E.L. and Guarnieri, F. 1999. A self‐consistent, microenvironment modulated screened coulomb potential approximation to calculate pH‐dependent electrostatic effects in proteins. Biophys. J. 75:3‐22. | |
| Neria, E., Fischer, S., and Karplus, M. 1996. Simulation of activation free energies in molecular systems. J. Chem. Phys. 105:1902‐1921. | |
| Nielsen, J.E. and Vriend, G. 2001. Optimizing the hydrogen‐bond network in Poisson‐Boltzmann equation‐based pK(a) calculations. Proteins 43:403‐412. | |
| Nielsen, J.E., Andersen, K.V., Honig, B., Hooft, R.W.W., Klebe, G., Vriend, G., and Wade, R.C. 1999. Improving macromolecular electrostatics calculations. Prot. Eng. 12:657‐662. | |
| Oberoi, H. and Allewell, N.M. 1993. Multigrid solution of the nonlinear Poisson‐Boltzmann Equation and calculation of titration curves. Biophys. J. 65:48‐55. | |
| Ondrechen, M.J., Clifton, J.G., and Ringe, D. 2001. THEMATICS: A simple computational predictor of enzyme structure from function. Proc. Natl. Acad. Sci. U.S.A. 98:12473‐12478. | |
| Richards, F.M. 1977. Areas, volumes, packing, and protein structure. Annu. Rev. Biophys. Bioeng. 6:151‐176. | |
| Schaefer, M., Van Vlijmen, H.W.T., and Karplus, M. 1998. Electrostatic contributions to molecular free energies in solution. In Advances In Protein Che |
推荐方法
移动版：