Kinetic evidence for an on-pathway intermediate in the folding of cytochrome c

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding-unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.

Kinetically trapped metastable intermediate of a disulfide-deficient mutant of the starch-binding domain of glucoamylase

Refolding of a thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase was investigated using differential scanning calorimetry, isothermal titration calorimetry, CD, and (1)H NMR. When the protein solution was rapidly cooled from a higher temperature, a kinetic intermediate was formed during refolding. The intermediate was unexpectedly stable compared with typical folding intermediates that have short half-lives. It was shown that this intermediate contained substantial secondary structure and tertiary packing and had the same binding ability with beta-cyclodextrin as the native state, suggesting that the intermediate is highly-ordered and native-like on the whole. These characteristics differ from those of partially folded intermediates such as molten globule states. Far-UV CD spectra showed that the secondary structure was once disrupted during the transition from the intermediate to the native state. These results suggest that the intermediate could be an off-pathway type, possibly a misfolded state, that has to undergo unfolding on its way to the native state.

Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

The beta alpha-repeat class of proteins, represented by the (beta alpha)(8) barrel and the alpha/beta/alpha sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of beta alpha-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Gō-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY.

NMR analysis of partially folded states and persistent structure in the alpha subunit of tryptophan synthase: implications for the equilibrium folding mechanism of a 29-kDa TIM barrel protein

Structural insights into the equilibrium folding mechanism of the alpha subunit of tryptophan synthase (alpha TS) from Escherichia coli, a (beta alpha)(8) TIM barrel protein, were obtained with a pair of complementary nuclear magnetic resonance (NMR) spectroscopic techniques. The secondary structures of rare high-energy partially folded states were probed by native-state hydrogen-exchange NMR analysis of main-chain amide hydrogens. 2D heteronuclear single quantum coherence NMR analysis of several (15)N-labeled nonpolar amino acids was used to probe the side chains involved in stabilizing a highly denatured intermediate that is devoid of secondary structure. The dynamic broadening of a subset of isoleucine and leucine side chains and the absence of protection against exchange showed that the highest energy folded state on the free-energy landscape is stabilized by a hydrophobic cluster lacking stable secondary structure. The core of this cluster, centered near the N-terminus of alpha TS, serves as a nucleus for the stabilization of what appears to be nonnative secondary structure in a marginally stable intermediate. The progressive decrease in protection against exchange from this nucleus toward both termini and from the N-termini to the C-termini of several beta-strands is best described by an ensemble of weakly coupled conformers. Comparison with previous data strongly suggests that this ensemble corresponds to a marginally stable off-pathway intermediate that arises in the first few milliseconds of folding and persists under equilibrium conditions. A second, more stable intermediate, which has an intact beta-barrel and a frayed alpha-helical shell, coexists with this marginally stable species. The conversion of the more stable intermediate to the native state of alpha TS entails the formation of a stable helical shell and completes the acquisition of the tertiary structure.

Plasticity of the protein folding landscape: switching between on- and off-pathway intermediates

Proteins may fold via parallel routes partitioned by the relative effect of solvent conditions on the relevant transition states. Thus, intermediates may or may not necessarily be obligatory species accumulated during the folding process, but rather kinetic traps due to the ruggedness of the folding landscape. Implicit in this view is the notion of plasticity of the folding pathway: proteins can be rerouted through the energy landscape by mutational, topological or solvent perturbations. Our work was specifically aimed to the experimental identification of a switch in the folding mechanism of a c-type cytochrome from the thermophilic bacterium Hydrogenobacter thermophilus (HT cyt c(552)) induced by acidic conditions. We present evidence that, by destabilizing the relevant transition state, the native state of HT cyt c(552) can be reached along alternative folding routes, which may involve an off-pathway intermediate.

Identification and characterization of protein folding intermediates

In order to understand the mechanism by which a polypeptide chain folds into its functionally active native state it is necessary to characterize in detail all the species accumulated along the pathway. The elusive nature of protein folding intermediates poses their identification and characterization as an extremely difficult task in the protein folding field. In the case of small single domain proteins, the direct measurement of the thermodynamics and structural parameters of protein folding intermediates has provided new insights on the nature of the forces involved in the stabilization of nascent protein structures. Here we summarize some of the experimental approaches aimed at the detection and characterization of folding intermediates along with a discussion of some general structural features emerging from these studies.

A unified mechanism for protein folding: predetermined pathways with optional errors

