C. Nick Pace

A Helix Propensity Scale Based on Experimental Studies of Peptides and Proteins

The average globular protein contains 30% alpha-helix, the most common type of secondary structure. Some amino acids occur more frequently in alpha-helices than others; this tendency is known as helix propensity. Here we derive a helix propensity scale for solvent-exposed residues in the middle positions of alpha-helices. The scale is based on measurements of helix propensity in 11 systems, including both proteins and peptides. Alanine has the highest helix propensity, and, excluding proline, glycine has the lowest, approximately 1 kcal/mol less favorable than alanine. Based on our analysis, the helix propensities of the amino acids are as follows (kcal/mol): Ala = 0, Leu = 0.21, Arg = 0.21, Met = 0.24, Lys = 0.26, Gln = 0.39, Glu = 0.40, Ile = 0.41, Trp = 0.49, Ser = 0.50, Tyr = 0. 53, Phe = 0.54, Val = 0.61, His = 0.61, Asn = 0.65, Thr = 0.66, Cys = 0.68, Asp = 0.69, and Gly = 1.

CONFORMATIONAL STABILITY OF GLOBULAR PROTEINS

The conformational stability of ribonuclease T1 has been measured as a function of the variables of most interest to biochemists: temperature, pH, salt concentration, disulfide-bond content and amino acid sequence. The results provide insight into the forces that stabilize globular proteins.

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.

Linear extrapolation method of analyzing solvent denaturation curves.

The two most common methods of measuring the conformational stability of a protein are differential scanning calorimetry and an analysis of solvent denaturation curves by using the linear extrapolation method. In this article, we trace the history of the linear extrapolation method, review how the method is used to measure protein stability, and then discuss some of the other important uses. Proteins 2000;Suppl 4:1–7.

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.

Contribution of the hydrophobic effect to globular protein stability

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.

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.