Folding of dihydrofolate reductase from Escherichia coli

The urea-induced equilibrium unfolding transition of dihydrofolate reductase from Escherichia coli was monitored by UV difference, circular dichroism (CD), and fluorescence spectroscopy. Each of these data sets were well described by a two-state unfolding model involving only native and unfolded forms. The free energy of folding in the absence of urea at pH 7.8, 15 degrees C is 6.13 +/- 0.36 kcal mol-1 by difference UV, 5.32 +/- 0.67 kcal mol-1 by CD, and 5.42 +/- 1.04 kcal mol-1 by fluorescence spectroscopy. The midpoints for the difference UV, CD, and fluorescence transitions are 3.12, 3.08, and 3.18 M urea, respectively. The near-coincidence of the unfolding transitions monitored by these three techniques also supports the assignment of a two-state model for the equilibrium results. Kinetic studies of the unfolding and refolding reactions show that the process is complex and therefore that additional species must be present. Unfolding jumps in the absence of potassium chloride revealed two slow phases which account for all of the amplitude predicted by equilibrium experiments. Unfolding in the presence of 400 mM KCl results in the selective loss of the slower phase, implying that there are two native forms present in equilibrium prior to unfolding. Five reactions were observed in refolding: two slow phases designated tau 1 and tau 2 that correspond to the slow phases in unfolding and three faster reactions designated tau 3, tau 4, and tau 5 that were followed by stopped-flow techniques. The kinetics of the recovery of the native form was monitored by following the binding of methotrexate, a tight-binding inhibitor of dihydrofolate reductase, at 380 nm.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Replacement of proline with valine does not remove an apparent proline isomerization-dependent folding event in CRABP I

Site-directed mutagenesis has frequently been used to replace proline with other amino acids in order to determine if proline isomerization is responsible for a slow phase during refolding. Replacement of Pro 85 with alanine in cellular retinoic acid binding protein I (CRABP-I) abolished the slowest refolding phase, suggesting that this phase is due to proline isomerization in the unfolded state. To further test this assumption, we mutated Pro 85 to valine, which is the conservative replacement in the two most closely related proteins in the family (cellular retinoic acid binding protein II and cellular retinol binding protein I). The mutant protein was about 1 kcal/mole more stable than wild type. Retinoic acid bound equally well to wild type and P85V-CRABP I, confirming the functional integrity of this mutation. The refolding and unfolding kinetics of the wild-type and mutant proteins were characterized by stopped flow fluorescence and circular dichroism. The mutant P85V protein refolded with three kinetic transitions, the same number as wild-type protein. This result conflicts with the P85A mutant, which lost the slowest refolding rate. The P85V mutation also lacked a kinetic unfolding intermediate found for wild-type protein. These data suggest that proline isomerization may not be responsible for the slowest folding phase of CRABP I. As such, the loss of a slow refolding phase upon mutation of a proline residue may not be diagnostic for proline isomerization effects on protein folding.     

Substitutions of prolines examine their role in kinetic trap formation of the caspase recruitment domain (CARD) of RICK

Caspase recruitment domains (CARDs) are small helical protein domains that adopt the Greek key fold. For the two CARDs studied to date, RICK-CARD and caspase-1-CARD (CP1-CARD), the proteins unfold by an apparent two-state process at equilibrium. However, the folding kinetics are complex for both proteins and may contain kinetically trapped species on the folding pathway. In the case of RICK-CARD, the time constants of the slow refolding phases are consistent with proline isomerism. RICK-CARD contains three prolines, P47 in turn 3, and P85 and P87. The latter two prolines constitute a nonconserved PxP motif in helix 6. To examine the role of the prolines in the complex folding kinetics of RICK-CARD, we generated seven proline-to-alanine mutants, including three single mutants, three double mutants, and one triple mutant. We examined the spectroscopic properties, equilibrium folding, binding to CP1-CARD, and folding kinetics. The results show that P85 is critical for maintaining the function of the protein and that all mutations decrease the stability. Results from single mixing and sequential mixing stopped-flow studies strongly suggest the presence of parallel folding pathways consisting of at least two unfolded populations. The mutations affect the distribution of the two unfolded species, thereby affecting the population that folds through each channel. The two conformations also are present in the triple mutant, demonstrating that interconversion between them is not due to prolyl isomerism. Overall, the data show that the complex folding pathway, especially formation of kinetically trapped species, is not due to prolyl isomerism.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