There is a fundamental conflict between two different views of how proteins fold. Kinetic experiments and theoretical calculations are often interpreted in terms of different population fractions folding through different intermediates in independent unrelated pathways (IUP model). However, detailed structural information indicates that all of the protein population folds through a sequence of intermediates predetermined by the foldon substructure of the target protein and a sequential stabilization principle. These contrary views can be resolved by a predetermined pathway--optional error (PPOE) hypothesis. The hypothesis is that any pathway intermediate can incorporate a chance misfolding error that blocks folding and must be reversed for productive folding to continue. Different fractions of the protein population will then block at different steps, populate different intermediates, and fold at different rates, giving the appearance of multiple unrelated pathways. A test of the hypothesis matches the two models against extensive kinetic folding results for hen lysozyme which have been widely cited in support of independent parallel pathways. The PPOE model succeeds with fewer fitting constants. The fitted PPOE reaction scheme leads to known folding behavior, whereas the IUP properties are contradicted by experiment. The appearance of a conflict with multipath theoretical models seems to be due to their different focus, namely on multitrack microscopic behavior versus cooperative macroscopic behavior. The integration of three well-documented principles in the PPOE model (cooperative foldons, sequential stabilization, optional errors) provides a unifying explanation for how proteins fold and why they fold in that way.

A PDZ domain recapitulates a unifying mechanism for protein folding

A unifying view has been recently proposed according to which the classical diffusion-collision and nucleation-condensation models may represent two extreme manifestations of an underlying common mechanism for the folding of small globular proteins. We report here the characterization of the folding process of the PDZ domain, a protein that recapitulates the three canonical steps involved in this unifying mechanism, namely: (i) the early formation of a weak nucleus that determines the native-like topology of a large portion of the structure, (ii) a global collapse of the entire polypeptide chain, and (iii) the consolidation of the remaining partially structured regions to achieve the native state conformation. These steps, which are clearly detectable in the PDZ domain investigated here, may be difficult to distinguish experimentally in other proteins, which would thus appear to follow one of the two limiting mechanisms. The analysis of the (un)folding kinetics for other three-state proteins (when available) appears consistent with the predictions ensuing from this unifying mechanism, thus providing a powerful validation of its general nature.

The folding pathway of T4 lysozyme: an on-pathway hidden folding intermediate

T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. In earlier studies, an amide hydrogen/deuterium exchange pulse-labeling experiment detected a stable submillisecond intermediate that accumulates before the rate-limiting transition state. It involves the formation of structures in both the N and C-terminal regions. However, a native-state hydrogen exchange experiment subsequently detected an equilibrium intermediate that only involves the formation of the C-terminal domain. Here, using stopped-flow circular dichroism and fluorescence, amide hydrogen exchange-folding competition, and protein engineering methods, we re-examined the folding pathway of T4-lysozyme. We found no evidence for the existence of a stable folding intermediate before the rate-limiting transition state at neutral pH. In addition, using native-state hydrogen exchange-directed protein engineering, we created a mimic of the equilibrium intermediate. We found that the intermediate mimic folds with the same rate as the wild-type protein, suggesting that the equilibrium intermediate is an on-pathway intermediate that exists after the rate-limiting transition state.

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.

Unveiling a hidden folding intermediate in c-type cytochromes by protein engineering

Several investigators have highlighted a correlation between the basic features of the folding process of a protein and its topology, which dictates the folding pathway. Within this conceptual framework we proposed that different members of the cytochrome c (cyt c) family share the same folding mechanism, involving a consensus partially structured state. Pseudomonas aeruginosa cyt c(551) (Pa cyt c(551)) folds via an apparent two-state mechanism through a high energy intermediate. Here we present kinetic evidence demonstrating that it is possible to switch its folding mechanism from two to three state, stabilizing the high energy intermediate by rational mutagenesis. Characterization of the folding kinetics of one single-site mutant of the Pa cyt c(551) (Phe(7) to Ala) indeed reveals an additional refolding phase and a fast unfolding process which are explained by the accumulation of a partially folded species. Further kinetic analysis highlights the presence of two parallel processes both leading to the native state, suggesting that the above mentioned species is a non obligatory on-pathway intermediate. Determination of the crystallographic structure of F7A shows the presence of an extended internal cavity, which hosts three "bound" water molecules and a H-bond in the N-terminal helix, which is shorter than in the wild type protein. These two features allow us to propose a detailed structural interpretation for the stabilization of the native and especially the intermediate states induced by a single crucial mutation. These results show how protein engineering, x-ray crystallography and state-of-the-art kinetics concur to unveil a folding intermediate and the structural determinants of its stability.

