Equilibrium and kinetic folding of an alpha-helical Greek key protein domain: caspase recruitment domain (CARD) of RICK

We have characterized the equilibrium and kinetic folding of a unique protein domain, caspase recruitment domain (CARD), of the RIP-like interacting CLARP kinase (RICK) (RICK-CARD), which adopts a alpha-helical Greek key fold. At equilibrium, the folding of RICK-CARD is well described by a two-state mechanism representing the native and unfolded ensembles. The protein is marginally stable, with a DeltaG(H)()2(O) of 3.0 +/- 0.15 kcal/mol and an m-value of 1.27 +/- 0.06 kcal mol(-1) M(-1) (30 mM Tris-HCl, pH 8, 1 mM DTT, 25 degrees C). While the m-value is constant, the protein stability decreases in the presence of moderate salt concentrations (below 200 mM) and then increases at higher salt concentrations. The results suggest that electrostatic interactions are stabilizing in the native protein, and the favorable Coulombic interactions are reduced at low ionic strength. Above 200 mM salt, the results are consistent with Hofmeister effects. The unfolding pathway of RICK-CARD is complex and contains at least three non-native conformations. The refolding pathway of RICK-CARD also is complex, and the data suggest that the unfolded protein folds via two intermediate conformations prior to reaching the native state. Overall, the data suggest the presence of kinetically trapped, or misfolded, species that are on-pathway both in refolding and in unfolding.     

NMR analysis of staphylococcal nuclease thermal quench refolding kinetics.

Thermally unfolded staphylococcal nuclease has been rapidly quenched to temperatures near 0 degree C and the refolding behavior examined using an NMR kinetic experiment. Unfolded protein, exhibiting random coil chemical shifts, persists following the quench and refolds in two distinct kinetic phases. A protein folding intermediate with a trans Lys 116-Pro 117 peptide bond is transiently overpopulated and relaxes to the predominantly cis native cis-trans equilibrium. The rate of trans-->cis isomerization in the native-like nuclease intermediate is approximately 100-fold faster than that observed in a Lys-Pro model peptide. The activation enthalpy of 20 kcal/mol observed for the nuclease Lys 116-Pro 117 peptide bond is comparable to that observed for other X-Pro isomerizations.     

A change in the apparent m value reveals a populated intermediate under equilibrium conditions in Escherichia coli ribonuclease HI

Experimental studies of protein stability often rely on the determination of an "m value", which describes the denaturant dependence of the free energy change between two states (DeltaG = DeltaG(H2O) - m[denaturant]). Changes in the m value accompanying site specific mutations are usually attributed to structural alterations in the native or unfolded ensemble. Here, we provide an example of significant reduction in the m value resulting from a subtle deviation in two-state behavior not detected by traditional methods. The protein that is studied is a variant of Escherchia coli RNase H in which three residues predicted to be involved in a partially buried salt bridge network were mutated to alanine (R46A, D102A, and D148A). Equilibrium denaturant profiles monitored by both fluorescence and circular dichroism appeared to be cooperative, and a two-state analysis yielded a DeltaG(UN) of approximately -3 kcal/mol with an m value of 1.4 kcal mol(-1) M(-1) (vs 2.3 for RNase H). Analysis of kinetic refolding experiments suggests that the system is actually three-state at equilibrium with an appreciable concentration of an intermediate state under low denaturant concentrations. The stability of the native state determined from a fit of these kinetic data is -6.7 kcal/mol, suggesting that the stability determined by traditional two-state equilibrium analysis is a gross underestimate. The only hint to this loss of two-state behavior was a decrease in the apparent m value, and the presence of the equilibrium intermediate was only identified by a kinetic analysis. Our work serves as a cautionary note; the possibility of a three-state system should be closely addressed before interpreting a change in the m value as a change in the native or unfolded state.     

Stability and cooperativity of individual tertiary contacts in RNA revealed through chemical denaturation

