Molecular collapse: the rate-limiting step in two-state cytochrome c folding

Biophysical and in-silico studies on the structure-function relationship of Brugia malayi protein disulfide isomerase

Human Lymphatic filariasis is caused by parasitic nematodes Wuchereria bancrofti, Brugia malayi, and Brugia timori. Protein disulfide isomerase (PDI), a redox-active enzyme, helps to form and isomerize the disulfide bonds, thereby acting as a chaperone. Such activity is essential for activating many essential enzymes and functional proteins. Brugia malayi protein disulfide isomerase (BmPDI) is crucial for parasite survival and an important drug target. Here, we used a combination of spectroscopic and computational analysis to study the structural and functional changes in the BmPDI during unfolding. Tryptophan fluorescence data revealed two well-separated transitions during the unfolding process, suggesting that the unfolding of the BmPDI is non-cooperative. The binding of the fluorescence probe 8-anilino-1-naphthalene sulfonic acid dye (ANS) validated the results obtained by the pH unfolding. The dynamics of molecular simulation performed at different pH conditions revealed the structural basis of BmPDI unfolding. Detailed analysis suggested that under different pH, both the global structure and the conformational dynamics of the active site residues were differentially altered. Our multiparametric study reveals the differential dynamics and collective motions of BmPDI unfolding, providing insights into its structure-function relationship.Communicated by Ramaswamy H. Sarma.

Factors That Control the Force Needed to Unfold a Membrane Protein in Silico Depend on the Mode of Denaturation

Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane-protein interactions had only moderate influences, except for the extraction of the last helix where the membrane-protein interactions had a stronger influence. The up-down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.

Folding and misfolding of potassium channel monomers during assembly and tetramerization

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τfold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization.

Dramatic Shape Changes Occur as Cytochrome c Folds

Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different

The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

Synergy in protein-osmolyte mixtures

Virtually all taxa use osmolytes to protect cells against biochemical stress. Osmolytes often occur in mixtures, such as the classical combination of urea with TMAO (trimethylamine N-oxide) in cartilaginous fish or the cocktail of at least six different osmolytes in the kidney. The concentration patterns of osmolyte mixtures found in vivo make it likely that synergy between them plays an important role. Using statistical mechanical n-component Kirkwood-Buff theory, we show from first principles that synergy in protein-osmolyte systems can arise from two separable sources: (1) mutual alteration of protein surface solvation and (2) effects mediated through bulk osmolyte chemical activities. We illustrate both effects in a four-component system with the experimental example of the unfolding of a notch ankyrin domain in urea-TMAO mixtures, which make urea a less effective denaturant and TMAO a more effective stabilizer. Protein surface effects are primarily responsible for this synergy. The specific patterns of surface solvation point to denatured state expansion as the main factor, as opposed to direct competition.

Morphological diversity and polymorphism of self-assembling collagen peptides controlled by length of hydrophobic domains

Synthetic collagen mimetic peptides are used to probe the role of hydrophobic forces in mediating protein self-assembly. Higher order association is an integral property of natural collagens, which assemble into fibers and meshes that comprise the extracellular matrix of connective tissues. The unique triple-helix fold fully exposes two-thirds of positions in the protein to solvent, providing ample opportunities for engineering interaction sites. Inclusion of just a few hydrophobic groups in a minimal peptide promotes a rich variety of self-assembly behaviors, resulting in hundred-nanometer to micron size nanodiscs and nanofibers. Morphology depends primarily on the length of hydrophobic domains. Peptide discs contain lipophilic domains capable of sequestering small hydrophobic dyes. Combining multiple peptide types result in composite structures of discs and fibers ranging from stars to plates-on-a-string. These systems provide valuable tools to shed insight into the fundamental principles underlying hydrophobicity-driven higher order protein association that will facilitate the design of self-assembling systems in biomaterials and nanomedical applications.

Kinetics study on the HIV-1 ectodomain protein quaternary structure formation reveals coupling of chain folding and self-assembly in the refolding cascade

