The early folding kinetics of apomyoglobin.

The folding pathway of apomyoglobin has been experimentally shown to have early kinetic intermediates involving the A, B, G, and H helices. The earliest detected kinetic events occur on a ns to micros time scale. We show that the early folding kinetics of apomyoglobin may be understood as the association of nascent helices through a network of diffusion-collision-coalescence steps G + H <--> GH + A <--> AGH + B <--> ABGH obtained by solving the diffusion-collision model in a chemical kinetics approximation. Our reproduction of the experimental results indicates that the model is a useful way to analyze folding data. One prediction from our fit is that the nascent A and H helices should be relatively more helix-like before coalescence than the other apomyoglobin helices.     

Stein and Moore Award address. The molten globule intermediate of apomyoglobin and the process of protein folding.

The molten globule model for the beginning of the folding process, which originated with Kuwajima's studies of alpha-lactalbumin (Kuwajima, K., 1989, Proteins Struct. Funct. Genet. 6, 87-103, and references therein), states that, for those proteins that exhibit equilibrium molten globule intermediates, the molten globule is a major kinetic intermediate near the start of the folding pathway. Pulsed hydrogen-deuterium exchange measurements confirm this model for apomyoglobin (Jennings, P.A. & Wright, P.E., in prep.). The energetics of the acid-induced unfolding transition, which have been determined by fitting a minimal three-state model (N<-->I<-->U; N = native, I = molten globule intermediate, U = unfolded) show that I is more stable than U at neutral pH (Barrick, D. & Baldwin, R.L., 1993, Biochemistry 32, in press), which provides an explanation for why I is formed from U at the start of folding. Hydrogen exchange rates measured by two-dimensional NMR for individual peptide NH protons, taken together with the CD spectrum of I, indicate that moderately stable helices are present in I at the locations of the A, G, and H helices of native myoglobin (Hughson, F.M., Wright, P.E., & Baldwin, R.L., 1990, Science 249, 1544-1548). Directed mutagnesis experiments indicate that the interactions between the A, G, and H helices in I are loose (Hughson, F.M., Barrick, D., & Baldwin, R.L., 1991, Biochemistry 30, 4113-4118), which can explain why I is formed rapidly from U at the start of folding.(ABSTRACT TRUNCATED AT 250 WORDS)     

Energy barriers, cooperativity, and hidden intermediates in the folding of small proteins

Current theoretical views of the folding process of small proteins (< approximately 100 amino acids) postulate that the landscape of potential mean force (PMF) for the formation of the native state has a funnel shape and that the free energy barrier to folding arises from the chain configurational entropy only. However, recent theoretical studies on the formation of hydrophobic clusters with explicit water suggest that a barrier should exist on the PMF of folding, consistent with the fact that protein folding generally involves a large positive activation enthalpy at room temperature. In addition, high-resolution structural studies of the hidden partially unfolded intermediates have revealed the existence of non-native interactions, suggesting that the correction of the non-native interactions during folding should also lead to barriers on PMF. To explore the effect of a PMF barrier on the folding behavior of proteins, we modified Zwanzig's model for protein folding with an uphill landscape of PMF for the formation of transition states. We found that the modified model for short peptide segments can satisfy the thermodynamic and kinetic criteria for an apparently two-state folding. Since the Levinthal paradox can be solved by a stepwise folding of short peptide segments, a landscape of PMF with a locally uphill search for the transition state and cooperative stabilization of folding intermediates/native state is able to explain the available experimental results for small proteins. We speculate that the existence of cooperative hidden folding intermediates in small proteins could be the consequence of the highly specific structures of the native state, which are selected by evolution to perform specific functions and fold in a biologically meaningful time scale.     

Diffusion-collision model study of misfolding in a four-helix bundle protein.

Proteins with complex folding kinetics will be susceptible to misfolding at some stage in the folding process. We simulate this problem by using the diffusion-collision model to study non-native kinetic intermediate misfolding in a four-helix bundle protein. We find a limit on the size of the pairwise hydrophobic area loss in non-native intermediates, such that burying above this limit creates long-lasting non-native kinetic intermediates that would disrupt folding and prevent formation of the native state. Our study of misfolding suggests a method for limiting the production of misfolded kinetic intermediates for helical proteins and could, perhaps, lead to more efficient production of proteins in bulk.     

