Differential scanning calorimetry of the unfolding of myosin subfragment 1, subfragment 2, and heavy meromyosin

The thermal unfolding of rabbit skeletal heavy meromyosin (HMM), myosin subfragment 1, and subfragment 2 has been studied by differential scanning calorimetry (DSC). Two distinct endotherms are observed in the DSC scan of heavy meromyosin. The first endotherm, with a Tm of 41 degrees C at pH 7.9 in 0.1 M KCl, is assigned to the unfolding of the subfragment 2 domain of HMM based on scans of isolated subfragment 2. The unfolding of the subfragment 2 domain is reversible both in the isolated form and in HMM. The unfolding of subfragment 2 in HMM can be fit as a single two-state transition with a delta Hvh and delta Hcal of 161 kcal/mol, indicating that subfragment 2 exists as a single domain in HMM. The unfolding of subfragment 2 is characterized by an extraordinarily large delta Cp of approximately 30,000 cal/(deg.mol). In the presence of nucleotides, the high-temperature HMM endotherm with a Tm of 48 degrees C shifts to higher temperature, indicating that this peak corresponds to the unfolding of the subfragment 1 domain. This assignment has been confirmed by comparison with isolated subfragment 1. The stabilizing effect of AMPPNP was significantly greater than that of ADP. The vanadate-trapped ADP species was slightly more stable than M.AMPPNP with a Tm at 58 degrees C. The unfolding of subfragment 1, both in the isolated form and in HMM, was irreversible. Only a single endotherm was noted in the DSC scans of the subfragment 1 domain of HMM and in freshly prepared subfragment 1 complexes.(ABSTRACT TRUNCATED AT 250 WORDS)     

The native state conformational ensemble of the SH3 domain from alpha-spectrin.

The folding/unfolding equilibrium of the alpha-spectrin SH3 domain has been measured by NMR-detected hydrogen/deuterium exchange and by differential scanning calorimetry. Protection factors against exchange have been obtained under native conditions for more than half of the residues in the domain. Most protected residues are located at the beta-strands, the short 3(10) helix, and part of the long RT loop, whereas the loops connecting secondary structure elements show no measurable protection. Apparent stability constants per residue and their corresponding Gibbs energies have been calculated from the exchange experiments. The most stable region of the SH3 domain is defined by the central portions of the beta-strands. The peptide binding region, on the other hand, is composed of a highly stable region (residues 53-57) and a highly unstable region, the loop between residues 34-41 (n-Src loop). All residues in the domain have apparent Gibbs energies lower than the global unfolding Gibbs energy measured by differential scanning calorimetry, indicating that under our experimental conditions the amide exchange of all residues in the SH3 domain occurs primarily via local unfolding reactions. A structure-based thermodynamic analysis has allowed us to predict correctly the thermodynamics of the global unfolding of the domain and to define the ensemble of conformational states that quantitatively accounts for the observed pattern of hydrogen exchange protection. These results demonstrate that under native conditions the SH3 domain needs to be considered as an ensemble of conformations and that the hydrogen exchange data obtained under those conditions cannot be interpreted by a two-state equilibrium. The observation that specific regions of a protein are able to undergo independent local folding/unfolding reactions indicates that under native conditions the scale of cooperative interactions is regional rather than global.     

A test of the linear extrapolation of unfolding free energy changes over an extended denaturant concentration range.

