Thermal stability of proteins does not correlate with the energy of intramolecular interactions

Small monomeric proteins from mesophilic and thermophilic organisms were studied. They have close structural and physical and chemical properties but vary in thermal stability. A thermodynamic analysis of heat unfolding was made and integral enthalpy of unfolding (DeltaH(unf)), heat capacity of hydration (DeltaC(p)(hyd)) and enthalpy of hydration (DeltaH(hyd)) and of the buried surface area (DeltaASA) of nonpolar and polar groups as well as the enthalpy of disruption of intramolecular interaction (DeltaH(int) in gas phase) at 298 K were determined. The absence of correlation between protein thermostability and energetic components suggests that regulatory mechanism of protein thermal stabilization has entropic nature.     

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.     

Effects of hydrated water on protein unfolding

The conformational stability of a protein in aqueous solution is described in terms of the thermodynamic properties such as unfolding Gibbs free energy, which is the difference in the free energy (Gibbs function) between the native and random conformations in solution. The properties are composed of two contributions, one from enthalpy due to intramolecular interactions among constituent atoms and chain entropy of the backbone and side chains, and the other from the hydrated water around a protein molecule. The hydration free energy and enthalpy at a given temperature for a protein of known three-dimensional structure can be calculated from the accessible surface areas of constituent atoms according to a method developed recently. Since the hydration free energy and enthalpy for random conformations are computed from those for an extended conformation, the thermodynamic properties of unfolding are evaluated quantitatively. The evaluated hydration properties for proteins of known transition temperature (Tm) and unfolding enthalpy (delta Hm) show an approximately linear dependence on the number of constituent heavy atoms. Since the unfolding free energy is zero at Tm, the enthalpy originating from interatomic interactions of a polypeptide chain and the chain entropy are evaluated from an experimental value of delta Hm and computed properties due to the hydrated water around the molecule at Tm. The chain enthalpy and entropy thus estimated are largely compensated by the hydration enthalpy and entropy, respectively, making the unfolding free energy and enthalpy relatively small. The computed temperature dependences of the unfolding free energy and enthalpy for RNase A, T4 lysozyme, and myoglobin showed a good agreement with the experimental ones.(ABSTRACT TRUNCATED AT 250 WORDS)     

The heat capacity of proteins

The heat capacity plays a major role in the determination of the energetics of protein folding and molecular recognition. As such, a better understanding of this thermodynamic parameter and its structural origin will provide new insights for the development of better molecular design strategies. In this paper we have analyzed the absolute heat capacity of proteins in different conformations. The results of these studies indicate that three major terms account for the absolute heat capacity of a protein: (1) one term that depends only on the primary or covalent structure of a protein and contains contributions from vibrational frequencies arising from the stretching and bending modes of each valence bond and internal rotations; (2) a term that contains the contributions of noncovalent interactions arising from secondary and tertiary structure; and (3) a term that contains the contributions of hydration. For a typical globular protein in solution the bulk of the heat capacity at 25°C is given by the covalent structure term (close to 85% of the total). The hydration term contributes about 15 and 40% to the total heat capacity of the native and unfolded states, respectively. The contribution of non-covalent structure to the total heat capacity of the native state is positive but very small and does not amount to more than 3% at 25°C. The change in heat capacity upon unfolding is primarily given by the increase in the hydration term (about 95%) and to a much lesser extent by the loss of noncovalent interactions (up to ～5%). It is demonstrated that a single universal mathematical function can be used to represent the partial molar heat capacity of the native and unfolded states of proteins in solution. This function can be experimentally written in terms of the molecular weight, the polar and apolar solvent accessible surface areas, and the total area buried from the solvent. This unique function accurately predicts the different magnitude and temperature dependences of the heat capacity of both the native and unfolded states, and therefore of the heat capacity changes associated with folding/unfolding transitions. © 1995 Wiley-Liss, Inc.     

The reversible two-state unfolding of a monocot mannose-binding lectin from garlic bulbs reveals the dominant role of the dimeric interface in its stabilization