Entry of HIV-1 into the target cell is mediated by the envelope glycoprotein consisting of noncovalently associated surface subunit gp120 and transmembrane subunit gp41. To form a functional gp41 complex, the protein undergoes hairpin formation and self-assembly. The fusion event can be inhibited by gp41-derived peptides at nanomolar concentration and is highly dependent on the time of addition, implying a role of folding kinetics on the inhibitory action. Oligomerization of the gp41 ectodomain was demonstrated by light scattering measurements. Kinetic study by stopped-flow fluorescence and absorption measurements (i) revealed a multistate folding pathway and stable intermediates; (ii) showed a dissection of fast and slow components for early and late stages of folding, respectively, with 3 orders of magnitude difference in the time scale; (iii) showed the slow process was attributed to misfolding and unzipping of the hairpin; and (iv) showed retardation of the native hairpin formation is assumed to lead to coupling of the correctly registered hairpin and self-assembly. This coupling allows the deduction on the time scale of intrachain folding (0.1-1 s) for the protein. The folding reaction was illustrated by a free energy profile to explain the temporal dichotomy of fast and slow steps of folding as well as effective inhibition by gp41-derived peptide.

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.

Ultrafast hydrogen exchange reveals specific structural events during the initial stages of folding of cytochrome c

Many proteins undergo a sharp decrease in chain dimensions during early stages of folding, prior to the rate-limiting step in folding. However, it remains unclear whether compact states are the result of specific folding events or a general hydrophobic collapse of the poly peptide chain driven by the change in solvent conditions. To address this fundamental question, we extended the temporal resolution of NMR-detected H/D exchange labeling experiments into the microsecond regime by adopting a microfluidics approach. By observing the competition between H/D exchange and folding as a function of labeling pH, coupled with direct measurement of exchange rates in the unfolded state, we were able to monitor hydrogen-bond formation for over 50 individual backbone NH groups within the initial 140 microseconds of folding of horse cytochrome c. Clusters of solvent-shielded amide protons were observed in two α-helical segments in the C-terminal half of the protein, while the N-terminal helix remained largely unstructured, suggesting that proximity in the primary structure is a major factor in promoting helix formation and association at early stages of folding, while the entropically more costly long-range contacts between the N- and C-terminal helices are established only during later stages. Our findings clearly indicate that the initial chain condensation in cytochrome c is driven by specific interactions among a subset of α-helical segments rather than a general hydrophobic collapse.

Simplified protein models: predicting folding pathways and structure using amino acid sequences

We demonstrate the ability of simultaneously determining a protein's folding pathway and structure using a properly formulated model without prior knowledge of the native structure. Our model employs a natural coordinate system for describing proteins and a search strategy inspired by the observation that real proteins fold in a sequential fashion by incrementally stabilizing nativelike substructures or "foldons." Comparable folding pathways and structures are obtained for the twelve proteins recently studied using atomistic molecular dynamics simulations [K. Lindorff-Larsen, S. Piana, R. O. Dror, D. E. Shaw, Science 334, 517 (2011)], with our calculations running several orders of magnitude faster. We find that nativelike propensities in the unfolded state do not necessarily determine the order of structure formation, a departure from a major conclusion of the molecular dynamics study. Instead, our results support a more expansive view wherein intrinsic local structural propensities may be enhanced or overridden in the folding process by environmental context. The success of our search strategy validates it as an expedient mechanism for folding both in silico and in vivo.

Interconnection of salt-induced hydrophobic compaction and secondary structure formation depends on solution conditions: revisiting early events of protein folding at single molecule resolution

What happens in the early stage of protein folding remains an interesting unsolved problem. Rapid kinetics measurements with cytochrome c using submillisecond continuous flow mixing devices suggest simultaneous formation of a compact collapsed state and secondary structure. These data seem to indicate that collapse formation is guided by specific short and long range interactions (heteropolymer collapse). A contrasting interpretation also has been proposed, which suggests that the collapse formation is rapid, nonspecific, and a trivial solvent related compaction, which could as well be observed by a homopolymer (homopolymer collapse). We address this controversy using fluorescence correlation spectroscopy (FCS), which enables us to monitor the salt-induced compaction accompanying collapse formation and the associated time constant directly at single molecule resolution. In addition, we follow the formation of secondary structure using far UV CD. The data presented here suggest that both these models (homopolymer and heteropolymer) could be applicable depending on the solution conditions. For example, the formation of secondary structure and compact state is not simultaneous in aqueous buffer. In aqueous buffer, formation of the compact state occurs through a two-state co-operative transition following heteropolymer formalism, whereas secondary structure formation takes place gradually. In contrast, in the presence of urea, a compaction of the protein radius occurs gradually over an extended range of salt concentration following homopolymer formalism. The salt-induced compaction and the formation of secondary structure take place simultaneously in the presence of urea.

