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.     

Role of the B helix in early folding events in apomyoglobin: evidence from site-directed mutagenesis for native-like long range interactions

The folding pathways of four mutants in which bulky hydrophobic residues in the B helix of apomyoglobin (ApoMb) are replaced by alanine (I28A, L29A, I30A, and L32A) have been analyzed using equilibrium and kinetic methods employing NMR, CD, fluorescence and mass spectrometry. Hydrogen exchange pulse-labeling followed by mass spectrometry reveals detectable intermediates in the kinetic folding pathways of each of these mutants. Comparison of the quench-flow data analyzed by NMR for the wild-type protein and the mutants showed that the substitutions I28A, L29A and L32A lead to destabilization of the B helix in the burst phase kinetic intermediate, relative to wild-type apomyoglobin. In contrast, the I30A mutation apparently has a slight stabilizing effect on the B helix in the burst phase intermediate; under weak labeling conditions, residues in the C helix region were also relatively stabilized in the mutant compared to the wild-type protein. This suggests that native-like helix B/helix C packing interactions occur in the folding intermediate. The L32A mutant showed significantly lower proton occupancies in the burst phase for several residues in the G helix, specifically F106, I107, E109 and A110, which are in close proximity to L32 in the X-ray structure of myoglobin, providing direct evidence that native-like helix B/helix G contacts are formed in the apomyoglobin burst phase intermediate. The L29A mutation resulted in an increase in burst phase proton occupancies for several residues in the E helix. Since these regions of the B and E helices are not in contact in the native myoglobin structure, these effects suggest the possibility of non-native B/E packing interactions in the kinetic intermediate. The differing effects of these B helix mutations on the apomyoglobin folding process suggests that each side-chain plays a different and important role in forming stable structure in the burst phase intermediate, and points to a role for both native-like and non-native contacts in stabilization of the folding intermediate.     

Multiple tryptophan probes reveal that ubiquitin folds via a late misfolded intermediate

Much of our understanding of protein folding mechanisms is derived from experiments using intrinsic fluorescence of natural or genetically inserted tryptophan (Trp) residues to monitor protein refolding and site-directed mutagenesis to determine the energetic role of amino acids in the native (N), intermediate (I) or transition (T) states. However, this strategy has limited use to study complex folding reactions because a single fluorescence probe may not detect all low-energy folding intermediates. To overcome this limitation, we suggest that protein refolding should be monitored with different solvent-exposed Trp probes. Here, we demonstrate the utility of this approach by investigating the controversial folding mechanism of ubiquitin (Ub) using Trp probes located at residue positions 1, 28, 45, 57, and 66. We first show that these Trp are structurally sensitive and minimally perturbing fluorescent probes for monitoring folding/unfolding of the protein. Using a conventional stopped-flow instrument, we show that ANS and Trp fluorescence detect two distinct transitions during the refolding of all five Trp mutants at low concentrations of denaturant: T₁, a denaturant-dependent transition and T₂, a slower transition, largely denaturant-independent. Surprisingly, some Trp mutants (Ub^M1W, Ub^S57W) display Trp fluorescence changes during T₁ that are distinct from the expected U → N transition suggesting that the denaturant-dependent refolding transition of Ub is not a U → N transition but represents the formation of a structurally distinct I-state (U → I). Alternatively, this U → I transition could be also clearly distinguished by using a combination of two Trp mutations Ub^F45W-T66W for which the two Trp probes that display fluorescence changes of opposite sign during T₁ and T₂ (Ub^F45W-T66W). Global fitting of the folding/unfolding kinetic parameters and additional folding-unfolding double-jump experiments performed on Ub^M1W, a mutant with enhanced fluorescence in the I-state, demonstrate that the I-state is stable, compact, misfolded, and on-pathway. These results illustrate how transient low-energy I-states can be characterized efficiently in complex refolding reactions using multiple Trp probes.     

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.     

Multiple folding mechanisms of protein ubiquitin