Reversible unfolding of bovine odorant binding protein induced by guanidinium hydrochloride at neutral pH

An analysis of the unfolding and refolding curves at equilibrium of dimeric bovine odorant binding protein (bOBP) has been performed. Unfolding induced by guanidinium chloride (GdnHCl) is completely reversible as far as structure and ligand binding capacity are concerned. The transition curves, as obtained by fluorescence and ellipticity measurements, are very similar and have the same protein concentration-independent midpoint (C1/2 approximately 2.6 M). This result implies a sequential, rather than a concerted, unfolding mechanism, with the involvement of an intermediate. However, since it has not been detected, this intermediate must be present in small amounts or have the same optical properties of either native or denatured protein. The thermodynamic best fit parameters, obtained according to a simple two-state model, are: deltaG degrees un,w = 5.0 +/- 0.6 kcal mol(-1), m = 1.9 +/- 0.2 kcal mol(-1) M(-1) and C1/2 = 2.6 +/- 0.1 M. The presence of the ligand dihydromyrcenol has a stabilising effect against unfolding by GdnHCl, with an extrapolated deltaG degrees un,w of 22.2 +/- 0.9 kcal mol(-1), a cooperative index of 3.2 +/- 0.3 and a midpoint of 4.6 +/- 0.4 M. The refolding curves, recorded after 24 h from dilution of denaturant are not yet at equilibrium: they show an apparently lower midpoint (C1/2 = 2.2 M), but tend to overlap the unfolding curve after several days. In contrast to chromatographic unfolding data, which fail to reveal the presence of folded intermediates, chromatographic refolding data as a function of time clearly show a rapid formation of folded monomers, followed by a slower step leading to folded dimers. Therefore, according to this result, we believe that the preferential unfolding/refolding mechanism is one in which dimer dissociation occurs before unfolding rather than the reverse.     

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).     

The equilibrium unfolding intermediate observed at pH 4 and its relationship with the kinetic folding intermediates in green fluorescent protein

Folding mechanisms of a variant of green fluorescent protein (F99S/M153T/V163A) were investigated by a wide variety of spectroscopic techniques. Equilibrium measurements on acid-induced denaturation of the protein monitored by chromophore and tryptophan fluorescence and small-angle X-ray scattering revealed that this protein accumulates at least two equilibrium intermediates, a native-like intermediate and an unfolding intermediate, the latter of which exhibits the characteristics of the molten globule state under moderately denaturing conditions at pH 4. To elucidate the role of the equilibrium unfolding intermediate in folding, a series of kinetic refolding experiments with various combinations of initial and final pH values, including pH 7.5 (the native condition), pH 4.0 (the moderately denaturing condition where the unfolding intermediate is accumulated), and pH 2.0 (the acid-denaturing condition) were carried out by monitoring chromophore and tryptophan fluorescence. Kinetic on-pathway intermediates were accumulated during the folding on the refolding reaction from pH 2.0 to 7.5. However, the signal change corresponding to the conversion from the acid-denatured to the kinetic intermediate states was significantly reduced on the refolding reaction from pH 4.0 to pH 7.5, whereas only the signal change corresponding to the above conversion was observed on the refolding reaction from pH 2.0 to pH 4.0. These results indicate that the equilibrium unfolding intermediate is composed of an ensemble of the folding intermediate species accumulated during the folding reaction, and thus support a hierarchical model of protein folding.     