Guanidine hydrochloride (GdnHCl) and thermally induced unfolding measurements on the oxidized form of Escherichia coli thioredoxin at pH 7 were combined for the purpose of assessing the functional dependence of unfolding free energy changes on denaturant concentration over an extended GdnHCl concentration range. Conventional analysis of GdnHCl unfolding exhibits a linear plot of unfolding delta G vs [GdnHCl] in the transition zone. In order to extend unfolding delta G measurements outside of that narrow concentration range, thermal unfolding measurements were performed using differential scanning calorimetry (DSC) in the presence of low to moderate concentrations of GdnHCl. The unfolding delta G values from the DSC measurements were corrected to 25 degrees C using the Gibbs-Helmholtz equation and mapped onto the delta G vs [GdnHCl] plot. The dependence of unfolding delta G on [GdnHCl] was found to be linear over the full denaturant concentration range, provided that the chloride ion concentration was kept at a threshold of greater than or equal to 1.5 M. In the DSC experiments performed in the presence of GdnHCl, chloride concentrations were maintained at 1.5 M by addition of appropriate amounts of NaCl. The linear extrapolation method (LEM) gives an unfolding free energy change in the absence of denaturant (delta G degrees N-U) in excellent agreement with the delta G determined by DSC measurement in 1.5 M NaCl. The various methods give a consensus unfolding delta G value of 8.0 kcal/mol at 25 degrees C in the absence of denaturant.(ABSTRACT TRUNCATED AT 250 WORDS)     

Calorimetric study of tryptophan synthase alpha-subunit and two mutant proteins

Heat-denaturation of tryptophan synthase alpha-subunit from E. coli and two mutant proteins (Glu 49 leads to Gln or Ser; called Gln 49 or Ser 49, respectively) has been studied by the scanning microcalorimetric method at various pH, in an attempt to elucidate the role of individual amino acid residues in the conformational stability of a protein. The partial specific heat capacity in the native state at 20 degrees, Cp20, has been found to be (0.43 +/- 0.02) cal . k-1 . g-1, the unfolding heat capacity change, delta dCp, (0.10 +/- 0.01) cal . K-1 . g-1, and the unfolding enthalpy value extrapolated to 110 degrees, delta dh110, (9.3 +/- 0.5) cal . g-1 for the three proteins. The value of Cp20 was larger than those found for "fully compact protein" and that of delta dh110 was smaller. Unfolding Gibbs energy, delta dG at 25 degrees for Wild-type, Gln 49, and Ser 49 were 5.8, 8.4, and 7.1 kcal . mol-1 at pH 9.3, respectively. Unfolding enthalpy, delta dH, of the three proteins seemed to be the same and equal to (23.2 +/- 1.2) kcal . mol-1 at 25 degrees. As a consequence of the same value of delta dH and the different value in delta dG, substantial differences in unfolding entropy, delta dS, were found for the three proteins. The values of delta dG for the three proteins at 25 degrees coincided with those from equilibrium methods of denaturation by guanidine hydrochloride.     

A thermodynamic study of the 434-repressor N-terminal domain and of its covalently linked dimers.

The isolated N-terminal 1-69 domain of the 434-phage repressor, R69, and its covalently linked (head-to-tail and tail-to-tail) dimers have been studied by differential scanning microcalorimetry (DSC) and CD. At neutral solvent conditions the R69 domain maintains its native structure, both in isolated form and within the dimers. The stability of the domain depends highly upon pH within the acidic range, thus at pH 2 and low ionic strength R69 is already partially unfolded at room temperature. The thermodynamic parameters of unfolding calculated from the DSC data are typical for small globular proteins. At neutral pH and moderate ionic strength, the domains of the dimers behave as two independent units with unfolding parameters similar to those of the isolated domain, which means that linking two R69 domains, either by a long peptide linker or by a designed C-terminal disulfide bridge, does not induce any cooperation between them.     

Stability of recombinant Lys25-ribonuclease T1

