Branching in the sequential folding pathway of cytochrome c

Genetically Encoded Fluorescent Probe for Detection of Heme-Induced Conformational Changes in Cytochrome c

Cytochrome c (Cytc) is a key redox protein for energy metabolism and apoptosis in cells. The activation of Cytc is composed of several steps, including its transfer to the mitochondrial membrane, binding to cytochrome c heme lyase (CCHL) and covalent attachment to heme. The spectroscopic methods are often applied to study the structural changes of Cytc. However, they require the isolation of Cytc from cells and have limited availability under physiological conditions. Despite recent studies to elucidate the tightly regulated folding mechanism of Cytc, the role of these events and their association with different conformational states remain elusive. Here, we provide a genetically encoded fluorescence method that allows monitoring of the conformational changes of Cytc upon binding to heme and CCHL. Cerulean and Venus fluorescent proteins attached at the N and C terminals of Cytc can be used to determine its unfolded, intermediate, and native states by measuring FRET amplitude. We found that the noncovalent interaction of heme in the absence of CCHL induced a shift in the FRET signal, indicating the formation of a partially folded state. The higher concentration of heme and coexpression of CCHL gave rise to the recovery of Cytc native structure. We also found that Cytc was weakly associated with CCHL in the absence of heme. As a result, a FRET-based fluorescence approach was demonstrated to elucidate the mechanism of heme-induced Cytc conformational changes with spatiotemporal resolution and can be applied to study its interaction with small molecules and other protein partners in living cells.

Initial Protein Unfolding Events in Ubiquitin, Cytochrome c and Myoglobin Are Revealed with the Use of 213 nm UVPD Coupled to IM-MS

The initial stages of protein unfolding may reflect the stability of the entire fold and can also reveal which parts of a protein can be perturbed, without restructuring the rest. In this work, we couple UVPD with activated ion mobility mass spectrometry to measure how three model proteins start to unfold. Ubiquitin, cytochrome c and myoglobin ions produced via nESI from salty solutions are subjected to UV irradiation pre-mobility separation; experiments are conducted with a range of source conditions which alter the conformation of the precursor ion as shown by the drift time profiles. For all three proteins, the compact structures result in less fragmentation than more extended structures which emerge following progressive in-source activation. Cleavage sites are found to differ between conformational ensembles, for example, for the dominant charge state of cytochrome c [M + 7H]⁷⁺, cleavage at Phe10, Thr19 and Val20 was only observed in activating conditions whilst cleavage at Ala43 is dramatically enhanced. Mapping the photo-cleaved fragments onto crystallographic structures provides insight into the local structural changes that occur as protein unfolding progresses, which is coupled to global restructuring observed in the drift time profiles. Graphical Abstract.

Interrupted Pressure-Jump NMR Experiments Reveal Resonances of On-Pathway Protein Folding Intermediate

Previous pressure-jump NMR experiments on a pressure-sensitized double mutant of ubiquitin showed evidence that its folding occurs via two parallel, comparably efficient pathways: a single barrier and a two-barrier pathway. An interrupted folding NMR experiment is introduced, where for a brief period the pressure is dropped to atmospheric conditions (1 bar), followed by a jump back to high pressure for signal detection. Conventional, forward sampling of the indirect dimension during the low-pressure period correlates the ¹⁵N or ¹³C' chemical shifts of the unfolded protein at 1 bar to the ¹H frequencies of both the unfolded and folded proteins at high pressure. Remarkably, sampling the data of the same experiment in the reverse direction yields the frequencies of proteins present at the end of the low-pressure interval, which include unfolded, intermediate, and folded species. Although the folding intermediate ¹⁵N shifts differ strongly from natively folded protein, its ¹³C' chemical shifts, which are more sensitive probes for secondary structure, closely match those of the folded protein and indicate that the folding intermediate must have a structure that is quite similar to the native state.

A folding nucleus and minimal ATP binding domain of Hsp70 identified by single-molecule force spectroscopy