Allium sativum agglutinin (ASAI) is a heterodimeric mannose-specific bulb lectin possessing two polypeptide chains of molecular mass 11.5 and 12.5 kDa. The thermal unfolding of ASAI, characterized by differential scanning calorimetry and circular dichroism, shows it to be highly reversible and can be defined as a two-state process in which the folded dimer is converted directly to the unfolded monomers (A2 if 2U). Its conformational stability has been determined as a function of temperature, GdnCl concentration, and pH using a combination of thermal and isothermal GdnCl-induced unfolding monitored by DSC, far-UV CD, and fluorescence, respectively. Analyses of these data yielded the heat capacity change upon unfolding (DeltaC(p) and also the temperature dependence of the thermodynamic parameters, namely, DeltaG, DeltaH, and DeltaS. The fit of the stability curve to the modified Gibbs-Helmholtz equation provides an estimate of the thermodynamic parameters DeltaH(g), DeltaS(g), and DeltaC(p) as 174.1 kcal x mol(-1), 0.512 kcal x mol(-1) x K(-1), and 3.41 kcal x mol(-1) x K(-1), respectively, at T(g) = 339.4 K. Also, the free energy of unfolding, DeltaG(s), at its temperature of maximum stability (T(s) = 293 K) is 13.13 kcal x mol(-1). Unlike most oligomeric proteins studied so far, the lectin shows excellent agreement between the experimentally determined DeltaC(p) (3.2 +/- 0.28 kcal x mol(-1) x K(-1)) and those evaluated from a calculation of its accessible surface area. This in turn suggests that the protein attains a completely unfolded state irrespective of the method of denaturation. The absence of any folding intermediates suggests the quaternary interactions to be the major contributor to the conformational stability of the protein, which correlates well with its X-ray structure. The small DeltaC(p) for the unfolding of ASAI reflects a relatively small, buried hydrophobic core in the folded dimeric protein.     

Elimination of the C-cap in ubiquitin - structure, dynamics and thermodynamic consequences

Single amino acid substitutions rarely produce substantial changes in protein structure. Here we show that substitution of the C-cap residue in the alpha-helix of ubiquitin with proline (34P variant) leads to dramatic structural changes. The resulting conformational perturbation extends over the last two turns of the alpha-helix and leads to enhanced flexibility for residues 27-37. Thermodynamic analysis of this ubiquitin variant using differential scanning calorimetry reveals that the thermal unfolding transition remains highly cooperative, exhibiting two-state behavior. Similarities with the wild type in the thermodynamic parameters (heat capacity change upon unfolding and m-value) of unfolding monitored by DSC and chemical denaturation suggests that the 34P variant has comparable buried surface area. The hydrophobic core of 34P variant is not packed as well as that of the wild type protein as manifested by a lower enthalpy of unfolding. The increased mobility of the polypeptide chain of this ubiquitin variant allows the transient opening of the hydrophobic core as evidenced by ANS binding. Taken together, these results suggest exceptional robustness of cooperativity in protein structures.     

Amino acid substitutions affecting protein dynamics in eglin C do not affect heat capacity change upon unfolding

The heat capacity change upon unfolding (deltaC(p)) is a thermodynamic parameter that defines the temperature dependence of the thermodynamic stability of proteins; however, physical basis of the heat capacity change is not completely understood. Although empirical surface area-based calculations can predict heat capacity changes reasonably well, accumulating evidence suggests that changes in hydration of those surfaces is not the only parameter contributing to the observed heat capacity changes upon unfolding. Because packing density in the protein interior is similar to that observed in organic crystals, we hypothesized that changes in protein dynamics resulting in increased rigidity of the protein structure might contribute to the observed heat capacity change upon unfolding. Using differential scanning calorimetry we characterized the thermodynamic behavior of a serine protease inhibitor eglin C and two eglin C variants with altered native state dynamics, as determined by NMR. We found no evidence of changes in deltaC(p) in either of the variants, suggesting that changes in rigidity do not contribute to the heat capacity change upon unfolding in this model system.     

Thermodynamics of unfolding of an integral membrane protein in mixed micelles

Quantitative studies of membrane protein folding and unfolding can be difficult because of difficulties with efficient refolding as well as a pronounced propensity to aggregate. However, mixed micelles, consisting of the anionic detergent sodium dodecyl sulfate and the nonionic detergent dodecyl maltoside facilitate reversible and quantitative unfolding and refolding. The 4-transmembrane helix protein DsbB from the inner membrane of Escherichia coli unfolds in mixed micelles according to a three-state mechanism involving an unfolding intermediate I. The temperature dependence of the kinetics of this reaction between 15 degrees and 45 degrees C supports that unfolding from I to the denatured state D is accompanied by a significant decrease in heat capacity. For water-soluble proteins, the heat capacity increases upon unfolding, and this is generally interpreted as the increased binding of water to the protein as it unfolds, exposing more surface area. The decrease in DsbB's heat capacity upon unfolding is confirmed by independent thermal scans. The decrease in heat capacity is not an artifact of the use of mixed micelles, since the water soluble protein S6 shows conventional heat-capacity changes in detergent. We speculate that it reflects the binding of SDS to parts of DsbB that are solvent-exposed in the native DM-bound state. This implies that the periplasmic loops of DsbB are relatively unstructured. This anomalous thermodynamic behavior has not been observed for beta-barrel membrane proteins, probably because they do not bind SDS so extensively. Thus the thermodynamic behavior of membrane proteins appears to be intimately connected to their detergent-binding properties.     

Oxidation of methionine residues in recombinant human interleukin-1 receptor antagonist: implications of conformational stability on protein oxidation kinetics

Oxidation of methionine residues is involved in several biochemical processes and in degradation of therapeutic proteins. The relationship between conformational stability and methionine oxidation in recombinant human interleukin-1 receptor antagonist (rhIL-1ra) was investigated to document how thermodynamics of unfolding affect methionine oxidation in proteins. Conformational stability of rhIL-1ra was monitored by equilibrium urea denaturation, and thermodynamic parameters of unfolding (DeltaGH2O, m, and Cm) were estimated at different temperatures. Methionine oxidation induced by hydrogen peroxide at varying temperatures was monitored during "coincubation" of rhIL-1ra with peptides mimicking specific regions of the reactive methionine residues in the protein. The coincubation study allowed estimation of oxidation rates in protein and peptide at each temperature from which normalized oxidation rate constants and activation energies were calculated. The rate constants for buried Met-11 in the protein were lower than for methionine in the peptide with an associated increase in activation energy. The rate constants and activation energy of solvent exposed methionines in protein and peptide were similar. The results showed that conformational stability, monitored using the Cm value, has an effect on oxidation rates of buried methionines. The rate constant of buried Met-11 correlated well with the Cm value but not DeltaGH2O. No correlation was observed for the oxidation rates of solvent-exposed methionines with any thermodynamic parameters of unfolding. The findings presented have implications in protein engineering, in design of accelerated stability studies for protein formulation development, and in understanding disease conditions involving protein oxidation.     

Probing the contribution of internal cavities to the volume change of protein unfolding under pressure.

The structural origin of the decrease in system volume upon protein denaturation by pressure has remained a puzzle for decades. This negative volume change upon unfolding is assumed to arise globally from more intimate interactions between the polypeptide chain and water, including electrostriction of buried charges that become exposed upon unfolding, hydration of the polypeptide backbone and amino acid side chains and elimination of packing defects and internal void volumes upon unfolding of the chain. However, the relative signs and magnitudes of each of these contributing factors have not been experimentally determined. Our laboratory has probed the fundamental basis for the volume change upon unfolding of staphylococcal nuclease (Snase) using variable solution conditions and point mutants of Snase (Royer CA et al., 1993, Biochemistry 32:5222-5232; Frye KJ et al., 1996, Biochemistry 35:10234-10239). Our prior results indicate that for Snase, neither electrostriction nor polar or nonpolar hydration contributes significantly to the value of the volume change of unfolding. In the present work, we investigate the pressure induced unfolding of three point mutants of Snase in which internal cavity size is altered. The experimentally determined volume changes of unfolding for the mutants suggest that loss of internal void volume upon unfolding represents the major contributing factor to the value of the volume change of Snase unfolding.     

Hydration state change of proteins upon unfolding in sugar solutions

Change in hydration number of proteins upon unfolding, Deltan, was obtained from the analysis of thermal unfolding behavior of proteins in various sugar solutions with water activity, a(W), varied. By applying the reciprocal form of Wyman-Tanford equation, Deltan was determined to be 133.9, 124.1, and 139.2 per protein molecule for ribonuclease A at pH=5.5, 4.2, and 2.8, respectively, 201.4 for lysozyme at pH=5.5, and 100.1 for alpha-chymotripnogen A at pH=2.0. Among the sugars tested, reducing sugars gave the lower apparent Deltan as compared with nonreducing sugars probably because of the direct interaction of reducing terminal with amino group of proteins at a high temperature. From the knowledge of Deltan, a new thermodynamic model for protein stability was proposed with explicit consideration for hydration state change of protein upon unfolding. From this model, the contribution of a(W) was proven to be always positive for stabilization of proteins and its effect is not negligible depending on Deltan and a(W).     

Applications of pressure perturbation calorimetry to study factors contributing to the volume changes upon protein unfolding

BackgroundPressure perturbation calorimetry (PPC) is a biophysical method that allows direct determination of the volume changes upon conformational transitions in macromolecules.Scope of this reviewThis review provides novel details of the use of PPC to analyze unfolding transitions in proteins. The emphasis is made on the data analysis as well as on the validation of different structural factors that define the volume changes upon unfolding. Four case studies are presented that show the application of these concepts to various protein systems.Major conclusionsThe major conclusions are:

• 1.Knowledge of the thermodynamic parameters for heat induced unfolding facilitates the analysis of the PPC profiles.
• 2.The changes in the thermal expansion coefficient upon unfolding appear to be temperature dependent.
• 3.Substitutions on the protein surface have negligible effects on the volume changes upon protein unfolding.
• 4.Structural plasticity of proteins defines the position dependent effect of amino acid substitutions of the residues buried in the native state.
• 5.Small proteins have positive volume changes upon unfolding which suggests difference in balance between the cavity/void volume in the native state and the hydration volume changes upon unfolding as compared to the large proteins that have negative volume changes.

General significanceThe information provided here gives a better understanding and deeper insight into the role played by various factors in defining the volume changes upon protein unfolding. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

A molecular dynamics simulation of SNase and its hydration shell at high temperature and high pressure

Temperature- and pressure-induced unfolding of staphylococcal nuclease (SNase) was studied by Royer, Winter et al. using a variety of experimental techniques (SAXS, FT-IR and fluorescence spectroscopy, DSC, PPC, densimetry). For a more detailed understanding of the underlying mechanistic processes of the different unfolding scenarios, we have carried out a series of molecular dynamics (MD) computer simulations on SNase. We investigated the initial changes of the structure of the protein upon application of pressure (up to 5 kbar) and discuss volumetric and structural differences between the native and pressure pre-denatured state. Additionally, we have obtained the compressibility of the protein and hydration water and compare these data with experimental results. As water plays a crucial role in determining the structure, dynamics and function of proteins, we undertook a detailed analysis of the structure of the interfacial water and the protein-solvent H-bond network as well. Moreover, we report here also MD results on the temperature-induced unfolding of SNase. The time evolution of the protein volume and solvent accessible surface area during thermal unfolding have been investigated, and we present a detailed discussion of the temperature-induced unfolding pathway of SNase in terms of secondary and tertiary structural changes.     

Comparison of the structures and stabilities of coiled‐coil proteins containing hexafluoroleucine and t‐butylalanine provides insight into the stabilizing effects of highly fluorinated amino acid side‐chains

Highly fluorinated analogs of hydrophobic amino acids are well known to increase the stability of proteins toward thermal unfolding and chemical denaturation, but there is very little data on the structural consequences of fluorination. We have determined the structures and folding energies of three variants of a de novo designed 4‐helix bundle protein whose hydrophobic cores contain either hexafluoroleucine (hFLeu) or t‐butylalanine (tBAla). Although the buried hydrophobic surface area is the same for all three proteins, the incorporation of tBAla causes a rearrangement of the core packing, resulting in the formation of a destabilizing hydrophobic cavity at the center of the protein. In contrast, incorporation of hFLeu, causes no changes in core packing with respect to the structure of the nonfluorinated parent protein which contains only leucine in the core. These results support the idea that fluorinated residues are especially effective at stabilizing proteins because they closely mimic the shape of the natural residues they replace while increasing buried hydrophobic surface area.     

Role of Protein Cavities on Unfolding Volume Change and on Internal Dynamics under Pressure

The effects of two single point cavity forming mutations, F110S and I7S, on the unfolding volume change (DeltaV(0)) of azurin from Pseudomonas aeruginosa and on the internal dynamics of the protein fold under pressure were probed by the fluorescence and phosphorescence emission of Trp-48, deeply buried in the compact hydrophobic core of the macromolecule. Pressure-induced unfolding, monitored by the shift of the center of mass of the fluorescence spectrum, showed that DeltaV(0) is in the range of 60-70 mL/mol, not significantly different between cavity mutants and compact azurin species such as the wild-type and the mutant C3A/C26A, in which the superficial disulphide has been removed. The lack of extra volume in F110S and I7S proves that the engineered cavities, 40 A(3) in I7S and 100 A(3) in F110S, are filled with water molecules. Changes in flexibility of the protein matrix around the chromophore were monitored by the intrinsic phosphorescence lifetime (tau(0)). The application of pressure in the predenaturation range initially decreases the internal flexibility of azurin, the trend eventually reverting on approaching unfolding. The main difference between compact folds, wild-type and C3A/C26A, and cavity mutants is that the inversion point is powered from approximately 3 kbar to 1.5 kbar for F110S and <0.1 kbar for I7S, meaning that in the latter species pressure-induced internal hydration dominates very early over any compaction of the globular fold resulting from the reduction of internal free volume. The similar response between wild-type and the significantly less-stable C3A/C26A mutant suggests that thermodynamic stability per se is not the dominant factor regulating pressure-induced internal hydration of proteins.     

Characteristic thermodynamic properties of hydrated water for 20 amino acid residues in globular proteins

Thermodynamic properties associated with hydrated water of proteins of known three-dimensional structure were computed and average values of hydration free energy, enthalpy, and heat capacity of unfolding for every amino acid residue were obtained. Each amino acid residue had characteristic values; in particular, the quantities for a side chain reflected the character of the amino acid, while those for the main chain were more or less the same except for glycine, alanine, and proline. The major contribution to the quantities was from the end group(s) of a side chain. The following interesting features were found. 1) The hydration quantity of unfolding derived from the native and extended conformations for a protein was approximately equal to the sum of the corresponding average quantities of component amino acid residues in the protein. 2) The profile of a quantity such as hydration free energy of unfolding along the sequence computed from the accessible surface areas of the native and extended conformations showed a strong correlation with the profile obtained by allocating the average value for the amino acid residue at every position on the sequence. The correlation coefficients between two profiles for unfolding quantities of hydration, i.e., free energy, enthalpy, heat capacity, and free energy of side chain are 0.72, 0.62, 0.80, and 0.75, respectively. Thus, every amino acid residue in the native conformation of a globular protein seems to be located in such a position that a thermodynamic quantity for each residue is approximately equal to its average value.     