The conformational stability of recombinant Lys25-ribonuclease T1 has been determined by differential scanning microcalorimetry (DSC), UV-monitored thermal denaturation measurements, and isothermal Gdn.HCl unfolding studies. Although rather different extrapolation procedures are involved in calculating the Gibbs free energy of stabilization, there is fair agreement between the delta G degrees values derived from the three different experimental techniques at pH 5, theta = 25 degrees C: DSC, 46.6 +/- 2.1 kJ/mol; UV melting curves, 48.7 +/- 5 kJ/mol; Gdn.HCl transition curves, 40.8 +/- 1.5 kJ/mol. Thermal unfolding of the enzyme is a reversible process, and the ratio of the van't Hoff and calorimetric enthalpy, delta HvH/delta Hcal, is 0.97 +/- 0.06. This result strongly suggests that the unfolding equilibrium of Lys25-ribonuclease T1 is adequately described by a simple two-state model. Upon unfolding the heat capacity increases by delta Cp degrees = 5.1 +/- 0.5 kJ/(mol.K). Similar values have been found for the unfolding of other small proteins. Surprisingly, this denaturational heat capacity change practically vanishes in the presence of moderate NaCl concentrations. The molecular origin of this effect is not clear; it is not observed to the same extent in the unfolding of bovine pancreatic ribonuclease A, which was employed in control experiments. NaCl stabilizes Lys25-ribonuclease T1. The transition temperature varies with NaCl activity in a manner that suggests two limiting binding equilibria to be operative. Below approximately 0.2 M NaCl activity unfolding is associated with dissociation of about one ion, whereas above that concentration about four ions are released in the unfolding reaction.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamic analysis of the folding of the streptococcal protein G IgG-binding domains B1 and B2: why small proteins tend to have high denaturation temperatures.

We have cloned, expressed, and characterized two naturally occurring variations of the IgG-binding domain of streptococcal protein G. The domain is a stable cooperative folding unit of 56 amino acids, which maintains a unique folded structure without disulfide cross-links or tight ligand binding. We have studied the thermodynamics of the unfolding reaction for the two versions of this domain, designated B1 and B2, which differ by six amino acids. They have denaturation temperatures of 87.5 degrees C and 79.4 degrees C, respectively at pH 5.4, as determined by differential scanning calorimetry. Thermodynamic state functions for the unfolding reaction (delta G, delta H, delta S, and delta Cp) have been determined and reveal several interesting insights into the behavior of very small proteins. First, though the B1 domain has a heat denaturation point close to 90 degrees C, it is not unusually stable at physiologically relevant temperatures (delta G = 25 kJ/mol at 37 degrees C). This behavior occurs because the stability profile (delta G vs temperature) is flat and shallow due to the small delta S and delta Cp for unfolding. Related to this point is the second observation that small changes in the free energy of unfolding of the B-domain due to mutation or change in solvent conditions lead to large shifts in the heat denaturation temperature. Third, the magnitude and relative contributions of hydrophobic vs nonhydrophobic forces (per amino acid residue) to the total free energy of folding of the B-domain are remarkably typical of other globular proteins of much larger size.     

Thermodynamics of unfolding for turkey ovomucoid third domain: thermal and chemical denaturation.

We have used thermal and chemical denaturation to characterize the thermodynamics of unfolding for turkey ovomucoid third domain (OMTKY3). Thermal denaturation was monitored spectroscopically at a number of wave-lengths and data were subjected to van't Hoff analysis; at pH 2.0, the midpoint of denaturation (Tm) occurs at 58.6 +/- 0.4 degrees C and the enthalpy of unfolding at this temperature (delta Hm) is 40.8 +/- 0.3 kcal/mol. When Tm was perturbed by varying pH and denaturant concentration, the resulting plots of delta Hm versus Tm yield a mean value of 590 +/- 120 cal/(mol.K) for the change in heat capacity upon unfolding (delta Cp). A global fit of the same data to an equation that includes the temperature dependence for the enthalpy of unfolding yielded a value of 640 +/- 110 cal/(mol.K). We also performed a variation of the linear extrapolation method described by Pace and Laurents, which is an independent method for determining delta Cp (Pace, C.N. & Laurents, D., 1989, Biochemistry 28, 2520-2525). First, OMTKY3 was thermally denatured in the presence of a variety of denaturant concentrations. Linear extrapolations were then made from isothermal slices through the transition region of the denaturation curves. When extrapolated free energies of unfolding (delta Gu) were plotted versus temperature, the resulting curve appeared linear; therefore, delta Cp could not be determined. However, the data for delta Gu versus denaturant concentration are linear over an extraordinarily wide range of concentrations. Moreover, extrapolated values of delta Gu in urea are identical to values measured directly.     