Monte Carlo Simulations of Proteins in Cages: Influence of Confinement on the Stability of Intermediate States

We study the folding of small proteins inside confining potentials. Proteins are described using an effective potential model that contains the Ramachandran angles as degrees of freedom and does not need any a priori information about the native state. Hydrogen bonds, dipole-dipole-, and hydrophobic interactions are taken explicitly into account. An interesting feature displayed by this potential is the presence of metastable intermediates between the unfolded and native states. We consider different types of confining potentials to describe proteins folding inside cages with repulsive or attractive walls. Using the Wang-Landau algorithm, we determine the density of states and analyze in detail the thermodynamical properties of the confined proteins for different sizes of the cages. We show that confinement dramatically reduces the phase space available to the protein and that the presence of intermediate states can be controlled by varying the properties of the confining potential. Cages with strongly attractive walls destabilize the intermediate states and lead to a two-state folding into a configuration that is less stable than the native structure. However, cages with slightly attractive walls enhance the stability of native structure and induce a folding process, which occurs through intermediate configurations.     

Early events during folding of wild-type staphylococcal nuclease and a single-tryptophan variant studied by ultrarapid mixing

A continuous-flow mixing device with a dead time of 100 micros coupled with intrinsic tryptophan and 1-anilinonaphthalene-8-sulfonate (ANS) fluorescence was used to monitor structure formation during early stages of the folding of staphylococcal nuclease (SNase). A variant with a unique tryptophan fluorophore in the N-terminal beta-barrel domain (Trp76 SNase) was obtained by replacing the single Trp140 in wild-type SNase with His in combination with Trp substitution of Phe76. A common background of P47G, P117G and H124L mutations was chosen in order to stabilize the protein and prevent accumulation of cis proline isomers under native conditions. In contrast to WT(*) SNase, which shows no changes in tryptophan fluorescence prior to the rate-limiting folding step ( approximately 100 ms), the F76W/W140H variant shows additional changes (enhancement) during an early folding phase with a time constant of 75 micros. Both proteins exhibit a major increase in ANS fluorescence and identical rates for this early folding event. These findings are consistent with the rapid accumulation of an ensemble of states containing a loosely packed hydrophobic core involving primarily the beta-barrel domain while the specific interactions in the alpha-helical domain involving Trp140 are formed only during the final stages of folding. The fact that both variants exhibit the same number of kinetic phases with very similar rates confirms that the folding mechanism is not perturbed by the F76W/W140H mutations. However, the Trp at position 76 reports on the rapid formation of a hydrophobic cluster in the N-terminal beta-sheet region while the wild-type Trp140 is silent during this early stage of folding. Quantitative modeling of the (un)folding kinetics and thermodynamics of these two proteins versus urea concentration revealed that the F76W/W140H mutation selectively destabilizes the native state relative to WT(*) SNase while the stability of transient intermediates remains unchanged, leading to accumulation of intermediates under equilibrium conditions at moderate denaturant concentrations.     

Structural changes and fluctuations of proteins. II. Analysis of the denaturation of globular proteins.

The statistical thermodynamic model of protein structure proposed in paper I is developed with special attention to the hydrophobic interaction. Calorimetric measurements of the thermal denaturation of five globular proteins, ribonuclease A, lysozyme, alpha-chymotrypsin, cytochrome c, and myoglobin, are quantitatively analyzed using the model. The thermodynamic parameters obtained by the least squares method reflect the global, average properties of proteins and are in good agreement with the expected values estimated from experimental and theoretical studies for model peptides. The average bond energy epsilon is well related to the tertiary structure of each protein. However, the difference in the parameters between different proteins is not observed for the cooperative energy ZJ and the chain entropy alpha. The individuality of a protein as far as its structural stability is concerned, is mainly reflected by the parameter gamma specifying the hydrophobic nature of a protein. The model is further applied in the analysis of several aspects of the structural stability of globular proteins. Denaturation induced by denaturants, salts, and pH are also explained by the model in a unified manner.     