Structural reorganization triggered by charging of Lys residues in the hydrophobic interior of a protein

Structural consequences of ionization of residues buried in the hydrophobic interior of proteins were examined systematically in 25 proteins with internal Lys residues. Crystal structures showed that the ionizable groups are buried. NMR spectroscopy showed that in 2 of 25 cases studied, the ionization of an internal Lys unfolded the protein globally. In five cases, the internal charge triggered localized changes in structure and dynamics, and in three cases, it promoted partial or local unfolding. Remarkably, in 15 proteins, the ionization of the internal Lys had no detectable structural consequences. Highly stable proteins appear to be inherently capable of withstanding the presence of charge in their hydrophobic interior, without the need for specialized structural adaptations. The extent of structural reorganization paralleled loosely with global thermodynamic stability, suggesting that structure-based pK(a) calculations for buried residues could be improved by calculation of thermodynamic stability and by enhanced conformational sampling.     

Thermodynamics of unfolding of the alpha-amylase inhibitor tendamistat. Correlations between accessible surface area and heat capacity.

Unfolding of the small alpha-amylase inhibitor tendamistat (74 residues, 2 disulfide bridges) has been characterized thermodynamically by high sensitivity scanning microcalorimetry. To link the stability parameters with structural information we use heat capacity group parameters and water accessible surface areas to calculate the change in heat capacity on unfolding of tendamistat. Our results show that both the group parameter and surface area approaches provide a reasonable, though not perfect, basis for delta Cp calculations. When using the experimentally determined temperature-independent heat capacity increase of 2.89 kJ mol-1 K-1 tendamistat exhibits convergence of thermodynamic parameters at about 140 degrees C, in agreement with recent predictions of the temperature at which the hydrophobic hydration is supposed to disappear. Despite the apparent support of this new view of the hydrophobic effect, there are inconsistencies in the interpretation of the thermodynamic parameters and these are addressed in the Discussion. The specific stability of tendamistat is similar to that of modified bovine pancreatic trypsin inhibitor, with only two of the native three disulfide bridges intact. This observation confirms our previous conclusion that disulfide bridges affect significantly the enthalpy and entropy of unfolding. The recent study by Doig & Williams provides additional convincing support for this conclusion. The predictive scheme proposed by these authors permits a fair estimate of the Gibbs free energy and enthalpy changes of these two proteins.