The folding pathways of large proteins are complex, with many of them requiring the aid of chaperones and others folding spontaneously. Along the folding pathways, partially folded intermediates are frequently populated; their role in the driving of the folding process is unclear. The structures of these intermediates are generally not amenable to high-resolution structural techniques because of their transient nature. Here we employed single-molecule force measurements to scrutinize the hierarchy of intermediate structures along the folding pathway of the nucleotide binding domain (NBD) of Escherichia coli Hsp70 DnaK. DnaK-NBD is a member of the sugar kinase superfamily that includes Hsp70s and the cytoskeletal protein actin. Using optical tweezers, a stable nucleotide-binding competent en route folding intermediate comprising lobe II residues (183-383) was identified as a critical checkpoint for productive folding. We obtained a structural snapshot of this folding intermediate that shows native-like conformation. To assess the fundamental role of folded lobe II for efficient folding, we turned our attention to yeast mitochondrial NBD, which does not fold without a dedicated chaperone. After replacing the yeast lobe II residues with stable E. coli lobe II, the obtained chimeric protein showed native-like ATPase activity and robust folding into the native state, even in the absence of chaperone. In summary, lobe II is a stable nucleotide-binding competent folding nucleus that is the key to time-efficient folding and possibly resembles a common ancestor domain. Our findings provide a conceptual framework for the folding pathways of other members of this protein superfamily.

Structural characterization of cardiolipin-driven activation of cytochrome c into a peroxidase and membrane perturbation

The interaction between cardiolipin (CL) and cytochrome c (cyt-c) results in a gain of function of peroxidase activity by cyt-c. Despite intensive research, disagreements on nature and molecular details of this interaction remain. In particular, it is still not known how the interaction triggers the onset of apoptosis. Enzymatic characterization of peroxidase activity has highlighted the need for a critical threshold concentration of CL, a finding of profound physiological relevance in vivo. Using solution NMR, fluorescence spectroscopy, and in silico modeling approaches we here confirm that full binding of cyt-c to the membrane requires a CL:cyt-c threshold ratio of 5:1. Among three binding sites, the simultaneous binding of two sites, at two opposing sides of the heme, provides a mechanism to open the heme crevice to substrates. This results in "productive binding" in which cyt-c then sequesters CL, inducing curvature in the membrane. Membrane perturbation along with lipid peroxidation, due to interactions of heme/CL acyl chains, initiates the next step in the apoptotic pathway of making the membrane leaky. The third CL binding site while allowing interaction with the membrane, does not cluster CL or induce subsequent events, making this interaction "unproductive".

Methionine ligand lability of homologous monoheme cytochromes c

Direct electrochemical analysis of adsorbed bacterial monoheme cytochromes c has revealed a phenomenological loss of the axial methionine when examined using pyrolytic "edge-plane" graphite (EPG) electrodes. While prior findings have reported that the Met-loss state may be quantitatively understood using the cytochrome c from Hydrogenobacter thermophilus as a model system, here we demonstrate that the formation of the Met-loss state upon EPG electrodes can be observed for a range of cytochrome orthologs. Through an electrochemical comparison of the wild-type proteins from organisms of varying growth temperature optima, we establish that Met-ligand losses at graphite surfaces have similar energetics to the "foldons" for known protein folding pathways. Furthermore, a downward shift in reduction potential to approximately -100 mV vs standard hydrogen electrode was observed, similar to that of the alkaline transition found in mitochondrial cytochromes c. Pourbaix diagrams for the Met-loss forms of each cytochrome, considered here in comparison to mutants where the Met-ligand has been substituted to His or Ala, suggest that the nature of the Met-loss state is distinct from either a His-/aquo- or a bis-His-ligated heme center, yet more closely matches the pKa values found for bis-His-ligated hemes., We find the propensity for adoption of the Met-loss state in bacterial monoheme cytochromes c scales with their overall thermal stability, though not with the specific stability of the Fe-Met bond.

The nature of protein folding pathways

How do proteins fold, and why do they fold in that way? This Perspective integrates earlier and more recent advances over the 50-y history of the protein folding problem, emphasizing unambiguously clear structural information. Experimental results show that, contrary to prior belief, proteins are multistate rather than two-state objects. They are composed of separately cooperative foldon building blocks that can be seen to repeatedly unfold and refold as units even under native conditions. Similarly, foldons are lost as units when proteins are destabilized to produce partially unfolded equilibrium molten globules. In kinetic folding, the inherently cooperative nature of foldons predisposes the thermally driven amino acid-level search to form an initial foldon and subsequent foldons in later assisted searches. The small size of foldon units, ∼ 20 residues, resolves the Levinthal time-scale search problem. These microscopic-level search processes can be identified with the disordered multitrack search envisioned in the "new view" model for protein folding. Emergent macroscopic foldon-foldon interactions then collectively provide the structural guidance and free energy bias for the ordered addition of foldons in a stepwise pathway that sequentially builds the native protein. These conclusions reconcile the seemingly opposed new view and defined pathway models; the two models account for different stages of the protein folding process. Additionally, these observations answer the "how" and the "why" questions. The protein folding pathway depends on the same foldon units and foldon-foldon interactions that construct the native structure.

