Protein stability curves

The stability curve of a protein is defined as the plot of the free energy of unfolding as a function of temperature. For most proteins the change in heat capacity on denaturation, or unfolding, is large but approximately constant. When unfolding is s two-state process, most of the salient features of the stability curves of proteins can be derived from this fact. A number of relations are obtained, including the special features of low-temperature denaturation, the properties of the maximum in stability, and the interrelationships of the characteristic temperatures of the protein. The paper closes with a formula that permits one to calculate small changes in stabilization free energy from changes in the melting temperature of the protein.     

A thermodynamic comparison of mesophilic and thermophilic ribonucleases H

The mechanisms by which thermophilic proteins attain their increased thermostability remain unclear, as usually the sequence and structure of these proteins are very similar to those of their mesophilic homologues. To gain insight into the basis of thermostability, we have determined protein stability curves describing the temperature dependence of the free energy of unfolding for two ribonucleases H, one from the mesophile Escherichia coli and one from the thermophile Thermus thermophilus. The circular dichroism signal was monitored as a function of temperature and guanidinium chloride concentration, and the resulting free energies of unfolding were fit to the Gibbs-Helmholtz equation to obtain a set of thermodynamic parameters for these proteins. Although the maximal stabilities for these proteins occur at similar temperatures, the heat capacity of unfolding for T. thermophilus RNase H is lower, resulting in a smaller temperature dependence of the free energy of unfolding and therefore a higher thermal melting temperature. In addition, the stabilities of these proteins are similar at the optimal growth temperatures for their respective organisms, suggesting that a balance of thermodynamic stability and flexibility is important for function.     

Role of entropy in protein thermostability: folding kinetics of a hyperthermophilic cold shock protein at high temperatures using 19F NMR

We used (19)F NMR to extend the temperature range accessible to detailed kinetic and equilibrium studies of a hyperthermophilic protein. Employing an optimized incorporation strategy, the small cold shock protein from the bacterium Thermotoga maritima (TmCsp) was labeled with 5-fluorotryptophan. Although chaotropically induced unfolding transitions revealed a significant decrease in the stabilization free energy upon fluorine labeling, the protein's kinetic folding mechanism is conserved. Temperature- and guanidinium chloride-dependent equilibrium unfolding transitions monitored by (19)F NMR agree well with the results from optical spectroscopy, and provide a stringent test of the two-state folding character of TmCsp. Folding and unfolding rate constants at high temperatures were determined from the (19)F NMR spectra close to the midpoint of thermal unfolding by global line shape analysis. In combination with results from stopped-flow experiments at lower temperatures, they show that the folding rate constant of TmCsp and its temperature dependence closely resemble those of its mesophilic homologue from Bacillus subtilis, BsCspB. However, the unfolding rate constant of TmCsp is two orders of magnitude lower over the entire temperature range that was investigated. Consequently, the difference in conformational stability between the two proteins is solely due to the unfolding rate constant over a wide temperature range. A thermodynamic analysis points to an important role of entropic factors in the stabilization of TmCsp relative to its mesophilic homologues.     

Thermodynamics of the temperature-induced unfolding of globular proteins.

The heat capacity, enthalpy, entropy, and Gibbs energy changes for the temperature-induced unfolding of 11 globular proteins of known three-dimensional structure have been obtained by microcalorimetric measurements. Their experimental values are compared to those we calculate from the change in solvent-accessible surface area between the native proteins and the extended polypeptide chain. We use proportionality coefficients for the transfer (hydration) of aliphatic, aromatic, and polar groups from gas phase to aqueous solution, we estimate vibrational effects, and we discuss the temperature dependence of each constituent of the thermodynamic functions. At 25 degrees C, stabilization of the native state of a globular protein is largely due to two favorable terms: the entropy of non-polar group hydration and the enthalpy of interactions within the protein. They compensate the unfavorable entropy change associated with these interactions (conformational entropy) and with vibrational effects. Due to the large heat capacity of nonpolar group hydration, its stabilizing contribution decreases quickly at higher temperatures, and the two unfavorable entropy terms take over, leading to temperature-induced unfolding.     

Thermodynamic stability and folding of proteins from hyperthermophilic organisms

