Folding kinetics of staphylococcal nuclease studied by tryptophan engineering and rapid mixing methods

Folding of Staphylococcal Nuclease Induced by Binding of Chemically Modified Substrate Analogues Sheds Light on Mechanisms of Coupled Folding/Binding Reactions

Several proteins have been shown to undergo a shift in the mechanism of ligand binding-induced folding from conformational selection (CS; folding precedes binding) to induced fit (IF; binding precedes folding) with increasing ligand concentration. In previous studies of the coupled folding/binding reaction of staphylococcal nuclease (SNase) in the presence of a substrate analogue, adenosine-3',5'-diphosphate (prAp), we found that the two phosphate groups make important energetic contributions toward stabilizing its complex with the native protein as well as transient conformational states encountered at high ligand concentrations favoring IF. However, the structural contributions of each phosphate group during the reaction remain unclear. To address this question, we relied on fluorescence, nuclear magnetic resonance (NMR), absorption, and isothermal titration calorimetry to study the effects of deletion of the phosphate groups of prAp on the kinetics of ligand-induced folding, using a strategy analogous to mutational ϕ-value analysis to interpret the results. Kinetic measurements over a wide range of ligand concentrations, together with structural characterization of a transient protein-ligand encounter complex using 2D NMR, indicated that, at high ligand concentrations favoring IF, (i) the 5'-phosphate group interacts weakly with denatured SNase during early stages of the reaction, resulting in loose docking of the two domains of SNase, and (ii) the 3'-phosphate group engages in some specific contacts with the polypeptide in the transition state prior to formation of the native SNase-prAp complex.

Highly Sensitive Tryptophan Fluorescence Probe for detecting Rhythmic Conformational changes of KaiC in the Cyanobacterial Circadian Clock System

KaiC, a core protein of the cyanobacterial circadian clock, consists of an N-terminal CI domain and a C-terminal CII domain, and assembles into a double-ring hexamer upon binding with ATP. KaiC rhythmically phosphorylates and dephosphorylates its own two adjacent residues Ser431 and Thr432 at the CII domain with a period of approximately 24h through assembly and disassembly with the other clock proteins, KaiA and/or KaiB. In this study, to understand how KaiC alters its conformation as the source of circadian rhythm, we investigated structural changes of an inner-radius side of the CII ring using time-resolved Trp fluorescence spectroscopy. A KaiC mutant harboring a Trp fluorescence probe at a position of 419 exhibited a robust circadian rhythm with little temperature sensitivity in the presence of KaiA and KaiB. Our fluorescence observations show a remarkable environmental change at the inner-radius side of the CII ring during circadian oscillation. Crystallographic analysis revealed that a side chain of Trp at the position of 419 was oriented toward a region undergoing a helix-coil transition, which is considered to be a key event to allosterically regulate the CI ring that plays a crucial role in determining the cycle period. The present study provides a dynamical insight into how KaiC generates circadian oscillation.

Advances in Mixer Design and Detection Methods for Kinetics Studies of Macromolecular Folding and Binding on the Microsecond Time Scale

Many important biological processes such as protein folding and ligand binding are too fast to be fully resolved using conventional stopped-flow techniques. Although advances in mixer design and detection methods have provided access to the microsecond time regime, there is room for improvement in terms of temporal resolution and sensitivity. To address this need, we developed a continuous-flow mixing instrument with a dead time of 12 to 27 µs (depending on solution viscosity) and enhanced sensitivity, sufficient for monitoring tryptophan or tyrosine fluorescence changes at fluorophore concentrations as low as 1 µM. Relying on commercially available laser microfabrication services, we obtained an integrated mixer/flow-cell assembly on a quartz chip, based on a cross-channel configuration with channel dimensions and geometry designed to minimize backpressure. By gradually increasing the width of the observation channel downstream from the mixing region, we are able to monitor a reaction progress time window ranging from ~10 µs out to ~3 ms. By combining a solid-state UV laser with a Galvano-mirror scanning strategy, we achieved highly efficient and uniform fluorescence excitation along the flow channel. Examples of applications, including refolding of acid-denatured cytochrome c triggered by a pH jump and binding of a peptide ligand to a PDZ domain, demonstrate the capability of the technique to resolve fluorescence changes down to the 10 µs time regime on modest amounts of reagents.