For proteins, understanding tertiary interactions involved in local versus global unfolding has become increasingly important for understanding the nature of the native state ensemble, the mechanisms of unfolding, and the stability of both the native and intermediate states in folding. In this work we have addressed related questions with respect to RNA structure by combining chemical denaturation and hydroxyl radical footprinting methods. We have determined unfolding isotherms for each of 26 discrete sites of protection located throughout the Tetrahymena thermophila group I ribozyme. The cooperativity of folding, m-value, and the free energy, DeltaG degrees N-U, associated with formation of each tertiary contact was determined by analysis of the isotherms. The DeltaG degrees N-U values measured in this study vary from 1.7 +/- 0.2 to 7. 6 +/- 1.2 kcal mol-1. Thus, the stability of these discrete tertiary contacts vary by almost 104. In addition, an intradomain contact and three interdomain contacts show high cooperativity (m-values of 1.1 +/- 0.2 to 1.7 +/- 0.3 kcal mol-1 M-1) indicating that these contacts exhibit global cooperatively in their folding behavior. This new approach to examining RNA stability provides an exciting comparison to our understanding of protein structure and folding mechanisms.     

Energetics and cooperativity of tertiary hydrogen bonds in RNA structure.

Tertiary interactions that allow RNA to fold into intricate three-dimensional structures are being identified, but little is known about the thermodynamics of individual interactions. Here we quantify the tertiary structure contributions of individual hydrogen bonds in a "ribose zipper" motif of the recently crystallized Tetrahymena group I intron P4-P6 domain. The 2'-hydroxyls of P4-P6 nucleotides C109/A184 and A183/G110 participate in forming the "teeth" of the zipper. These four nucleotides were substituted in all combinations with their 2'-deoxy and (separately) 2'-methoxy analogues, and thermodynamic effects on the tertiary folding DeltaG degrees ' were assayed by the Mg2+ dependence of electrophoretic mobility in nondenaturing gels. The 2'-deoxy series showed a consistent trend with an average contribution to the tertiary folding DeltaG degrees' of -0.4 to -0.5 kcal/mol per hydrogen bond. Contributions were approximately additive, reflecting no cooperativity among the hydrogen bonds. Each "tooth" of the ribose zipper (comprising two hydrogen bonds) thus contributes about -1.0 kcal/mol to the tertiary folding DeltaG degrees'. Single 2'-methoxy substitutions destabilized folding by approximately 1 kcal/mol, but the trend reversed with multiple 2'-methoxy substitutions; the folding DeltaG degrees' for the quadruple 2'-methoxy derivative was approximately unchanged relative to wild-type. On the basis of these data and on temperature-gradient gel results, we conclude that entropically favorable hydrophobic interactions balance enthalpically unfavorable hydrogen bond deletions and steric clashes for multiple 2'-methoxy substitutions. Because many of the 2'-deoxy derivatives no longer have the characteristic hydrogen-bond patterns of the ribose zipper motif but simply have individual long-range ribose-base or ribose-ribose hydrogen bonds, we speculate that the energetic value of -0.4 to -0.5 kcal/mol per tertiary hydrogen bond may be more generally applicable to RNA folding.     

Contribution of hydrogen bonding to protein stability estimated from isotope effects

An unresolved issue in structural biology concerns the relative contribution of H bonds to protein stability. We use the small molecules 4-acetamidobenzoic acid and N-acetylanthranilic acid as model compounds to relate the energetic contribution from hydrogen bonds (H bonds) to the deuterium/hydrogen amide isotope effect. N-Acetylanthranilic acid models carbonyl-amide H bonds formed during protein folding; 4-acetamidobenzoic acid models the unfolded state in which the amide H bonds to water. NMR is used to measure shifts in the pK(a) of the ionizable carboxyl group when the amides of the compounds are either protonated or deuterated. From the pK(a) shift, we obtain a quantitative scale factor: SF = partial partial differential(DeltaG(HB))/partial partial differential(RT ln Phi), where DeltaG(HB) is the change in free energy of an H bond upon isotope substitution and Phi is the fractionation factor. Isotope effect data also are reported for a small globular protein, lambda repressor, using the "C(m) experiment". The protein's isotope effect, which reports on the shape of the energy well, is converted to H-bonding free energy by applying the scale factor. We estimate that amide-related H bonds (amide-carbonyl and amide-water) contribute favorably to protein stability by approximately 30-50 kcal/mol in lambda repressor, GCN4 coiled coil, and cytochrome c but unfavorably by approximately 6 kcal/mol in ubiquitin. The results indicate that H-bond strength varies from one protein to another and presumably at different sites within the same protein.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

Thermodynamics of the alpha-helix-coil transition of amphipathic peptides in a membrane environment: implications for the peptide-membrane binding equilibrium