The nature of the rate-limiting steps in the refolding of the cofactor-dependent protein aspartate aminotransferase

The refolding of mitochondrial aspartate aminotransferase (mAAT; EC 2.6.1.1) has been studied following unfolding in 6 m guanidine hydrochloride for different periods of time. Whereas reactivation of equilibrium-unfolded mAAT is sigmoidal, reactivation of the short term unfolded protein displays a double exponential behavior consistent with the presence of fast and slow refolding species. The amplitude of the fast phase decreases with increasing unfolding times (k approximately 0.75 min(-1) at 20 degrees C) and becomes undetectable at equilibrium unfolding. According to hydrogen exchange and stopped-flow intrinsic fluorescence data, unfolding of mAAT appears to be complete in less than 10 s, but hydrolysis of the Schiff base linking the coenzyme pyridoxal 5'-phosphate (PLP) to the polypeptide is much slower (k approximately 0.08 min(-1)). This implies the existence in short term unfolded samples of unfolded species with PLP still attached. However, since the disappearance of the fast refolding phase is about 10-fold faster than the release of PLP, the fast refolding phase does not correspond to folding of the coenzyme-containing molecules. The fast refolding phase disappears more rapidly in the pyridoxamine and apoenzyme forms of mAAT, both of which lack covalently attached cofactor. Thus, bound PLP increases the kinetic stability of the fast refolding unfolding intermediates. Conversion between fast and slow folding forms also takes place in an early folding intermediate. The presence of cyclophilin has no effect on the reactivation of either equilibrium or short term unfolded mAAT. These results suggest that proline isomerization may not be the only factor determining the slow refolding of this cofactor-dependent protein.     

荧光研究表明去辅基天青蛋白突变体M121L至少存在两种不同的构象

以往对绿脓杆菌去辅基天青蛋白变性机制的研究认为它经历了一个复杂的反应过程,相比之下,锌离子替代的天青蛋白的变性符合简单的二态模型。以脲为变性剂对去辅基天青蛋白突变体M121L的变性过程进行了研究。结果表明,虽然稳态条件下突变体的变性／复性符合二态模型,但其动力学过程复杂,并可用溶液中存在着两种可以相互转化的构象的变性／复性来解释。天然态N1去折叠的速度快,其重折叠的速度也快,N1的折叠机制可用存在着折叠途径上的快速折叠中间体模型来描述;天然态N2的去折叠速度慢,其重折叠主要是首先生成天然态N1,然后再缓慢地转化成N2。添加Zn^2 能够把两种构象整合成一种构象,相应地,Zn^2 替代的天青蛋白突变体的变性过程简化为单指数过程。对该突变体的研究加深了对天青蛋白去折叠机制的理解。     

Folding mechanism of an ankyrin repeat protein: scaffold and active site formation of human CDK inhibitor p19(INK4d)

The p19(INK4d) protein consists of five ankyrin repeats (ANK) and controls the human cell cycle by inhibiting the cyclin D-dependent kinases (CDK) 4 and 6. We investigated the folding of p19(INK4d) by urea-induced unfolding transitions, kinetic analyses of unfolding and refolding, including double-mixing experiments and a special assay for folding intermediates. Folding is a sequential two-step reaction via a hyperfluorescent on-pathway intermediate. This intermediate is present under all conditions, during unfolding, refolding and at equilibrium. The folding mechanism was confirmed by a quantitative global fit of a consistent set of equilibrium and kinetic data revealing the thermodynamics and intrinsic folding rates of the different states. Surprisingly, the N<-->I transition is much faster compared to the I<-->U transition. The urea-dependence of the intrinsic folding rates causes population of the intermediate at equilibrium close to the transition midpoint. NMR detected hydrogen/deuterium exchange and the analysis of truncated variants showed that the C-terminal repeats ANK3-5 are already folded in the on-pathway intermediate, whereas the N-terminal repeats 1 and 2 are not folded. We suggest that during refolding, repeats ANK3-ANK5 first form the scaffold for the subsequent assembly of repeats ANK1 and ANK2. The binding function of p19(INK4d) resides in the latter repeats. We propose that the graded stability and the facile unfolding of repeats 1 and 2 is a prerequisite for the down-regulation of the inhibitory activity of p19(INK4d) during the cell-cycle.     

