Equilibrium unfolding of dimeric and engineered monomeric forms of lambda Cro (F58W) repressor and the effect of added salts: evidence for the formation of folded monomer induced by sodium perchlorate

The equilibrium unfolding transitions of Cro repressor variants, dimeric variant Cro F58W and monomer Cro K56[DGEVK]F58W, have been studied by urea and guanidine hydrochloride to probe the folding mechanism. The unfolding transitions of a dimeric variant are well described by a two state process involving native dimer and unfolded monomer with a free energy of unfolding, DeltaG(0,un)(0), of approximately 10-11 kcal/mol. The midpoint of transition curves is dependent on total protein concentration and DeltaG(0,un)(0) is independent of protein concentration, as expected for this model. Unfolding of Cro monomer is well described by the standard two state model. The stability of both forms of protein increases in the presence of salt but decreases with the decrease in pH. Because of the suggested importance of a N2<-->2F dimerization process in DNA binding, we have also studied the effect of sodium perchlorate, containing the chaotropic perchlorate anion, on the conformational transition of Cro dimer by CD, fluorescence and NMR (in addition to urea and guanidine hydrochloride) in an attempt both to characterize the thermodynamics of the process and to identify conditions that lead to an increase in the population of the folded monomers. Data suggest that sodium perchlorate stabilizes the protein at low concentration (<1.5 M) and destabilizes the protein at higher perchlorate concentration with the formation of a "significantly folded" monomer. The tryptophan residue in the "significantly folded" monomer induced by perchlorate is more exposed to the solvent than in native dimer.     

Isolation and in vitro characterization of temperature-sensitive mutants of the bacteriophage f1 gene V protein

In vivo selections were used to isolate 43 temperature-sensitive gene V mutants of the bacteriophage f1 from a collection of mutants constructed by saturation mutagenesis of the gene. The sites of temperature-sensitive substitutions are found in both the beta-sheets and the turns of the protein, and some sites are exposed to the solvent while others are not. Thirteen of the variant proteins were purified and characterized to evaluate their free energy changes upon unfolding and their affinities for single-stranded DNA, and eight were tested for their tendencies to aggregate at 42 degrees C. Each of the three temperature-sensitive mutants at buried sites and six of ten at surface sites had free energy changes of unfolding substantially lower (less stabilizing) than the wild-type at 25 degrees C. A seventh mutant at a surface site had a substantially altered unfolding transition and its free energy of unfolding was not estimated. The affinities of the mutant proteins for single-stranded DNA varied considerably, but two mutants at a surface site, Lys69, had much weaker binding to single-stranded DNA than any of the other mutants, while two mutants at another surface site, Glu30, had the highest DNA-binding affinities. The wild-type gene V protein is stable at 42 degrees C, but six of the eight mutants tested aggregated within a few minutes and the remaining two aggregated within 30 minutes at this temperature. Overall, each of the temperature-sensitive proteins tested had a tendency to aggregate at 42 degrees C, and most also had either a low free energy of unfolding (at 25 degrees C), or weak DNA binding. We suggest that any of these properties can lead to a temperature-sensitive gene V phenotype.     

Equilibrium and kinetic measurements of the conformational transition of thioredoxin in urea

Addition of urea to solutions of Escherichia coli thioredoxin results in a cooperative unfolding of the protein centered at 6.7 M urea at 25 degrees C and 5.1 M urea at 2 degrees C and neutral pH as judged by changes in tryptophan fluorescence emission, far-ultraviolet circular dichroism, and exclusion chromatography. Kinetic profiles of changes in tryptophan fluorescence emission intensity were analyzed following either manual or stopped-flow mixing to initiate unfolding or refolding. Unfolding of the native protein occurs in a single kinetic phase whose time constant is markedly dependent on urea concentration. Refolding of the urea-denatured protein occurs in a multiplicity of kinetic phases whose time constants and fractional amplitudes are also dependent upon urea concentration. Urea gradient gel electrophoretic and exclusion chromatographic measurements suggest the transient accumulation of at least one and likely two compact nativelike intermediate conformations during refolding. Simulations of both electrophoretic and chromatographic results suggest that the intermediate conformations are generated by the concerted action of the middle and fast refolding phases.     