Amphipathic alpha-helices are the membrane binding motif in many proteins. The corresponding peptides are often random coil in solution but are folded into an alpha-helix upon interaction with the membrane. The energetics of this ubiquitous folding process are still a matter of conjecture. Here, we present a new method to quantitatively analyze the thermodynamics of peptide folding at the membrane interface. We have systematically varied the helix content of a given amphipathic peptide when bound to the membrane and have correlated the thermodynamic binding parameters determined by isothermal titration calorimetry with the alpha-helix content obtained by circular dichroism spectroscopy. The peptides investigated were the antibiotic magainin 2 amide and three analogs in which two adjacent amino acid residues were substituted by their d-enantiomers. The thermodynamic parameters controlling the alpha-helix formation were found to be linearly related to the helicity of the membrane-bound peptides. Helix formation at the membrane surface is characterized by an enthalpy change of DeltaH(helix) approximately -0.7 kcal/mol per residue, an entropy change of DeltaS(helix) approximately -1.9 cal/molK residue and a free energy change of DeltaG(helix)=-0.14 kcal/mol residue. Helix formation is a strong driving force of peptide insertion into the membrane and accounts for about 50 % of the free energy of binding. An increase in temperature entails an unfolding of the membrane-bound helix. The temperature dependence can be described with the Zimm-Bragg theory and the enthalpy of unfolding agrees with that deduced from isothermal titration calorimetry.     

Folding mechanism of FIS, the intertwined, dimeric factor for inversion stimulation

FIS, the factor for inversion stimulation, from Escherichia coli and other enteric bacteria, is an interwined alpha-helical homodimer. Size exclusion chromatography and static light scattering measurements demonstrated that FIS is predominately a stable dimer at the concentrations (1-10 microM monomer) and buffer conditions employed in this study. The folding and unfolding of FIS were studied with both equilibrium and kinetic methods by circular dichroism using urea and guanidinium chloride (GdmCl) as the perturbants. The equilibrium folding is reversible and well-described by a two-state folding model, with stabilities at 10 degrees C of 15.2 kcal mol(-1) in urea and 13.5 kcal mol(-1) in GdmCl. The kinetic data are consistent with a two-step folding reaction where the two unfolded monomers associate to a dimeric intermediate within the mixing time for the stopped-flow instrument (<5 ms), and a slower, subsequent folding of the dimeric intermediate to the native dimer. Fits of the burst phase amplitudes as a function of denaturant showed that the free energy for the formation of the dimeric intermediate constitutes the majority of the stability of the folding (9.6 kcal mol(-1) in urea and 10.5 kcal mol(-1) in GdmCl). Folding-to-unfolding double jump kinetic experiments were also performed to monitor the formation of native dimer as a function of folding delay times. The data here demonstrate that the dimeric intermediate is obligatory and on-pathway. The folding mechanism of FIS, when compared to other intertwined, alpha-helical, homodimers, suggests that a transient kinetic dimeric intermediate may be a common feature of the folding of intertwined, segment-swapped, alpha-helical dimers.     

Folding of a de novo designed native-like four-helix bundle protein

The folding of a model native-like dimeric four-helix bundle protein, (alpha(2))(2), was investigated using guanidine hydrochloride, hydrostatic pressure, and low temperature. Unfolding by guanidine hydrochloride followed by circular dichroism and intrinsic fluorescence spectroscopy revealed a highly cooperative transition between the native-like and unfolded states, with free energy of unfolding determined from CD data, DeltaG(unf) = 14.3 +/- 0.8 kcal/mol. However, CD and intrinsic fluorescence data were not superimposable, indicating the presence of an intermediate state during the folding transition. To stabilize the folding intermediate, we used hydrostatic pressure and low temperature. In both cases, dissociation of the dimeric native-like (alpha(2))(2) into folded monomers (alpha(2)) was observed. van't Hoff analysis of the low temperature experiments, assuming a two-state dimer 171-monomer transition, yielded a free energy of dissociation of (alpha(2))(2) of DeltaG(diss) = 11.4 +/- 0.4 kcal/mol, in good agreement with the free energy determined from pressure dissociation experiments (DeltaG(diss) = 10.5 +/- 0.1 kcal/mol). Binding of the hydrophobic fluorescent probe 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) to the pressure- and cold-dissociated states of (alpha(2))(2) indicated the existence of molten-globule monomers. In conclusion, we demonstrate that the folding pathway of (alpha(2))(2) can be described by a three-state transition including a monomeric molten globule-like state.     

A general approach for detecting folding intermediates from steady-state and time-resolved fluorescence of single-tryptophan-containing proteins