Life grows almost everywhere on earth, including in extreme environments and under harsh conditions. Organisms adapted to high temperatures are called thermophiles (growth temperature 45-75 degrees C) and hyperthermophiles (growth temperature >or= 80 degrees C). Proteins from such organisms usually show extreme thermal stability, despite having folded structures very similar to their mesostable counterparts. Here, we summarize the current data on thermodynamic and kinetic folding/unfolding behaviors of proteins from hyperthermophilic microorganisms. In contrast to thermostable proteins, rather few (i.e. less than 20) hyperthermostable proteins have been thoroughly characterized in terms of their in vitro folding processes and their thermodynamic stability profiles. Examples that will be discussed include co-chaperonin proteins, iron-sulfur-cluster proteins, and DNA-binding proteins from hyperthermophilic bacteria (i.e. Aquifex and Theromotoga) and archea (e.g. Pyrococcus, Thermococcus, Methanothermus and Sulfolobus). Despite the small set of studied systems, it is clear that super-slow protein unfolding is a dominant strategy to allow these proteins to function at extreme temperatures.     

Enthalpic and entropic contributions mediate the role of disulfide bonds on the conformational stability of interleukin-4

The role of disulfide bridges in the structure, stability, and folding pathways of proteins has been the subject of wide interest in the fields of protein design and engineering. However, the relative importance of entropic and enthalpic contributions for the stabilization of proteins provided by disulfides is not always clear. Here, we perform a detailed analysis of the role of disulfides in the conformational stability of human Interleukin-4 (IL4), a four-helix bundle protein. In order to evaluate the contribution of two out of the three disulfides to the structure and stability of IL4, two IL4 mutants, C3T-IL4 and C24T-IL4, were used. NMR and ANS binding experiments were compatible with altered dynamics and an increase of the nonpolar solvent-accessible surface area of the folded state of the mutant proteins. Chemical and thermal unfolding experiments followed by fluorescence and circular dichroism revealed that both mutant proteins have lower conformational stability than the wild-type protein. Transition temperatures of unfolding decreased 14 degrees C for C3T-IL4 and 10 degrees C for C24T-IL4, when compared to WT-IL4, and the conformational stability, at 25 degrees C, decreased 4.9 kcal/mol for C3T-IL4 and 3.2 kcal/mol for C24T-IL4. Interestingly, both the enthalpy and the entropy of unfolding, at the transition temperature, decreased in the mutant proteins. Moreover, a smaller change in heat capacity of unfolding was also observed for the mutants. Thus, disulfide bridges in IL4 play a critical role in maintaining the thermodynamic stability and core packing of the helix bundle.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Thermal stability of wild type and disulfide bridge containing mutant of poplar plastocyanin

A comparative study of the thermal stability of wild type poplar plastocyanin and of a mutant form containing a disulfide bridge between residues 21 and 25 was performed using differential scanning calorimetry and optical spectroscopic techniques. For wild type plastocyanin the transition temperature, determined from the calorimetric profiles, is 62.7 degrees C at the scan rate of 60 degrees C/h, whereas for the mutant it is reduced to 58.0 degrees C. In both cases, the endothermic peak is followed by an exothermic one at higher temperatures. The unfolding process monitored by optical absorption at 596 nm also reveals a reduced thermal stability of the mutated plastocyanin compared to the wild type protein, with transition temperatures of 54.8 and 58.0 degrees C, respectively. For both proteins, the denaturation process was found to be irreversible and dependent on the scan rate preventing the thermodynamic analysis of the unfolding process. In parallel, small conformational changes between wild type and mutant plastocyanin emerge from fluorescence spectroscopy measurements. Here, a difference in the interaction of the two proteins between the microenvironment surrounding the fluorophores and the solvent was proposed. The destabilization observed in the disulfide containing mutant of plastocyanin suggests that the double mutation, Ile21Cys and Glu25Cys, introduces strain into the protein which offsets the stabilizing effect expected from the formation of a covalent crosslink.     

Thermodynamics and mechanism of cutinase stabilization by trehalose

Trehalose has been widely used to stabilize cellular structures such as membranes and proteins. The effect of trehalose on the stability of the enzyme cutinase was studied. Thermal unfolding of cutinase reveals that trehalose delays thermal unfolding, thus increasing the temperature at the midpoint of unfolding by 7.2 degrees . Despite this stabilizing effect, trehalose also favors pathways that lead to irreversible denaturation. Stopped-flow kinetics of cutinase folding and unfolding was measured and temperature was introduced as experimental variable to assess the mechanism and thermodynamics of protein stabilization by trehalose. The main stabilizing effect of trehalose was to delay the rate constant of the unfolding of an intermediate. A full thermodynamic analysis of this step has revealed that trehalose induces the phenomenon of entropy-enthalpy compensation, but the enthalpic contribution increases more significantly leading to a net stabilizing effect that slows down unfolding of the intermediate. Regarding the molecular mechanism of stabilization, trehalose increases the compactness of the unfolded state. The conformational space accessible to the unfolded state decreases in the presence of trehalose when the unfolded state acquires residual native interactions that channel the folding of the protein. This residual structure results into less hydrophobic groups being newly exposed upon unfolding, as less water molecules are immobilized upon unfolding.     

