Population of on-pathway intermediates in the folding of ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.     

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.     

Folding kinetics of the S100A11 protein dimer studied by time-resolved electrospray mass spectrometry and pulsed hydrogen-deuterium exchange

This study reports the application of electrospray ionization (ESI) mass spectrometry (MS) with on-line rapid mixing for millisecond time-resolved studies of the refolding and assembly of a dimeric protein complex. Acid denaturation of S100A11 disrupts the native homodimeric protein structure. Circular dichroism and HSQC nuclear magnetic resonance measurements reveal that the monomeric subunits unfold to a moderate degree but retain a significant helicity and some tertiary structural elements. Following a rapid change in solution conditions to a slightly basic pH, the native protein reassembles with an effective rate constant of 6 s(-)(1). The ESI charge state distributions measured during the reaction suggest the presence of three kinetic species, namely, a relatively unfolded monomer (M(U)), a more tightly folded monomeric reaction intermediate (M(F)), and dimeric S100A11. These three forms exhibit distinct calcium binding properties, with very low metal loading levels for M(U), up to two calcium ions for M(F), and up to four for the dimer. Surprisingly, on-line pulsed hydrogen-deuterium exchange (HDX) reveals that each of the monomeric forms of the protein comprises two subspecies that can be distinguished on the basis of their isotope exchange levels. As the reaction proceeds, the more extensively labeled species are depleted. The exponential nature of the measured intensity-time profiles implies that the rate-determining step of the overall process is a unimolecular event. The kinetics are consistent with a sequential folding and assembly mechanism involving two increasingly nativelike monomeric intermediates en route to the native S100A11 dimer.     

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.     

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.     

Distinguishing between two-state and three-state models for ubiquitin folding

Conflicting results exist regarding whether the folding of mammalian ubiquitin at 25 degrees C is a simple, two-state kinetic process or a more complex, three-state process with a defined kinetic intermediate. We have measured folding rate constants up to about 1000 s(-1) using conventional rapid mixing methods in single-jump, double-jump, and continuous-flow modes. The linear dependence of folding rates on denaturant concentration and the lack of an unaccounted "burst-phase" change for the fluorescence signal indicate that a two-state folding model is adequate to describe the folding pathway. This behavior also is seen for folding in the presence of the stabilizing additives 0.23 M sodium sulfate and 1 M sodium chloride. These results stress the need for caution in interpreting deviations from ideal two-state "chevron" behavior when folding is heterogeneous or folding rate constants are near the detection limit.     

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.     

Refolding of a small all beta-sheet protein proceeds with accumulation of kinetic intermediates

The refolding kinetics of Cobrotoxin (CBTX), a small all beta-sheet protein is investigated using a variety of biophysical techniques including quenched-flow hydrogen-deuterium (H/D) exchange in conjunction with two-dimensional NMR spectroscopy. Urea-induced equilibrium unfolding of CBTX follows a two-state mechanism with no distinct intermediates. The protein is observed to fold very rapidly within 250 ms. Both the refolding and the unfolding limbs of the chevron plot of CBTX show a prominent curvature suggesting the accumulation of kinetic intermediates. Quenched-flow H/D exchange data suggest the presence of a broad continuum of kinetic intermediates between the unfolded and native states of the protein. Comparison of the native state hydrogen exchange data and the results of the quenched-flow H/D exchange experiments, reveals that the residues constituting the folding core of CBTX are not a subset of the slow exchange core. To our knowledge, this is the first report wherein the refolding of a small all beta-sheet protein is shown to be a multi-step process involving the accumulation of kinetic intermediates.     

Structure is lost incrementally during the unfolding of barstar

Coincidental equilibrium unfolding transitions observed by multiple structural probes are taken to justify the modeling of protein unfolding as a two-state, N <==> U, cooperative process. However, for many of the large number of proteins that undergo apparently two-state equilibrium unfolding reactions, folding intermediates are detected in kinetic experiments. The small protein barstar is one such protein. Here the two-state model for equilibrium unfolding has been critically evaluated in barstar by estimating the intramolecular distance distribution by time-resolved fluorescence resonance energy transfer (TR-FRET) methods, in which fluorescence decay kinetics are analyzed by the maximum entropy method (MEM). Using a mutant form of barstar containing only Trp 53 as the fluorescence donor and a thionitrobenzoic acid moiety attached to Cys 82 as the fluorescence acceptor, the distance between the donor and acceptor has been shown to increase incrementally with increasing denaturant concentration. Although other probes, such as circular dichroism and fluorescence intensity, suggest that the labeled protein undergoes two-state equilibrium unfolding, the TR-FRET probe clearly indicates multistate equilibrium unfolding. Native protein expands progressively through a continuum of native-like forms that achieve the dimensions of a molten globule, whose heterogeneity increases with increasing denaturant concentration and which appears to be separated from the unfolded ensemble by a free energy barrier.     

Relationship between the native-state hydrogen exchange and folding pathways of a four-helix bundle protein