Energetics and kinetics of substrate analog-coupled staphylococcal nuclease folding revealed by a statistical mechanical approach

Protein conformational changes associated with ligand binding, especially those involving intrinsically disordered proteins, are mediated by tightly coupled intra- and intermolecular events. Such reactions are often discussed in terms of two limiting kinetic mechanisms, conformational selection (CS), where folding precedes binding, and induced fit (IF), where binding precedes folding. It has been shown that coupled folding/binding reactions can proceed along both CS and IF pathways with the flux ratio depending on conditions such as ligand concentration. However, the structural and energetic basis of such complex reactions remains poorly understood. Therefore, we used experimental, theoretical, and computational approaches to explore structural and energetic aspects of the coupled-folding/binding reaction of staphylococcal nuclease in the presence of the substrate analog adenosine-3',5'-diphosphate. Optically monitored equilibrium and kinetic data, combined with a statistical mechanical model, gave deeper insight into the relative importance of specific and Coulombic protein-ligand interactions in governing the reaction mechanism. We also investigated structural aspects of the reaction at the residue level using NMR and all-atom replica-permutation molecular dynamics simulations. Both approaches yielded clear evidence for accumulation of a transient protein-ligand encounter complex early in the reaction under IF-dominant conditions. Quantitative analysis of the equilibrium/kinetic folding revealed that the ligand-dependent CS-to-IF shift resulted from stabilization of the compact transition state primarily by weakly ligand-dependent Coulombic interactions with smaller contributions from specific binding energies. At a more macroscopic level, the CS-to-IF shift was represented as a displacement of the reaction "route" on the free energy surface, which was consistent with a flux analysis.

Lessons from pressure denaturation of proteins

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

Evolutionarily Conserved Pattern of Interactions in a Protein Revealed by Local Thermal Expansion Properties

The way in which the network of intramolecular interactions determines the cooperative folding and conformational dynamics of a protein remains poorly understood. High-pressure NMR spectroscopy is uniquely suited to examine this problem because it combines the site-specific resolution of the NMR experiments with the local character of pressure perturbations. Here we report on the temperature dependence of the site-specific volumetric properties of various forms of staphylococcal nuclease (SNase), including three variants with engineered internal cavities, as measured with high-pressure NMR spectroscopy. The strong temperature dependence of pressure-induced unfolding arises from poorly understood differences in thermal expansion between the folded and unfolded states. A significant inverse correlation was observed between the global thermal expansion of the folded proteins and the number of strong intramolecular hydrogen bonds, as determined by the temperature coefficient of the backbone amide chemical shifts. Comparison of the identity of these strong H-bonds with the co-evolution of pairs of residues in the SNase protein family suggests that the architecture of the interactions detected in the NMR experiments could be linked to a functional aspect of the protein. Moreover, the temperature dependence of the residue-specific volume changes of unfolding yielded residue-specific differences in expansivity and revealed how mutations impact intramolecular interaction patterns. These results show that intramolecular interactions in the folded states of proteins impose constraints against thermal expansion and that, hence, knowledge of site-specific thermal expansivity offers insight into the patterns of strong intramolecular interactions and other local determinants of protein stability, cooperativity, and potentially also of function.

Experimental evolution of a green fluorescent protein composed of 19 unique amino acids without tryptophan

