Folding of the alphaII-spectrin SH3 domain under physiological salt conditions

The SH3 domain has often been used as a model for protein folding due to its typical two-state behaviour. However, recent experimental data at low pH as well as molecular dynamic simulations have indicated that the folding process of SH3 probably is more complicated, and may involve intermediate states. Using both kinetic and equilibrium measurements we have obtained evidence that under native-like conditions the folding of the spectrin SH3 domain does not follow a classic two-state behaviour. The curvature we observed in the Chevron plots is a strong indication of a non-linear activation energy relationship due to the presence of high-energy intermediates. In addition, circular dichroism measurements indicated that refolding after thermal denaturation did not follow the same pattern as thermal unfolding but rather implied less cooperativity and that the refolding transition increased with increasing protein concentration. Further, NMR experiments indicated that upon refolding the SH3 domain gave rise to more than one conformation. Therefore, our results suggest that the folding of the SH3 domain of αII-spectrin does not follow a classical two-state process under high-salt conditions and neutral pH. Heterogeneous folding pathways, which can include folding intermediates as well as misfolded intermediates, might give a more reasonable insight into the folding behaviour of the αII-spectrin SH3 domain.     

Determining denaturation midpoints in multiprobe equilibrium protein folding experiments

Multiprobe equilibrium unfolding experiments in the downhill regime (i.e., maximal barrier < 3 RT) can resolve the folding process with atomic resolution [ Munoz ( 2002) Int. J. Quantum Chem. 90, 1522 -1528] . Such information is extracted from hundreds of heterogeneous atomic equilibrium unfolding curves, which are characterized according to their denaturation midpoint (e.g., T m for thermal denaturation). Using statistical methods, we analyze T m accuracy when determined from the extremum of the derivative of the unfolding curve and from two-state fits under different sets of simulated experimental conditions. We develop simple procedures to discriminate between real unfolding heterogeneity at the atomic level and experimental uncertainty in the single T m of conventional two-state folding. We apply these procedures to the recently published multiprobe NMR experiments of BBL [ Sadqi et al. ( 2006) Nature 442, 317 -321 ] and conclude that for the 122 single transition atomic unfolding curves reported for this protein the mean T m accuracy is better than 1.8 K for both methods, compared to the 60 K spread in T m determined experimentally. Importantly, we also find that when the pre- or posttransition baseline is incomplete, the two-state fits systematically drift the estimated T m value toward the center of the experimental range. Therefore, the reported 60 K T m spread in BBL is in fact a lower limit. The derivative method is significantly less sensitive to this problem and thus is a better choice for multiprobe experiments with a broad T m distribution. The results we obtain in this work lay the foundations for the quantitative analysis of future multiprobe unfolding experiments in fast-folding proteins.     

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.     

T-jump infrared study of the folding mechanism of coiled-coil GCN4-p1

Partially folded intermediates have been frequently observed in equilibrium and kinetic protein folding studies. However, folding intermediates that exist at the native side of the rate-limiting step are rather difficult to study because they often evade detection by conventional folding kinetic methods. Here, we demonstrated that a laser-induced temperature-jump method can potentially be used to identify the existence of such post-transition or hidden intermediates. Specifically, we studied two cross-linked variants of GCN4-p1 coiled-coil. The GCN4 leucine zipper has been studied extensively and most of these studies have regarded it as a two-state folder. Our static circular dichroism and infrared data also indicate that the thermal unfolding of these two monomeric coiled-coils can be adequately described by an apparent two-state model. However, their temperature-jump-induced relaxation kinetics exhibit non-monoexponential behavior, dependent upon sequence and temperature. Taken together, our results support a folding mechanism wherein at least one folding intermediate populates behind the main rate-limiting step.     

Folding of a designed simple ankyrin repeat protein

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.     

Detection and characterization of partially unfolded oligomers of the SH3 domain of alpha-spectrin

For the purpose of equilibrium and kinetic folding-unfolding studies, the SH3 domain of alpha-spectrin (spc-SH3) has long been considered a classic two-state folding protein. In this work we have indeed observed that the thermal unfolding curves of spc-SH3 measured at pH 3.0 by differential scanning calorimetry, circular dichroism, and NMR follow apparently the two-state model when each unfolding profile is considered individually. Nevertheless, we have found that protein concentration has a marked effect upon the thermal unfolding profiles. This effect cannot be properly explained in terms of the two-state unfolding model and can only be interpreted in terms of the accumulation of intermediate associated states in equilibrium with the monomeric native and unfolded states. By chemical cross-linking and pulsed-field gradient NMR diffusion experiments we have been able to confirm the existence of associated states formed during spc-SH3 unfolding. A three-state model, in which a dimeric intermediate state is assumed to be significantly populated, provides the simplest interpretation of the whole set of thermal unfolding data and affords a satisfactory explanation for the concentration effects observed. Whereas at low concentrations the population of the associated intermediate state is negligible and the unfolding process consequently takes place in a two-state fashion, at concentrations above approximately 0.5 mM the population of the intermediate state becomes significant at temperatures between 45 degrees C and 80 degrees C and reaches up to 50% at the largest concentration investigated. The thermodynamic properties of the intermediate state implied by this analysis fall in between those of the unfolded state and the native ones, indicating a considerably disordered conformation, which appears to be stabilized by oligomerization.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.     

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.     