Based on the C(alpha) Go-type model, the folding kinetics and mechanisms of protein ubiquitin with mixed alpha/beta topology are studied by molecular dynamics simulations. The relaxation kinetics shows that there are three phases, namely the major phase, the intermediate phase and the slowest minor phase. The existence of these three phases are relevant to the phenomenon found in experiments. According to our simulations, the folding at high temperatures around the folding transition temperature T(f) is of a two-state process, and the folding nucleus is consisted of contacts between the front end of alpha-helix and the turn(4). The folding at low temperature (approximately T = 0.8) is also studied, where an A-state like structure is found lying on the major folding pathway. The appearance of this structure is related to the stability of the first part (residue 1-51) of protein ubiquitin. As the temperature decreases, the formation of secondary structures, tertiary structures and collapse of the protein are found to be decoupled gradually and the folding mechanism changes from the nucleation-condensation to the diffusion-collision. This feature indicates a unifying common folding mechanism for proteins. The intermediate phase is also studied and is found to represent a folding process via a long-lived intermediate state which is stabilized by strong interactions between the beta(1) and the beta(5) strand. These strong interactions are important for the function of protein ubiquitin as a molecular chaperone. Thus the intermediate phase is assumed as a byproduct of the requirement of protein function. In addition, the validity of the current Go-model is also investigated, and a lower limited temperature for protein ubiquitin T(limit) = 0.8 is proposed. At temperatures higher than this value, the kinetic traps due to glass dynamics cannot be significantly populated and the intermediate states can be reliably identified although there is slight chevron rollover in the folding rates. At temperature lower than T(limit), however, the traps due to glass dynamics become dominant and may be mistaken for real intermediate states. This limitation of valid temperature range prevents us to reveal the burst phase intermediate in the major folding phase since it might only be stabilized at temperatures lower than T(limit), according to experiments. Our works show that caution must be taken when studying low-temperature intermediate states by using the C(alpha) Go-models.     

Stabilization of the ribosomal protein S6 by trehalose is counterbalanced by the formation of a putative off-pathway species

The effect of trehalose on folding and stability of the small ribosomal protein S6 was studied. Non-disruptive point mutations distributed along the protein structure were analyzed to characterize the stabilizing effect of trehalose and map the folding pathway of S6. On average, the stability of the wild-type and S6 mutants increases by 3 kcal/mol M trehalose. Despite the non-specific thermodynamic stabilization mechanism, trehalose particularly stabilizes the less destabilized mutants. Folding/unfolding kinetics shows clearly that trehalose induces the collapse of the unfolded state to an off-pathway intermediate with non-native diffuse contacts. This state is similar to the collapsed state induced by high concentrations of stabilizing salts, as previously reported. Although it leads to the accumulation of this off-pathway intermediate, trehalose does not change the compactness of the transition state ensemble. Furthermore, the productive folding pathway of S6 is not affected by trehalose as shown by a Phi-value analysis. The unfolded state ensemble of S6 should be more compact in the presence of trehalose and therefore destabilized due to decreased conformational entropy. Increased compaction of the unfolded state ensemble might also occur for more stable mutants of S6, thus explaining the synergistic effect of trehalose and point mutations on protein stabilization.     

Mechanism of the kinetically-controlled folding reaction of subtilisin

Like many secreted proteases, subtilisin is kinetically stable in the mature form but unable to fold without assistance from its prodomain. The existence of high kinetic barriers to folding challenges many widely accepted ideas, namely, the thermodynamic determination of native structure and the sufficiency of thermodynamic stability to determine a pathway. The purpose of this article is to elucidate the physical nature of the kinetic barriers to subtilisin folding and to show how the prodomain overcomes these barriers. To address these questions, we have studied the bimolecular folding reaction of the subtilisin prodomain and a series of subtilisin mutants, which were designed to explore the steps in the folding reaction. Our analysis shows that inordinately slow folding of the mature form of subtilisin results from the accrued effects of two slow and sequential processes: (1) the formation of an unstable and topologically challenged intermediate and (2) the proline-limited isomerization of the intermediate to the native state. The low stability of nascent folding intermediates results in part from subtilisin's high dependence on metal binding for stability. Native subtilisin is thermodynamically unstable in the absence of bound metals. Because the two metal binding sites are formed late in folding, however, they contribute little to the stability of folding intermediates. The formation of productive folding intermediates is further hindered by the topological challenge of forming a left-handed crossover connection between beta-strands S2 and S3. This connection is critical to propagate the folding reaction. In the presence of the prodomain, folding proceeds through one major intermediate, which is stabilized by prodomain binding, independent of metal concentration and proline isomerization state. The prodomain also catalyzes the late proline isomerizations needed to form metal site B. Rate-limiting proline isomerization is common in protein folding, but its effect in slowing subtilisin folding is amplified because of the instability of the intermediate and an apparent need for simultaneous isomerization of multiple prolines in order to create metal site B. Thus, the kinetically controlled folding reaction of subtilisin, although unusual, is explained by the accrued effects of events found in other proteins.     

