Determining the Conformational Stability of a Protein from Urea and Thermal Unfolding Curves.

This article contains protocols for determining the conformational stability of a globular protein from either urea or thermal unfolding curves. Circular dichroism is the optical spectroscopic technique most commonly used to monitor protein unfolding. These protocols describe how to analyze data from an unfolding curve to obtain the thermodynamic parameters necessary to calculate conformational stability, and how to determine differences in stability between protein variants. Curr. Protoc. Protein Sci. 71:28.4.1-28.4.14. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Determining protein conformational stability from urea-induced unfolding curves Support Protocol 1: Preparing a urea stock solution Support Protocol 2: Analyzing urea unfolding curves Basic Protocol 2: Determining the conformational stability of a protein from thermal unfolding curves Support Protocol 3: Analyzing thermal unfolding curves Support Protocol 4: Determining differences in conformational stability for protein variants.

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.

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.

Determining the conformational stability of a protein from urea and thermal unfolding curves.

This unit contains basic protocols for determining the conformational stability of a globular protein from either urea or thermal unfolding curves. Circular dichroism is the optical spectroscopic technique most commonly used to monitor protein unfolding. The protocols describe how to analyze data from an unfolding curve to obtain the thermodynamic parameters necessary to calculate conformational stability, and how to determine differences in stability between protein variants.

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.

Urea denatured state ensembles contain extensive secondary structure that is increased in hydrophobic proteins.

The goal of this article is to gain a better understanding of the denatured state ensemble (DSE) of proteins through an experimental and computational study of their denaturation by urea. Proteins unfold to different extents in urea and the most hydrophobic proteins have the most compact DSE and contain almost as much secondary structure as folded proteins. Proteins that unfold to the greatest extent near pH 7 still contain substantial amounts of secondary structure. At low pH, the DSE expands due to charge-charge interactions and when the net charge per residue is high, most of the secondary structure is disrupted. The proteins in the DSE appear to contain substantial amounts of polyproline II conformation at high urea concentrations. In all cases considered, including staph nuclease, the extent of unfolding by urea can be accounted for using the data and approach developed in the laboratory of Wayne Bolen (Auton et al., Proc Natl Acad Sci 2007; 104:15317-15323).

Increasing protein stability by improving beta-turns.

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-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 beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-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 beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-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 beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.

Determining the conformational stability of a protein using urea denaturation curves.

The stability of globular proteins is an important factor in determining their usefulness in basic research and medicine. A number of environmental factors contribute to the conformational stability of a protein, including pH, temperature, and ionic strength. In addition, variants of proteins may show remarkable differences in stability from their wild-type form. In this chapter, we describe the method and analysis of urea denaturation curves to determine the conformational stability of a protein. This involves relatively simple experiments that can be done in a typical biochemistry laboratory, especially when using ordinary spectroscopic techniques to follow unfolding.

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.

Solvent denaturation of proteins and interpretations of the m value.

The stability of globular proteins is important in medicine, proteomics, and basic research. The conformational stability of the folded state can be determined experimentally by analyzing urea, guanidinium chloride, and thermal denaturation curves. Solvent denaturation curves in particular may give useful information about a protein such as the existence of domains or the presence of stable folding intermediates. The linear extrapolation method (LEM) for analyzing solvent denaturation curves gives the parameter m, which is a measure of the dependence of ΔG on denaturant concentration. There is much recent interest in the m value as it relates to the change in accessible surface area of a protein when it unfolds and what it may reveal about the denatured states of proteins.

Tryptophan fluorescence reveals the presence of long-range interactions in the denatured state of ribonuclease Sa.

Characterizing the denatured state ensemble is crucial to understanding protein stability and the mechanism of protein folding. The aim of this research was to see if fluorescence could be used to gain new information on the denatured state ensemble. Ribonuclease Sa (RNase Sa) contains no Trp residues. We made five variants of RNase Sa by adding Trp residues at locations where they are found in other members of the microbial ribonuclease family. To better understand the protein denatured state, we also studied the fluorescence properties of the following peptides: N-acetyl-Trp-amide (NATA), N-acetyl-Ala-Trp-Ala-amide (AWA), N-acetyl-Ala-Ala-Trp-Ala-Ala-amide (AAWAA), and the five pentapeptides with the same sequence as the Trp substitution sites in RNase Sa. The major conclusions are: 1), the wavelength of maximum fluorescence intensity, lambda(max), does not differ significantly for the peptides and the denatured proteins; 2), the fluorescence intensity at lambda(max), I(F), differs significantly for the five Trp containing variants of RNase Sa; 3), the I(F) differences for the denatured proteins are mirrored in the peptides, showing that the short-range effects giving rise to the I(F) differences in the peptides are also present in the proteins; 4) the I(F) values for the denatured proteins are more than 30% greater than for the peptides, showing the presence of long-range effects in the proteins; 5), fluorescence quenching of Trp by acrylamide and iodide is more than 50% greater in the peptides than in the denatured proteins, showing that long-range effects limit the accessibility of the quenchers to the Trp side chains in the proteins; and 6), these results show that nonlocal effects in the denatured states of proteins influence Trp fluorescence and accessibility significantly.

Peptide sequence and conformation strongly influence tryptophan fluorescence.

This article probes the denatured state ensemble of ribonuclease Sa (RNase Sa) using fluorescence. To interpret the results obtained with RNase Sa, it is essential that we gain a better understanding of the fluorescence properties of tryptophan (Trp) in peptides. We describe studies of N-acetyl-L-tryptophanamide (NATA), a tripeptide: AWA, and six pentapeptides: AAWAA, WVSGT, GYWHE, HEWTV, EAWQE, and DYWTG. The latter five peptides have the same sequence as those surrounding the Trp residues studied in RNase Sa. The fluorescence emission spectra, the fluorescence lifetimes, and the fluorescence quenching by acrylamide and iodide were measured in concentrated solutions of urea and guanidine hydrochloride. Excited-state electron transfer from the indole ring of Trp to the carbonyl groups of peptide bonds is thought to be the most important mechanism for intramolecular quenching of Trp fluorescence. We find the maximum fluorescence intensities vary from 49,000 for NATA with two carbonyls, to 24,400 for AWA with four carbonyls, to 28,500 for AAWAA with six carbonyls. This suggests that the four carbonyls of AWA are better able to quench Trp fluorescence than the six carbonyls of AAWAA, and this must reflect a difference in the conformations of the peptides. For the pentapeptides, EAWQE has a fluorescence intensity that is more than 50% greater than DYWTG, showing that the amino acid sequence influences the fluorescence intensity either directly through side-chain quenching and/or indirectly through an influence on the conformational ensemble of the peptides. Our results show that peptides are generally better models for the Trp residues in proteins than NATA. Finally, our results emphasize that we have much to learn about Trp fluorescence even in simple compounds.

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins.

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.

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 degrees C are: 3.67+/-0.03 (alpha-carboxyl), 3.67+/-0.04 (Asp), 4.25+/-0.05 (Glu), 6.54+/-0.04 (His), 8.00+/-0.03 (alpha-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); alpha-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.

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 = epsilon x l x c, where e is the molar absorption coefficient and l is the cell path length. This unit provides protocols for calculation of epsilon for a folded or unfolded protein, making use of the average epsilon 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 epsilon 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.