The H2A-H2B dimeric kinetic intermediate is stabilized by widespread hydrophobic burial with few fully native interactions

The H2A-H2B histone heterodimer folds via monomeric and dimeric kinetic intermediates. Within ∼5 ms, the H2A and H2B polypeptides associate in a nearly diffusion limited reaction to form a dimeric ensemble, denoted I₂ and I₂*, the latter being a subpopulation characterized by a higher content of nonnative structure (NNS). The I₂ ensemble folds to the native heterodimer, N₂, through an observable, first-order kinetic phase. To determine the regions of structure in the I₂ ensemble, we characterized 26 Ala mutants of buried hydrophobic residues, spanning the three helices of the canonical histone folds of H2A and H2B and the H2B C-terminal helix. All but one targeted residue contributed significantly to the stability of I₂, the transition state and N₂; however, only residues in the hydrophobic core of the dimer interface perturbed the I₂* population. Destabilization of I₂* correlated with slower folding rates, implying that NNS is not a kinetic trap but rather accelerates folding. The pattern of Φ values indicated that residues forming intramolecular interactions in the peripheral helices contributed similar stability to I₂ and N₂, but residues involved in intermolecular interactions in the hydrophobic core are only partially folded in I₂. These findings suggest a dimerize-then-rearrange model. Residues throughout the histone fold contribute to the stability of I₂, but after the rapid dimerization reaction, the hydrophobic core of the dimer interface has few fully native interactions. In the transition state leading to N₂, more native-like interactions are developed and nonnative interactions are rearranged.

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.

Role of protein stabilizers on the conformation of the unfolded state of cytochrome c and its early folding kinetics: investigation at single molecular resolution

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r(H)) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r(H) of 29 A) over an extended one (r(H) of 40 A), which forms in their absence. Also, the time constant of a kinetic component (tau(R)) of about 30 micros has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.

Mutational studies uncover non-native structure in the dimeric kinetic intermediate of the H2A-H2B heterodimer

The folding pathway of the histone H2A-H2B heterodimer minimally includes an on-pathway, dimeric, burst-phase intermediate, I(2). The partially folded H2A and H2B monomers populated at equilibrium were characterized as potential monomeric kinetic intermediates. Folding kinetics were compared for initiation from isolated, folded monomers and the heterodimer unfolded in 4 M urea. The observed rates were virtually identical above 0.4 M urea, exhibiting a log-linear relationship on the final denaturant concentration. Below approximately 0.4 M urea (concentrations inaccessible from the 4-M urea unfolded state), a rollover in the rates was observed; this suggests that a component of the I(2) ensemble contains non-native structure that rearranges/isomerizes to a more native-like species. The contribution of helix propensity to the stability of the I(2) ensemble was assessed with a set of H2A-H2B mutants containing Ala and Gly replacements at nine sites, focusing mainly on the long, central alpha2 helix. Equilibrium and kinetic folding/unfolding data were collected to determine the effects of the mutations on the stability of I(2) and the transition state between I(2) and N(2). This limited mutational study indicated that residues in the alpha2 helices of H2A and H2B as well as alpha1 of H2B and both the C-terminus of alpha3 and the short alphaC helix of H2A contribute to the stability of the I(2) burst-phase species. Interestingly, at least eight of the nine targeted residues stabilize I(2) by interactions that are non-native to some extent. Given that destabilizing I(2) and these non-native interactions does not accelerate folding, it is concluded that the native and non-native structures present in the I(2) ensemble enable efficient folding of H2A-H2B.

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.

Irreversible thermal denaturation of cytochrome C studied by electrospray mass spectrometry