Kinetic folding mechanism of PDZ2 from PTP-BL

PDZ domains represent a large family of protein-interaction modules associated with a variety of unrelated proteins with different functions. We report a complete characterization of the kinetic folding mechanism of a fluorescent variant of PDZ2 from PTP-BL, investigated under a variety of different experimental conditions. For this purpose, we engineered a fluorescent variant of this protein Y43W (called pseudo-wild-type, pWT43). The results suggest the presence of a high-energy intermediate in the folding of PDZ2, as revealed by a pronounced non-linear dependence of the unfolding rate constant on denaturant concentration. Such an intermediate may or may not be detectable depending on the experimental conditions, giving rise to apparent two-state folding under stabilizing conditions (e.g. in the presence of sodium sulfate). Interestingly, even under these conditions, three-state folding can be restored by selectively destabilizing the native-like rate-limiting barrier by one specific mutation (V44A). Finally, we show that data taken on pWT43 under different experimental conditions (e.g. different pH values from 2.1 to 8.0 or in the presence of a stabilizing salt) and also data on a site-directed conservative mutant can be rationalized in terms of a simple reaction scheme involving a single set of intermediates and transition states.

A protein folding pathway with multiple folding intermediates at atomic resolution

Using native-state hydrogen-exchange-directed protein engineering and multidimensional NMR, we determined the high-resolution structure (rms deviation, 1.1 angstroms) for an intermediate of the four-helix bundle protein: Rd-apocytochrome b562. The intermediate has the N-terminal helix and a part of the C-terminal helix unfolded. In earlier studies, we also solved the structures of two other folding intermediates for the same protein: one with the N-terminal helix alone unfolded and the other with a reorganized hydrophobic core. Together, these structures provide a description of a protein folding pathway with multiple intermediates at atomic resolution. The two general features for the intermediates are (i) native-like backbone topology and (ii) nonnative side-chain interactions. These results have implications for important issues in protein folding studies, including large-scale conformation search, -value analysis, and computer simulations.

Protein misfolding: optional barriers, misfolded intermediates, and pathway heterogeneity

To investigate the character and role of misfolded intermediates in protein folding, a recombinant cytochrome c without the normally blocking histidine to heme misligation was studied. Folding remains heterogeneous as in the wild-type protein. Half of the population folds relatively rapidly to the native state in a two-state manner. The other half collapses (fluorescence quenching) and forms a full complement of helix (CD) with the same rate and denaturant dependence as the fast folding fraction but then is blocked and reaches the native structure (695nm absorbance) much more slowly. The factors that transiently block folding are not intrinsic to the folding process but depend on ambient conditions, including protein aggregation (f(concentration)), N terminus to heme misligation (f(pH)), and proline mis-isomerization (f(U state equilibration time)). The misfolded intermediate populated by the slowly folding fraction was characterized by hydrogen exchange pulse labeling. It is very advanced with all of the native-like elements fairly stably formed but not the final Met80-S to heme iron ligation, similar to a previously studied molten globule form induced by low pH. To complete final native state acquisition, some small back unfolding is required (error repair) but the misfolded intermediate does not revisit the U state before proceeding to N. These properties show that the intermediate is a normal on-pathway form that contains, in addition, adventitious misfolding errors that transiently block its forward progress. Related observations for other proteins (partially misfolded intermediates, pathway heterogeneity) might be similarly explained in terms of the optional insertion of error-dependent barriers into a classical folding pathway.

Effect of Axial Coordination on the Kinetics of Assembly and Folding of the Two Halves of Horse Heart Cytochrome c