At some stage of evolution, genes of organisms may have encoded proteins that were synthesized using fewer than 20 unique amino acids. Similar to evolution of the natural 19-amino-acid proteins GroEL/ES, proteins composed of 19 unique amino acids would have been able to evolve by accumulating beneficial mutations within the 19-amino-acid repertoire encoded in an ancestral genetic code. Because Trp is thought to be the last amino acid included in the canonical 20-amino-acid repertoire, this late stage of protein evolution could be mimicked by experimental evolution of 19-amino-acid proteins without tryptophan (Trp). To further understand the evolution of proteins, we tried to mimic the evolution of a 19-amino-acid protein involving the accumulation of beneficial mutations using directed evolution by random mutagenesis on the whole targeted gene sequence. We created active 19-amino-acid green fluorescent proteins (GFPs) without Trp from a poorly fluorescent 19-amino-acid mutant, S1-W57F, by using directed evolution with two rounds of mutagenesis and selection. The N105I and S205T mutations showed beneficial effects on the S1-W57F mutant. When these two mutations were combined on S1-W57F, we observed an additive effect on the fluorescence intensity. In contrast, these mutations showed no clear improvement individually or in combination on GFPS1, which is the parental GFP mutant composed of 20 amino acids. Our results provide an additional example for the experimental evolution of 19-amino-acid proteins without Trp, and would help understand the mechanisms underlying the evolution of 19-amino-acid proteins. (236 words).

Nonuniform chain collapse during early stages of staphylococcal nuclease folding detected by fluorescence resonance energy transfer and ultrarapid mixing methods

The development of tertiary structure during folding of staphylococcal nuclease (SNase) was studied by time-resolved fluorescence resonance energy transfer measured using continuous- and stopped-flow techniques. Variants of this two-domain protein containing intradomain and interdomain fluorescence donor/acceptor pairs (Trp and Cys-linked fluorophore or quencher) were prepared to probe the intradomain and interdomain structural evolution accompanying SNase folding. The intra-domain donor/acceptor pairs are within the β-barrel domain (Trp27/Cys64 and Trp27/Cys97) and the interdomain pair is between the α-helical domain and the β-barrel domain (Trp140/Cys64). Time-resolved energy transfer efficiency accompanying folding and unfolding at different urea concentrations was measured over a time range from 30 μs to ≈ 10 s. Information on average donor/acceptor distances at different stages of the folding process was obtained by using a quantitative kinetic modeling approach. The average distance for the donor/acceptor pairs in the β-barrel domain decreases to nearly native values whereas that of the interdomain donor/acceptor pairs remains unchanged in the earliest intermediate (<500 μs of refolding). This indicates a rapid nonuniform collapse resulting in an ensemble of heterogeneous conformations in which the central region of the β-barrel domain is well developed while the C-terminal α-helical domain remains disordered. The distance between Trp140 and Cys64 decreases to native values on the 100-ms time scale, indicating that the α-helical domain docks onto the preformed β-barrel at a late stage of the folding. In addition, the unfolded state is found to be more compact under native conditions, suggesting that changes in solvent conditions may induce a nonspecific hydrophobic collapse.

Remodeling of the folding free energy landscape of staphylococcal nuclease by cavity-creating mutations

The folding of staphylococcal nuclease (SNase) is known to proceed via a major intermediate in which the central OB subdomain is folded and the C-terminal helical subdomain is disordered. To identify the structural and energetic determinants of this folding free energy landscape, we have examined in detail, using high-pressure NMR, the consequences of cavity creating mutations in each of the two subdomains of an ultrastable SNase, Δ+PHS. The stabilizing mutations of Δ+PHS enhanced the population of the major folding intermediate. Cavity creation in two different regions of the Δ+PHS reference protein, despite equivalent effects on global stability, had very distinct consequences on the complexity of the folding free energy landscape. The L125A substitution in the C-terminal helix of Δ+PHS slightly suppressed the major intermediate and promoted an additional excited state involving disorder in the N-terminus, but otherwise decreased landscape heterogeneity with respect to the Δ+PHS background protein. The I92A substitution, located in the hydrophobic OB-fold core, had a much more profound effect, resulting in a significant increase in the number of intermediate states and implicating the entire protein structure. Denaturant (GuHCl) had very subtle and specific effects on the landscape, suppressing some states and favoring others, depending upon the mutational context. These results demonstrate that disrupting interactions in a region of the protein with highly cooperative, unfrustrated folding has very profound effects on the roughness of the folding landscape, whereas the effects are less pronounced for an energetically equivalent substitution in an already frustrated region.

Cavities determine the pressure unfolding of proteins