Ruggedness in the Free Energy Landscape Dictates Misfolding of the Prion Protein

Experimental determination of the key features of the free energy landscapes of proteins, which dictate their adeptness to fold correctly, or propensity to misfold and aggregate and which are modulated upon a change from physiological to aggregation-prone conditions, is a difficult challenge. In this study, sub-millisecond kinetic measurements of the folding and unfolding of the mouse prion protein reveal how the free energy landscape becomes more complex upon a shift from physiological (pH 7) to aggregation-prone (pH 4) conditions. Folding and unfolding utilize the same single pathway at pH 7, but at pH 4, folding occurs on a pathway distinct from the unfolding pathway. Moreover, the kinetics of both folding and unfolding at pH 4 depend not only on the final conditions but also on the conditions under which the processes are initiated. Unfolding can be made to switch to occur on the folding pathway by varying the initial conditions. Folding and unfolding pathways appear to occupy different regions of the free energy landscape, which are separated by large free energy barriers that change with a change in the initial conditions. These barriers direct unfolding of the native protein to proceed via an aggregation-prone intermediate previously identified to initiate the misfolding of the mouse prion protein at low pH, thus identifying a plausible mechanism by which the ruggedness of the free energy landscape of a protein may modulate its aggregation propensity.     

Unfolding by temperature and guanidine hydrochloride of chicken pancreatic polypeptide

Thermal unfolding of chicken pancreatic polypeptide at two different concentrations was studied at various pH values. The thermal stability was higher at higher protein concentrations. The transition temperatures at two different protein concentrations changed with pH in parallel and decreased by about 30 degrees C on lowering pH from 5 to 2. The results on the thermal unfolding were analyzed by assuming that the dimerization constant is independent of pH, that the thermal unfolding occurs only after the pancreatic polypeptide dimers dissociated into the monomers, and that one ionizable group participates in the acid unfolding of the monomer. The free energy change for the unfolding of the pancreatic polypeptide monomer was estimated to be 1.4 kcal/mol. The unfolding of pancreatic polypeptide by guanidine hydrochloride at pH 6.0 and 25 degrees C was also studied. The stability to guanidine hydrochloride was higher at higher protein concentrations.     

The V108M mutation decreases the structural stability of catechol O-methyltransferase

The human gene for catechol O-methyltransferase has a common single-nucleotide polymorphism that results in substitution of methionine (M) for valine (V) 108 in the soluble form of the enzyme (s-COMT). 108M s-COMT loses enzymatic activity more rapidly than 108V s-COMT at physiological temperature, and the 108M allele has been associated with increased risk of breast cancer and several neuropsychiatric disorders. We used circular dichroism (CD), dynamic light scattering, and fluorescence spectroscopy to examine how the 108V/M polymorphism affects the stability of the purified, recombinant protein to heat and guanidine hydrochloride (GuHCl). COMT contains two tryptophan residues, W143 and W38Y, which are located in loops that border the S-adenosylmethionine (SAM) and catechol binding sites. We therefore also studied the single-tryptophan mutants W38Y and W143Y in order to dissect the contributions of the individual tryptophans to the fluorescence signals. The 108V and 108M proteins differed in the stability of both the tertiary structure surrounding the active site, as probed by the fluorescence yields and emission spectra, and their global secondary structure as reflected by CD. With either probe, the midpoint of the thermal transition of 108M s-COMT was 5 to 7 degrees C lower than that of 108V s-COMT, and the free energy of unfolding at 25 degrees C was smaller by about 0.4 kcal/mol. 108M s-COMT also was more prone to aggregation or partial unfolding to a form with an increased radius of hydration at 37 degrees C. The co-substrate SAM stabilized the secondary structure of both 108V and 108M s-COMT. W143 dominates the tryptophan fluorescence of the folded protein and accounts for most of the decrease in fluorescence that accompanies unfolding by GuHCl. While replacing either tryptophan by tyrosine was mildly destabilizing, the lower stability of the 108M variant was retained in all cases.     

