Protein pKa calculations
From Wikipedia, the free encyclopedia
In computational biology, protein pKa calculations are used to estimate the pKa values of amino acids as they exist within proteins. These calculations complement the pKa values reported for amino acids in their free state, and are used frequently within the fields of molecular modeling, structural bioinformatics, and computational biology.
Contents |
[edit] Amino acid pKa values
pKa values of amino acid side chains play an important role in defining the pH-dependent characteristics of a protein. The pH-dependence of the activity displayed by enzymes and the pH-dependence of protein stability, for example, are properties that are determined by the pKa values of amino acid side chains.
The pKa values of an amino acid side chain in solution is typically inferred from the pKa values of model compounds (compounds that are similar to the side chains of amino acids). (See Amino acid for the pKa values of all amino acid side chains inferred in such a way.) The table below lists the model pKa values that are normally used in a protein pKa calculation.
| Amino Acid | pKa |
|---|---|
| Asp | 4.0 |
| Glu | 4.4 |
| Arg | 12.0 |
| Lys | 9.0 |
| His | 6.08 |
| Cys | 8.28 (-SH) |
| Tyr | 9.84 |
[edit] pKa calculation methods
There are several methods available for calculating protein pKa values:
[edit] Using the Poisson-Boltzmann equation
Some methods are based on solutions to the Poisson-Boltzmann equation (PBE), often referred to as FDPB-based methods (FDPB is for "finite difference Poisson-Boltzmann"). The PBE is a modification of Poisson's equation that incorporates a description of the effect of solvent ions on the electrostatic field around a molecule.
FDPB-based methods calculate the change in the pKa value of an amino acid side chain when that side chain is moved from a hypothetical fully solvated state to its position in the protein. To perform such a calculation, one needs theoretical methods that can calculate the effect of the protein interior on a pKa value, and knowledge of the pKa values of amino acid side chains in their fully solvated state.
[edit] Theoretical methods
The pKa values of amino acid side chains in their fully solvated state are often inferred from comparisons with pKa values of so-called model compounds.
[edit] Empirical methods
A set of empirical rules relating the protein structure to the pKa values of ionizable residues have been developed by Li, Robertson, and Jensen. These rules form the basis for the web-accessible program called PROPKA for rapid predictions of pKa values.
[edit] Molecular dynamics (MD)-based methods
Molecular dynamics methods of calculating pKa values involve computationally measuring the free energy difference between the protonated and deprotonated forms of the molecule. This free energy differences is measured using methods such as free energy perturbation, thermodynamic integration and the Bennett acceptance ratio.
Molecular dynamics is typically a much more computationally expensive way to predict pKa's than using the Poisson Boltzmann equation. Furthermore, it is currently much less accurate. This is because currently used molecular force fields do not take polarisability into account. Polarisability is an important property for protonation energies.
[edit] Using pH titration curve
The pH value at half the equivalence point is equal to the pKa. Thus when the pH of a solution is at the pKa of an acid HA, half will be present in the deprotonated (A-) form and half will be protonated (HA). This is immediately evident from the Henderson-Hasselbalch equation.
In proteins there are often many ionazable sites in the same molecule, which makes it much more difficult to interpret experimental evidence.
 |
This article is uncategorized.
Please help improve this article by adding it to one or more categories, so it may be associated with related articles (how?). Please remove this tag after categorizing. |