It has been known for nearly 100 years that pressure unfolds proteins, yet the physical basis of this effect is not understood. Unfolding by pressure implies that the molar volume of the unfolded state of a protein is smaller than that of the folded state. This decrease in volume has been proposed to arise from differences between the density of bulk water and water associated with the protein, from pressure-dependent changes in the structure of bulk water, from the loss of internal cavities in the folded states of proteins, or from some combination of these three factors. Here, using 10 cavity-containing variants of staphylococcal nuclease, we demonstrate that pressure unfolds proteins primarily as a result of cavities that are present in the folded state and absent in the unfolded one. High-pressure NMR spectroscopy and simulations constrained by the NMR data were used to describe structural and energetic details of the folding landscape of staphylococcal nuclease that are usually inaccessible with existing experimental approaches using harsher denaturants. Besides solving a 100-year-old conundrum concerning the detailed structural origins of pressure unfolding of proteins, these studies illustrate the promise of pressure perturbation as a unique tool for examining the roles of packing, conformational fluctuations, and water penetration as determinants of solution properties of proteins, and for detecting folding intermediates and other structural details of protein-folding landscapes that are invisible to standard experimental approaches.

The role of the local environment of engineered Tyr to Trp substitutions for probing the denaturation mechanism of FIS

Factor for inversion stimulation (FIS), a 98-residue homodimeric protein, does not contain tryptophan (Trp) residues but has four tyrosine (Tyr) residues located at positions 38, 51, 69, and 95. The equilibrium denaturation of a P61A mutant of FIS appears to occur via a three-state (N(2) ⇆ I(2) ⇆ 2U) process involving a dimeric intermediate (I(2)). Although it was suggested that this intermediate had a denatured C-terminus, direct evidence was lacking. Therefore, three FIS double mutants, P61A/Y38W, P61A/Y69W, and P61A/Y95W were made, and their denaturation was monitored by circular dichroism and Trp fluorescence. Surprisingly, the P61A/Y38W mutant best monitored the N(2) ⇆ I(2) transition, even though Trp38 is buried within the dimer removed from the C-terminus. In addition, although Trp69 is located on the protein surface, the P61A/Y69W FIS mutant exhibited clearly biphasic denaturation curves. In contrast, P61A/Y95W FIS was the least effective in decoupling the two transitions, exhibiting a monophasic fluorescence transition with modest concentration-dependence. When considering the local environment of the Trp residues and the effect of each mutation on protein stability, these results not only confirm that P61A FIS denatures via a dimeric intermediate involving a disrupted C-terminus but also suggest the occurrence of conformational changes near Tyr38. Thus, the P61A mutation appears to compromise the denaturation cooperativity of FIS by failing to propagate stability to those regions involved mostly in intramolecular interactions. Furthermore, our results highlight the challenge of anticipating the optimal location to engineer a Trp residue for investigating the denaturation mechanism of even small proteins.

Probing the folding intermediate of Bacillus subtilis RNase P protein by nuclear magnetic resonance

Protein folding intermediates are often imperative for overall folding processes and consequent biological functions. However, the low population and transient nature of the intermediate states often hinder their biochemical and biophysical characterization. Previous studies have demonstrated that Bacillus subtilis ribonuclease P protein (P protein) is conformationally heterogeneous and folds with multiphasic kinetics, indicating the presence of an equilibrium and kinetic intermediate in its folding mechanism. In this study, nuclear magnetic resonance (NMR) spectroscopy was used to study the ensemble corresponding to this intermediate (I). The results indicate that the N-terminal and C-terminal helical regions are mostly unfolded in I. 1H−15N heteronuclear single-quantum coherence NMR spectra collected as a function of pH suggest that the protonation of His 22 may play a major role in the energetics of the equilibria among the unfolded, intermediate, and folded state ensembles of P protein. NMR paramagnetic relaxation enhancement experiments were also used to locate the small anion binding sites in both the intermediate and folded ensembles. The results for the folded protein are consistent with the previously modeled binding regions. These structural insights suggest a possible role for I in the RNase P holoenzyme assembly process.