Structural re-arrangement and peroxidase activation of cytochrome c by anionic analogues of vitamin E, tocopherol succinate and tocopherol phosphate

Cytochrome c is a multifunctional hemoprotein in the mitochondrial intermembrane space whereby its participation in electron shuttling between respiratory complexes III and IV is alternative to its role in apoptosis as a peroxidase activated by interaction with cardiolipin (CL), and resulting in selective CL peroxidation. The switch from electron transfer to peroxidase function requires partial unfolding of the protein upon binding of CL, whose specific features combine negative charges of the two phosphate groups with four hydrophobic fatty acid residues. Assuming that other endogenous small molecule ligands with a hydrophobic chain and a negatively charged functionality may activate cytochrome c into a peroxidase, we investigated two hydrophobic anionic analogues of vitamin E, α-tocopherol succinate (α-TOS) and α-tocopherol phosphate (α-TOP), as potential inducers of peroxidase activity of cytochrome c. NMR studies and computational modeling indicate that they interact with cytochrome c at similar sites previously proposed for CL. Absorption spectroscopy showed that both analogues effectively disrupt the Fe-S(Met(80)) bond associated with unfolding of cytochrome c. We found that α-TOS and α-TOP stimulate peroxidase activity of cytochrome c. Enhanced peroxidase activity was also observed in isolated rat liver mitochondria incubated with α-TOS and tBOOH. A mitochondria-targeted derivative of TOS, triphenylphosphonium-TOS (mito-VES), was more efficient in inducing H2O2-dependent apoptosis in mouse embryonic cytochrome c(+/+) cells than in cytochrome c(-/-) cells. Essential for execution of the apoptotic program peroxidase activation of cytochrome c by α-TOS may contribute to its known anti-cancer pharmacological activity.

Mapping polymerization and allostery of hemoglobin S using point mutations

Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J. Mol. Biol. 2013, 425, 647-661), now applicable to systems with multiple coupled binding sites, such as hemoglobin. First, we predict the relative stabilities of substates and microstates of hemoglobin, which are determined primarily by entropy within our model. Next, we predict the impact of 866 annotated mutations on hemoglobin's oxygen binding equilibrium. We then discuss a subset of 30 mutations that occur in the presence of the sickle cell mutation and whose effects on polymerization have been measured. Seven of these HbS mutations occur in three predicted druggable binding pockets that might be exploited to directly inhibit polymerization; one of these binding pockets is not apparent in the crystal structure, but only in structures generated by AllosMod. For the 30 mutations, we predict that mutation-induced conformational changes within a single tetramer tend not to significantly impact polymerization; instead, these mutations more likely impact polymerization by directly perturbing a polymerization interface. Finally, our analysis of allostery allows us to hypothesize why hemoglobin evolved to have multiple subunits and a persistent low frequency sickle cell mutation.

Investigations of the low-frequency spectral density of cytochrome c upon equilibrium unfolding

The equilibrium unfolding process of ferric horse heart cytochrome c (cyt c), induced by guanidinium hydrochloride (GdHCl), was studied using UV-vis absorption spectroscopy, resonance Raman spectroscopy, and vibrational coherence spectroscopy (VCS). The unfolding process was successfully fit using a three-state model which included the fully folded (N) and unfolded (U) states, along with an intermediate (I) assigned to a Lys bound heme. The VCS spectra revealed for the first time several low-frequency heme modes that are sensitive to cyt c unfolding: γ(a) (~50 cm(-1)), γ(b) (~80 cm(-1)), γ(c) (~100 cm(-1)), and ν(s)(His-Fe-His) at 205 cm(-1). These out-of-plane modes have potential functional relevance and are activated by protein-induced heme distortions. The free energies for the N-I and the I-U transitions at pH 7.0 and 20 °C were found to be 4.6 kcal/M and 11.6 kcal/M, respectively. Imidazole was also introduced to replace the methionine ligand so the unfolding can be modeled as a two-state system. The intensity of the mode γ(b)~80 cm(-1) remains nearly constant during the unfolding process, while the amplitudes of the other low frequency modes track with spectral changes observed at higher frequency. This confirms that the heme deformation changes are coupled to the protein tertiary structural changes that take place upon unfolding. These studies also reveal that damping of the coherent oscillations depends sensitively on the coupling between heme and the surrounding water solvent.

Becoming a peroxidase: cardiolipin-induced unfolding of cytochrome c