Sequential barriers and an obligatory metastable intermediate define the apparent two-state folding pathway of the ubiquitin-like PB1 domain of NBR1

The 90-residue N-terminal Phox and Bem1p (PB1) domain of NBR1 forms an α/β ubiquitin-like fold. Kinetic analysis using stopped-flow fluorescence reveals two-state kinetics; however, nonlinear effects in the denaturant dependence of the unfolding data demonstrate changes in the position of the rate-limiting barrier along the folding coordinate as the folding conditions change. The kinetics of wt-PB1 and several mutants show that this curvature is consistent with a single-pathway mechanism involving sequential transition states (TS1 and TS2) separated by a transiently populated high-energy intermediate, rather than movement of the transition state on a broad energy plateau. We show that the two transition states within the sequential model represent structurally and thermodynamically distinct species. TS1 is a collapsed state (α_TS1 = 0.71) with a large enthalpic barrier to formation that is rate-limiting under conditions that strongly favour folding. TS2 is highly native-like (α_TS2 = 0.93) and represents a late entropic barrier to formation of the native state. In support of the sequential transition state mechanism, we show that the G62A helix 2 substitution stabilises TS1 and the intermediate to such an extent that the latter becomes significantly populated, leading to the observation of a fast kinetic phase representing the initial U → I transition, with TS2 (α_TS2 = 0.87) becoming rate-limiting. The folding rate is not retarded by populating an intermediate, which would be expected for a misfold state, but is accelerated, suggesting that the I state is productive and on-pathway. The results show that the apparent two-state folding of the wt-PB1 domain occurs along a well-defined pathway involving structurally and thermodynamically distinct sequential transition states and an obligatory metastable intermediate that represents a productive local minimum in the energy landscape that increases the efficiency of barrier crossing through favourable effects on the entropy of activation.     

Experimental comparison of energy landscape features of ubiquitin family proteins

Small ubiquitin-related modifiers (SUMO1 and SUMO2) are ubiquitin family proteins, structurally similar to ubiquitin, differing in terms of their amino acid sequence and functions. Therefore, they provide a great platform for investigating sequence-structure-stability-function relationship. Here, we used chemical denaturation in comparing the folding-unfolding pathways of the SUMO proteins with their structural homologue ubiquitin (UF45W-pseudo wild-type [WT] tryptophan variant) with structurally analogous tryptophan mutations (SUMO1 [S1F66W], SUMO2 [S2F62W]). Equilibrium denaturation studies report that ubiquitin is the most stable protein among the three. The observed denaturant-dependent folding rates of SUMOs are much lower than ubiquitin and primarily exhibit a two-state folding pathway unlike ubiquitin, which has a kinetic folding intermediate. We hypothesize that, as SUMO proteins start off as slow folders, they avoid stabilizing their folding intermediates and the presence of which might further slow-down their folding rates. The denaturant-dependent unfolding of ubiquitin is the fastest, followed by SUMO2, and slowest for SUMO1. However, the spontaneous unfolding rate constant is the lowest for ubiquitin (~40 times), and similar for SUMOs. This correlation between thermodynamic stability and kinetic stability is achieved by having different unfolding transition state positions with respect to the solvent-accessible surface area, as quantified by the Tanford β _u values: ubiquitin (0.42) > SUMO2 (0.20) > SUMO1 (0.16). The results presented here highlight the unique energy landscape features which help in optimizing the folding-unfolding rates within a structurally homologous protein family.     

Low energy pathways and non-native interactions: the influence of artificial disulfide bridges on the mechanism of folding

