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.

Toward a molecular understanding of protein solubility: increased negative surface charge correlates with increased solubility.

Protein solubility is a problem for many protein chemists, including structural biologists and developers of protein pharmaceuticals. Knowledge about how intrinsic factors influence solubility is limited due to the difficulty of obtaining quantitative solubility measurements. Solubility measurements in buffer alone are difficult to reproduce, because gels or supersaturated solutions often form, making it impossible to determine solubility values for many proteins. Protein precipitants can be used to obtain comparative solubility measurements and, in some cases, estimations of solubility in buffer alone. Protein precipitants fall into three broad classes: salts, long-chain polymers, and organic solvents. Here, we compare the use of representatives from two classes of precipitants, ammonium sulfate and polyethylene glycol 8000, by measuring the solubility of seven proteins. We find that increased negative surface charge correlates strongly with increased protein solubility and may be due to strong binding of water by the acidic amino acids. We also find that the solubility results obtained for the two different precipitants agree closely with each other, suggesting that the two precipitants probe similar properties that are relevant to solubility in buffer alone.

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.

Hydrogen bonding of beta-turn structure is stabilized in D(2)O.

The lower critical solution temperature (LCST) of elastin-like polypeptides (ELPs) was investigated as a function of ELP chain length and guest residue chemistry. These measurements were made in both D(2)O and H(2)O. Differences in the LCST values with heavy and light water were correlated with secondary structure formation of the polypeptide chains. Such structural information was obtained by circular dichroism and infrared measurements. Additional thermodynamic data were obtained by differential scanning calorimetry. It was found that there is a greater change in the LCST value between H(2)O and D(2)O for those polypeptides which form the greatest amount of beta-turn/beta-aggregate structure. Moreover, these same molecules were the least hydrophobic ELPs. Therefore, hydrogen bonding rather than hydrophobicity was the key factor in the stabilization of the collapsed state of ELPs in D(2)O compared with H(2)O.

Carbodiimide EDC induces cross-links that stabilize RNase A C-dimer against dissociation: EDC adducts can affect protein net charge, conformation, and activity.

RNase A self-associates under certain conditions to form a series of domain-swapped oligomers. These oligomers show high catalytic activity against double-stranded RNA and striking antitumor actions that are lacking in the monomer. However, the dissociation of these metastable oligomers limits their therapeutic potential. Here, a widely used conjugating agent, 1-ethyl-3-(3-dimethylaminoisopropyl) carbodiimide (EDC), has been used to induce the formation of amide bonds between carboxylate and amine groups of different subunits of the RNase A C-dimer. A cross-linked C-dimer which does not dissociate was isolated and was found have augmented enzymatic activity toward double-stranded RNA relative to the unmodified C-dimer. Characterization using chromatography, electrophoresis, mass spectrometry, and NMR spectroscopy revealed that the EDC-treated C-dimer retains its structure and contains one to three novel amide bonds. Moreover, both the EDC-treated C-dimer and EDC-treated RNase A monomer were found to carry an increased number of positive charges (about 6 ± 2 charges per subunit). These additional positive charges are presumably due to adduct formation with EDC, which neutralizes a negatively charged carboxylate group and couples it to a positively charged tertiary amine. The increased net positive charge endowed by EDC adducts likely contributes to the heightened cleavage of double-stranded RNA of the EDC-treated monomer and EDC-treated C-dimer. Further evidence for EDC adduct formation is provided by the reaction of EDC with a dipeptide Ac-Asp-Ala-NH(2) monitored by NMR spectroscopy and mass spectrometry. To determine if EDC adduct formation with proteins is common and how this affects protein net charge, conformation, and activity, four well-characterized proteins, ribonuclease Sa, hen lysozyme, carbonic anhydrase, and hemoglobin, were incubated with EDC and the products were characterized. EDC formed adducts with all these proteins, as judged by mass spectrometry and electrophoresis. Moreover, all suffered conformational changes ranging from slight structural modifications in the case of lysozyme, to denaturation for hemoglobin as measured by NMR spectroscopy and enzyme assays. We conclude that EDC adduct formation with proteins can affect their net charge, conformation, and enzymatic activity.

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.

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.

RNase-induced apoptosis: fate of calcium-activated potassium channels.