Interactions of cytochrome c (cyt c) with a unique mitochondrial glycerophospholipid cardiolipin (CL) are relevant for the protein's function in oxidative phosphorylation and apoptosis. Binding to CL-containing membranes promotes cyt c unfolding and dramatically enhances the protein's peroxidase activity, which is critical in early stages of apoptosis. We have employed a collection of seven dansyl variants of horse heart cyt c to probe the sequence of steps in this functional transformation. Kinetic measurements have unraveled four distinct processes during CL-induced cyt c unfolding: rapid protein binding to CL liposomes; rearrangements of protein substructures with small unfolding energies; partial insertion of the protein into the lipid bilayer; and extensive protein restructuring leading to "open" extended structures. While early rearrangements depend on a hierarchy of foldons in the native structure, the later process of large-scale unfolding is influenced by protein interactions with the membrane surface. The opening of the cyt c structure exposes the heme group, which enhances the protein's peroxidase activity and also frees the C-terminal helix to aid in the translocation of the protein through CL membranes.

Stepwise protein folding at near amino acid resolution by hydrogen exchange and mass spectrometry

The kinetic folding of ribonuclease H was studied by hydrogen exchange (HX) pulse labeling with analysis by an advanced fragment separation mass spectrometry technology. The results show that folding proceeds through distinct intermediates in a stepwise pathway that sequentially incorporates cooperative native-like structural elements to build the native protein. Each step is seen as a concerted transition of one or more segments from an HX-unprotected to an HX-protected state. Deconvolution of the data to near amino acid resolution shows that each step corresponds to the folding of a secondary structural element of the native protein, termed a "foldon." Each folded segment is retained through subsequent steps of foldon addition, revealing a stepwise buildup of the native structure via a single dominant pathway. Analysis of the pertinent literature suggests that this model is consistent with experimental results for many proteins and some current theoretical results. Two biophysical principles appear to dictate this behavior. The principle of cooperativity determines the central role of native-like foldon units. An interaction principle termed "sequential stabilization" based on native-like interfoldon interactions orders the pathway.

Protein dynamics governed by interfaces of high polarity and low packing density

The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.

Methionine ligand lability in bacterial monoheme cytochromes c: an electrochemical study

The direct electrochemical analysis of adsorbed redox active proteins has proven to be a powerful technique in biophysical chemistry, frequently making use of the electrode material pyrolytic "edge-plane" graphite. However, many heme-bearing proteins such as cytochromes c have been also examined systematically at alkanethiol-modified gold surfaces, and previously we reported the characterization of the redox properties of a series of bacterial cytochromes c in a side-by-side comparison of carbon and gold electrode materials. In our prior findings, we reported an unanticipated, low potential (E(m) ∼ -100 mV vs SHE) redox couple that could be analogously observed when a variety of monoheme cytochromes c are adsorbed onto carbon-based electrodes. Here we demonstrate that our prior phenomological data can be understood quantitatively in the loss of the methionine ligand of the heme iron, using the cytochrome c from Hydrogenbacter thermophilum as a model system. Through the comparison of wild-type protein with M61H and M61A mutants, in direct electrochemical analyses conducted as a function of temperature and exogenous ligand concentration, we are able to show that Met-ligated cytochromes c have a propensity to lose their Met ligand at graphite surfaces, and that energetics of this process (6.3 ± 0.2 kJ/mol) is similar to the energies associated with "foldons" of known protein folding pathways.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.

The folding energy landscape and free energy excitations of cytochrome c