The hydrogen exchange behavior of a four-helix bundle protein in low concentrations of denaturant reveals some partially unfolded forms that are significantly more stable than the fully unfolded state. Kinetic folding of the protein, however, is apparently two-state in the absence of the accumulation of early folding intermediates. The partially unfolded forms are either as folded as or more folded than the rate-limiting transition state and appear to represent the major intermediates in a folding and unfolding reaction. These results are consistent with the suggestion that partially unfolded intermediates may form after the rate-limiting step for small proteins with apparent two-state folding kinetics.     

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.     

Insights into the stability of native and partially folded states of ubiquitin: effects of cosolvents and denaturants on the thermodynamics of protein folding

The thermodynamics of the native<-->A state and native<-->unfolded transitions for ubiquitin have been characterized in detail using the denaturants methanol and guanidinium chloride (Gdn.HCl) both separately and in combination. Gdn.HCl destabilizes the partially folded alcohol-induced A state such that the effects of alcoholic solvents on the native<-->unfolded transition can be investigated directly via a two-state model. The combined denaturing effects of methanol and Gdn.HCl appear to conform to a simple additive model. We show that ubiquitin folds and unfolds cooperatively in all cases, forming the same "native" state; however, the thermodynamics of the N<-->U transition change dramatically between alcoholic and Gdn.HCl solutions, with folding in aqueous methanol associated with large negative enthalpy and entropy terms at 298 K with a gradual falloff in DeltaC(p) at higher methanol concentrations, as previously reported for the N<-->A transition (Woolfson, D. N., Cooper, A., Harding, M. M., Williams, D. H., and Evans, P. A. (1993) J. Mol. Biol. 229, 502-511.). Both the N<-->U and the N<-->A transitions are enthalpy driven to a similar extent. We conclude that under these conditions van der Waals interactions in the packing of the nonpolar protein core, which is common to both the N<-->U and the N<-->A transitions, appear to drive folding in the absence of entropic effects associated with release of ordered solvent (hydrophobic effect). Solvent transfer studies of hydrocarbons into alcoholic solvents, with and without Gdn.HCl, are consistent with a large enthalpic driving force for burial of a nonpolar surface, with a linear dependence of protein stability (DeltaG(N)(<-->)(U)) on cosolvent concentration reflected in a similar linear dependence of hydrocarbon solubility. The data demonstrate that the hydrophobic effect is not a prerequisite for specific stabilization of the native state or the A state and that van der Waals packing of the nonpolar core appears to be the dominant factor in stabilization of the native state.     

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

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

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

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

The nature of protein folding pathways: The classical versus the new view

Summary Pulsed hydrogen exchange and other studies of the kinetic refolding pathways of several small proteins have established that folding intermediates with native-like secondary structures are well populated, but these studies have also shown that the folding kinetics are not well synchronized. Older studies of the kinetics of formation of the native protein, monitored by optical probes, indicate that the folding kinetics should be synchronized. The model commonly used in these studies is the simple sequential model, which postulates a unique folding pathway with defined and sequential intermediates. Theories of the folding process and Monte Carlo simulations of folding suggest that neither the folding pathway nor the set of folding intermediates is unique, and that folding intermediates accumulate because of kinetic traps caused by partial misfolding. Recent experiments with cytochrome c lend support to this new view of folding pathways. These different views of the folding process are discussed. Misfolding and consequent slowing down of the folding process as a result of cis-trans isomerization about prolyl peptide bonds in the unfolded protein are well known; isomerization occurs before refolding is initiated. The occurrence of equilibrium intermediates on the kinetic folding pathways of some proteins, such as -lactalbumin and apomyoglobin, argues that these intermediates are not caused by kinetic traps but rather are stable intermediates under certain conditions, and this conclusion is consistent with a sequential model of folding. Folding reactions with successive kinetic intermediates, in which late intermediates are more highly folded than early intermediates, indicate that folding is hierarchical. New experiments that test the predictions of the classical and the new views are needed.     

Pathway diversity and concertedness in protein folding: an ab-initio approach

Making use of an ab-initio folding simulator, we generate in vitro pathways leading to the native fold in moderate size single- domain proteins. The assessment of pathway diversity is not biased by any a priori information on the native fold. We focus on two study cases, hyperthermophile variant of protein G domain (1gb4) and ubiquitin (1ubi), with the same topology but different context dependence in their native folds. We demonstrate that a quenching of structural fluctuations is achieved once the proteins find a stationary plateau maximizing the number of highly protected hydrogen bonds. This enables us to identify the folding nucleus and show that folding does not become expeditious until a concerted event takes place generating a topology able to prevent water attack on a maximal number of hydrogen bonds. This result is consistent with the standard nucleation mechanism postulated for two-state folders. Pathway diversity is correlated with the extent of conflict between local structural propensity and large-scale context, rather than with contact order: In highly context-dependent proteins, the success of folding cannot rely on a single fortuitous event in which local propensity is overruled by large-scale effects. We predict mutational Pi values on individual pathways, compute ensemble averages and predict extent of surface burial and percentage of hydrogen bonding on each component of the transition state ensemble, thus deconvoluting individual folding-route contributions to the averaged two-state kinetic picture. Our predicted kinetic isotopic effects find experimental support and lead to further probes. Finally, the molecular redesign potentiality of the method, aimed at increasing folding expediency, is explored.     