Nonlocal interactions are responsible for tertiary structure formation in staphylococcal nuclease

Rapid molecular collapse mediated by nonlocal interactions is believed to be a crucial event for protein folding. To investigate the role of nonlocal interactions in tertiary structure formation, we performed a nonlocal interaction substitution mutation analysis on staphylococcal nuclease (SNase). Y54 and I139 of wild-type (WT) SNase and Delta140-149 were substituted by cysteine to form intramolecular disulfide bonds, respectively called WT-SS and Delta140-149-SS. Under physiological conditions, the reduced form of Delta140-149-SS appears to assume a denatured structure; in contrast, the oxidized form of Delta140-149-SS forms a native-like structure. From this result, we conclude that the C-terminal region participates in a nonlocal interaction that is indispensable for the native structure. Although the oxidized form of WT-SS assumes a more compact denatured structure under acidic conditions than the WT, the kinetic measurements reveal that the refolding reactions of both the reduced and oxidized forms of WT-SS are similar to those of the WT, suggesting that an intact nonlocal interaction is established within the dead time (22 ms). On the basis of these results, we propose that the native nonlocal contact established at the early stage of the folding process facilitates further secondary structure formation.

A natural missing link between activated and downhill protein folding scenarios

We propose protein PTB1 : 4W as a good candidate for engineering into a downhill folder. PTB1 : 4W has a probe-dependent thermal unfolding curve and sub-millisecond T-jump relaxation kinetics on more than one time scale. Its refolding rate in denaturant is a non-linear function of denaturant concentration (curved chevron plot). Yet at high denaturant concentration its unfolding is probe-independent, and the folding kinetics can be fitted to a single exponential decay. The domain appears to fold via a mechanism between downhill folding and activated folding over several small barriers, and when denaturant is added, one of these barriers greatly increases and simplifies the observed folding to apparent two-state kinetics. We predict the simplest free energy function consistent with the thermal denaturation and kinetics experiments by using the singular value Smoluchowski dynamics (SVSD) model. PTB1 : 4W is a natural 'missing link' between downhill and activated folding. We suggest mutations that could move the protein into the downhill folding limit.

Slow formation of aggregation-resistant beta-sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding-related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large beta-helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all-beta-sheet protein allows detailed analysis of the formation of beta-sheet structure in larger proteins. Using a combination of fluorescence and far-UV circular dichroism spectroscopy, we show that the pertactin beta-helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and beta-sheet-rich topology, pertactin refolding is reversible and not complicated by off-pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate-limiting step. Furthermore, site-specific labeling experiments indicate that the beta-helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, beta-sheet-rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, beta-sheet-rich refolding intermediates. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Equilibrium and kinetic approaches for studying oligomeric protein folding

This chapter describes the approaches and considerations necessary for extension of current protein folding methods to the equilibrium and kinetic reactions of oligomeric proteins, using dimers as the primary example. Spectroscopic and transport methods to monitor folding and unfolding transitions are summarized. The data collection and analyses to determine protein stability and kinetic folding mechanisms are discussed in the context of the additional dimension of complexity that arises in higher order folding processes, compared to first order monomeric proteins. As a case study to illustrate the data analysis process, equilibrium, and kinetic data are presented for SmtB, a homodimeric DNA-binding protein from Synechococcus PCC7942.

Comparison of Two ESI MS Based H/D Exchange Methods for Extracting Protein Folding Energies