The covalently bound heme cofactor plays a dominant role in the folding of cytochrome c. Because of the complicated inorganic chemistry of the heme, some might consider the folding of cytochrome c to be a special case, following principles different from those used to describe the folding of proteins without cofactors. Recent investigations, however, demonstrate that common models describing folding for many proteins work well for cytochrome c when heme is explicitly introduced, generally providing results that agree with experimental observations. In this Account, we first discuss results from simple native structure-based models. These models include attractive interactions between nonadjacent residues only if they are present in the crystal structure at pH 7. Because attractive nonnative contacts are not included in native structure-based models, their energy landscapes can be described as "perfectly funneled". In other words, native structure-based models are energetically guided towards the native state and contain no energetic traps that would hinder folding. Energetic traps are denoted sources of "frustration", which cause specific transient intermediates to be populated. Native structure-based models do, however, include repulsion between residues due to excluded volume. Nonenergetic traps can therefore exist if the chain, which cannot cross over itself, must partially unfold so that folding can proceed. The ability of native structure-based models to capture this kind of motion is partly responsible for their successful predictions of folding pathways for many types of proteins. Models without frustration describe the sequence of folding events for cytochrome c well (as inferred from hydrogen-exchange experiments), thereby justifying their use as a starting point. At low pH, the experimentally observed folding sequence of cytochrome c deviates from that at pH 7 and from models with perfectly funneled energy landscapes. Here, alternate folding pathways are a result of "chemical frustration". This frustration arises because some regions of the protein are destabilized more than others due to the heterogeneous distribution of titratable residues that are protonated at low pH. Beginning with native structure-based terms, we construct more complex models by adding chemical frustration. These more complex models only modestly perturb the energy landscape, which remains, overall, well funneled. These perturbed models can accurately describe how alternative folding pathways are used at low pH. At alkaline pH, cytochrome c populates distinctly different structural ensembles. For instance, lysine residues are deprotonated and compete for the heme ligation site. The same models that can describe folding at low pH also predict well the structures and relative stabilities of intermediates populated at alkaline pH. The success of models based on funneled energy landscapes suggest that cytochrome c folding is driven primarily by native contacts. The presence of heme appears to add chemical complexity to the folding process, but it does not require fundamental modification of the general principles used to describe folding. Moreover, its added complexity provides a valuable means of probing the folding energy landscape in greater detail than is possible with simpler systems.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

Kinetic barriers and the role of topology in protein and RNA folding

This review compares the folding behavior of proteins and RNAs. Topics covered include the role of topology in the determination of folding rates, major folding events including collapse, properties of denatured states, pathway heterogeneity, and the influence of the mode of initiation on the folding pathway.

Protein folding: independent unrelated pathways or predetermined pathway with optional errors

The observation of heterogeneous protein folding kinetics has been widely interpreted in terms of multiple independent unrelated pathways (IUP model), both experimentally and in theoretical calculations. However, direct structural information on folding intermediates and their properties now indicates that all of a protein population folds through essentially the same stepwise pathway, determined by cooperative native-like foldon units and the way that the foldons fit together in the native protein. It is essential to decide between these fundamentally different folding mechanisms. This article shows, contrary to previous supposition, that the heterogeneous folding kinetics observed for the staphylococcal nuclease protein (SNase) does not require alternative parallel pathways. SNase folding kinetics can be fit equally well by a single predetermined pathway that allows for optional misfolding errors, which are known to occur ubiquitously in protein folding. Structural, kinetic, and thermodynamic information for the folding intermediates and pathways of many proteins is consistent with the predetermined pathway-optional error (PPOE) model but contrary to the properties implied in IUP models.

The foldon substructure of staphylococcal nuclease

To search for submolecular foldon units, the spontaneous reversible unfolding and refolding of staphylococcal nuclease under native conditions was studied by a kinetic native-state hydrogen exchange (HX) method. As for other proteins, it appears that staphylococcal nuclease is designed as an assembly of well-integrated foldon units that may define steps in its folding pathway and may regulate some other functional properties. The HX results identify 34 amide hydrogens that exchange with solvent hydrogens under native conditions by way of large transient unfolding reactions. The HX data for each hydrogen measure the equilibrium stability (Delta G(HX)) and the kinetic unfolding and refolding rates (k(op) and k(cl)) of the unfolding reaction that exposes it to exchange. These parameters separate the 34 identified residues into three distinct HX groupings. Two correspond to clearly defined structural units in the native protein, termed the blue and red foldons. The remaining HX grouping contains residues, not well separated by their HX parameters alone, that represent two other distinct structural units in the native protein, termed the green and yellow foldons. Among these four sets, a last unfolding foldon (blue) unfolds with a rate constant of 6 x 10(-6) s(-1) and free energy equal to the protein's global stability (10.0 kcal/mol). It represents part of the beta-barrel, including mutually H-bonding residues in the beta 4 and beta 5 strands, a part of the beta 3 strand that H-bonds to beta 5, and residues at the N-terminus of the alpha2 helix that is capped by beta 5. A second foldon (green), which unfolds and refolds more rapidly and at slightly lower free energy, includes residues that define the rest of the native alpha2 helix and its C-terminal cap. A third foldon (yellow) defines the mutually H-bonded beta1-beta2-beta 3 meander, completing the native beta-barrel, plus an adjacent part of the alpha1 helix. A final foldon (red) includes residues on remaining segments that are distant in sequence but nearly adjacent in the native protein. Although the structure of the partially unfolded forms closely mimics the native organization, four residues indicate the presence of some nonnative misfolding interactions. Because the unfolding parameters of many other residues are not determined, it seems likely that the concerted foldon units are more extensive than is shown by the 34 residues actually observed.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.