Four versions of a beta-sheet protein (CD2.d1) have been made, each with a single artificial disulfide bond inserted into hairpin structures. Folding kinetics of reduced and oxidized forms shows bridge position strongly influences its effect on the folding reaction. Bridging residues 58 and 62 does not affect the rapidly formed intermediate (I) or rate-limiting transition (t) state, whereas bridging 33 and 38, or 31 and 41, lowers the t-state energy, with the latter having the stronger influence. Bridging residues 79 and 90 stabilizes both I- and t-states. To assess additivity in the energetic effects of these bridges, four double-bridge variants have also been made. All show precise additivity of overall stability, with two showing additivity when ground states and the rate-limiting t-state are assessed, i.e. no measurable change in the folding mechanism occurs. However, combining 31-41 and 79-90 bridges produces a molecule that folds through a different pathway, with a much more stable intermediate than expected and a much higher t-state barrier. This is explained by the artificial introduction of stabilizing, non-native contacts in the I-state. More surprisingly, for another double-bridge version (58-62 and 79-90) both I- and t-states are less stable than expected, showing that conformational constraints introduced by the two bridges prevent formation of non-native contacts that would otherwise stabilize the I- and t-states, thereby lowering the energy of the folding landscape in the wild-type (unbridged) molecule. We conclude that the lowest energy path for folding has I- and t-state structures that are stabilized by non-native interactions.     

Aggregation events occur prior to stable intermediate formation during refolding of interleukin 1beta

A point mutation, lysine 97 --> isoleucine (K97I), in a surface loop in the beta-sheet protein interleukin 1beta (IL-1beta), exhibits increased levels of inclusion body (IB) formation relative to the wild-type protein (WT) when expressed in Escherichia coli. Despite the common observation that less stable proteins are often found in IBs, K97I is more stable than WT. We examined the folding pathway of the mutant and wild-type proteins at pH 6.5 and 25 degrees C with manual-mixing and stopped-flow optical spectroscopy to determine whether changes in the properties of transiently populated species in vitro correlate with the observation of increased aggregation in vivo. The refolding reactions of the WT and K97I proteins are both described by three exponential processes. Two exponential processes characterize fast events (0.1-1.0 s) in folding while the third exponential process correlates with a slow (70 s) single pathway to and from the native state. The K97I replacement affects the earlier steps in the refolding pathway. Aggregation, absent in the WT refolding reaction, occurs in K97I above a critical protein concentration of 18 microM. This observation is consistent with an initial nucleation step mediating protein aggregation. Stopped-flow kinetic studies of the K97I aggregation process demonstrate that K97I aggregates most rapidly during the earliest refolding times, when unfolded protein conformers remain highly populated and the concentration of folding intermediates is low. Folding and aggregation studies together support a model in which the formation of stable folding intermediates afford protection against further K97I aggregation.     

A rationally designed monomeric variant of anthranilate phosphoribosyltransferase from Sulfolobus solfataricus is as active as the dimeric wild-type enzyme but less thermostable

The anthranilate phosphoribosyltransferase from Sulfolobus solfataricus (ssAnPRT) forms a homodimer with a hydrophobic subunit interface. To elucidate the role of oligomerisation for catalytic activity and thermal stability of the enzyme, we loosened the dimer by replacing two apolar interface residues with negatively charged residues (mutations I36E and M47D). The purified double mutant I36E+M47D formed a monomer with wild-type catalytic activity but reduced thermal stability. The single mutants I36E and M47D were present in a monomer-dimer equilibrium with dissociation constants of about 1 μM and 20 μM, respectively, which were calculated from the concentration-dependence of their heat inactivation kinetics. The monomeric form of M47D, which is populated at low subunit concentrations, was as thermolabile as monomeric I36E+M47D. Likewise, the dimeric form of I36E, which was populated at high subunit concentrations, was as thermostable as dimeric wild-type ssAnPRT. These findings show that the increased stability of wild-type ssAnPRT compared to the I36E+M47D double mutant is not caused by the amino acid exchanges per se but by the higher intrinsic stability of the dimer compared to the monomer. In accordance with the negligible effect of the mutations on catalytic activity and stability, the X-ray structure of M47D contains only minor local perturbations at the dimer interface. We conclude that the monomeric double mutant resembles the individual wild-type subunits, and that ssAnPRT is a dimer for stability but not for activity reasons.     

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.     

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.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Contribution of individual disulfide bonds to the oxidative folding of ribonuclease A