Comparison of thermal properties of bovine spectrin and fodrin

Thermal properties of bovine brain fodrin have been studied by circular dichroism and electron spin resonance and compared to those of bovine erythrocyte spectrin. Protein unfolding was induced either by urea or by a combination of heat and urea. The denaturation profiles of the two proteins, as measured by the changes in ellipticity at 222 nm as a function of temperature, were very similar but fodrin denaturation occurred at both higher temperatures and higher urea concentrations. In the absence of urea the thermal transition of spectrin was centered at 51 degrees C and that of fodrin at 54.5 degrees C. Proteins were also labeled with a maleimide analog spin probe. Spin-labeled fodrin showed a thermal transition similar to that of spectrin but centered at 46 degrees C instead of 42 degrees C. These findings indicated a close structural similarity of the two proteins but a slightly higher conformational stability of fodrin to both heat and urea.     

Refinement of noncalorimetric determination of the change in heat capacity, DeltaC(p), of protein unfolding and validation across a wide temperature range

The change in heat capacity, DeltaC(p), on protein unfolding has been usually determined by calorimetry. A noncalorimetric method which employs the Gibbs-Helmholtz relationship to determine DeltaC(p) has seen some use. Generally, in this method the free energy change on unfolding of the protein is determined at a variety of temperatures and the temperature at which DeltaG is zero, T(m), and change in enthalpy at T(m) are determined by thermal denaturation and DeltaC(p) is then calculated using the Gibbs-Helmholtz equation. We show here that an abbreviated method with stability determinations at just two temperatures gives values of DeltaC(p) consistent with values from free energy change on unfolding determination at a much wider range of temperatures. Further, even the free energy change on unfolding from a single solvent denaturation at the proper temperature, when coupled with the melting temperature, T(m), and the van't Hoff enthalpy, DeltaH(vH), from a thermal denaturation, gives a reasonable estimate of DeltaC(p), albeit with greater uncertainty than solvent denaturations at two temperatures. We also find that nonlinear regression of the Gibbs-Helmholtz equation as a function of stability and temperature while simultaneously fitting DeltaC(p), T(m), and DeltaH(vH) gives values for the last two parameters that are in excellent agreement with experimental values.     

Investigation of charged polymer influence on green fluorescent protein thermal stability

Methods of stabilization and formulation of proteins are important in both biopharmaceutical and biocatalysis industries. Polymers are often used as modifiers of characteristics of biological macromolecules to improve the biochemical activity and stability of proteins or drug bioavailability. Green fluorescent protein (GFP) shows remarkable structural stability and high fluorescence; its stability can be directly related to its fluorescence output, among other characteristics. GFP is stable under increasing temperatures, and its thermal denaturation is highly reproducible. Relative thermal stability was undertaken by incubation of GFP at varying temperatures and GFP fluorescence was used as a reporter for unfolding. At 80°C, DEAE-dextran did not have any effect on GFP fluorescence, indicating that it does not confer stability.     

Thermal stability of proteins analyzed by temperature scanning pH-stat titration

Based on the fact that pH changes occur during the thermal unfolding of a protein, a pH-stat titrimetric procedure is described for the analysis of thermal stability. In all cases the agreement with other stability measurements was good, including a correlation with activity loss in enzymes. A model for the titration curves, assuming first-order denaturation kinetics, linear temperature increase, and validity of the Arrhenius equation, has been proposed and analyzed. Thus, thermodynamic constants can be calculated from tritration curves, or transition temperatures estimated if the Arrhenius constants are known. The equipment consists of a pH-stat, a programmable heating unit, and a temperature measuring/recording system. Analysis can be done quickly and on partially purified solutions, provided the buffer capacity is low, using about 20 mg protein/10 ml sample. The effects of pH, Ca²⁺ ions, substrate, chemical modification, etc., on thermal stability are conveniently analyzed up to about 90°C.     