This work uses electrospray ionization mass spectrometry (ESI-MS) in conjunction with hydrogen/deuterium exchange (HDX) and optical spectroscopy for characterizing the solution-phase properties of cytochrome c (cyt c) after heat exposure. Previous work demonstrated that heating results in irreversible denaturation for a subpopulation of proteins in the sample. However, that study did not investigate the physical reasons underlying this interesting effect. Here we report that the formation of oxidative modifications at elevated temperature plays a key role for the observed behavior. Tryptic digestion followed by tandem mass spectrometry is used to identify individual oxidation sites. Trp59 and Met80 are among the modified amino acids. In native cyt c both of these residues are buried deep within the protein structure, such that covalent modifications would be expected to be particularly disruptive. ESI-MS analysis after heat exposure results in a bimodal charge-state distribution. Oxidized protein appears predominantly in charge states around 11+, whereas a considerably lower degree of oxidation is observed for the 7+ and 8+ peaks. This finding confirms that different oxidation levels are associated with different solution-phase conformations. HDX measurements for different charge states are complicated by peak distortions arising from oxygen adduction. Nonetheless, comparison with simulated peak shapes clearly shows that the HDX properties are different for high- and low-charge states, confirming that interconversion between unfolded and folded conformers is blocked in solution. In addition to oxidation, partial aggregation upon heat exposure likely contributes to the formation of irreversibly denatured protein.

Early events, kinetic intermediates and the mechanism of protein folding in cytochrome C

Kinetic studies of the early events in cytochrome c folding are reviewed with a focus on the evidence for folding intermediates on the submillisecond timescale. Evidence from time-resolved absorption, circular dichroism, magnetic circular dichroism, fluorescence energy and electron transfer, small-angle X-ray scattering and amide hydrogen exchange studies on the t ≤ 1 ms timescale reveals a picture of cytochrome c folding that starts with the ~ 1-μs conformational diffusion dynamics of the unfolded chains. A fractional population of the unfolded chains collapses on the 1 – 100 μs timescale to a compact intermediate I_C containing some native-like secondary structure. Although the existence and nature of I_C as a discrete folding intermediate remains controversial, there is extensive high time-resolution kinetic evidence for the rapid formation of I_C as a true intermediate, i.e., a metastable state separated from the unfolded state by a discrete free energy barrier. Final folding to the native state takes place on millisecond and longer timescales, depending on the presence of kinetic traps such as heme misligation and proline mis-isomerization. The high folding rates observed in equilibrium molten globule models suggest that I_C may be a productive folding intermediate. Whether it is an obligatory step on the pathway to the high free energy barrier associated with millisecond timescale folding to the native state, however, remains to be determined.

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an alpha-helix and two beta-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 A). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.

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.

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.

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.

Folding and assembly of large macromolecular complexes monitored by hydrogen-deuterium exchange and mass spectrometry

Recent advances in protein mass spectrometry (MS) have enabled determinations of hydrogen deuterium exchange (HDX) in large macromolecular complexes. HDX-MS became a valuable tool to follow protein folding, assembly and aggregation. The methodology has a wide range of applications in biotechnology ranging from quality control for over-expressed proteins and their complexes to screening of potential ligands and inhibitors. This review provides an introduction to protein folding and assembly followed by the principles of HDX and MS detection, and concludes with selected examples of applications that might be of interest to the biotechnology community.

Unique fluorophores in the dimeric archaeal histones hMfB and hPyA1 reveal the impact of nonnative structure in a monomeric kinetic intermediate