Investigation of the folding pathway of the TEM-1 β-lactamase

The TEM-1 β-lactamase is a globular protein containing 12 proline residues. The folding mechanism of this enzyme was investigated by kinetic and equilibrium experiments with the help of fluorescence spectroscopy and circular dichroism. The equilibrium denaturation of the protein induced by guanidine hydrochloride occurs in two discrete steps, indicating the existence of a thermodynamically stable intermediate state. Thisstate is 5.2 ± 0.4 kcal/mol less stable than the native conformation and 5.7 ± 0.2 kcal/mol more stable than the fully denaturedprotein. This intermediate state exhibits a high content of native secondary structure elements but is devoid of specific tertiary organization; its relation to the “molten globule” is discussed. Refolding kinetic experimentsrevealed the existence of a transient intermediate conformation between thethermodynamically stable intermediate and the native protein. This transient intermediate appears rapidly during the folding reaction. It exhibits a secondary structure content very similar to that of the native protein and has also recovered a significant amount of tertiary organisation. The final refolding step of the TEM-1 β-lactamase, leading to the native enzyme, is dominated by two major slow kinetic phases which probablyreflect a very complex process kinetically limited by proline cis/transisomerization. © 1995 Wiley-Liss, Inc.     

Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: a comparative study of the folding reactions of alpha-lactalbumin and lysozyme

The refolding kinetics of alpha-lactalbumin at different concentrations of guanidine hydrochloride have been investigated by means of kinetic circular dichroism and stopped-flow absorption measurements. The refolding reaction consists of at least two stages, the instantaneous accumulation of the transient intermediate that has peptide secondary structure and the subsequent slow process associated with formation of tertiary structure. The transient intermediate is compared with the well-characterized equilibrium intermediate observed during the denaturant-induced unfolding. Stabilities of the secondary structures against the denaturant, affinities for Ca2+, and tryptophan absorption properties of the transient and equilibrium intermediates were investigated. In all of these respects, the transient intermediate is identical with the equilibrium one, demonstrating the validity of the use of the equilibrium intermediate as a model of the folding intermediate. Essentially the same transient intermediate was also detected in the folding of lysozyme, the protein known to be homologous to alpha-lactalbumin but whose equilibrium unfolding is represented as a two-state reaction. The stability and cooperativity of the secondary structure of the intermediate of lysozyme are compared with those of alpha-lactalbumin. The results show that the protein folding occurring via the intermediate is not limited to the proteins that show equilibrium intermediates. Although the unfolding equilibria of most proteins are well approximated as a two-state reaction, the two-state hypothesis may not be applicable to the folding reaction under the native condition. Two models of protein folding, intermediate-controlled folding model and multiple-pathway folding model, which are different in view of the role of the intermediate in determining the pathway of folding, are also discussed.     

Opposing effects of phosphoenolpyruvate and pyruvate with Mg(2+) on the conformational stability and dimerization of phosphotransferase enzyme I from Escherichia coli

