This example was prepared in collaboration with David Sept for a biomolecular modeling class.
OverviewAlthough pKa calculations may not seem like the most routine applications for demonstrating continuum electrostatics concepts, they have important scientific and teaching value. From a scientific standpoint, pKa values are important determinants of biomolecular (particularly enzymatic) function and can be used to assess functional activity and identify active sites. From a pedagogical perspective, pKa calculations require all of the important continuum electrostatics concepts and therefore serve as a "selfcontained" introduction to solvation and binding energies.
This tutorial covers PoissonBoltzmann methods for determining biomolecule pKa values. Other methods for pKa and titration state determination are given in the PDB2PQR examples. Comparison with PDB2PQR results can provide hours of additional entertainment.
This is a very brief introduction to the concepts behind biomolecular pKas and titration states. More information can be obtained from most biochemistry or biophysics textbooks as well as some of the original articles on pKa evaluation:
Recall that the acid dissociation constant K_{a} describes the dissocation of an acid into its components
where the equilibrium constant is defined related to the activities of the species:
which, under conditions of "ideality" (never realized in a biological system...), can be replaced with concentrations to give
You should also recall that chemical equilibrium constants can be related to free energies by
However, chemists found it easier to use quantities in base10, therefore this energetic quantity is often referred to as a pKa
Amino acid model pKa valuesIn many calculations, pKa values are assigned based on model values for amino acid side chains to mimic the reaction
http://apbs.wustl.edu/MediaWiki/images/8/81/Aciddissociationequilibriummodel.png
for an "isolated" amino acid in solution. Some of these model pKa values are given in the following table:
This data is from Nielsen JE, Vriend G. Optimizing the hydrogenbond network in PoissonBoltzmann equationbased pKa calculations. Proteins. 43 (4), 40312, 2001. As we'll see in the next section, these model values provide the basis for calculating pKa values in proteins.
pKa values in proteinsThe role of the model pKa values introduced in the previous sections is to move all the chemical (bond making and breaking) complexity of protonation into the model values. In particular, pKa values in proteins are calculated as pertubations of the model compounds according to the free energy cycles shown below. The pKa for the amino acid in the context of the protein is given by the free energy cycle:
We are interested in determining the unknown \Delta_a G_{\mathrm{HA}} from the known model \Delta_a G_{\mathrm{HA,model}}, and the unknown \Delta_{\mathrm{xfer}} G_{\mathrm{HA}} and \Delta_{\mathrm{xfer}} G_{\mathrm{A}^} according to:
or
In general, the quantities \Delta_{\mathrm{xfer}} G_{\mathrm{HA}} and \Delta_{\mathrm{xfer}} G_{\mathrm{A}^} are obtained from a computational approach. Nearly any free energy calculation method could be used to obtain these energies according to a scheme where the (de)solvation energies of the charged and uncharged amino acids are calculated according to:
Continuum electrostatics methods for pKa calculations in proteinsAlthough nearly any free energy method could be used to evaluate the energies of transferring the protonated and unprotonated amino acids from solution into the protein environment, continuum electrostatics offer a (usually) satisfying compromise between accuracy and computational efficiency.
The transfer free energies to be calculated, \Delta_{\mathrm{xfer}} G_{\mathrm{HA}} and \Delta_{\mathrm{xfer}} G_{\mathrm{A}^}, can be determined from PoissonBoltzmann (PB) energies. In particular, these energies can be calculated as effective “binding energy calculations” similar to those covered in the Binding energies section:
where
Note that, as with binding energies, \Delta_{\mathrm{xfer}} G_{X} can be evaluated two ways:
Both of these methods can be combined, via free energy cycles, to give the desired . However, when care is taken to use the same grids and conformations for all calculations, the direct method using total electrostatic energies is usually the most efficient.
Note that none of the methods discussed above have explicitly allowed for changes in titration state of other groups in the protein during protonation/deprotonation of the acid group of interest. Additionally, none of these methods explicitly provide for conformational changes in the protein coupled to protonation/deprotonation. As such, we are not true pKa values with this method. Instead, we are calculating socalled "intrinsic pKa" values.
ExamplesWe currently provide a lysozyme intrinsic pKa example to illustrate these principles.