In this report, the model proteins staphylococcal nuclease and ubiquitin were used to test the applicability of two new hydrogen/deuterium exchange (HX) electrospray ionization mass spectrometry (ESI-MS) methods for estimating protein folding energies. Both methods use the H/D exchange of globally protected amide protons (amide protons which are buried in the hydrophobic core) to elucidate protein folding energies. One method is a kinetic-based method and the other is equilibrium-based. The first method, the HX ESI-MS kinetic-based approach is conceptually identical to SUPREX (stability of unpurified proteins from rates of H/D exchange) method but is based on ESI-MS rather than MALDI-MS (matrix assisted laser desorption mass spectrometry). This method employs the time-dependence of H/D exchange using various denaturant concentrations to extract folding energies. Like SUPREX, this approach requires the assumption of EX2 exchange kinetics. The second method, which we call a protein equilibrium population snapshot (PEPS) by HX ESI-MS uses data collected only for a single time point (usually the shortest possible) to obtain a snapshot of the open and closed populations of the protein. The PEPS approach requires few assumptions in the derivation of the equations used for calculation of the folding energies. The extraction of folding energies from mass spectral data is simple and straightforward. The PEPS method is applicable for proteins that follow either EX1 or EX2 HX mechanisms. In our experiments the kinetic-based method produced less accurate ΔG(H(2)O) and m(GdHCl) values for wild-type staphylococcal nuclease and mutants undergoing H/D exchange by EX1, as would be expected. Better results were obtained for ubiquitin which undergoes HX by an EX2 mechanism. Using the PEPS method we obtained ΔG(H(2)O) and m(GdHCl) values that were in good agreement with literature values for both staphylococcal nuclease (EX1) and ubiquitin (EX2). We also show that the observation of straight lines in linear extrapolation method (LEM) plots is not a reliable indicator of the validity of the data obtained using the LEM approach.

Evaluation of folding co-operativity of a chimeric protein based on the molecular recognition between polyproline ligands and SH3 domains

In previous work, we designed a chimeric protein, named SPCp41, to evaluate the thermodynamics of the interaction between SH3 domains and proline-rich ligands by combining thermal unfolding measurements and mutagenesis. Here, we have investigated the energetic integrity of the chain extension corresponding to the ligand sequence into the native structure, since the opposite will produce changes in the folding mechanism of the SH3 domain that may give rise to undesirable contributions to the thermodynamic parameters. We have analysed the folding-unfolding kinetics under standard conditions (50 mM phosphate pH 7). Kinetic evolutions are well described by a bi-exponential where, on top of the main kinetic phase, a low-populated slower phase appears as a consequence of cis-trans isomerisation of Pro39, as demonstrated by the influence of prolyl isomerases and by mutational analysis. There is also a burst phase possibly due to a productive formation of some helical ensembles. The main evolution, accounting for the true folding kinetics of SPCp41, can be considered as a two-state process, where the folding transition state produces essentially the same picture shown by the circular permutant S19-P20s (the 'nucleus' of the design) and the ligand will dock at the latter stages of the two-state process. Thus, all conclusions argue in favour of the effectiveness of SPCp41 to study energetic, dynamic and structural aspects of SH3-ligand interactions.

Aromatic residues engineered into the beta-turn nucleation site of ubiquitin lead to a complex folding landscape, non-native side-chain interactions, and kinetic traps

The fast folding of small proteins is likely to be the product of evolutionary pressures that balance the search for native-like contacts in the transition state with the minimum number of stable non-native interactions that could lead to partially folded states prone to aggregation and amyloid formation. We have investigated the effects of non-native interactions on the folding landscape of yeast ubiquitin by introducing aromatic substitutions into the beta-turn region of the N-terminal beta-hairpin, using both the native G-bulged type I turn sequence (TXTGK) as well as an engineered 2:2 XNGK type I' turn sequence. The N-terminal beta-hairpin is a recognized folding nucleation site in ubiquitin. The folding kinetics for wt-Ub (TLTGK) and the type I' turn mutant (TNGK) reveal only a weakly populated intermediate, however, substitution with X = Phe or Trp in either context results in a high propensity to form a stable compact intermediate where the initial U-->I collapse is visible as a distinct kinetic phase. The introduction of Trp into either of the two host turn sequences results in either complex multiphase kinetics with the possibility of parallel folding pathways, or formation of a highly compact I-state stabilized by non-native interactions that must unfold before refolding. Sequence substitutions with aromatic residues within a localized beta-turn capable of forming non-native hydrophobic contacts in both the native state and partially folded states has the undesirable consequence that folding is frustrated by the formation of stable compact intermediates that evolutionary pressures at the sequence level may have largely eliminated.

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.