During denaturant-induced equilibrium (un)folding of wild-type apoflavodoxin from Azotobacter vinelandii, a molten globule-like folding intermediate is formed. This wild-type protein contains three tryptophans. In this study, we use a general approach to analyze time-resolved fluorescence and steady-state fluorescence data that are obtained upon denaturant-induced unfolding of a single-tryptophan-containing variant of apoflavodoxin [i.e., W74/F128/F167 (WFF) apoflavodoxin]. The experimental data are assembled in matrices, and subsequent singular-value decomposition of these matrices (i.e., based on either steady-state or time-resolved fluorescence data) shows the presence of three significant, and independent, components. Consequently, to further analyze the denaturation trajectories, we use a three-state protein folding model in which a folding intermediate and native and unfolded protein molecules take part. Using a global analysis procedure, we determine the relative concentrations of the species involved and show that the stability of WFF apoflavodoxin against global unfolding is ～4.1 kcal/mol. Analysis of time-resolved anisotropy data of WFF apoflavodoxin unfolding reveals the remarkable observation that W74 is equally well fixed within both the native protein and the molten globule-like folding intermediate. Slight differences between the direct environments of W74 in the folding intermediate and native protein cause different rotameric populations of the indole in both folding species as fluorescence lifetime analysis reveals. Importantly, thermodynamic analyses of the spectral denaturation trajectories of the double-tryptophan-containing protein variants WWF apoflavodoxin and WFW apoflavodoxin show that these variants are significantly more stable (5.9 kcal/mol and 6.8 kcal/mol, respectively) than WFF apoflavodoxin (4.1 kcal/mol) Hence, tryptophan residues contribute considerably to the 10.5 kcal/mol thermodynamic stability of native wild-type apoflavodoxin.     

Interleukin-1 beta folding between pH 5 and 7: experimental evidence for three-state folding behavior and robust transition state positions late in folding

The thermodynamic stability and folding kinetics of the all beta-sheet protein interleukin-1beta were measured between 0 and 4 M GdmCl concentrations and pH 5-7. Native interleukin-1beta undergoes a 3.5 kcal/mol decrease in thermodynamic stability, Delta, as pH is increased from 5 to 7. The native state parameter m(NU), measuring protein destabilization/[GdmCl], remains constant between pH 5 and 7, indicating that the solvent-exposed surface area difference between the native state and unfolded ensemble is unchanged across this pH range. Similarly, pH changes between 5 and 7 decrease only the thermodynamic stability, DeltaG(H)2(O), and not the m-values, of the kinetic intermediate and transition states. This finding is shown to be consistent with transition state configurations which continue to be the high-energy configurations of the transition state in the face of changing stability conditions. A three-state folding mechanism U right arrow over left arrow I right arrow over left arrow N is shown to be sufficient in characterizing IL-1beta folding under all conditions studied. The m-values of refolding transitions are much larger than the m-values of unfolding transitions, indicating that that the fast, T(2) (U right arrow over left arrow I), and slow, T(1) (I right arrow over left arrow N), transition states are highly similar to the intermediate I and native state N, respectively. Many of the folding properties of interleukin-1beta are shared among other members of the beta-trefoil protein family, although clear differences can exist.     

Effects of tryptophan microenvironment, soluble domain, and vesicle size on the thermodynamics of membrane protein folding: lessons from the transmembrane protein OmpA

Refolding curves of the integral membrane protein outer membrane protein A (OmpA) were measured to determine the conformational stabilities of this model system for membrane protein folding. Wild-type OmpA exhibits a free energy of unfolding (DeltaG degrees H2O) of 10.5 kcal/mol. Mutants, containing a single tryptophan residue at the native positions 7, 15, 57, 102, or 143, are less stable than wild-type OmpA, with DeltaG degrees H2O values of 6.7, 4.8, 2.4, 4.7, and 2.8 kcal/mol, respectively. The trend observed here is discussed in terms of noncovalent interactions, including aromatic interactions and hydrogen bonding. The effect of the soluble tail on the conformational stability of the transmembrane domain of OmpA was also investigated via truncated single-Trp mutants; DeltaG degrees H2O values for four of the five truncated mutants are greater by >2.7 kcal/mol relative to the full-length versions, suggesting that the absence of the soluble domain may destabilize the unfolded transmembrane domain. Finally, dynamic light scattering experiments were performed to measure the effects of urea and protein on vesicle size and stability. Urea concentrations greater than 1 M cause an increase in vesicle size, and these diameters are unaltered in the presence of protein. These dynamic light scattering results complement the fluorescence studies and illustrate the important effects of vesicle size on protein conformational stability.     