Decoding Structural Properties of a Partially Unfolded Protein Substrate: En Route to Chaperone Binding

Many proteins comprising of complex topologies require molecular chaperones to achieve their unique three-dimensional folded structure. The E.coli chaperone, GroEL binds with a large number of unfolded and partially folded proteins, to facilitate proper folding and prevent misfolding and aggregation. Although the major structural components of GroEL are well defined, scaffolds of the non-native substrates that determine chaperone-mediated folding have been difficult to recognize. Here we performed all-atomistic and replica-exchange molecular dynamics simulations to dissect non-native ensemble of an obligate GroEL folder, DapA. Thermodynamics analyses of unfolding simulations revealed populated intermediates with distinct structural characteristics. We found that surface exposed hydrophobic patches are significantly increased, primarily contributed from native and non-native β-sheet elements. We validate the structural properties of these conformers using experimental data, including circular dichroism (CD), 1-anilinonaphthalene-8-sulfonic acid (ANS) binding measurements and previously reported hydrogen-deutrium exchange coupled to mass spectrometry (HDX-MS). Further, we constructed network graphs to elucidate long-range intra-protein connectivity of native and intermediate topologies, demonstrating regions that serve as central “hubs”. Overall, our results implicate that genomic variations (or mutations) in the distinct regions of protein structures might disrupt these topological signatures disabling chaperone-mediated folding, leading to formation of aggregates.     

OneG: a computational tool for predicting cryptic intermediates in the unfolding kinetics of proteins under native conditions

Understanding the relationships between conformations of proteins and their stabilities is one key to address the protein folding paradigm. The free energy change (ΔG) of unfolding reactions of proteins is measured by traditional denaturation methods and native hydrogen-deuterium (H/D) exchange methods. However, the free energy of unfolding (ΔG(U)) and the free energy of exchange (ΔG(HX)) of proteins are not in good agreement, though the experimental conditions of both methods are well matching to each other. The anomaly is due to any one or combinations of the following reasons: (i) effects of cis-trans proline isomerisation under equilibrium unfolding reactions of proteins (ii) inappropriateness in accounting the baselines of melting curves (iii) presence of cryptic intermediates, which may elude the melting curve analysis and (iv) existence of higher energy metastable states in the H/D exchange reactions of proteins. Herein, we have developed a novel computational tool, OneG, which accounts the discrepancy between ΔG(U) and ΔG(HX) of proteins by systematically accounting all the four factors mentioned above. The program is fully automated and requires four inputs: three-dimensional structures of proteins, ΔG(U), ΔG(U)(*) and residue-specific ΔG(HX) determined under EX2-exchange conditions in the absence of denaturants. The robustness of the program has been validated using experimental data available for proteins such as cytochrome c and apocytochrome b(562) and the data analyses revealed that cryptic intermediates of the proteins detected by the experimental methods and the cryptic intermediates predicted by the OneG for those proteins were in good agreement. Furthermore, using OneG, we have shown possible existence of cryptic intermediates and metastable states in the unfolding pathways of cardiotoxin III and cobrotoxin, respectively, which are homologous proteins. The unique application of the program to map the unfolding pathways of proteins under native conditions have been brought into fore and the program is publicly available at http://sblab.sastra.edu/oneg.html.     

Conformational specificity of the chaperonin GroEL for the compact folding intermediates of alpha-lactalbumin.