Presence of the cofactor speeds up folding of Desulfovibrio desulfuricans flavodoxin

Flavodoxin is an alpha/beta protein with a noncovalently bound flavin-mononucleotide (FMN) cofactor. The apo-protein adopts a structure identical to that of the holo-form, although there is more dynamics in the FMN-binding loops. The equilibrium unfolding processes of Azotobacter vinelandii apo-flavodoxin, and Desulfovibrio desulfuricans ATCC strain 27774 apo- and holo-flavodoxins involve rather stable intermediates. In contrast, we here show that both holo- and apo-forms of flavodoxin from D. desulfuricans ATCC strain 29577 (75% sequence similarity with the strain 27774 protein) unfold in two-state equilibrium processes. Moreover, the FMN cofactor remains bound to the unfolded holo-protein. The folding and unfolding kinetics for holo-flavodoxin exhibit two-state behavior, albeit an additional slower phase is present at very low denaturant concentrations. The extrapolated folding time in water for holo-flavodoxin, approximately 280 microsec, is in excellent agreement with that predicted from the protein's native-state topology. Unlike the holo-protein behavior, the folding and unfolding reactions for apo-flavodoxin are best described by two kinetic phases, with rates differing approximately 15-fold, suggesting the presence of a kinetic intermediate. Both folding phases for apo-flavodoxin are orders of magnitude slower (40- and 530-fold, respectively) than that for the holo-protein. We conclude that polypeptide-cofactor interactions in the unfolded state of D. desulfuricans strain 29577 flavodoxin alter the kinetic-folding path towards two-state and speed up the folding reaction.     

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.     

Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.     

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.     

Modulation of the multistate folding of designed TPR proteins through intrinsic and extrinsic factors

Tetratricopeptide repeats (TPRs) are a class of all alpha-helical repeat proteins that are comprised of 34-aa helix-turn-helix motifs. These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. Several studies on designed TPR proteins have shown that this modular structure is reflected in their folding, that is, modular multistate folding is observed as opposed to two-state folding. Here we show that TPR multistate folding can be suppressed to approximate two-state folding through modulation of intrinsic stability or extrinsic environmental variables. This modulation was investigated by comparing the thermodynamic unfolding under differing buffer regimes of two distinct series of consensus-designed TPR proteins, which possess different intrinsic stabilities. A total of nine proteins of differing sizes and differing consensus TPR motifs were each thermally and chemically denatured and their unfolding monitored using differential scanning calorimetry (DSC) and CD/fluorescence, respectively. Analyses of both the DSC and chemical denaturation data show that reducing the total stability of each protein and repeat units leads to observable two-state unfolding. These data highlight the intimate link between global and intrinsic repeat stability that governs whether folding proceeds by an observably two-state mechanism, or whether partial unfolding yields stable intermediate structures which retain sufficient stability to be populated at equilibrium.     

Role of entropy in protein thermostability: folding kinetics of a hyperthermophilic cold shock protein at high temperatures using 19F NMR

We used (19)F NMR to extend the temperature range accessible to detailed kinetic and equilibrium studies of a hyperthermophilic protein. Employing an optimized incorporation strategy, the small cold shock protein from the bacterium Thermotoga maritima (TmCsp) was labeled with 5-fluorotryptophan. Although chaotropically induced unfolding transitions revealed a significant decrease in the stabilization free energy upon fluorine labeling, the protein's kinetic folding mechanism is conserved. Temperature- and guanidinium chloride-dependent equilibrium unfolding transitions monitored by (19)F NMR agree well with the results from optical spectroscopy, and provide a stringent test of the two-state folding character of TmCsp. Folding and unfolding rate constants at high temperatures were determined from the (19)F NMR spectra close to the midpoint of thermal unfolding by global line shape analysis. In combination with results from stopped-flow experiments at lower temperatures, they show that the folding rate constant of TmCsp and its temperature dependence closely resemble those of its mesophilic homologue from Bacillus subtilis, BsCspB. However, the unfolding rate constant of TmCsp is two orders of magnitude lower over the entire temperature range that was investigated. Consequently, the difference in conformational stability between the two proteins is solely due to the unfolding rate constant over a wide temperature range. A thermodynamic analysis points to an important role of entropic factors in the stabilization of TmCsp relative to its mesophilic homologues.     

Guanidine-induced equilibrium unfolding of a homo-hexameric enzyme 4-oxalocrotonate tautomerase (4-OT)