Thermodynamics of amide hydrogen bond formation in polar and apolar solvents

We present the initial findings of a theoretical study of hydrogen bond formation between two formamide molecules in water and in carbon tetrachloride. These systems were chosen as the simplest models for secondary structure formation in the polar environment near the protein surface and the apolar environment of the protein interior. We have employed thermodynamic simulation methods to obtain absolute binding free energies and free energy profiles for the formation of peptide hydrogen bonds in the two solvents. We find that the amide hydrogen bond is stable by 8.4 kcal/mol in CCl4, and by 0.3 kcal/mol in water. Our results indicate also that the hydrogen-bonded dimer is 2.2 kcal/mol more stable in water than it is in CCl4. We compare our results with those from experiment, and discuss their use in interpreting mechanisms of protein folding.     

Simulation of folding of a small alpha-helical protein in atomistic detail using worldwide-distributed computing

By employing thousands of PCs and new worldwide-distributed computing techniques, we have simulated in atomistic detail the folding of a fast-folding 36-residue alpha-helical protein from the villin headpiece. The total simulated time exceeds 300 micros, orders of magnitude more than previous simulations of a molecule of this size. Starting from an extended state, we obtained an ensemble of folded structures, which is on average 1.7A and 1.9A away from the native state in C(alpha) distance-based root-mean-square deviation (dRMS) and C(beta) dRMS sense, respectively. The folding mechanism of villin is most consistent with the hydrophobic collapse view of folding: the molecule collapses non-specifically very quickly ( approximately 20ns), which greatly reduces the size of the conformational space that needs to be explored in search of the native state. The conformational search in the collapsed state appears to be rate-limited by the formation of the aromatic core: in a significant fraction of our simulations, the C-terminal phenylalanine residue packs improperly with the rest of the hydrophobic core. We suggest that the breaking of this interaction may be the rate-determining step in the course of folding. On the basis of our simulations we estimate the folding rate of villin to be approximately 5micros. By analyzing the average features of the folded ensemble obtained by simulation, we see that the mean folded structure is more similar to the native fold than any individual folded structure. This finding highlights the need for simulating ensembles of molecules and averaging the results in an experiment-like fashion if meaningful comparison between simulation and experiment is to be attempted. Moreover, our results demonstrate that (1) the computational methodology exists to simulate the multi-microsecond regime using distributed computing and (2) that potential sets used to describe interatomic interactions may be sufficiently accurate to reach the folded state, at least for small proteins. We conclude with a comparison between our results and current protein-folding theory.     

SecA folds via a dimeric intermediate

Though many proteins in the cell are large and multimeric, their folding has not been extensively studied. We have chosen SecA as a folding model because it is a large, homodimeric protein (monomer molecular mass of 102 kDa) with multiple folding domains. SecA is the ATPase for the Sec-dependent preprotein translocase of many bacteria. SecA is a soluble protein that can penetrate into the membrane during preprotein translocation. Because SecA may partially unfold prior to its insertion into the membrane, studies of its stability and folding pathway are important for understanding how it functions in vivo. Kinetic folding transitions in the presence of urea were monitored using circular dichroism and tryptophan fluorescence, while equilibrium folding transitions were monitored using the same techniques as well as a fluorescent ATP analogue. The reversible equilibrium folding transition exhibited a plateau, indicating the presence of an intermediate. Based on the data presented here, we propose a three-state model, N(2) if I(2) if 2U, where the native protein unfolds to a dimeric intermediate which then dissociates into two unfolded monomers. The SecA dimer was determined to have an overall stability (DeltaG) of -22.5 kcal/mol. We also investigated the stability of SecA using analytical ultracentrifugation equilibrium and velocity sedimentation, which again indicated that native or refolded SecA was a stable dimer. The rate-limiting step in the folding pathway was conversion of the dimeric intermediate to the native dimer. Unfolding of native, dimeric SecA was slow with a relaxation time in H(2)O of 3.3 x 10(4) s. Since SecA is a stable dimer, dissociation to monomeric subunits during translocation is unlikely.     

Distributions in protein conformation space: implications for structure prediction and entropy