The chaperonin GroEL binds unfolded polypeptides, preventing aggregation, and then mediates their folding in an ATP-dependent process. To understand the structural features in non-native polypeptides recognized by GroEL, we have used alpha-lactalbumin (alpha LA) as a model substrate. alpha LA (14.2 kDa) is stabilized by four disulfide bonds and a bound Ca2+ ion, offering the possibility of trapping partially folded disulfide intermediates between the native and the fully unfolded state. The conformers of alpha LA with high affinity for GroEL are compact, containing up to three disulfide bonds, and have significant secondary structure, but lack stable tertiary structure and expose hydrophobic surfaces. Complex formation requires almost the complete alpha LA sequence and is strongly dependent on salts that stabilize hydrophobic interactions. Unfolding of alpha LA to an extended state as well as the burial of hydrophobic surface upon formation of ordered tertiary structure prevent the binding to GroEL. Interestingly, GroEL interacts only with a specific subset of the many partially folded disulfide intermediates of alpha LA and thus may influence in vitro the kinetics of the folding pathways that lead to disulfide bonds with native combinations. We conclude that the chaperonin interacts with the hydrophobic surfaces exposed by proteins in a flexible compact intermediate or molten globule state.     

Quantitative analysis of the kinetics of denaturation and renaturation of barstar in the folding transition zone.

The fluorescence-monitored kinetics of folding and unfolding of barstar by guanidine hydrochloride (GdnHCl) in the folding transition zone, at pH 7, 25 degrees C, have been quantitatively analyzed using a 3-state mechanism: U(S)<-->UF<-->N. U(S) and UF are slow-refolding and fast-refolding unfolded forms of barstar, and N is the native protein. U(S) and UF probably differ in possessing trans and cis conformations, respectively, of the Tyr 47-Pro 48 bond. The 3-state model could be used because the kinetics of folding and unfolding of barstar show 2 phases, a fast phase and a slow phase, and because the relative amplitudes of the 2 phases depend only on the final refolding conditions and not on the initial conditions. Analysis of the observed kinetics according to the 3-state model yields the values of the 4 microscopic rate constants that describe the transitions between the 3 states at different concentrations of GdnHCl. The value of the equilibrium unfolded ratio U(S):UF (K21) and the values of the rate constants of the U(S)-->UF and UF-->U(S) reactions, k12 and k21, respectively, are shown to be independent of the concentration of GdnHCl. K21 has a value of 2.1 +/- 0.1, and k12 and k21 have values of 5.3 x 10(-3) s-1 and 11.2 x 10(-3) s-1, respectively. Double-jump experiments that monitor reactions that are silent to fluorescence monitoring were used to confirm the values of K21, k12, and k21 obtained from the 3-state analysis and thereby the validity of the 3-state model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamics and kinetics of protein folding: an evolutionary perspective

This article appeals to an evolutionary model which postulates that primordial proteins were described by small polypeptide chains which (i) lack disulfide bridges, and (ii) display slow folding rates with multi-state kinetics, to determine relations between structural properties of proteins and their folding kinetics. We parameterize the energy landscape of proteins in terms of thermodynamic activation variables. The model studies evolutionary changes in these thermodynamic parameters, and we invoke relations between these activation variables and structural properties of the protein to predict the following correspondence between protein structure and folding kinetics. 1. Proteins with inter- and intra-chain disulfide bridges: large variability in both folding rates and stability of intermediates, multi-state kinetics. 2. Proteins which lack inter and intra-chain disulfide bridges. 2.1 Single-domain chains: fast folding rates; unstable intermediates; two-state kinetics. 2.2 Multi-domain monomers: intermediate rates; metastable intermediates; multi-state kinetics. 2.3 Multi-domain oligomers: slow rates; metastable intermediates; multi-state kinetics. The evolutionary model thus provides a kinetic characterization of one important subfamily of proteins which we describe by the following properties: Folding dynamics of single-domain proteins which lack disulfide bridges are described by two-state kinetics. Folding rate of this class of proteins is positively correlated with the thermodynamic stability of the folded state.     

Secondary structure length as a determinant of folding rate of proteins with two- and three-state kinetics

We present a simple method for determining the folding rates of two- and three-state proteins from the number of residues in their secondary structures (secondary structure length). The method is based on the hypothesis that two- and three-state foldings share a common pattern. Three-state proteins first condense into metastable intermediates, subsequent forming of alpha-helices, turns, and beta-sheets at slow rate-limiting step. The folding rate of such proteins anticorrelate with the length of these beta-secondary structures. It is also assumed that in two-state folding, rapidly folded alpha-helices and turns may facilitate formation of fleeting unobservable intermediates and thus show two-state behavior. There is an inverse relationship between the folding rate and the length of beta-sheets and loops. Our study achieves 94.0 and 88.1% correlations with folding rates determined experimentally for 21 three- and 38 two-state proteins, respectively, suggesting that protein-folding rates are determined by the secondary structure length. The kinetic kinds are selected on the basis of a competitive formation of hydrophobic collapse and alpha-structure in early intermediates.     