Nonideality and protein thermal denaturation

We studied the thermal denaturation of eglin c by using CD spectropolarimetry and differential scanning calorimetry (DSC). At low protein concentrations, denaturation is consistent with the classical two-state model. At concentrations greater than several hundred microM, however, the calorimetric enthalpy and the midpoint transition temperature increase with increasing protein concentration. These observations suggested the presence of intermediates and/or native state aggregation. However, the transitions are symmetric, suggesting that intermediates are absent, the DSC data do not fit models that include aggregation, and analytical ultracentrifugation (AUC) data show that native eglin c is monomeric. Instead, the AUC data show that eglin c solutions are nonideal. Analysis of the AUC data gives a second virial coefficient that is close to values calculated from theory and the DSC data are consistent with the behavior expected for nonideal solutions. We conclude that the concentration dependence is caused by differential nonideality of the native and denatured states. The nondeality arises from the high charge of the protein at acid pH and is exacerbated by low buffer concentrations. Our conclusion may explain differences between van't Hoff and calorimetric denaturation enthalpies observed for other proteins whose behavior is otherwise consistent with the classical two-state model.     

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.     

The time-dependent distribution of phosphorylated intermediates in native sarcoplasmic reticulum Ca2+-ATPase from skeletal muscle is not compatible with a linear kinetic model

Quenched-flow mixing was used to characterize the kinetic behavior of the intermediate reactions of the skeletal muscle sarcoplasmic reticulum (SR) Ca-ATPase (SERCA1) at 2 and 21 degrees C. At 2 degrees C, phosphorylation of SR Ca-ATPase with 100 microM ATP labeled one-half of the catalytic sites with a biphasic time dependence [Mahaney, J. E., Froehlich, J. P., and Thomas, D. D. (1995) Biochemistry 34, 4864-4879]. Chasing the phosphoenzyme (EP) with 1.66 mM ADP 10 ms after the start of phosphorylation revealed mostly ADP-insensitive E2P (95% of EP(total)), consistent with its rapid formation from ADP-sensitive E1P. The consecutive relationship of the phosphorylated intermediates predicts a decrease in the proportion of E1P ([E1P]/[EP(total)]) with increasing phosphorylation time. Instead, after 10 ms the proportion of E1P increased and that of E2P decreased until they reached a constant 1:1 stoichiometry ([E1P]:[E2P] approximately 1). At 21 degrees C, phosphorylation displayed a transient overshoot associated with an inorganic phosphate (P(i)) burst, reflecting increased turnover of E2P at the higher temperature. The P(i) burst exceeded the decay of the EP overshoot, suggesting that rephosphorylation of the enzyme occurs before the recycling step (E2 --> E1). This behavior and the reversed order of accumulation of phosphorylated intermediates at 2 degrees C are not compatible with the conventional linear consecutive reaction mechanism: E1 + ATP --> E1.ATP --> E1P + ADP --> E2P --> E2.P(i) --> E1 + P(i). Solubilization of the Ca-ATPase into monomers using the nonionic detergent C(12)E(8) gave a pattern of phosphorylation in which E1P and E2P behave like consecutive intermediates. Kinetic modeling of the C(12)E(8)-solubilized SR Ca-ATPase showed that it behaves according to the conventional Ca-ATPase reaction mechanism, consistent with monomeric catalytic function. We conclude that the nonconforming features of native SERCA1 arise from oligomeric protein conformational interactions that constrain the subunits to a staggered or out-of-phase mode of operation.     

Pathway Diversity and Concertedness in Protein Folding: An ab-initio Approach

Abstract

Making use of an ab-initio folding simulator, we generate in vitro pathways leading to the native fold in moderate size single-domain proteins. The assessment of pathway diversity is not biased by any a-priori information on the native fold. We focus on two study cases, hyperthermophile variant of protein G domain (1gb4) and ubiquitin (1ubi), with the same topology but different context dependence in their native folds. We demonstrate that a quenching of structural fluctuations is achieved once the proteins find a stationary plateau maximizing the number of highly protected hydrogen bonds. This enables us to identify the folding nucleus and show that folding does not become expeditious until a concerted event takes place generating a topology able to prevent water attack on a maximal number of hydrogen bonds. This result is consistent with the standard nucleation mechanism postulated for two-state folders. Pathway diversity is correlated with the extent of conflict between local structural propensity and large-scale context, rather than with contact order: In highly context-dependent proteins, the success of folding cannot rely on a single fortuitous event in which local propensity is overruled by large-scale effects. We predict mutational Φ values on individual pathways, compute ensemble averages and predict extent of surface burial and percentage of hydrogen bonding on each component of the transition state ensemble, thus deconvoluting individual folding-route contributions to the averaged two-state kinetic picture. Our predicted kinetic isotopic effects find experimental support and lead to further probes. Finally, the molecular redesign potentiality of the method, aimed at increasing folding expediency, is explored.