Characterization of the stabilizing effect of point mutations of pyruvate oxidase from Lactobacillus plantarum: protection of the native state by modulating coenzyme binding and subunit interaction.

Point mutations in the gene of pyruvate oxidase from Lactobacillus plantarum, with proline residue 178 changed to serine, serine 188 to asparagine, and alanine 458 to valine, as well as a combination of the three single point mutations, lead to a significant functional stabilization of the protein. The enzyme is a tetrameric flavoprotein with tightly bound cofactors, FAD, TPP, and divalent metal ions. Thus, stabilization may be achieved either at the level of tertiary or quaternary interactions, or by enhanced cofactor binding. In order to discriminate between these alternatives, unfolding, dissociation, and cofactor binding of the mutant proteins were analyzed. The point mutations do not affect the secondary and tertiary structure, as determined by circular dichroism and protein fluorescence. Similarly, the amino acid substitutions neither modulate the enzymatic properties of the mutant proteins nor do they stabilize the structural stability of the apoenzymes. This holds true for both the local and the global structure with unfolding transitions around 2.5 M and 5 M urea, respectively. On the other hand, deactivation of the holoenzyme (by urea or temperature) is significantly decreased. The most important stabilizing effect is caused by the Ala-Val exchange in the C-terminal domain of the molecule. Its contribution is close to the value observed for the triple mutant, which exhibits maximum stability, with a shift in the thermal transition of ca. 10 degrees C. The effects of the point mutations on FAD binding and subunit association are interconnected. Because FAD binding is linked to oligomerization, the stability of the mutant apoenzyme-FAD complexes is increased. Accordingly, mutants with maximum apparent FAD binding exhibit maximum stability. Analysis of the quaternary structure of the mutant enzymes in the absence and in the presence of coenzymes gives clear evidence that both improved ligand binding and subunit interactions contribute to the observed thermal stabilization.     

On the thermal unfolding character of globular proteins

A theoretical model is presented to study the stepwise thermal unfolding of globular proteins using the stabilizing/destabilizing characters of amino acid residues in protein crystals. A multiple regression relation connecting the melting temperature and the amounts of stabilizing and destabilizing groups of residues in a protein, when used for the thermal behavior of peptide segments, provides reliable results on the stepwise unfolding nature of the protein. In ribonuclease A, the shell residues 16–22 are predicted to unfold earlier in the temperature range 30–45°C; the -sheet structures undergo thermal denaturation as a single cooperative unit and there is evidence indicating the segment 106–118 as a nucleation site. In ribonuclease S, the S-peptide unfolds earlier than S-protein. The predicted average and the range of melting temperatures, and the folding pathways of a set of globular proteins, agree very well with the experimental results. The results obtained in the present study indicate that (i) most of the nucleation parts possess high relative thermal stability, (ii) the unfolded state retains some residual structure, and (iii) some segments undergo gradual and overlapping thermal denaturation.     

Quantitation of protein-protein interactions by thermal stability shift analysis

Thermal stability shift analysis is a powerful method for examining binding interactions in proteins. We demonstrate that under certain circumstances, protein-protein interactions can be quantitated by monitoring shifts in thermal stability using thermodynamic models and data analysis methods presented in this work. This method relies on the determination of protein stabilities from thermal unfolding experiments using fluorescent dyes such as SYPRO Orange that report on protein denaturation. Data collection is rapid and straightforward using readily available real-time polymerase chain reaction instrumentation. We present an approach for the analysis of the unfolding transitions corresponding to each partner to extract the affinity of the interaction between the proteins. This method does not require the construction of a titration series that brackets the dissociation constant. In thermal shift experiments, protein stability data are obtained at different temperatures according to the affinity- and concentration-dependent shifts in unfolding transition midpoints. Treatment of the temperature dependence of affinity is, therefore, intrinsic to this method and is developed in this study. We used the interaction between maltose-binding protein (MBP) and a thermostable synthetic ankyrin repeat protein (Off7) as an experimental test case because their unfolding transitions overlap minimally. We found that MBP is significantly stabilized by Off7. High experimental throughput is enabled by sample parallelization, and the ability to extract quantitative binding information at a single partner concentration. In a single experiment, we were able to quantify the affinities of a series of alanine mutants, covering a wide range of affinities (～ 100 nM to ～ 100 μM).     

Thermodynamics elucidation of the structural stability of a thermophilic protein