The eight cysteine residues of ribonuclease A form four disulfide bonds in the native protein. We have analyzed the folding of three double RNase A mutants (C65A/C72A, C58A/C110A, and C26A/C84A, lacking the C65-C72, C58-C110, and C26-C84 disulfide bonds, respectively) and two single mutants (C110A and C26A), in which a single cysteine is replaced with an alanine and the paired cysteine is present in the reduced form. The folding of these mutants was carried out in the presence of oxidized and reduced glutathione, which constitute the main redox agents present within the ER. The use of mass spectrometry in the analysis of the folding processes allowed us (i) to follow the formation of intermediates and thus the pathway of folding of the RNase A mutants, (ii) to quantitate the intermediates that formed, and (iii) to compare the rates of formation of intermediates. By comparison of the folding kinetics of the mutants with that of wild-type RNase A, the contribution of each disulfide bond to the folding process has been evaluated. In particular, we have found that the folding of the C65A/C72A mutant occurs on the same time scale as that of the wild-type protein, thus suggesting that the removal of the C65-C72 disulfide bond has no effect on the kinetics of RNase A folding. Conversely, the C58A/C110A and C26A/C84A mutants fold much more slowly than the wild-type protein. The removal of the C58-C110 and C26-C84 disulfide bonds has a dramatic effect on the kinetics of RNase A folding. Results described in this paper provide specific information about conformational folding events in the regions involving the mutated cysteine residues, thus contributing to a better understanding of the complex mechanism of oxidative folding.     

The Plastic Energy Landscape of Protein Folding: A TRIANGULAR FOLDING MECHANISM WITH AN EQUILIBRIUM INTERMEDIATE FOR A SMALL PROTEIN DOMAIN*

Protein domains usually fold without or with only transiently populated intermediates, possibly to avoid misfolding, which could result in amyloidogenic disease. Whether observed intermediates are productive and obligatory species on the folding reaction pathway or dispensable by-products is a matter of debate. Here, we solved the crystal structure of a small protein domain, SAP97 PDZ2 I342W C378A, and determined its folding pathway. The presence of a folding intermediate was demonstrated both by single and double-mixing kinetic experiments using urea-induced (un)folding as well as ligand-induced folding. This protein domain was found to fold via a triangular scheme, where the folding intermediate could be either on- or off-pathway, depending on the experimental conditions. Furthermore, we found that the intermediate was present at equilibrium, which is rarely seen in folding reactions of small protein domains. The folding mechanism observed here illustrates the roughness and plasticity of the protein folding energy landscape, where several routes may be employed to reach the native state. The results also reconcile the folding mechanisms of topological variants within the PDZ domain family.     

Unfolding and folding kinetics of amyotrophic lateral sclerosis-associated mutant Cu,Zn superoxide dismutases

More than 110 mutations in dimeric, Cu,Zn superoxide dismutase (SOD) have been linked to the fatal neurodegenerative disease, amyotrophic lateral sclerosis (ALS). In both human patients and mouse model studies, protein misfolding has been implicated in disease pathogenesis. A central step in understanding the misfolding/aggregation mechanism of this protein is the elucidation of the folding pathway of SOD. Here we report a systematic analyses of unfolding and folding kinetics using single- and double-jump experiments as well as measurements as a function of guanidium chloride, protein, and metal concentration for fully metallated (holo) pseudo wild-type and ALS-associated mutant (E100G, G93R, G93A, and metal binding mutants G85R and H46R) SODs. The kinetic mechanism for holo SODs involves native dimer, monomer intermediate, and unfolded monomer, with variable metal dissociation from the monomeric states depending on solution conditions. The effects of the ALS mutations on the kinetics of the holoproteins in guanidium chloride are markedly different from those observed previously for acid-induced unfolding and for the unmetallated (apo) forms of the proteins. The mutations decrease the stability of holo SOD mainly by increasing unfolding rates, which is particularly pronounced for the metal-binding mutants, and have relatively smaller effects on the observed folding kinetics. Mutations also seem to favour increased formation of a Zn-free monomer intermediate, which has been implicated in the formation of toxic aggregates. The results reveal the kinetic basis for the extremely high stability of wild-type holo SOD and the possible consequences of kinetic changes for disease.     

The folding pathway of ubiquitin from all-atom molecular dynamics simulations

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.