By considering how polymer structures are distributed in conformation space, we show that it is possible to quantify the difficulty of structural prediction and to provide a measure of progress for prediction calculations. The critical issue is the probability that a conformation is found within a specified distance of another conformer. We address this question by constructing a cumulative distribution function (CDF) for the average probability from observations about its limiting behavior at small displacements and numerical simulations of polyalanine chains. We can use the CDF to estimate the likelihood that a structure prediction is better than random chance. For example, the chance of randomly predicting the native backbone structure of a 150-amino-acid protein to low resolution, say within 6 A, is 10(-14). A high-resolution structural prediction, say to 2 A, is immensely more difficult (10(-57)). With additional assumptions, the CDF yields the conformational entropy of protein folding from native-state coordinate variance. Or, using values of the conformational entropy change on folding, we can estimate the native state's conformational span. For example, for a 150-mer protein, equilibrium alpha-carbon displacements in the native ensemble would be 0.3-0.5 A based on T Delta S of 1.42 kcal/(mol residue).     

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.     

Multistate equilibrium unfolding of Escherichia coli dihydrofolate reductase: thermodynamic and spectroscopic description of the native, intermediate, and unfolded ensembles

The thermodynamic and spectroscopic properties of a cysteine-free variant of Escherichia coli dihydrofolate reductase (AS-DHFR) were investigated using the combined effects of urea and temperature as denaturing agents. Circular dichroism (CD), absorption, and fluorescence spectra were recorded during temperature-induced unfolding at different urea concentrations and during urea-induced unfolding at different temperatures. The first three vectors obtained by singular-value decomposition of each set of unfolding spectra were incorporated into a global analysis of a unique thermodynamic model. Although individual unfolding profiles can be described as a two-state process, a simultaneous fit of 99 vectors requires a three-state model as the minimal scheme to describe the unfolding reaction along both perturbation axes. The model, which involves native (N), intermediate (I), and unfolded (U) states, predicts a maximum apparent stability, DeltaG degrees (NU), of 6 kcal mol(-)(1) at 15 degrees C, an apparent m(NU) value of 2 kcal mol(-)(1) M(-)(1), and an apparent heat capacity change, DeltaC(p)()(-NU), of 2.5 kcal mol(-)(1) K(-)(1). The intermediate species has a maximum stability of approximately 2 kcal mol(-)(1) and a compactness closer to that of the native than to that of the unfolded state. The population of the intermediate is maximal ( approximately 70%) around 50 degrees C and falls below the limits of detection of > or =2 M urea or at temperatures of <35 or >65 degrees C. The fluorescence properties of the equilibrium intermediate resemble those of a transient intermediate detected during refolding from the urea-denatured state, suggesting that a tryptophan-containing hydrophobic cluster in the adenosine-binding domain plays a key role in both the equilibrium and kinetic reactions. The CD spectroscopic properties of the native state reveal the presence of two principal isoforms that differ in ligand binding affinities and in the packing of the adenosine-binding domain. The relative populations of these species change slightly with temperature and do not depend on the urea concentration, implying that the two native isoforms are well-structured and compact. Global analysis of data from multiple spectroscopic probes and several methods of unfolding is a powerful tool for revealing structural and thermodynamic properties of partially and fully folded forms of DHFR.     

Dimeric procaspase-3 unfolds via a four-state equilibrium process.

We have examined the folding and assembly of a catalytically inactive mutant of procaspase-3, a homodimeric protein that belongs to the caspase family of proteases. The caspase family, and especially caspase-3, is integral to apoptosis. The equilibrium unfolding data demonstrate a plateau between 3 and 5 M urea, consistent with an apparent three-state unfolding process. However, the midpoint of the second transition as well as the amplitude of the plateau are dependent on the protein concentration. Overall, the data are well described by a four-state equilibrium model in which the native dimer undergoes an isomeration to a dimeric intermediate, and the dimeric intermediate dissociates to a monomeric intermediate, which then unfolds. By fitting the four-state model to the experimental data, we have determined the free energy change for the first step of unfolding to be 8.3 +/- 1.3 kcal/mol. The free energy change for the dissociation of the dimeric folding intermediate to two monomeric intermediates is 10.5 +/- 1 kcal/mol. The third step in the unfolding mechanism represents the complete unfolding of the monomeric intermediate, with a free energy change of 7.0 +/- 0.5 kcal/mol. These results show two important points. First, dimerization of procaspase-3 occurs as a result of the association of two monomeric folding intermediates, demonstrating that procaspase-3 dimerization is a folding event. Second, the stability of the dimer contributes significantly to the conformational free energy of the protein (18.8 of 25.8 kcal/mol).