The kinetics of the assembly of two complementary fragments of oxidized horse heart cytochrome c (cyt c), namely the heme-containing fragment-(1-56) and the fragment-(57-104), have been characterized at different pH values. At neutral pH the fragment-(1-56) is hexacoordinated and has two histidines axially ligated to the heme-Fe(III) (Santucci, R., Fiorucci, L., Sinibaldi, F., Polizio, F., Desideri, A., and Ascoli, F. (2000) Arch. Biochem. Biophys. 379, 331-336), thus mimicking what occurs in the folding intermediate of cyt c. The kinetics of the formation of the complex between the two fragments are characterized at pH 7.0 by a slow rate constant that is independent of the concentration of the reactants; conversely, at a low pH the kinetics are much faster and depend on the concentration of the fragments. This behavior suggests that the rate-limiting step observed in the recombination process of the fragments at neutral pH (that leads to the final coordination of Met-80) has to be ascribed to the detachment of the "misligated" histidine. Thus, the faster recombination rate at a low pH can be related to the fact that histidine is protonated and not able to coordinate to the metal. Furthermore, the independence of the rate constant on the concentration of the reactants observed at pH 7.0 can be accounted for by the occurrence of a conformational transition, which takes place immediately after the two fragments collapse together, likely simulating what induces the detachment of the misligated histidine during cytochrome folding.

The 40s ?-loop plays a critical role in the stability and the alkaline conformational transition of cytochrome c

The structural and redox properties of a non-covalent complex reconstituted upon mixing two non-contiguous fragments of horse cytochrome c, the residues 1-38 heme-containing N-fragment with the residues 57-104 C-fragment, have been investigated. With respect to native cyt c, the complex lacks a segment of 18 residues, corresponding, in the native protein, to an omega (Omega)-loop region. The fragment complex shows compact structure, native-like alpha-helix content but a less rigid atomic packing and reduced stability with respect to the native protein. Structural heterogeneity is observed at pH 7.0, involving formation of an axially misligated low-spin species and consequent partial displacement of Met80 from the sixth coordination position of the heme-iron. Spectroscopic data suggest that a lysine (located in the Met80-containing loop, namely Lys72, Lys73, or Lys79) replaces the methionine residue. The residues 1-38/57-104 fragment complex shows an unusual biphasic alkaline titration characterized by a low (p K(a1)=6.72) and a high p K(a)-associated state transition (p K(a2)=8.56); this behavior differs from that of native cyt c, which shows a monophasic alkaline transition (p K(a)=8.9). The data indicate that the 40s Omega-loop plays an important role in the stability of cyt c and in ensuring a correct alkaline conformational transition of the protein.

How Cytochrome c Folds, and Why: Submolecular Foldon Units and their Stepwise Sequential Stabilization

Native state hydrogen exchange experiments have shown that the cytochrome c (Cyt c) protein consists of five cooperative folding-unfolding units, called foldons. These are named, in the order of increasing unfolding free energy, the nested-Yellow, Red, Yellow, Green, and Blue foldons. Previous results suggest that these units unfold in a stepwise sequential way so that each higher energy partially unfolded form includes all of the previously unfolded lower free energy units. If this is so, then selectively destabilizing any given foldon should equally destabilize each subsequent unfolding step above it in the unfolding ladder but leave the lower ones before it unaffected. To perform this test, we introduced the mutation Glu62Gly, which deletes a salt link in the Yellow unit and destabilizes the protein by 0.8 kcal/mol. Native state hydrogen exchange and other experiments show that the stability of the Yellow unit and the states above it in the free energy ladder are destabilized by about the same amount while the lower lying states are unaffected. These results help to confirm the sequential stepwise nature of the Cyt c unfolding pathway and therefore a similar refolding pathway. The steps in the pathway are dictated by the concerted folding-unfolding property of the individual unit foldons; the order of steps is determined by the sequential stabilization of progressively added foldons in the native context. Much related information for Cyt c strongly conforms with this mechanism. Its generality is supported by available information for other proteins.

Demonstration of a low-energy on-pathway intermediate in a fast-folding protein by kinetics, protein engineering, and simulation

It is controversial whether fast-folding proteins can form productive on-pathway intermediates that are more stable than the denatured state because noncovalent intermediates are usually evanescent. Here, we apply the classical criteria for the existence of intermediates: namely, the intermediates form and react rapidly enough to be on pathway and they can be isolated and characterized. The folding of the 71-residue, mainly alpha-helical FF domain from human HYPA/FBP11 fulfills these classical criteria, as was found for Im7. The FF domain folds in two phases, one on the micros and the other on the ms time scale. An engineered mutant folds only to a partly folded state, with some 20-40% of the native helical content. The kinetic properties of the mutant are identical to those found for the fast phase of the wild-type protein, and it is likely that the mutant folds just to the intermediate state. A full kinetic analysis of the folding of wild-type protein, using the amplitudes of its native and denatured states and the observed values for the mutant, rules out an off-pathway scheme but fits an on-pathway scheme, with a low energy intermediate that is modeled by the mutant. The experimental proof benchmarks a molecular dynamics method that identifies an obligatory intermediate observed in multiple simulations. The conformational space defining this intermediate is visited several times in the simulations, leading to high populations consistent with the presence of a low energy intermediate.