Stability of a homo-dimeric Ca(2+)-binding member of the beta gamma-crystallin superfamily: DSC measurements on spherulin 3a from Physarum polycephalum.

Spherulin 3a (S3a) from Physarum polycephalum represents the only known single-domain member of the superfamily of beta gamma eye-lens crystallins. It shares the typical two Greek-key motif and is stabilized by dimerization and Ca(2+)-binding. The temperature and denaturant-induced unfolding of S3a in the absence and in the presence of Ca2+ were investigated by differential scanning calorimetry and fluorescence spectroscopy. To accomplish reversibility without chemical modification of the protein during thermal denaturation, the only cysteine residue (Cys4) was substituted by serine; apart from that, the protein was destabilized by adding 0.5-1.8 M guanidinium chloride (GdmCl). The Cys4Ser mutant was found to be indistinguishable from natural S3a. The equilibrium unfolding transitions obey the two-state model according to N2-->2 U, allowing thermodynamic parameters to be determined by linear extrapolation to zero GdmCl concentration. The corresponding transition temperatures TM for the Ca(2+)-free and Ca(2+)-loaded protein were found to be 65 and 85 degrees C, the enthalpy changes delta Hcal, 800 and 1280 kJ/mol(dimer), respectively. The strong dependencies of TM and delta Hcal on the GdmCl concentration allow the molar heat capacity change delta Cp to be determined. As a result, delta Cp = 18 kJ/(K mol(dimer)) was calculated independent of Ca2+. No significant differences were obtained between the free energy delta G degree calculated from delta Hcal and TM, and extrapolated from the stability curves in the presence of different amounts of denaturant. The free energy derived from thermal unfolding was confirmed by the spectral results obtained from GdmCl-induced equilibrium transitions at different temperatures for the Ca(2+)-free or the Ca(2+)-loaded protein, respectively. Within the limits of error, the delta G degree values extrapolated from the transitions of chemical denaturation to zero denaturant concentration are identical with the calorimetric results.     

Denaturant m values and heat capacity changes: relation to changes in accessible surface areas of protein unfolding.

Denaturant m values, the dependence of the free energy of unfolding on denaturant concentration, have been collected for a large set of proteins. The m value correlates very strongly with the amount of protein surface exposed to solvent upon unfolding, with linear correlation coefficients of R = 0.84 for urea and R = 0.87 for guanidine hydrochloride. These correlations improve to R = 0.90 when the effect of disulfide bonds on the accessible area of the unfolded protein is included. A similar dependence on accessible surface area has been found previously for the heat capacity change (delta Cp), which is confirmed here for our set of proteins. Denaturant m values and heat capacity changes also correlate well with each other. For proteins that undergo a simple two-state unfolding mechanism, the amount of surface exposed to solvent upon unfolding is a main structural determinant for both m values and delta Cp.     

A comparative study of the unfolding thermodynamics of vertebrate metmyoglobins

Differential scanning microcalorimetry (DSC) of horse, rat, opossum, raccoon, carp, and armadillo metmyoglobins at alkaline pH gave data that fit the two-state unfolding model well. Monte Carlo studies were used to assess the impact of truncating DSC scans on the reliability of the calculated results when aggregation exotherms overlapped the unfolding endotherm at the high-temperature end of the scan. The DSC estimates for the conformational free energy at pH 8 and 298 K are compared to earlier results from isothermal acid and guanidinium chloride unfolding. Stability estimates at pH 8 for these six metmyoglobins obtained by DSC experiments do not agree with free energy estimates at pH 8 from guanidinium chloride unfolding. This is true for all three models used to extrapolate the free energy change to 0 M guanidinium chloride. Among these six myoglobins, significant variation appears in the temperature at which the myoglobin is half-unfolded, in the change in heat capacity upon unfolding, and in the change in enthalpy at 310 K. Calculations made with the hydrophobic model for protein folding [Baldwin, R.L. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 8069] suggest that a sizable variation exists for that portion of the unfolding enthalpy change assigned to forces other than the hydrophobic effect.     