Homodimeric archaeal histones and heterodimeric eukaryotic histones share a conserved structure but fold through different kinetic mechanisms, with a correlation between faster folding/association rates and the population of kinetic intermediates. Wild-type hMfB (from Methanothermus fervidus) has no intrinsic fluorophores; Met35, which is Tyr in hyperthermophilic archaeal histones such as hPyA1 (from Pyrococcus strain GB-3A), was mutated to Tyr and Trp. Two Tyr-to-Trp mutants of hPyA1 were also characterized. All fluorophores were introduced into the long, central alpha-helix of the histone fold. Far-UV circular dichroism (CD) indicated that the fluorophores did not significantly alter the helical content of the histones. The equilibrium unfolding transitions of the histone variants were two-state, reversible processes, with DeltaG degrees (H2O) values within 1 kcal/mol of the wild-type dimers. The hPyA1 Trp variants fold by two-state kinetic mechanisms like wild-type hPyA1, but with increased folding and unfolding rates, suggesting that the mutated residues (Tyr-32 and Tyr-36) contribute to transition state structure. Like wild-type hMfB, M35Y and M35W hMfB fold by a three-state mechanism, with a stopped-flow CD burst-phase monomeric intermediate. The M35 mutants populate monomeric intermediates with increased secondary structure and stability but exhibit decreased folding rates; this suggests that nonnative interactions occur from burial of the hydrophobic Tyr and Trp residues in this kinetic intermediate. These results implicate the long central helix as a key component of the structure in the kinetic monomeric intermediates of hMfB as well as the dimerization transition state in the folding of hPyA1.

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.

Folding of a universal ribozyme: the ribonuclease P RNA

Ribonuclease P is among the first ribozymes discovered, and is the only ubiquitously occurring ribozyme besides the ribosome. The bacterial RNase P RNA is catalytically active without its protein subunit and has been studied for over two decades as a model system for RNA catalysis, structure and folding. This review focuses on the thermodynamic, kinetic and structural frameworks derived from the folding studies of bacterial RNase P RNA.

Strange kinetic phase in the extremely early folding process of beta-lactoglobulin