The connection between the action of microbial RNases and Ca2+-activated K+ (KCa) channels was investigated in human embryo kidney cells HEKhSK4 artificially expressing the channels. These channels protected HEKhSK4 cells from apoptosis induced by binase and 5K charge reversal mutant of RNase Sa. After the first 24h, potassium current increased without increase in intracellular Ca2+, and mitochondrial potential remained high. After 72 h, the concentration of calcium increased and mitochondria lost their potential. Whole-cell recordings of membrane currents through KCa channels in RNase-treated cells demonstrated a biphasic pattern: initially their activity in cell population increased, peaked at 24h, and then gradually decreased. In each individual cell we observed either an increase of the amplitude of KCa current, or a complete shutdown of the channels. The activity of KCa channels could be restored by removing RNases from the media. Based on this pattern and especially its timing, we hypothesize that toxic RNases downregulate KCa channels at the level of transcription or translation. Our results indicate that new anticancer agents could be created on the basis of microbial RNases targeting KCa channels.

Amino acid contribution to protein solubility: Asp, Glu, and Ser contribute more favorably than the other hydrophilic amino acids in RNase Sa.

Poor protein solubility is a common problem in high-resolution structural studies, formulation of protein pharmaceuticals, and biochemical characterization of proteins. One popular strategy to improve protein solubility is to use site-directed mutagenesis to make hydrophobic to hydrophilic mutations on the protein surface. However, a systematic investigation of the relative contributions of all 20 amino acids to protein solubility has not been done. Here, 20 variants at the completely solvent-exposed position 76 of ribonuclease (RNase) Sa are made to compare the contributions of each amino acid. Stability measurements were also made for these variants, which occur at the i+1 position of a type II beta-turn. Solubility measurements in ammonium sulfate solutions were made at high positive net charge, low net charge, and high negative net charge. Surprisingly, there was a wide range of contributions to protein solubility even among the hydrophilic amino acids. The results suggest that aspartic acid, glutamic acid, and serine contribute significantly more favorably than the other hydrophilic amino acids especially at high net charge. Therefore, to increase protein solubility, asparagine, glutamine, or threonine should be replaced with aspartic acid, glutamic acid or serine.

Increasing protein conformational stability by optimizing beta-turn sequence.

Protein conformational stability is an important concern in many fields. Here, we describe a strategy for significantly increasing conformational stability by optimizing beta-turn sequence. Proline and glycine residues are statistically preferred at several beta-turn positions, presumably because their unique side-chains contribute favorably to conformational stability in certain beta-turn positions. However, beta-turn sequences often deviate from preferred proline or preferred glycine. Therefore, our strategy involves replacing non-proline and non-glycine beta-turn residues with preferred proline or preferred glycine residues. Here, we develop guidelines for selecting appropriate mutations, and present results for five mutations (S31P, S42G, S48P, T76P, and Q77G) that significantly increase the conformational stability of RNase Sa. The increases in stability ranged from 0.7 kcal/mol to 1.3 kcal/mol. The strategy was successful in overlapping or isolated beta-turns, at buried (up to 50%) or completely exposed sites, and at relatively flexible or inflexible sites. Considering the significant number of beta-turn residues in every globular protein and the frequent deviation of beta-turn sequences from preferred proline and preferred glycine residues, this simple, efficient strategy will be useful for increasing the conformational stability of proteins.

NMR study and molecular dynamics simulations of optimized beta-hairpin fragments of protein G.

The stability and structure of several beta-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the beta-turn conformation. Finally, we show that the effects of combining these single mutations are additive, suggesting that independent stabilizing interactions can be isolated and evaluated in a simple model system. Our results indicate that the structure and stability of this beta-hairpin peptide can be modulated in numerous ways and thus contributes toward a more complete understanding of this important model beta-hairpin as well as to the folding and stability of larger peptides and proteins.

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.

Single-molecule electrophoresis of beta-hairpin peptides by electrical recordings and Langevin dynamics simulations.

We used single-channel electrical recordings and Langevin molecular dynamics simulations to explore the electrophoretic translocation of various beta-hairpin peptides across the staphylococcal alpha-hemolysin (alphaHL) protein pore at single-molecule resolution. The beta-hairpin peptides, which varied in their folding properties, corresponded to the C terminal residues of the B1 domain of protein G. The translocation time was strongly dependent on the electric force and was correlated with the folding features of the beta-hairpin peptides. Highly unfolded peptides entered the pore in an extended conformation, resulting in fast single-file translocation events. In contrast, the translocation of the folded beta-hairpin peptides occurred more slowly. In this case, the beta-hairpin peptides traversed the alphaHL pore in a misfolded or fully folded conformation. This study demonstrates that the interaction between a polypeptide and a beta-barrel protein pore is dependent on the folding features of the polypeptide.