Analysis of the two-state behavior of the thermal unfolding serum retinol binding protein containing a single retinol ligand.

Through the use of CD and DSC, the thermal unfolding of holo serum retinol binding protein containing a single, tightly bound retinol ligand was studied at pH 7.4. The DSC endotherm of the holoprotein ([retinol]/[protein] = 1) was asymmetric about the transition temperature of 78 degrees C. Using changes in ellipticity at 230 nm, the thermal unfolding curve was also asymmetric about the inflection point centered near 78 degrees C. van't Hoff enthalpies were determined by three means and compared to the calorimetric enthalpy (delta Hcal) of 200 kcal/mol. A van't Hoff enthalpy of 190 kcal/mol was determined from the dependence of transition temperature on the concentration of the ligand-bound protein. This value agreed well with the van't Hoff enthalpies found from fits of the DSC (delta HvH = 184 kcal/mol) and spectroscopic (delta HvH = 181 kcal/mol) curves to a two-state thermodynamic model that included ligand dissociation (NR in equilibrium with U+R, where NR is the native holoprotein, U is the unfolded apoprotein, and R is retinol). Poor agreement was obtained with a two-state model that ignored ligand dissociation (N in equilibrium with U). Furthermore, the NR in equilibrium with U+R model accounted for the asymmetry in both CD and DSC transitions and yielded a much improved fit of the data over the N in equilibrium with U model. From these considerations and simulations on other equilibrium models, it is suggested that the NR in equilibrium with U+R model is the simplest model that describes the thermal unfolding of this ligand-bound protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Influence of transition rates and scan rate on kinetic simulations of differential scanning calorimetry profiles of reversible and irreversible protein denaturation.

The thermodynamic parameters characterizing protein folding can be obtained directly using differential scanning calorimetry (DSC). They are meaningful only for reversible unfolding at equilibrium, which holds for small globular proteins; however, the unfolding or denaturation of most large, multidomain or multisubunit proteins is either partially or totally irreversible. The simplest kinetic model describing partially irreversible denaturation requires three states: Formula [see text] We obtain numerical solutions for N, U, and D as a function of temperature for this model and derive profiles of excess specific heat (Cp) in terms of the reduced variables v/ki and k1/k3, where v is the scan rate. The three-state model reduces to the two-state reversible or irreversible models for very large or very small values of k1/k3, respectively. The apparent transition temperature (Tapp) is always reduced by the irreversible step (U-->D). For all values of k3, Tapp is independent of v/k1 at sufficiently slow scan rates, even when denaturation is highly irreversible, but increases identically for all models at fast scan rates in which case the excess specific heat profile is determined by the rate of unfolding. Accurate values of delta H and delta S can be obtained for the reversible step only when k1 is more than 2000-50,000 times greater than k3. In principle, approximate values for the ratio k1/k3 can be obtained from plots of fraction unfolded vs fraction irreversibly denatured as a function of temperature; however, the fraction irreversibly denatured is difficult to measure accurately by DSC alone.(ABSTRACT TRUNCATED AT 250 WORDS)     

A new method for determining the heat capacity change for protein folding

In order to use results from calorimetry or thermal unfolding curves to estimate the free energy change for protein unfolding at 25 degrees C, it is necessary to know the change in heat capacity for unfolding, delta Cp. We describe a new method for measuring delta Cp which is based on results from urea and thermal unfolding curves but does not require a calorimeter. We find that delta Cp = 1650 +/- 200 cal/(deg.mol) for the unfolding of ribonuclease T1 and that delta Cp = 2200 +/- 300 cal/(deg.mol) for the unfolding of ribonuclease A.     

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.     

Basic pancreatic trypsin inhibitor has unusual thermodynamic stability parameters