The activity of enzyme I (EI), the first protein in the bacterial PEP:sugar phosphotransferase system, is regulated by a monomer-dimer equilibrium where a Mg(2+)-dependent autophosphorylation by PEP requires the homodimer. Using inactive EI(H189A), in which alanine is substituted for the active-site His189, substrate-binding effects can be separated from those of phosphorylation. Whereas 1 mM PEP (with 2 mM Mg(2+)) strongly promotes dimerization of EI(H189A) at pH 7.5 and 20 degrees C, 5 mM pyruvate (with 2 mM Mg(2+)) has the opposite effect. A correlation between the coupling of N- and C-terminal domain unfolding, measured by differential scanning calorimetry, and the dimerization constant for EI, determined by sedimentation equilibrium, is observed. That is, when the coupling between N- and C-terminal domain unfolding produced by 0.2 or 1.0 mM PEP and 2 mM Mg(2+) is inhibited by 5 mM pyruvate, the dimerization constant for EI(H189A) decreases from > 10(8) to < 5 x 10(5) or 3 x 10(7) M(-1), respectively. Incubation of the wild-type, dephospho-enzyme I with the transition-state analog phosphonopyruvate and 2 mM Mg(2+) also increases domain coupling and the dimerization constant approximately 42-fold. With 2 mM Mg(2+) at 15-25 degrees C and pH 7.5, PEP has been found to bind to one site/monomer of EI(H189A) with K(A)' approximately 10(6) M(-1) (deltaG' = -8.05 +/- 0.05 kcal/mole and deltaH = +3.9 kcal/mole at 20 degrees C); deltaC(p) = -0.33 kcal K(-1) mole(-1). The binding of PEP to EI(H189A) is synergistic with that of Mg(2+). Thus, physiological concentrations of PEP and Mg(2+) increase, whereas pyruvate and Mg(2+) decrease the amount of dimeric, active, dephospho-enzyme I.     

Effects of sulphate and urea on the stability and reversible unfolding of beta-lactamase from Staphylococcus aureus. Implications for the folding pathway of beta-lactamase

The reversible denaturation by urea of beta-lactamase from Staphylococcus aureus was followed in the presence and absence of ammonium sulphate by circular dichroism studies, difference absorption spectroscopy and measurement of enzyme activity. The multiple unfolding and refolding transitions demonstrate the existence of a thermodynamically stable state of intermediate conformation in equilibrium with the native (N) and fully unfolded (U) states. Its physical properties show that it is identical to the state H found on denaturation by guanidinium chloride. State H is 10.1 (+/-1.5) kJ mol-1 less stable than the native state and 10.1 (+/-1.6) kJ mol-1 more stable than the unfolded state. Ammonium sulphate shifts both the N in equilibrium H and H in equilibrium U transitions to concentrations of urea higher by 5.3 M per mole of sulphate. It has markedly different effects on the thermodynamic stabilities of states N and H, making delta G'N-H, O and delta G'H-U, O more negative by 41 kJ mol and 20 kJ mole, respectively, per mole of ammonium sulphate. The change in equilibrium constant for the N-H transition is reflected almost exclusively in a dramatic change of the unfolding rate constant, which is decreased by a factor of 10(11) on addition of 1.4 M-sulphate. The presence of the substrate benzyl penicillin has little effect on the equilibria or kinetics of the N-H transition. The results are discussed in terms of the nature of the N-H transition and of the ordering of intermediate states on the folding pathway.     

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.     

Folding and stability of trp aporepressor from Escherichia coli

Equilibrium and kinetic studies of the urea-induced unfolding of trp aporepressor from Escherichia coli were performed to probe the folding mechanism of this intertwined, dimeric protein. The equilibrium unfolding transitions at pH 7.6 and 25 degrees C monitored by difference absorbance, fluorescence, and circular dichroism spectroscopy are coincident within experimental error. All three transitions are well described by a two-state model involving the native dimer and the unfolded monomer; the free energy of folding in the absence of denaturant and under standard-state conditions is estimated to be 23.3 +/- 0.9 kcal/mol of dimer. The midpoint of the equilibrium unfolding transition increases with increasing protein concentration in the manner expected from the law of mass action for the two-state model. We find no evidence for stable folding intermediates. Kinetic studies reveal that unfolding is governed by a single first-order reaction whose relaxation time decreases exponentially with increasing urea concentration and also decreases with increasing protein concentration in the transition zone. Refolding involves at least three phases that depend on both the protein concentration and the final urea concentration in a complex manner. The relaxation time of the slowest of these refolding phases is identical with that for the single phase in unfolding in the transition zone, consistent with the results expected for a reaction that is kinetically reversible. The two faster refolding phases are presumed to arise from slow isomerization reactions in the unfolded form and reflect parallel folding channels.     