Deciphering the Hidden Informational Content of Protein Sequences: FOLDABILITY OF PROINSULIN HINGES ON A FLEXIBLE ARM THAT IS DISPENSABLE IN THE MATURE HORMONE*

Protein sequences encode both structure and foldability. Whereas the interrelationship of sequence and structure has been extensively investigated, the origins of folding efficiency are enigmatic. We demonstrate that the folding of proinsulin requires a flexible N-terminal hydrophobic residue that is dispensable for the structure, activity, and stability of the mature hormone. This residue (Phe^B1 in placental mammals) is variably positioned within crystal structures and exhibits ¹H NMR motional narrowing in solution. Despite such flexibility, its deletion impaired insulin chain combination and led in cell culture to formation of non-native disulfide isomers with impaired secretion of the variant proinsulin. Cellular folding and secretion were maintained by hydrophobic substitutions at B1 but markedly perturbed by polar or charged side chains. We propose that, during folding, a hydrophobic side chain at B1 anchors transient long-range interactions by a flexible N-terminal arm (residues B1–B8) to mediate kinetic or thermodynamic partitioning among disulfide intermediates. Evidence for the overall contribution of the arm to folding was obtained by alanine scanning mutagenesis. Together, our findings demonstrate that efficient folding of proinsulin requires N-terminal sequences that are dispensable in the native state. Such arm-dependent folding can be abrogated by mutations associated with β-cell dysfunction and neonatal diabetes mellitus.     

An obligatory intermediate in the folding pathway of cytochrome c552 from Hydrogenobacter thermophilus

The folding mechanism of many proteins involves the population of partially organized structures en route to the native state. Identification and characterization of these intermediates is particularly difficult, as they are often only transiently populated and may play different mechanistic roles, being either on-pathway productive species or off-pathway kinetic traps. Following different spectroscopic probes, and employing state-of-the-art kinetic analysis, we present evidence that the folding mechanism of the thermostable cytochrome c552 from Hydrogenobacter thermophilus does involve the presence of an elusive, yet compact, on-pathway intermediate. Characterization of the folding mechanism of this cytochrome c is particularly interesting for the purpose of comparative folding studies, because H. thermophilus cytochrome c552 shares high sequence identity and structural homology with its homologue from the mesophilic bacterium Pseudomonas aeruginosa cytochrome c551, which refolds through a broad energy barrier without the accumulation of intermediates. Analysis of the folding kinetics and correlation with the three-dimensional structure add new evidence for the validity of a consensus folding mechanism in the cytochrome c family.     

RNA and protein folding: common themes and variations

Visualizing the navigation of an ensemble of unfolded molecules through the bumpy energy landscape in search of the native state gives a pictorial view of biomolecular folding. This picture, when combined with concepts in polymer theory, provides a unified theory of RNA and protein folding. Just as for proteins, the major folding free energy barrier for RNA scales sublinearly with the number of nucleotides, which allows us to extract the elusive prefactor for RNA folding. Several folding scenarios can be anticipated by considering variations in the energy landscape that depend on sequence, native topology, and external conditions. RNA and protein folding mechanism can be described by the kinetic partitioning mechanism (KPM) according to which a fraction (Phi) of molecules reaches the native state directly, whereas the remaining fraction gets kinetically trapped in metastable conformations. For two-state folders Phi approximately 1. Molecular chaperones are recruited to assist protein folding whenever Phi is small. We show that the iterative annealing mechanism, introduced to describe chaperonin-mediated folding, can be generalized to understand protein-assisted RNA folding. The major differences between the folding of proteins and RNA arise in the early stages of folding. For RNA, folding can only begin after the polyelectrolyte problem is solved, whereas protein collapse requires burial of hydrophobic residues. Cross-fertilization of ideas between the two fields should lead to an understanding of how RNA and proteins solve their folding problems.     