The structural stability of the protein, phycocyanin isolated from two strains of cyanophyta, Synechococcus lividus (thermophile) and Phormidium luridum (mesophile), are investigated by comparative thermal and denaturant unfolding, using differential scanning calorimetry, visible absorption spectrophotometry, and circular dichroism. The thermophilic protein exhibits a much higher temperature and enthalpy of unfolding from the native to the denatured state. The concentration of urea at half-completion of thermal unfolding is essentially the same between the thermophilic and mesophilic proteins; in contrast, the corresponding temperature and the enthalpy of thermal unfolding are much higher for the thermophilic protein. In addition, the concentration of urea at which the non-thermal (denaturant) unfolding of protein is half-completed, as detected by either circular dichroism or absorption spectroscopy, is significantly higher in the thermophilic protein, while the apparent free energy of unfolding only shows a moderate difference between the two proteins. The distinct differences in the enthalpy of thermal unfolding and the free energy of denaturant unfolding are interpreted in terms of a significant entropy change associated with the unfolding of these proteins. This entropy contribution is much higher in the thermophilic protein, and may be derived from its more rigid overall structure that possesses higher internal hydrophobicity and stronger internal packing.     

Some thermodynamic implications for the thermostability of proteins

An analysis of the thermodynamics of protein stability reveals a general tendency for proteins that denature at higher temperatures to have greater free energies of maximal stability. To a reasonable approximation, the temperature of maximal stability for the set of globular, water-soluble proteins surveyed by Robertson and Murphy occurs at T* approximately 283K, independent of the heat denaturation temperature, T(m). This observation indicates, at least for these proteins, that thermostability tends to be achieved through elevation of the stability curve rather than by broadening or through a horizontal shift to higher temperatures. The relationship between the free energy of maximal stability and the temperature of heat denaturation is such that an increase in maximal stability of approximately 0.008 kJ/mole/residue is, on average, associated with a 1 degrees C increase in T(m). An estimate of the energetic consequences of thermal expansion suggests that these effects may contribute significantly to the destabilization of the native state of proteins with increasing temperature.     

Conformational stability of HPr: the histidine-containing phosphocarrier protein from Bacillus subtilis.

The conformational stability of the histidine-containing phosphocarrier protein (HPr) from Bacillus subtilis has been determined using a combination of thermal unfolding and solvent denaturation experiments. The urea-induced denaturation of HPr was monitored spectroscopically at fixed temperatures and thermal unfolding was performed in the presence of fixed concentrations of urea. These data were analyzed in several different ways to afford a measure of the cardinal parameters (delta Hg, Tg, delta Sg, and delta Cp) that describe the thermodynamics of folding for HPr. The method of Pace and Laurents (Pace CN, Laurents DV, 1989, Biochemistry 28:2520-2525) was used to estimate delta Cp as was a global analysis of the thermal- and urea-induced unfolding data. Each method used to analyze the data gives a similar value for delta Cp (1,170 +/- 50 cal mol-1K-1). Despite the high melting temperature for HPr (Tg = 73.5 degrees C), the maximum stability of the protein, which occurs at 26 degrees C, is quite modest (delta Gs = 4.2 kcal mol-1). In the presence of moderate concentrations of urea, HPr exhibits cold denaturation, and thus a complete stability curve for HPr, including a measure of delta Cp, can be achieved using the method of Chen and Schellman (Chen B, Schellman JA, 1989, Biochemistry 28:685-691). A comparison of the different methods for the analysis of solvent denaturation curves is provided and the effects of urea on the thermal stability of this small globular protein are discussed. The methods presented will be of general utility in the characterization of the stability curve for many small proteins.     

Relationship of protein flexibility to thermostability

Thermostability of proteins arises from the simultaneous effect of several forces, which in fact lead to decreased flexibility of the polypeptide chain. This is verified by flexibility indices, which are derived from normalized B-values of individual amino acids in several refined three-dimensional structures. Flexibility indices show that overall flexibility is reduced when thermostability is increased. Protein molecules require both flexibility and rigidity to function, but the higher the temperature optimum and stability the more rigid is the structure needed to compensate for increased thermal fluctuations. Flexibilities of proteins performing the same catalytic activity seem to be about the same at their temperature optima, but the more rigid thermostable proteins reach the flexibility of thermolabile proteins at higher temperatures. In several proteins such as allosteric enzymes, some local sites of flexibility are highly conserved. The relevance of reduced flexibility to overall stability of proteins is also discussed. Flexibility indices and profiles can be used in the design of more stable proteins by site-directed mutagenesis.