4-Oxalocrotonate tautomerase (4-OT) is a bacterial enzyme that is comprised of 6 identical 62 amino acid subunits. The 4-OT enzyme is an attractive model system in which to study the interrelationship between protein folding, subunit assembly, and catalytic function. Here we report on the GuHCl-induced equilibrium unfolding properties of wild-type 4-OT using catalytic activity measurements and using far-UV circular dichroism (CD) spectroscopy. We demonstrate that the unfolding of wild-type 4-OT in 50 mM phosphate buffers containing 6 M GuHCl is reversible at pHs 6.0, 7.4, and 8.5; and we find that there is both an enzyme concentration dependence and a pH dependence to the equilibrium unfolding properties of 4-OT. Our data suggests that the GuHCl-induced unfolding of 4-OT in 50 mM phosphate buffer at pH 8.5 can be modeled as a two-state process involving folded hexamer and unfolded monomer. On the basis of this model, we determined a free-energy value for the unfolding of 4-OT at pH 8.5 to be 68.7 +/- 3.2 kcal/mol under standard state conditions (1 M hexamer). In 50 mM phosphate buffers at pHs 6.0 and 7.4, only the catalytic activity denaturation curves are consistent with a two-state folding mechanism. At the lower pHs the far-UV-CD transitions are not well described by a two-state model. Our results at pHs 6.0 and 7.4 suggest that intermediate state(s) are populated in the equilibrium unfolding reaction at these lower pHs and that these intermediate state(s) have some helical content but no measurable catalytic activity.     

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.     

Kinetically robust monomeric protein from a hyperthermophile

Equilibrium and kinetic studies were carried out under denaturation conditions to clarify the energetic features of the high stability of a monomeric protein, ribonuclease HII, from a hyperthermophile, Thermococcus kodakaraensis (Tk-RNase HII). Guanidine hydrochloride (GdnHCl)-induced unfolding and refolding were measured with circular dichroism at 220 nm, and heat-induced denaturation was studied with differential scanning calorimetry. Both GdnHCl- and heat-induced denaturation are very reversible. It was difficult to obtain the equilibrated unfolding curve of Tk-RNase HII below 40 degrees C, because of the remarkably slow unfolding. The two-state unfolding and refolding reactions attained equilibrium at 50 degrees C after 2 weeks. The Gibbs energy change of GdnHCl-induced unfolding (DeltaG(H(2)O)) at 50 degrees C was 43.6 kJ mol(-1). The denaturation temperature in the DSC measurement shifted as a function of the scan rate; the denaturation temperature at a scan rate of 90 degrees C h(-1) was higher than at a scan rate of 5 degrees C h(-1). The unfolding and refolding kinetics of Tk-RNase HII were approximated as a first-order reaction. The ln k(u) and ln k(r) values depended linearly on the denaturant concentration between 10 and 50 degrees C. The DeltaG(H(2)O) value obtained from the rate constant in water using the two-state model at 50 degrees C, 44.5 kJ mol(-1), was coincident with that from the equilibrium study, 43.6 kJ mol(-1), suggesting the two-state folding of Tk-RNase HII. The values for the rate constant in water of the unfolding for Tk-RNase HII were much smaller than those of E. coli RNase HI and Thermus thermophilus RNase HI, which has a denaturation temperature similar to that of Tk-RNase HII. In contrast, little difference was observed in the refolding rates among these proteins. These results indicate that the stabilization mechanism of monomeric protein from a hyperthermophile, Tk-RNase HII, with reversible two-state folding is characterized by remarkably slow unfolding.     

Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9

The helical bacterial immunity proteins Im7 and Im9 have been shown to fold via kinetic mechanisms of differing complexity, despite having 60 % sequence identity. At pH 7.0 and 10 degrees C, Im7 folds in a three-state mechanism involving an on-pathway intermediate, while Im9 folds in an apparent two-state transition. In order to examine the folding mechanisms of these proteins in more detail, the folding kinetics of both Im7 and Im9 (at 10 degrees C in 0.4 M sodium sulphate) have been examined as a function of pH. Kinetic modelling of the folding and unfolding data for Im7 between pH 5.0 and 8.0 shows that the on-pathway intermediate is stabilised by more acidic conditions, whilst the native state is destabilised. The opposing effect of pH on the stability of these states results in a significant population of the intermediate at equilibrium at pH 6.0 and below. At pH 7.0, the folding and unfolding kinetics for Im9 can be fitted adequately by a two-state model, in accord with previous results. However, under acidic conditions there is a clear change of slope in the plot of the logarithm of the folding rate constant versus denaturant concentration, consistent with the population of one or more intermediate(s) early during folding. The kinetic data for Im9 at these pH values can be fitted to a three-state model, where the intermediate ensemble is stabilised and the native state destabilised as the pH is reduced, rationalising previous results that showed that an intermediate is not observed experimentally at pH 7.0. The data suggest that intermediate formation is a general step in immunity protein folding and demonstrate that it is necessary to explore a wide range of refolding conditions in order to show that intermediates do not form in the folding of other small, single-domain proteins.     

The kinetic pathway of folding of barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.