Intimate View of a Kinetic Protein Folding Intermediate: Residue-resolved Structure, Interactions, Stability, Folding and Unfolding Rates, Homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.

Exploring the cytochrome c folding mechanism: cytochrome c552 from thermus thermophilus folds through an on-pathway intermediate

Understanding the role of partially folded intermediate states in the folding mechanism of a protein is a crucial yet very difficult problem. We exploited a kinetic approach to demonstrate that a transient intermediate of a thermostable member of the widely studied cytochrome c family (cytochrome c552 from Thermus thermophilus) is indeed on-pathway. This is the first clear indication of an obligatory intermediate in the folding mechanism of a cytochrome c. The fluorescence properties of this intermediate demonstrate that the relative position of the heme and of the only tryptophan residue cannot correspond to their native orientation. Based on an analysis of the three-dimensional structure of cytochrome c552, we propose an interpretation of the data which explains the residual fluorescence of the intermediate and is consistent with the established role played by some conserved interhelical interactions in the folding of other members of this family. A limited set of topologically conserved contacts may guide the folding of evolutionary distant cytochromes c through the same partially structured state, which, however, can play different kinetic roles, acting either as an intermediate or a transition state.

Folding units govern the cytochrome c alkaline transition

The alkaline transition of cytochrome c is a model for protein structural switching in which the normal heme ligand is replaced by another group. Stopped flow data following a jump to high pH detect two slow kinetic phases, suggesting two rate-limiting structure changes. Results described here indicate that these events are controlled by the same structural unfolding reactions that account for the first two steps in the reversible unfolding pathway of cytochrome c. These and other results show that the cooperative folding-unfolding behavior of protein foldons can account for a variety of functional activities in addition to determining folding pathways.

Cooperative omega loops in cytochrome c: role in folding and function

Hydrogen exchange experiments under slow exchange conditions show that an omega loop in cytochrome c (residues 40-57) acts as a cooperative unfolding/refolding unit under native conditions. This unit behavior accounts for an initial step on the unfolding pathway, a final step in refolding, and a number of other structural, functional and evolutionary properties.

Hidden intermediates and levinthal paradox in the folding of small proteins

It has long been suggested that existence of partially folded intermediates may be essential for proteins to fold in a biologically meaningful time scale. Although partially folded intermediates have been commonly observed in larger proteins, they are generally not detectable in the kinetic folding of smaller proteins (approximately 100 amino acids or less). Recent native-state hydrogen exchange studies suggest that partially folded intermediates may exist behind the rate-limiting transition state in small proteins and evade detection by conventional kinetic methods.

Kinetics of cytochrome C folding: atomically detailed simulations