Lessons in stability from thermophilic proteins.

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.

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.

Hydrogen bonding increases packing density in the protein interior.

The contribution of hydrogen bonds and the burial of polar groups to protein stability is a controversial subject. Theoretical studies suggest that burying polar groups in the protein interior makes an unfavorable contribution to the stability, but experimental studies show that burying polar groups, especially those that are hydrogen bonded, contributes favorably to protein stability. Understanding the factors that are not properly accounted for by the theoretical models would improve the models so that they more accurately describe experimental results. It has been suggested that hydrogen bonds may contribute to protein stability, in part, by increasing packing density in the protein interior, and thereby increasing the contribution of van der Waals interactions to protein stability. To investigate the influence of hydrogen bonds on packing density, we analyzed 687 crystal structures and determined the volume of buried polar groups as a function of their extent of hydrogen bonding. Our findings show that peptide groups and polar side chains that form hydrogen bonds occupy a smaller volume than the same groups when they do not form hydrogen bonds. For example, peptide groups in which both polar groups are hydrogen bonded occupy a volume, on average, 5.2 A3 less than a peptide group that is not hydrogen bonded.

Terminal ion pairs stabilize the second beta-hairpin of the B1 domain of protein G.

The effects of terminal ion pairs on the stability of a beta-hairpin peptide corresponding to the C-terminal residues of the B1 domain of protein G were determined using thermal unfolding as monitored by nuclear magnetic resonance and circular dichroism spectroscopy. Molecular dynamics (MD) simulations were also performed to examine the effect of ion pairs on the structures. Eight peptides were studied including the wild type (G41) and the N-terminal modified sequences that had the first residue deleted (E42), replaced with a Lys (K41), or extended by an additional Gly (G40). Acetylated variants were made to examine the effect of removing the positive N-terminal charge on beta-hairpin stability. The rank in stability determined experimentally is K41 > E42 approximately G41 approximately G40 > Ac-K41 > Ac-E42 approximately Ac-G41 > Ac-G40. The Tm of the K41 peptide is 12 degrees C higher than G41, while the Tm values for the acetylated peptides are less than their unacetylated forms by more than 15 degrees C. NOE cross-peaks between side-chain methylene groups at the N- and C-termini and larger CalphaH shifts compared to random values are seen for K41. The addition of 20% methanol increases the stability in K41 and G41. The MD studies complement these results by showing that the charged N-terminus is important to stability. The type of ion pair observed varies with peptide, and when formed the simulations show that the ion pair can prevent fraying of the beta-strands through electrostatic and hydrophobic contacts. Therefore, introducing favorable electrostatic interactions at the N- and C-termini can substantially enhance beta-hairpin stability and help define the structure.

The HPr proteins from the thermophile Bacillus stearothermophilus can form domain-swapped dimers.

The study of proteins from extremophilic organisms continues to generate interest in the field of protein folding because paradigms explaining the enhanced stability of these proteins still elude us and such studies have the potential to further our knowledge of the forces stabilizing proteins. We have undertaken such a study with our model protein HPr from a mesophile, Bacillus subtilis, and a thermophile, Bacillus stearothermophilus. We report here the high-resolution structures of the wild-type HPr protein from the thermophile and a variant, F29W. The variant proved to crystallize in two forms: a monomeric form with a structure very similar to the wild-type protein as well as a domain-swapped dimer. Interestingly, the structure of the domain-swapped dimer for HPr is very different from that observed for a homologous protein, Crh, from B.subtilis. The existence of a domain-swapped dimer has implications for amyloid formation and is consistent with recent results showing that the HPr proteins can form amyloid fibrils. We also characterized the conformational stability of the thermophilic HPr proteins using thermal and solvent denaturation methods and have used the high-resolution structures in an attempt to explain the differences in stability between the different HPr proteins. Finally, we present a detailed analysis of the solution properties of the HPr proteins using a variety of biochemical and biophysical methods.

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa.

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.