Unfolding of the regulatory subunit of cAMP-dependent protein kinase I.

The unfolding of the recombinant regulatory subunit of cAMP-dependent protein kinase I was followed by monitoring the intrinsic protein fluorescence. Unfolding proceeds in at least two stages. First, the quenching of fluorescence due to cAMP binding is abolished at relatively low levels of urea (less than 2 M) and is observed as an increase in intensity at 340 nm. The high-affinity binding of cAMP is retained in 3 M urea even though the quenching is lost. The second stage of unfolding, presumably representing unfolding of the polypeptide chain, is seen as a shift in lambda max from 340 to 353 nm. The midpoint concentration, Cm, for this process is 5.0 M. Cyclic AMP binding activity is lost at a half-maximal urea concentration of 3.5 M and precedes the shift in lambda max. Unfolding of the protein in the presence of urea was fully reversible; furthermore, the presence of excess levels of cAMP stabilized the regulatory subunit. A free energy value (delta GDH2O) of 7.1 +/- 0.2 kcal/mol was calculated for the native form of the protein when denaturation was induced with either urea or guanidine hydrochloride. Iodide quenching of tryptophan fluorescence was used to elucidate the number of tryptophan residues accessible during various stages of the unfolding process. In the native cAMP-bound form of the regulatory subunit, only one of the three tryptophans in the regulatory subunit is quenched by iodide while more than two tryptophans can be quenched with iodide in the presence of 3 M urea.     

Thermodynamic stabilities of yeast iso-1-cytochromes c having amino acid substitutions for lysine 32

Iso-1-cytochromes c having lysine 32 replaced by leucine, glutamine, tyrosine, and tryptophan were prepared from strains of bakers' yeast, Saccharomyces cerevisiae, and chemically blocked at cysteine 107 with methyl methanethiolsulfonate to prevent dimerization. These modified ferricytochromes c were guanidine denatured, and the unfolding thermodynamics were determined by circular dichroism and fluorescence measurements. Thermal unfolding was also monitored by absorbance measurements. The guanidine denaturation midpoints for the altered proteins are smaller than the wild type, while the orders of stability from unfolding free energy changes are: Lys-32 (wild type) approximately Leu-32 approximately Gln-32 (circular dichroism), greater than Gln-32 (fluorescence) greater than Tyr-32 approximately Trp-32. Midpoints and differences in free energy changes for thermal unfolding parallel the fluorescence free energy changes for guanidine-induced unfolding. Thus, the blocked Leu-32 and Lys-32 proteins are equally stable with respect to both chemical and thermal denaturation. The reported data indicate that single replacements may significantly modify protein stability, and that substitution for an evolutionarily retained residue in normal cytochrome c structures does not always destabilize the protein. In addition, in vitro thermal stabilities approximately correlate with in vivo specific activities.     

Comparative analyses of the conformational stability of a hyperthermophilic protein and its mesophilic counterpart.

Comparison of the conformational stability of an O(6)-methylguanine-DNA methyltransferase (MGMT) from the hyperthermophilic archaeon Thermococcus kodakaraensis strain KOD1 (Tk-MGMT), and its mesophilic counterpart C-terminal Ada protein from Escherichia coli (Ec-AdaC) was performed in order to obtain information about the relationship between thermal stability and other factors, such as thermodynamic parameters, thermodynamic stability and other unfolding conditions. Tk-MGMT unfolded at Tm = 98.6 degrees C, which was 54.8 degrees C higher than the unfolding temperature of Ec-AdaC. The maximum free energy (DeltaG(max)) of the proteins were different; the value of Tk-MGMT (42.9 kJ.mol-1 at 29.5 degrees C) was 2.6 times higher than that of Ec-AdaC (16.6 kJ.mol-1 at 7.4 degrees C). The high conformational stability of Tk-MGMT was attributed to a 1.6-fold higher enthalpy value than that of Ec-AdaC. In addition, the DeltaG(max) temperature of Tk-MGMT was considerably higher (by 22.1 degrees C). The apparent heat capacity of denaturation (DeltaC(p)) of Tk-MGMT was 0.7-fold lower than that of Ec-AdaC. These three synergistic effects, increasing DeltaGmax, shifted DeltaG vs. temperature curve, and low DeltaC(p), give Tk-MGMT its thermal stability. Unfolding profiles of the two proteins, tested with four alcohols and three denaturants, showed that Tk-MGMT possessed higher stability than Ec-AdaC in all conditions studied. These results indicate that the high stability of Tk-MGMT gives resistance to chemical unfolding, in addition to thermal unfolding.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Thermal versus guanidine-induced unfolding of ubiquitin. An analysis in terms of the contributions from charge-charge interactions to protein stability.