The vast range of time scales (from nanoseconds to seconds) during protein folding is a challenge for experiments and computations. To make concrete predictions on folding mechanisms, atomically detailed simulations of protein folding, using potentials derived from chemical physics principles, are desired. However, due to their computational complexity, straightforward molecular dynamics simulations of protein folding are impossible today. An alternative algorithm is used that makes it possible to compute approximate atomically detailed long time trajectories (the Stochastic Difference Equation in Length). This algorithm is used to compute 26 atomically detailed folding trajectories of cytochrome c (a millisecond process). The early collapse of the protein chain (with marginal formation of secondary structure), and the earlier formation of the N and C helices (compare to the 60's helix) are consistent with the experiment. The existence of an energy barrier upon entry to the molten globule is examined as well. In addition to (favorable) comparison to experiments, we show that non-native contacts drive the formation of the molten globule. In contrast to popular folding models, the non-native contacts do not form off-pathway kinetic traps in cytochrome c.

Submolecular cooperativity produces multi-state protein unfolding and refolding

Hydrogen exchange experiments show that cytochrome c and other proteins under native conditions reversibly unfold in a multi-step manner. The step from one intermediate to the next is determined by the intrinsically cooperative nature of secondary structural elements, which is retained in the native protein. Folding uses the same pathway in the reverse direction, moving from the unfolded to the native state through relatively discrete intermediate forms by the sequential addition of native-like secondary structural units.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

Cytochrome c folding pathway: kinetic native-state hydrogen exchange

Native-state hydrogen exchange experiments under EX1 conditions can distinguish partially unfolded intermediates by their formation rates and identify the amide hydrogens exposed and protected in each. Results obtained define a cytochrome c intermediate seen only poorly before and place it early on the major unfolding pathway. Four distinct unfolding steps are found to be kinetically ordered in the same pathway sequence inferred before.

Structural restoration of inactive recombinant fish growth hormones by chemical chaperonin and solvent restraint approaches

Recombinant proteins may undergo conformational distortion, leading to aggregation and loss of function, when they are expressed in heterologous systems. The structural and functional restoration of such inactive proteins is highly desirable. We have over-expressed recombinant growth hormones from the fish ayu (Plecoglossus altivelis) and yellow grouper (Epinephelus awoara) by a pET expression system. Both recombinant proteins accumulate as insoluble form in Escherichia coli. We refolded these inactive proteins into the active form using a stepwise refolding process with a dilute denaturing agent as a steric blocker and chemical chaperonin. Optical characterization showed that stable folding intermediates with a helical conformation can be detected in the molten globule state. Moreover, the function of restored recombinant growth hormones was demonstrated by its ability to stimulate proliferation in zebrafish liver cells.

Structural transitions in neutral and charged proteins in vacuo

In vacuo proteins provide a simple laboratory to explore the roles of sequence, temperature, charge state, and initial configuration in protein folding. Moreover, by the very absence of solvent, the study of anhydrous proteins in vacuo will also help us to understand specific environmental effects. From the experimental viewpoint, these systems are now beginning to be characterized at low resolution. Molecular dynamics (MD) simulations, in combination with tools for protein shape analysis, can complement experiments and provide further insights on the folding-unfolding transitions of these proteins. We review some aspects of this issue by using the results from a detailed MD study of hen egg-white lysozyme. For lysozyme ions, unfolding can be triggered by Coulombic repulsion. In neutral lysozyme, unfolding can be induced by centrifugal forces and also by weakening the monomer-monomer interaction. In both cases, the resulting unfolded transients can be used as initial configurations for relaxation dynamics. All trajectories are analyzed in terms of global molecular shape features of the backbone, including its anisometry and chain entanglement complexity. This strategy allows us to quantify separately the degree of polymer collapse and the evolution of large-scale folding features. Using these last two notions, we discuss some basic questions regarding the nature of the accessible paths associated with unfolding from, and refolding into, compact conformers.

Folding of horse cytochrome c in the reduced state

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10 degrees C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded DeltaG(H2O)=18.8(+/-1.45) kcal mol(-1), and C(m)=5.1(+/-0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of -164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U<-->N) of 80-90 % of unfolded molecules with a time constant in the range 190-550 micros estimated by linear extrapolation and model calculations. The remaining 10-20 % of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I<-->U<-->N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10-20 % of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pK(a) value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is approximately 300(+/-125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I<-->U<-->N scheme have yielded the microscopic rate constants.

Folding and signaling share the same pathway in a photoreceptor

The photoreceptor photoactive yellow protein (PYP) was used as a model system to study receptor activation and protein folding. Refolding was studied by stopped-flow absorbance spectroscopy for PYP with either a trans or a cis chromophore. Chromophore trans to cis isomerization, the mechanism of light detection by PYP, greatly affects the protein folding process. When the cis chromophore is present, refolding from the unfolded state proceeds through the putative signaling state of PYP as an on-pathway intermediate. In addition, moderate denaturant concentrations result in the specific unfolding of the signaling state of PYP. Thus, the signaling state is common to the pathways of folding and signaling. This result provides an avenue for the study of protein folding. We demonstrate how this approach can be used to establish whether a folding intermediate is on-pathway or off-pathway. The results also reveal transient partial unfolding as a molecular mechanism for signaling.

Snapshots of protein folding. A study on the multiple transition state pathway of cytochrome c(551) from Pseudomonas aeruginosa

Cytochrome c(551) (cyt c(551)) from Pseudomonas aeruginosa is a small protein (82 residues) that folds via a three-state pathway with the accumulation in the microsecond time-range of a compact collapsed intermediate. The presence of a single His residue, at position 16, permits the study of the refolding at pH 7.0 in the absence of miscoordination events. Here, we report on folding kinetics in the millisecond time-range as a function of urea under different pH conditions. Analysis of this process (over-and-above proline cis-trans isomerization) at pH 7.0, suggests the existence of a multiple transition state pathway in which we postulate three transition states. Taking advantage of site-directed mutagenesis we propose that the first "unfolded-like" transition state (t(1)) originates from the electrostatic properties of the collapsed state, while the second transition state (t(2)) involves the interaction between the N and C-terminal helices and is stabilized by the salt bridge between Lys10 and Glu70 ( approximately 1 kcal mol(-1)). Our results suggest that, contrary to other cytochromes c, the roll-over effect observed for cyt c(551) at low denaturant concentration can be interpreted in terms of a broad energy barrier without population of any intermediates. The third and more "native-like" transition state (M) can be associated with the breaking/formation of the Fe(3+)-Met61 bond. This strong interaction is stabilized by the hydrogen bond between Trp56 and heme propionate 17 (HP-17) as suggested by the increase in the unfolding rate at high denaturant concentration of the Trp56Phe site-directed mutant.

Protein folding intermediates and pathways studied by hydrogen exchange

In order to solve the immensely difficult protein-folding problem, it will be necessary to characterize the barriers that slow folding and the intermediate structures that promote it. Although protein-folding intermediates are not accessible to the usual structural studies, hydrogen exchange (HX) methods have been able to detect and characterize intermediates in both kinetic and equilibrium modes--as transient kinetic folding intermediates on a subsecond time scale, as labile equilibrium molten globule intermediates under destabilizing conditions, and as infinitesimally populated intermediates in the high free-energy folding landscape under native conditions. Available results consistently indicate that protein-folding landscapes are dominated by a small number of discrete, metastable, native-like partially unfolded forms (PUFs). The PUFs appear to be produced, one from another, by the unfolding and refolding of the protein's intrinsically cooperative secondary structural elements, which can spontaneously create stepwise unfolding and refolding pathways. Kinetic experiments identify three kinds of barrier processes: (a) an initial intrinsic search-nucleation-collapse process that prepares the chain for intermediate formation by pinning it into a condensed coarsely native-like topology; (b) smaller search-dependent barriers that put the secondary structural units into place; and (c) optional error-dependent misfold-reorganization barriers that can cause slow folding, intermediate accumulation, and folding heterogeneity. These conclusions provide a coherent explanation for the grossly disparate folding behavior of different globular proteins in terms of distinct folding pathways.

An amino acid code for protein folding

Direct structural information obtained for many proteins supports the following conclusions. The amino acid sequences of proteins can stabilize not only the final native state but also a small set of discrete partially folded native-like intermediates. Intermediates are formed in steps that use as units the cooperative secondary structural elements of the native protein. Earlier intermediates guide the addition of subsequent units in a process of sequential stabilization mediated by native-like tertiary interactions. The resulting stepwise self-assembly process automatically constructs a folding pathway, whether linear or branched. These conclusions are drawn mainly from hydrogen exchange-based methods, which can depict the structure of infinitesimally populated folding intermediates at equilibrium and kinetic intermediates with subsecond lifetimes. Other kinetic studies show that the polypeptide chain enters the folding pathway after an initial free-energy-uphill conformational search. The search culminates by finding a native-like topology that can support forward (native-like) folding in a free-energy-downhill manner. This condition automatically defines an initial transition state, the search for which sets the maximum possible (two-state) folding rate. It also extends the sequential stabilization strategy, which depends on a native-like context, to the first step in the folding process. Thus the native structure naturally generates its own folding pathway. The same amino acid code that translates into the final equilibrium native structure-by virtue of propensities, patterning, secondary structural cueing, and tertiary context-also produces its kinetic accessibility.

Integrity of thermus thermophilus cytochrome c552 synthesized by Escherichia coli cells expressing the host-specific cytochrome c maturation genes, ccmABCDEFGH: biochemical, spectral, and structural characterization of the recombinant protein

We describe the design of Escherichia coli cells that synthesize a structurally perfect, recombinant cytochrome c from the Thermus thermophilus cytochrome c552 gene. Key features are (1) construction of a plasmid-borne, chimeric cycA gene encoding an Escherichia coli-compatible, N-terminal signal sequence (MetLysIleSerIleTyrAlaThrLeu AlaAlaLeuSerLeuAlaLeuProAlaGlyAla) followed by the amino acid sequence of mature Thermus cytochrome c552; and (2) coexpression of the chimeric cycA gene with plasmid-borne, host-specific cytochrome c maturation genes (ccmABCDEFGH). Approximately 1 mg of purified protein is obtained from 1 L of culture medium. The recombinant protein, cytochrome rsC552, and native cytochrome c552 have identical redox potentials and are equally active as electron transfer substrates toward cytochrome ba3, a Thermus heme-copper oxidase. Native and recombinant cytochromes c were compared and found to be identical using circular dichroism, optical absorption, resonance Raman, and 500 MHz 1H-NMR spectroscopies. The 1.7 A resolution X-ray crystallographic structure of the recombinant protein was determined and is indistinguishable from that reported for the native protein (Than, ME, Hof P, Huber R, Bourenkov GP, Bartunik HD, Buse G, Soulimane T, 1997, J Mol Biol 271:629-644). This approach may be generally useful for expression of alien cytochrome c genes in E. coli.

Cytochrome c(553), a small heme protein that lacks misligation in its unfolded state, folds with rapid two-state kinetics

Cytochrome c(553) (cyt c(553)) from Desulfovibrio vulgaris is a small helical heme protein that displays apparent two-state equilibrium-unfolding behavior. The covalently attached heme is low-spin, ligated by Met and His residues, in the native state but becomes high-spin upon unfolding at pH 7. Here, we show that in contrast to other c-type heme proteins, where misligations in the unfolded states are prominent, cyt c(553) refolding kinetics at pH 7 proceeds rapidly without detectable intermediates. The extrapolated folding rate constant in water for oxidized cyt c(553) matches exactly that predicted from the cyt c(553) native-state topology: 5300 s(-1 )(experimental) versus 5020 s(-1) (predicted). We therefore conclude that the presence of the oxidized cofactor does not affect the intrinsic formation speed of the cyt c(553 )structural motif.

An essential intermediate in the folding of dihydrofolate reductase

The folding of Escherichia coli dihydrofolate reductase was examined at pH 7.8 and 15 degrees C by using stopped-flow fluorescence and absorbance spectroscopies. The formation of a highly fluorescent intermediate occurs with relaxation times ranging between 142 and 343 msec, whereas stopped-flow absorbance spectroscopy using methotrexate binding assays shows a distinct lag phase during these time frames for the native state. The lag in absorbance kinetics and the lack of fast-track folding events indicate that the formation of this ensemble of intermediates is an obligatory step in the folding reaction.

Coupled kinetic traps in cytochrome c folding: His-heme misligation and proline isomerization

The effect of His-heme misligation on folding has been investigated for a triple mutant of yeast iso-2 cytochrome c (N26H,H33N,H39K iso-2). The variant contains a single misligating His residue at position 26, a location at which His residues are found in several cytochrome c homologues, including horse, tuna, and yeast iso-1. The amplitude for fast phase folding exhibits a strong initial pH dependence. For GdnHCl unfolded protein at an initial pH<5, the observed refolding at final pH 6 is dominated by a fast phase (tau(2f)=20 ms, alpha(2f)=90 %) that represents folding in the absence of misligation. For unfolded protein at initial pH 6, folding at final pH 6 occurs in a fast phase of reduced amplitude (alpha(2f) approximately 20 %) but the same rate (tau(2f)=20 ms), and in two slower phases (tau(m)=6-8 seconds, alpha(m) approximately 45 %; and tau(1b)=16-20 seconds, alpha(1b) approximately 35 %). Double jump experiments show that the initial pH dependence of the folding amplitudes results from a slow pH-dependent equilibrium between fast and slow folding species present in the unfolded protein. The slow equilibrium arises from coupling of the His protonation equilibrium to His-heme misligation and proline isomerization. Specifically, Pro25 is predominantly in trans in the unligated low-pH unfolded protein, but is constrained in a non-native cis isomerization state by His26-heme misligation near neutral pH. Refolding from the misligated unfolded form proceeds slowly due to the large energetic barrier required for proline isomerization and displacement of the misligated His26-heme ligand.

Kinetic evidence of an on-pathway intermediate in the folding of lysozyme

By means of a kinetic test, it was demonstrated that one of the folding intermediates (Ialpha) of hen lysozyme with alpha-domain folded and beta-domain unfolded is on the folding pathway under the classical definition. Ialpha folds to the native (N) state directly (unfolded (U) <==> Ialpha <==> N) without having to unfold to U and then refold to N through alternative folding pathways as in Ialpha <==> U <==> N.