A continuous-wave probed laser-induced temperature jump system was constructed and applied to monitor the changes in tryptophan fluorescence of the beta-lactoglobulin during its folding; the kinetic phases were traced from 300 ns to 10 ms after a temperature jump. Notably, an early phase with typical squeezed-exponential characteristics, [exp[-(kt)(beta)], beta>1.0], was observed around several tens of microseconds after the temperature jump, which is actually the earliest phase ever observed for beta-lactoglobulin. This process can be explained by conformational shift occurring within the unfolded ensemble (U-->U'), which is followed by the non-native intermediate (I) formation of this protein.

Calculation of the free energy and cooperativity of protein folding

Calculation of the free energy of protein folding and delineation of its pre-organization are of foremost importance for understanding, predicting and designing biological macromolecules. Here, we introduce an energy smoothing variant of parallel tempering replica exchange Monte Carlo (REMS) that allows for efficient configurational sampling of flexible solutes under the conditions of molecular hydration. Its usage to calculate the thermal stability of a model globular protein, Trp cage TC5b, achieves excellent agreement with experimental measurements. We find that the stability of TC5b is attained through the coupled formation of local and non-local interactions. Remarkably, many of these structures persist at high temperature, concomitant with the origin of native-like configurations and mesostates in an otherwise macroscopically disordered unfolded state. Graph manifold learning reveals that the conversion of these mesostates to the native state is structurally heterogeneous, and that the cooperativity of their formation is encoded largely by the unfolded state ensemble. In all, these studies establish the extent of thermodynamic and structural pre-organization of folding of this model globular protein, and achieve the calculation of macromolecular stability ab initio, as required for ab initio structure prediction, genome annotation, and drug design.

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.

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.

Aggregation as the basis for complex behaviour of cutinase in different denaturants

We have previously described the complexity of the folding of the lipolytic enzyme cutinase from F. solani pisi in guanidinium chloride. Here we extend the refolding analysis by refolding from the pH-denatured state and analyze the folding behaviour in the presence of the weaker denaturant urea and the stronger denaturant guanidinium thiocyanate. In urea there is excellent consistency between equilibrium and kinetic data, and the intermediate accumulating at low denaturant concentrations is off-pathway. However, in GdmCl, refolding rates, and consequently the stability of the native state, vary significantly depending on whether refolding takes place from the pH- or GdmCl-denatured state, possibly due to transient formation of aggregates during folding from the GdmCl-denatured state. In GdmSCN, stability is reduced by several kcal/mol with significant aggregation in the unfolding transition region. The basis for the large variation in folding behaviour may be the denaturants' differential ability to support formation of exposed hydrophobic regions and consequent changes in aggregative properties during refolding.

Fully reduced ribonuclease A does not expand at high denaturant concentration or temperature

The dimensions of a denatured protein, fully reduced ribonuclease A (r-RNase A), have been measured using synchrotron-based small angle X-ray scattering. The radius of gyration, 34-35 A, is unchanged from 0-6 M guanidinium chloride and from 20-90 degrees C at pH 2.5, and agrees with the known scaling behavior for a multitude of chemically denatured states. The polypeptide is behaving as a statistical coil in the non-interacting, high-temperature limit.

Equilibrium unfolding of DLC8 monomer by urea and guanidine hydrochloride: Distinctive global and residue level features

We present circular dichroism (CD), steady state fluorescence and multidimensional NMR investigations on the equilibrium unfolding of monomeric dynein light chain protein (DLC8) by urea and guanidine hydrochloride (GdnHCl). Quantitative analysis of the CD and fluorescence denaturation curves reveals that urea unfolding is a two-state process, whereas guanidine unfolding is more complex. NMR investigations in the native state and in the near native states created by low denaturant concentrations enabled residue level characterization of the early structural and dynamic perturbations by the two denaturants. Firstly, (15)N transverse relaxation rates in the native state indicate that the regions around N10, Q27, the loop between beta2 and beta4 strands, and K87 at the C-terminal are potential unfolding initiation sites in the protein. Amide and (15)N chemical shift perturbations indicate different accessibilities of the residues along the chain and help identify locations of the early perturbations by the two denaturants. Guanidine and urea are seen to interact at several sites some of which are different in the two cases. Notable among the common interaction site is that around K87 which is in close proximity to W54 on the protein structure, but the interaction modes of the two denaturants are different. The secondary chemical shifts indicate that the structural perturbation by 1M urea is small, compared to that by guanidine which is more encompassing over the length of the chain. The probable (phi, psi) changes at the individual residues have been calculated using the TALOS algorithm. It appears that the helices in the protein are significantly perturbed by guanidine. Further, comparison of the spectral density functions of the native and the two near native states in the two denaturants implicate greater loosening of the structure by guanidine as compared to that by urea, even though the structures are still in the native state ensemble. These differences in the early perturbations of the native state structure and dynamics by the two denaturants might direct the protein along different pathways, as the unfolding progresses on further increasing the denaturant concentration.

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.

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.

Detection of a hidden folding intermediate in the focal adhesion target domain: Implications for its function and folding

The focal adhesion target (FAT) domain of focal adhesion kinase has a four-helix bundle structure. Based on a hydrogen exchange-constrained computer simulation study and some indirect experimental results, it has been suggested that a partially unfolded state of the FAT domain with the N-terminal helix unfolded plays an important role in its biological function. Here, using a native-state hydrogen exchange method, we directly detected an intermediate with the N-terminal helix unfolded in a mutant (Y925E) of the FAT domain. In addition, kinetic folding studies on the FAT domain suggest that this intermediate exists on the native side of the rate-limiting transition state for folding. These results provide more direct evidence of the existence of the proposed intermediate and help to understand the folding mechanism of small single domain proteins.

Small Proteins Fold Through Transition States With Native-like Topologies

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.

Secondary structure determines protein topology

Using a test set of 13 small, compact proteins, we demonstrate that a remarkably simple protocol can capture native topology from secondary structure information alone, in the absence of long-range interactions. It has been a long-standing open question whether such information is sufficient to determine a protein's fold. Indeed, even the far simpler problem of reconstructing the three-dimensional structure of a protein from its exact backbone torsion angles has remained a difficult challenge owing to the small, but cumulative, deviations from ideality in backbone planarity, which, if ignored, cause large errors in structure. As a familiar example, a small change in an elbow angle causes a large displacement at the end of your arm; the longer the arm, the larger the displacement. Here, correct secondary structure assignments (alpha-helix, beta-strand, beta-turn, polyproline II, coil) were used to constrain polypeptide backbone chains devoid of side chains, and the most stable folded conformations were determined, using Monte Carlo simulation. Just three terms were used to assess stability: molecular compaction, steric exclusion, and hydrogen bonding. For nine of the 13 proteins, this protocol restricts the main chain to a surprisingly small number of energetically favorable topologies, with the native one prominent among them.