Ferricytochrome c. Refolding and the methionine 80-sulfur-iron linkage

The refolding of urea-denatured horse heart ferricytochrome c in the presence of imidazole, 0.5 M, pH 7.0, has been examined using stopped-flow and equilibrium measurements at 407.5 nm. Thermodynamically, imidazole-cytochrome c folds and unfolds via a single transition with [urea]1/2 of 5.9 M. Kinetically, the refolding is a triphasic process: (i) a slow, urea-independent phase, time constant of 22 +/- 6 s, and an amplitude of 10-13%; (ii) an intermediate reaction, with a slightly positive urea-dependent rate constant, average time constant of 150 ms; and (iii) a fast phase with negative urea dependence of the rate constant from 4-6 M urea and positive dependence above the 6 M concentration, with the largest time constant, 25 +/- 6 ms, at 5.8 M urea, the midpoint of the transition. The amplitudes of the intermediate and the fast phases exhibit inverse dependence on the final urea concentrations, favoring the intermediate form at higher concentrations, while maintaining an almost constant sum of the two amplitudes throughout the range. The temperature dependence of the three apparent rate constants for the refolding from denatured base-line to midpoint of the transition, 9 to 6.03 M urea, yields linear Arrhenius plots with activation energies of 14, 19, and 23 +/- 3 kcal/mol for the slow, intermediate, and rapid reactions, respectively. These findings show that the slow reaction, time constant in decaseconds , does not require, directly or indirectly, the coordination of Met-80-S to heme iron. The formation of this linkage during the folding of the urea-denatured protein in the absence of extrinsic ligand, however, does alter the course of the refolding process. From a comparison of the proposed mechanisms and of the kinetic parameters for the folding of urea-denatured and of guanidine hydrochloride-denatured ferricytochrome c, it has been suggested that the two systems are distinct in detail, although both systems exhibit the slow, decasecond process.     

Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape

The role of native contact topology in the folding of a TIM barrel model based on the alpha-subunit of tryptophan synthase (alphaTS) from Salmonella typhimurium (Protein Data Bank structure 1BKS) was studied using both equilibrium and kinetic simulations. Equilibrium simulations of alphaTS reveal the population of two intermediate ensembles, I1 and I2, during unfolding/refolding at the folding temperature, Tf = 335 K. Equilibrium intermediate I1 demonstrates discrete structure in regions alpha0-beta6 whereas intermediate I2 is a loose ensemble of states with N-terminal structure varying from at least beta1-beta3 (denoted I2A) to alpha0-beta4 at most (denoted I2B). The structures of I1 and I2 match well with the two intermediate states detected in equilibrium folding experiments of Escherichia coli alphaTS. Kinetic folding simulations of alphaTS reveal the sequential population of four intermediate ensembles, I120Q, I200Q, I300Q, and I360Q, during refolding. Kinetic intermediates I120Q, I200Q, and I300Q are highly similar to equilibrium alphaTS intermediates I2A, I2B, and I1, respectively, consistent with kinetic experiments on alphaTS from E. coli. A small population (approximately 10%) of kinetic trajectories are trapped in the I120Q intermediate ensemble and require a slow and complete unfolding step to properly refold. Both the on-pathway and off-pathway I120Q intermediates show structure in beta1-beta3, which is also strikingly consistent with kinetic folding experiments of alphaTS. In the off-pathway intermediate I(120Q), helix alpha2 is wrapped in a nonnative chiral arrangement around strand beta3, sterically preventing the subsequent folding step between beta3 and beta4. These results demonstrate the success of combining kinetic and equilibrium simulations of minimalist protein models to explore TIM barrel folding and the folding of other large proteins.