8-anilino-1-naphthalene sulfonic acid (ANS) induces folding of acid unfolded cytochrome c to molten globule state as a result of electrostatic interactions.

Hydrophobic interaction of 8-anilino-1-naphthalene sulfonic acid (ANS) with proteins is one of the widely used methods for characterizing/detecting partially folded states of proteins. We have carried out a systematic investigation on the effect of ANS, a charged hydrophobic fluorescent dye, on structural properties of acid-unfolded horse heart cytochrome c at pH 2.0 by a combination of optical methods and electrospray ionization mass spectroscopy (ESI MS). ANS was found to induce, a secondary structure similar to native protein and quenching of fluorescence of tryptophan residue, in the acid-unfolded protein. However, the tertiary structure was found to be disrupted thus indicating that ANS stabilizes a molten globule state in acid-unfolded protein. To understand the mechanism of ANS-induced folding of acid-unfolded cytochrome c, comparative ESI MS, soret absorption, and tryptophan fluorescence studies using nile red, a neutral hydrophobic dye, and ANS were carried out. These studies suggested that, at low pH, electrostatic interactions between negatively charged ANS molecules and positively charged amino acid residues present in acid-unfolded cytochrome c are probably responsible for ANS-induced folding of acid-unfolded protein to partially folded compact state or molten globule state. This is the first experimental demonstration of ANS induced folding of unfolded protein and puts to question the usefulness of ANS for characterization/determination of partially folded intermediates of proteins observed under low pH conditions.     

Protein folding stabilities are a major determinant of oxidation rates for buried methionine residues

The oxidation of protein-bound methionines to form methionine sulfoxides has a broad range of biological ramifications, making it important to delineate factors that influence methionine oxidation rates within a given protein. This is especially important for biopharmaceuticals, where oxidation can lead to deactivation and degradation. Previously, neighboring residue effects and solvent accessibility have been shown to impact the susceptibility of methionine residues to oxidation. In this study, we provide proteome-wide evidence that oxidation rates of buried methionine residues are also strongly influenced by the thermodynamic folding stability of proteins. We surveyed the Escherichia coli proteome using several proteomic methodologies and globally measured oxidation rates of methionine residues in the presence and absence of tertiary structure, as well as the folding stabilities of methionine-containing domains. These data indicated that buried methionines have a wide range of protection factors against oxidation that correlate strongly with folding stabilities. Consistent with this, we show that in comparison to E. coli, the proteome of the thermophile Thermus thermophilus is significantly more stable and thus more resistant to methionine oxidation. To demonstrate the utility of this correlation, we used native methionine oxidation rates to survey the folding stabilities of E. coli and T. thermophilus proteomes at various temperatures and propose a model that relates the temperature dependence of the folding stabilities of these two species to their optimal growth temperatures. Overall, these results indicate that oxidation rates of buried methionines from the native state of proteins can be used as a metric of folding stability.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Modeling the role of disulfide bonds in protein folding: Entropic barriers and pathways

The role of disulfide bonds in directing protein folding is studied using lattice models. We find that the stability and the specificity of the disulfide bond interactions play quite different roles in the folding process: Under some conditions, the stability decreases the overall rate of folding; the specificity, however, by yielding a simpler connectivity of intermediates, always increases the rate of folding. This conclusion is intimately related to the selection mechanism entailed by entropic driving forces, such as the loop formation probability, and entropic barriers separating the native and the many native-like metastable states. The folding time is found to be a minimum for a certain range of the effective disulfide bond interaction. Examination of a model, which allows for the formation of disulfide bonded intermediates, suggests that folding proceeds via a threestage multiple pathways kinetics. We show that there are pathways to the native state involving only native-like intermediates, as well as those that are mediated by nonnative intermediates. These findings are interpreted in terms of the appropriate energy landscape describing the barriers connecting low energy conformations. The consistency of our conclusions with several experimental studies is also discussed. © 1995 Wiley-Liss, Inc.