We have characterized the guanidine-induced unfolding of both yeast and bovine ubiquitin at 25 degrees C and in the acidic pH range on the basis of fluorescence and circular dichroism measurements. Unfolding Gibbs energy changes calculated by linear extrapolation from high guanidine unfolding data are found to depend very weakly on pH. A simple explanation for this result involves the two following assumptions: (1) charged atoms of ionizable groups are exposed to the solvent in native ubiquitin (as supported by accessible surface area calculations), and Gibbs energy contributions associated with charge desolvation upon folding (a source of pK shifts) are small; (2) charge-charge interactions (another source of pK shifts upon folding) are screened out in concentrated guanidinium chloride solutions. We have also characterized the thermal unfolding of both proteins using differential scanning calorimetry. Unfolding Gibbs energy changes calculated from the calorimetric data do depend strongly on pH, a result that we attribute to the pH dependence of charge-charge interactions (not eliminated in the absence of guanidine). In fact, we find good agreement between the difference between the two series of experimental unfolding Gibbs energy changes (determined from high guanidine unfolding data by linear extrapolation and from thermal denaturation data in the absence of guanidine) and the theoretical estimates of the contribution from charge-charge interactions to the Gibbs energy change for ubiquitin unfolding obtained by using the solvent-accessibility-corrected Tanford-Kirkwood model, together with the Bashford-Karplus (reduced-set-of-sites) approximation. This contribution is found to be stabilizing at neutral pH, because most charged groups on the native protein interact mainly with groups of the opposite charge, a fact that, together with the absence of large charge-desolvation contributions, may explain the high stability of ubiquitin at neutral pH. In general, our analysis suggests the possibility of enhancing protein thermal stability by adequately redesigning the distribution of solvent-exposed, charged residues on the native protein surface.     

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40 MPa amplitude (5 ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500 MPa, between 30 and 50 degrees C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50 degrees C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113-Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect.     

Slow unfolding and refolding kinetics of the mesophilic Rop wild-type protein in the transition range.

We describe the guanidinium hydrochloride induced folding kinetics of the four-helix-bundle protein Rop wild-type (wt) under equilibrium conditions at three temperatures. The choice of appropriate denaturant conditions inside the transition range permitted, in combination with equilibrium transition curves, the determination of both unfolding and refolding rate constants. The ratio of the rate constants at zero denaturant concentration provided equilibrium constants and standard free energy changes that are in good agreement with values obtained in previous differential scanning calorimetry studies. The DeltaG0D values for 19, 25 and 40 degrees C calculated from the present kinetic studies are, respectively, 66.8, 70.8 and 57.2 kJ.mol-1. The unfolding reactions are extremely slow under these conditions. Equilibrium was reached only after 18, 12 and 6 days at 19, 25 and 40 degrees C. These results demonstrate that for Rop wt high stability correlates with slow folding kinetics.     

Probing folding and fluorescence quenching in human gammaD crystallin Greek key domains using triple tryptophan mutant proteins

Human gammaD crystallin (HgammaD-Crys), a major component of the human eye lens, is a 173-residue, primarily beta-sheet protein, associated with juvenile and mature-onset cataracts. HgammaD-Crys has four tryptophans, with two in each of the homologous Greek key domains, which are conserved throughout the gamma-crystallin family. HgammaD-Crys exhibits native-state fluorescence quenching, despite the absence of ligands or cofactors. The tryptophan absorption and fluorescence quenching may influence the lens response to ultraviolet light or the protection of the retina from ambient ultraviolet damage. To provide fluorescence reporters for each quadrant of the protein, triple mutants, each containing three tryptophan-to-phenylalanine substitutions and one native tryptophan, have been constructed and expressed. Trp 42-only and Trp 130-only exhibited fluorescence quenching between the native and denatured states typical of globular proteins, whereas Trp 68-only and Trp 156-only retained the anomalous quenching pattern of wild-type HgammaD-Crys. The three-dimensional structure of HgammaD-Crys shows Tyr/Tyr/His aromatic cages surrounding Trp 68 and Trp 156 that may be the source of the native-state quenching. During equilibrium refolding/unfolding at 37 degrees C, the tryptophan fluorescence signals indicated that domain I (W42-only and W68-only) unfolded at lower concentrations of GdnHCl than domain II (W130-only and W156-only). Kinetic analysis of both the unfolding and refolding of the triple-mutant tryptophan proteins identified an intermediate along the HgammaD-Crys folding pathway with domain I unfolded and domain II intact. This species is a candidate for the partially folded intermediate in the in vitro aggregation pathway of HgammaD-Crys.     

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.     

Osmolyte effects on kinetics of FKBP12 C22A folding coupled with prolyl isomerization

Unfolding and refolding kinetics of human FKBP12 C22A were monitored by fluorescence emission over a wide range of urea concentration in the presence and absence of protecting osmolytes glycerol, proline, sarcosine and trimethylamine-N-oxide (TMAO). Unfolding is well described by a mono-exponential process, while refolding required a minimum of two exponentials for an adequate fit throughout the urea concentration range considered. The bi-exponential behavior resulted from complex coupling between protein folding, and prolyl isomerization in the denatured state in which the urea-dependent rate constant for folding was greater than, equal to, and less than the rate constants for prolyl isomerization within the urea concentration range of zero to five molar. Amplitudes and the observed folding and unfolding rate constants were fitted to a reversible three-state model composed of two sequential steps involving the native state and a folding-competent denatured species thermodynamically linked to a folding-incompetent denatured species. Excellent agreement between thermodynamic parameters for FKBP12 C22A folding calculated from the kinetic parameters and those obtained directly from equilibrium denaturation assays provides strong support for the applicability of the mechanism, and provides evidence that FKBP12 C22A folding/unfolding is two-state, with prolyl isomer heterogeneity in the denatured ensemble. Despite the chemical diversity of the protecting osmolytes, they all exhibit the same kinetic behavior of increasing the rate constant of folding and decreasing the rate constant for unfolding. Osmolyte effects on folding/unfolding kinetics are readily explained in terms of principles established in understanding osmolyte effects on protein stability. These principles involve the osmophobic effect, which raises the Gibbs energy of the denatured state due to exposure of peptide backbone, thereby increasing the folding rate. This effect also plays a key role in decreasing the unfolding rate when, as is often the case, the activated complex exposes more backbone than is exposed in the native state.     

The role of Trp-82 in the folding of intestinal fatty acid binding protein

Multiple phases have been observed during the folding and unfolding of intestinal fatty acid binding protein (WT-IFABP) by stopped-flow fluorescence. Site-directed mutagenesis has been used to examine the role of each of the two tryptophans of this protein in these processes. The unfolding and refolding kinetics of the mutant protein containing only tryptophan 82 (W6Y-IFABP) showed that the tryptophan at this location was critical to the fluorescence signal changes observed throughout the unfolding reaction and early in the refolding reaction. However, the kinetic patterns of the mutant protein containing only tryptophan 6 (W82Y-IFABP) indicated that the tryptophan at this location participated in the fluorescence signal changes observed early in the unfolding reaction and late in the refolding reaction. Together, these data suggest that native-like structure was formed first in the vicinity of tryptophan 82, near the center of the hydrophobic core of this beta-sheet protein, prior to formation of native-like structure in the periphery of the protein.     

Equilibrium unfolding of Bombyx mori glycyl-tRNA synthetase

Unfolding of Bombyx mori glycyl-tRNA synthetase was examined by multiple spectroscopic techniques. Tryptophan fluorescence of wild type enzyme and an N-terminally truncated form (N55) increased at low concentrations of urea or guanidine-HCl followed by a reduction in intensity at intermediate denaturant concentrations; a transition at higher denaturant was detected as decreased fluorescence intensity and a red-shifted emission. Solute quenching of fluorescence indicated that tryptophans become progressively solvent-exposed during unfolding. Wild type enzyme had stronger negative CD bands between 220 and 230 nm than the mutant, indicative of greater alpha-helical content. Urea or guanidine-HCl caused a reduction in ellipticity at 222 nm at low denaturant concentration with the wild type enzyme, a transition that is absent in the mutant; both enzymes exhibited a cooperative transition at higher denaturant concentrations. Both enzymes dissociate to monomers in 1.5 m urea. Unfolding of wild type enzyme is described by a multistate unfolding and a parallel two state unfolding; the two-state component is absent in the mutant. Changes in spectral properties associated with unfolding were largely reversible after dilution to low denaturant. Unfolding of glycyl-tRNA synthetase is complex with a native state, a native-like monomer, partially unfolded states, and the unfolded state.     

Conformational stability of ribonuclease T1 determined by hydrogen-deuterium exchange.

The hydrogen-deuterium exchange kinetics of 37 backbone amide residues in RNase T1 have been monitored at 25, 40, 45, and 50 degrees C at pD 5.6 and at 40 and 45 degrees C at pD 6.6. The hydrogen exchange rate constants of the hydrogen-bonded residues varied over eight orders of magnitude at 25 degrees C with 13 residues showing exchange rates consistent with exchange occurring as a result of global unfolding. These residues are located in strands 2-4 of the central beta-pleated sheet. The residues located in the alpha-helix and the remaining strands of the beta-sheet exhibited exchange behaviors consistent with exchange occurring due to local structural fluctuations. For several residues at 25 degrees C, the global free energy change calculated from the hydrogen exchange data was over 2 kcal/mol greater than the free energy of unfolding determined from urea denaturation experiments. The number of residues showing this unexpected behavior was found to increase with temperature. This apparent inconsistency can be explained quantitatively if the cis-trans isomerization of the two cis prolines, Pro-39 and Pro-55, is taken into account. The cis-trans isomerization equilibrium calculated from kinetic data indicates the free energy of the unfolded state will be 2.6 kcal/mol higher at 25 degrees C when the two prolines are cis rather than trans (Mayr LM, Odefey CO, Schutkowski M, Schmid FX. 1996. Kinetic analysis of the unfolding and refolding of ribonuclease T1 by a stopped-flow double-mixing technique. Biochemistry 35: 5550-5561). The hydrogen exchange results are consistent with the most slowly exchanging hydrogens exchanging from a globally higher free energy unfolded state in which Pro-55 and Pro-39 are still predominantly in the cis conformation. When the conformational stabilities determined by hydrogen exchange are corrected for the proline isomerization equilibrium, the results are in excellent agreement with those from an analysis of urea denaturation curves.     

Low-temperature unfolding of a mutant of phage T4 lysozyme. 1. Equilibrium studies

The mutant protein I3C-C97/C54T of phage T4 lysozyme is free of sulfhydryl groups and has a genetically engineered disulfide bridge between positions 3 and 97 (Perry & Wetzel, 1986). This protein has a maximum stability at 12 degrees C in 3 M guanidinium chloride and undergoes reversible high- and low-temperature melting at 28 and -3 degrees C, respectively, in this medium. The free energy of stabilization of the protein has been studied over a range of temperature that includes both melting transitions. The stability curve fits a constant delta Cp model over the entire range, permitting an unusually complete determination of the thermodynamic parameters of the protein and demonstrating that the low-temperature unfolded form of the protein may be interpreted as an extrapolation with constant delta Cp of the high-temperature unfolded form. The free energy of unfolding is a linear function of guanidinium concentration within experimental error which permits a rough estimate of the stability of the protein at low temperatures and of the differential interaction of the unfolded protein with guanidinium chloride. These equilibrium studies provide a basis for the interpretation of the kinetic studies reported in the following paper.