Gerald R. Grimsley

Protein ionizable groups: pK values and their contribution to protein stability and solubility.

The structure, stability, solubility, and function of proteins depend on their net charge and on the ionization state of the individual residues. Consequently, biochemists are interested in the pK values of the ionizable groups in proteins and how these pK values depend on their environment. We review what has been learned about pK values of ionizable groups in proteins from experimental studies and discuss the important contributions they make to protein stability and solubility.

A summary of the measured pK values of the ionizable groups in folded proteins

We tabulated 541 measured pK values reported in the literature for the Asp, Glu, His, Cys, Tyr, and Lys side chains, and the C and N termini of 78 folded proteins. The majority of these values are for the Asp, Glu, and His side chains. The average pK values are Asp 3.5 ± 1.2 (139); Glu 4.2 ± 0.9 (153); His 6.6 ± 1.0 (131); Cys 6.8 ± 2.7 (25); Tyr 10.3 ± 1.2 (20); Lys 10.5 ± 1.1 (35); C‐terminus 3.3 ± 0.8 (22) and N‐terminus 7.7 ± 0.5 (16). We compare these results with the measured pK values of these groups in alanine pentapeptides, and comment on our overall findings.

pK values of the ionizable groups of proteins

We have used potentiometric titrations to measure the pK values of the ionizable groups of proteins in alanine pentapeptides with appropriately blocked termini. These pentapeptides provide an improved model for the pK values of the ionizable groups in proteins. Our pK values determined in 0.1 M KCl at 25°C are: 3.67±0.03 (α‐carboxyl), 3.67±0.04 (Asp), 4.25±0.05 (Glu), 6.54±0.04 (His), 8.00±0.03 (α‐amino), 8.55±0.03 (Cys), 9.84±0.11 (Tyr), and 10.40±0.08 (Lys). The pK values of some groups differ from the Nozaki and Tanford (N&T) pK values often used in the literature: Asp (3.67 this work vs. 4.0 N&T); His (6.54 this work vs. 6.3 N&T); α‐amino (8.00 this work vs. 7.5 N&T); Cys (8.55 this work vs. 9.5 N&T); and Tyr (9.84 this work vs. 9.6 N&T). Our pK values will be useful to those who study pK perturbations in folded and unfolded proteins, and to those who use theory to gain a better understanding of the factors that determine the pK values of the ionizable groups of proteins.

The effect of net charge on the solubility, activity, and stability of ribonuclease Sa

The net charge and isoelectric pH (pI) of a protein depend on the content of ionizable groups and their pK values. Ribonuclease Sa (RNase Sa) is an acidic protein with a pI = 3.5 that contains no Lys residues. By replacing Asp and Glu residues on the surface of RNase Sa with Lys residues, we have created a 3K variant (D1K, D17K, E41K) with a pI = 6.4 and a 5K variant (3K + D25K, E74K) with a pI = 10.2. We show that pI values estimated using pK values based on model compound data can be in error by >1 pH unit, and suggest how the estimation can be improved. For RNase Sa and the 3K and 5K variants, the solubility, activity, and stability have been measured as a function of pH. We find that the pH of minimum solubility varies with the pI of the protein, but that the pH of maximum activity and the pH of maximum stability do not.

Contribution of hydrophobic interactions to protein stability.

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI = 3.5) and a basic variant (pI = 10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulombs Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulombs Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.

Forces stabilizing proteins.

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site‐directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH2− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non‐polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.

Spectrophotometric determination of protein concentration.

The concentration of a purified protein in solution is most conveniently and accurately measured using absorbance spectroscopy. The absorbance, A, is a linear function of the molar concentration, C, according to the Beer‐Lambert law: A = ɛ × l × c, where e is the molar absorption coefficient and l is the cell path length. This unit provides protocols for calculation of ɛ for a folded or unfolded protein, making use of the average ɛ values for the three contributing chromophores in proteins (the side chains of Trp, Tyr, and Cys). A basic protocol describes how to measure the concentration of a protein using the calculated ɛ and the Beer‐Lambert law. A sensitive method is provided for measuring the concentration of proteins that contain few if any tryptophan or tyrosine residues, and a simple method is provided for estimating total protein concentration in crude extracts.

Increasing protein stability by improving beta-turns.

Our goal was to gain a better understanding of how protein stability can be increased by improving β‐turns. We studied 22 β‐turns in nine proteins with 66–370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some β‐turn positions. We studied: Cold shock protein B (CspB), Histidine‐containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase α‐subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 ± 0.3; Range = 0.3–1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II β‐turns and at the i position in Type II β‐turns. Other turn positions can sometimes be used if the φ angle is near −60° for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in β‐turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in β‐turns that could be replaced by Gly to increase protein stability. Improving β‐turns by substituting Pro residues is a generally useful way of increasing protein stability. Proteins 2009.

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.