The stability parameters delta Gst, delta Hst and delta Sst of native basic pancreatic trypsin inhibitor (BPTI) have been characterized by microcalorimetric unfolding studies in various buffer solutions, at different pH values and in the presence of guanidine hydrochloride. The unfolding enthalpy of BPTI, in contrast ot other globular proteins, exhibits a very small dependence on temperature, which results in a characteristic different temperature dependence of the Gibbs energy of stabilization. BPTI has a very high specific Gibbs energy of stabilization, which renders the slow exchange rates of amide protons understandable. Comparison of the unfolding entropy of BPTI at 110 degrees C with corresponding values of other proteins, revealed that the delta S values of BPTI are lower by 2.9 J/(K X residue). This lower value of the unfolding entropy is in good agreement with predictions of a theoretical study by Poland & Scheraga (1965) where the influence of crosslinks on the configurational entropy has been studied. Additionally, we were able to calculate an interaction enthalpy per site of -5.6 kJ/mol based on the measurements of unfolding of BPTI in 6 M-guanidine hydrochloride.     

Thermal denaturation of iso-1-cytochrome c variants: comparison with solvent denaturation.

Thermal denaturation studies as a function of pH were carried out on wild-type iso-1-cytochrome c and three variants of this protein at the solvent-exposed position 73 of the sequence. By examining the enthalpy and Tm at various pH values, the heat capacity increment (delta Cp), which is dominated by the degree of change in nonpolar hydration upon protein unfolding, was found for the wild type where lysine 73 is normally present and for three variants. For the Trp 73 variant, the delta Cp value (1.15 +/- 0.17 kcal/mol K) decreased slightly relative to wild-type iso-1-cytochrome c (1.40 +/- 0.06 kcal/mol K), while for the Ile 73 (1.65 +/- 0.07 kcal/mol K) and the Val 73 (1.50 +/- 0.06 kcal/mol K) variants, delta Cp increased slightly. In previous studies, the Trp 73, Ile 73, and Val 73 variants have been shown to have decreased m-values in guanidine hydrochloride denaturations relative to the wild-type protein (Hermann L, Bowler BE, Dong A, Caughey WS. 1995. The effects of hydrophilic to hydrophobic surface mutations on the denatured state of iso-1-cytochrome c: Investigation of aliphatic residues. Biochemistry 34:3040-3047). Both the m-value and delta Cp are related to the change in solvent exposure upon unfolding and other investigators have shown a correlation exists between these two parameters. However, for this subset of variants of iso-1-cytochrome c, a lack of correlation exists which implies that there may be basic differences between the guanidine hydrochloride and thermal denaturations of this protein. Spectroscopic data are consistent with different denatured states for thermal and guanidine hydrochloride unfolding. The different response of m-values and delta Cp for these variants will be discussed in this context.     

Effect of the central disulfide bond on the unfolding behavior of elongation factor Ts homodimer from Thermus thermophilus

Functionally active elongation factor Ts (EF-Ts) from Thermus thermophilus forms a homodimer. The dimerization interface of EF-Ts is composed of two antiparallel beta-sheets that can be connected by an intermolecular disulfide bond. The stability of EF-Ts from T. thermophilus in the presence and absence of the intermolecular disulfide bond was studied by differential scanning calorimetry and circular dichroism. The ratio of the van't Hoff and calorimetric enthalpies, delta H(vH)/delta H(cal), indicates that EF-Ts undergoes thermal unfolding as a dimer independently of the presence or absence of the disulfide bond. This can be concluded from (1) the presence of residual secondary structure above the thermal transition temperature, (2) the absence of concentration dependence, which would be expected for dissociation of the dimer prior to unfolding of the monomers, and (3) a relatively low heat capacity change (delta Cp) upon unfolding. The retained dimeric structure of the thermally denatured state allowed for the determination of the effect of the intermolecular disulfide bond on the conformational stability of EF-Ts, which is deltadelta G(S-S,SH HS) = 10.5 kJ/mol per monomer at 72.5 degrees C. The possible physiological implications of the dimeric EF-Ts structure and of the intersubunit disulfide bond for the extreme conformational stability of proteins in thermophiles are discussed.