Parallel protein-unfolding pathways revealed and mapped

Theoretical studies of protein folding suggest that multiple folding pathways should exist, but there is little experimental evidence to support this. Here we demonstrate changes in the flux between different transition states on parallel folding pathways, resulting in unprecedented upward curvature in the denaturant-dependent unfolding kinetics of a beta-sandwich protein. As denaturant concentration increases, the highly compact transition state of one pathway becomes destabilized and the dominant flux of protein molecules shifts toward another pathway with a less structured transition state. Furthermore, point mutations alter the relative accessibility of the pathways, allowing the structure of two transition states on separate, direct folding pathways to be mapped by systematic Phi-value analysis. It has been suggested that pathways with diffuse rather than localized transition states are evolutionarily selected to prevent misfolding, and indeed we find that the transition state favored at high concentrations of denaturant is more polarized than the physiologically relevant one.     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

A New Computational Approach for Mechanical Folding Kinetics of RNA Hairpins

Based on an ensemble of kinetically accessible conformations, we propose a new analytical model for RNA folding kinetics. The model gives populational kinetics, kinetic rates, transition states, and pathways from the rate matrix. Applications of the new kinetic model to mechanical folding of RNA hairpins such as trans-activation-responsive RNA reveal distinct kinetic behaviors in different force regimes, from zero force to forces much stronger than the critical force for the folding-unfolding transition. In the absence of force or a low force, folding can be initiated (nucleated) at any position by forming the first base stack and there exist many pathways for the folding process. In contrast, for a higher force, the folding/unfolding would predominantly proceed along a single zipping/unzipping pathway. Studies for different hairpin-forming sequences indicate that depending on the nucleotide sequence, a kinetic intermediate can emerge in the low force regime but disappear in high force regime, and a new kinetic intermediate, which is absent in the low and high force regimes, can emerge in the medium force range. Variations of the force lead to changes in folding cooperativity and rate-limiting steps. The predicted network of pathways for trans-activation-responsive RNA suggests two parallel dominant pathways. The rate-limiting folding steps (at f = 8 pN) are the formation of specific basepairs that are 2-4 basepairs away from the loop. At a higher force (f = 11 pN), the folding rate is controlled by the formation of the bulge loop. The predicted rates and transition states are in good agreement with the experimental data for a broad force regime.     

Thermal unfolding simulations of a multimeric protein--transition state and unfolding pathways

The folding of an oligomeric protein poses an extra challenge to the folding problem because the protein not only has to fold correctly; it has to avoid nonproductive aggregation. We have carried out over 100 molecular dynamics simulations using an implicit solvation model at different temperatures to study the unfolding of one of the smallest known tetramers, p53 tetramerization domain (p53tet). We found that unfolding started with disruption of the native tetrameric hydrophobic core. The transition state for the tetramer to dimer transition was characterized as a diverse ensemble of different structures using Phi value analysis in quantitative agreement with experimental data. Despite the diversity, the ensemble was still native-like with common features such as partially exposed tetramer hydrophobic core and shifts in the dimer-dimer arrangements. After passing the transition state, the secondary and tertiary structures continued to unfold until the primary dimers broke free. The free dimer had little secondary structure left and the final free monomers were random-coil like. Both the transition states and the unfolding pathways from these trajectories were very diverse, in agreement with the new view of protein folding. The multiple simulations showed that the folding of p53tet is a mixture of the framework and nucleation-condensation mechanisms and the folding is coupled to the complex formation. We have also calculated the entropy and effective energy for the different states along the unfolding pathway and found that the tetramerization is stabilized by hydrophobic interactions.     

Microscopic reversibility of protein folding in molecular dynamics simulations of the engrailed homeodomain

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions or, in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultrafast unfolding and folding protein at its melting temperature to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 micros of simulation in water, we found six transient denaturing events near the melting temperature (323 and 330 K) and an additional refolding event following a previously identified unfolding event at a high temperature (373 K). In each case, unfolding and refolding transition state ensembles were identified, and they agreed well with experiment on the basis of a comparison of S and Phi values. On the basis of several structural properties, these 13 transition state ensembles agreed very well with each other and with four previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other but the paths deviated.     

GroEL Can Unfold Late Intermediates Populated on the Folding Pathways of Monellin

The modulation of the folding mechanism of the small protein single-chain monellin (MNEI) by the Escherichia coli chaperone GroEL has been studied. In the absence of the chaperone, the folding of monellin occurs via three parallel routes. When folding is initiated in the presence of a saturating concentration of GroEL, only 50-60% of monellin molecules fold completely. The remaining 40-50% of the monellin molecules remain bound to the GroEL and are released only upon addition of ATP. It is shown that the basic folding mechanism of monellin is not altered by the presence of GroEL, but that it occurs via only one of the three available routes when folding is initiated in the presence of saturating concentrations of GroEL. Two pathways become nonoperational because GroEL binds very tightly to early intermediates that populate these pathways in a manner that makes the GroEL-bound intermediates incompetent to fold. This accounts for the monellin molecules that remain GroEL-bound at the end of the folding reaction. The third pathway remains operational because the GroEL-bound early intermediate on this pathway is folding-competent, suggesting that this early intermediate binds to GroEL in a manner that is different from that of the binding of the early intermediates on the other two pathways. It appears, therefore, that the same protein can bind GroEL in more than one way. The modulation of the folding energy landscape of monellin by GroEL occurs because GroEL binds folding intermediates on parallel folding pathways, in different ways, and with different affinities. Moreover, when GroEL is added to refolding monellin at different times after commencement of refolding, the unfolding of two late kinetic intermediates on two of the three folding pathways can be observed. It appears that the unfolding of late folding intermediates is enabled by a thermodynamic coupling mechanism, wherein GroEL binds more tightly to an early intermediate than to a late intermediate on a folding pathway, with preferential binding energy being larger than the stability of the late intermediate. Hence, it is shown that GroEL can inadvertently and passively cause, through its ability to bind different folding intermediates differentially, the unfolding of late productive intermediates on folding pathways, and that its unfolding action is not restricted solely to misfolded or kinetically trapped intermediates.     

Multiple folding pathways of the SH3 domain

Experimental observations suggest that proteins follow different folding pathways under different environmental conditions. We perform molecular dynamics simulations of a model of the c-Crk SH3 domain over a broad range of temperatures, and identify distinct pathways in the folding transition. We determine the kinetic partition temperature-the temperature for which the c-Crk SH3 domain undergoes a rapid folding transition with minimal kinetic barriers-and observe that below this temperature the model protein may undergo a folding transition by multiple folding pathways via only one or two intermediates. Our findings suggest the hypothesis that the SH3 domain, a protein fold for which only two-state folding kinetics was observed in previous experiments, may exhibit intermediate states under conditions that strongly stabilize the native state.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

Titration Properties and Thermodynamics of the Transition State for Folding: Comparison of Two-state and Multi-state Folding Pathways

CI2 folds and unfolds as a single cooperative unit by simple two-state kinetics, which enables the properties of the transition state to be measured from both the forward and backward rate constants. We have examined how the free energy of the transition state for the folding of chymotrypsin inhibitor 2 (CI2) changes with pH and temperature. In addition to the standard thermodynamic quantities, we have measured the overall acid-titration properties of the transition state and its heat capacity relative to both the denatured and native states. We were able to determine the latter by a method analogous to a well-established procedure for measuring the change in heat capacity for equilibrium unfolding: the enthalpy of activation of unfolding at different values of acid pH were plotted against the average temperature of each determination. Our results show that the transition state of CI2 has lost most of the electrostatic and van der Waals' interactions that are found in the native state, but it remains compact and this prevents water molecules from entering some parts of the hydrophobic core. The properties of the transition state of CI2 are then compared with the major folding transition state of the larger protein barnase, which folds by a multi-state mechanism, with the accumulation of a partly structured intermediate (D^physor I). CI2 folds from a largely unstructured denatured state under physiological conditionsviaa transition state which is compact but relatively uniformly unstructured, with tertiary and secondary structure being formed in parallel. We term this an expanded pathway. Conversely, barnase folds from a largely structured denatured state in which elements of structure are well formed through a transition state that has islands of folded elements of structure. We term this a compact pathway. These two pathways may correspond to the two extreme ends of a continuous spectrum of protein folding mechanisms. Although the properties of the two transition states are very different, the activation barrier for folding (D^phys→3 ) is very similar for both proteins.     

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering,Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.     

Is protein unfolding the reverse of protein folding? A lattice simulation analysis.

Simulations and experiments that monitor protein unfolding under denaturing conditions are commonly employed to study the mechanism by which a protein folds to its native state in a physiological environment. Due to the differences in conditions and the complexity of the reaction, unfolding is not necessarily the reverse of folding. To assess the relevance of temperature initiated unfolding studies to the folding problem, we compare the folding and unfolding of a 125-residue protein model by Monte Carlo dynamics at two temperatures; the lower one corresponds to the range used in T -jump experiments and the higher one to the range used in unfolding simulations of all-atom models. The trajectories that lead from the native state to the denatured state at these elevated temperatures are less diverse than those observed in the folding simulations. At the lower temperature, the system unfolds through a mandatory intermediate that corresponds to a local free energy minimum. At the higher temperature, no such intermediate is observed, but a similar pathway is followed. The structures contributing to the unfolding pathways resemble most closely those that make up the "fast track" of folding. The transition state for unfolding at the lower temperature (above Tm) is determined and is found to be more structured than the transition state for folding below the melting temperature. This shift towards the native state is consistent with the Hammond postulate. The implications for unfolding simulations of higher resolution models and for unfolding experiments of proteins are discussed.     

Interpretation of protein folding psi values

The structural characterization of transition states is essential for understanding the mechanism of protein folding. Analyzing the effect of mutations on protein stability and folding kinetics in phi-value analysis is commonly used to gain information about the presence of side-chain interactions in transition states. Recently, specific binding of ligands to engineered binding sites was applied to monitor the formation of local structures in transition states (psi analysis). A surprising result from psi analysis was the presence of parallel folding pathways in all reported studies and a major discrepancy between phi and psi values measured in the same protein. Here, we show that psi values cannot be analyzed in the same way as other rate-equilibrium free energy relationships due to the involvement of bimolecular reactions that may have different dissociation constants for the native, unfolded and transition state. As a consequence, psi values reflect the relative binding energy (kappa) of the transition state only for the extreme values of kappa=0 or kappa=1. In all other cases, non-linear rate-equilibrium free-energy relationships (Leffler plots) are observed. This apparently indicates the presence of parallel folding pathways even if folding occurs over a single homogeneous transition state. Consequently, the results from Leffler plots do not yield information about the structural properties of the transition state. This explains the lack of agreement between results from psi analysis and other methods used to characterize protein folding transition states. We further show that the same considerations apply for the analysis of the effect of pH on protein folding.     

The Folding Pathway of the Antibody VL Domain

Antibodies are modular proteins consisting of domains that exhibit a β-sandwich structure, the so-called immunoglobulin fold. Despite structural similarity, differences in folding and stability exist between different domains. In particular, the variable domain of the light chain V_L is unusual as it is associated with misfolding diseases, including the pathologic assembly of the protein into fibrillar structures. Here, we have analysed the folding pathway of a V_L domain with a view to determine features that may influence the relationship between productive folding and fibril formation. The V_L domain from MAK33 (murine monoclonal antibody of the subtype κ/IgG1) has not previously been associated with fibrillisation but is shown here to be capable of forming fibrils. The folding pathway of this V_L domain is complex, involving two intermediates in different pathways. An obligatory early molten globule-like intermediate with secondary structure but only loose tertiary interactions is inferred. The native state can then be formed directly from this intermediate in a phase that can be accelerated by the addition of prolyl isomerases. However, an alternative pathway involving a second, more native-like intermediate is also significantly populated. Thus, the protein can reach the native state via two distinct folding pathways. Comparisons to the folding pathways of other antibody domains reveal similarities in the folding pathways; however, in detail, the folding of the V_L domain is striking, with two intermediates populated on different branches of the folding pathway, one of which could provide an entry point for molecules diverted into the amyloid pathway.     

The family feud: do proteins with similar structures fold via the same pathway?

Theoretical and experimental studies of protein folding have suggested that the topology of the native state may be the most important factor determining the folding pathway of a protein, independent of its specific amino acid sequence. To test this concept, many experimental studies have been carried out with the aim of comparing the folding pathways of proteins that possess similar tertiary structures, but divergent sequences. Many of these studies focus on quantitative comparisons of folding transition state structures, as determined by Phi(f) value analysis of folding kinetic data. In some of these studies, folding transition state structures are found to be highly conserved, whereas in others they are not. We conclude that folds displaying more conserved transition state structures may have the most restricted number of possible folding pathways and that folds displaying low transition state structural conservation possess many potential pathways for reaching the native state.     

Oligomeric States along the Folding Pathways of β2-Microglobulin: Kinetics,Thermodynamics, and Structure

The transition of proteins from their soluble functional state to amyloid fibrils and aggregates is associated with the onset of several human diseases. Protein aggregation often requires some structural reshaping and the subsequent formation of intermolecular contacts. Therefore, the study of the conformation of excited protein states and their ability to form oligomers is of primary importance for understanding the molecular basis of amyloid fibril formation. Here, we investigated the oligomerization processes that occur along the folding of the amyloidogenic human protein β2-microglobulin. The combination of real-time two-dimensional NMR data with real-time small-angle X-ray scattering measurements allowed us to derive thermodynamic and kinetic information on protein oligomerization of different conformational states populated along the folding pathways. In particular, we could demonstrate that a long-lived folding intermediate (I-state) has a higher propensity to oligomerize compared to the native state. Our data agree well with a simple five-state kinetic model that involves only monomeric and dimeric species. The dimers have an elongated shape with the dimerization interface located at the apical side of β2-microglobulin close to Pro32, the residue that has a trans conformation in the I-state and a cis conformation in the native (N) state. Our experimental data suggest that partial unfolding in the apical half of the protein close to Pro32 leads to an excited state conformation with enhanced propensity for oligomerization. This excited state becomes more populated in the transient I-state due to the destabilization of the native conformation by the trans-Pro32 configuration.     

Protein Folding Failure Sets High-Temperature Limit on Growth of Phage P22 in Salmonella enterica Serovar Typhimurium

The high-temperature limit for growth of microorganisms differs greatly depending on their species and habitat. The importance of an organism's ability to manage thermal stress is reflected in the ubiquitous distribution of the heat shock chaperones. Although many chaperones function to reduce protein folding defects, it has been difficult to identify the specific protein folding pathways that set the high-temperature limit of growth for a given microorganism. We have investigated this for a simple system, phage P22 infection of Salmonella enterica serovar Typhimurium. Production of infectious particles exhibited a broad maximum of 150 phage per cell when host cells were grown at between 30 and 39°C in minimal medium. Production of infectious phage declined sharply in the range of 40 to 41°C, and at 42°C, production had fallen to less than 1% of the maximum rate. The host cells maintained optimal division rates at these temperatures. The decrease in phage infectivity was steeper than the loss of physical particles, suggesting that noninfectious particles were formed at higher temperatures. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis revealed a decrease in the tailspike adhesins assembled on phage particles purified from cultures incubated at higher temperatures. The infectivity of these particles was restored by in vitro incubation with soluble tailspike trimers. Examination of tailspike folding and assembly in lysates of phage-infected cells confirmed that the fraction of polypeptide chains able to reach the native state in vivo decreased with increasing temperature, indicating a thermal folding defect rather than a particle assembly defect. Thus, we believe that the folding pathway of the tailspike adhesin sets the high-temperature limit for P22 formation in Salmonella serovar Typhimurium.     

Kinetic studies of the folding of heterodimeric monellin: evidence for switching between alternative parallel pathways

Determining whether or not a protein uses multiple pathways to fold is an important goal in protein folding studies. When multiple pathways are present, defined by transition states that differ in their compactness and structure but not significantly in energy, they may manifest themselves by causing the dependence on denaturant concentration of the logarithm of the observed rate constant of folding to have an upward curvature. In this study, the folding mechanism of heterodimeric monellin [double-chain monellin (dcMN)] has been studied over a range of protein and guanidine hydrochloride (GdnHCl) concentrations, using the intrinsic tryptophan fluorescence of the protein as the probe for the folding reaction. Refolding is shown to occur in multiple kinetic phases. In the first stage of refolding, which is silent to any change in intrinsic fluorescence, the two chains of monellin bind to one another to form an encounter complex. Interrupted folding experiments show that the initial encounter complex folds to native dcMN via two folding routes. A productive folding intermediate population is identified on one route but not on both of these routes. Two intermediate subpopulations appear to form in a fast kinetic phase, and native dcMN forms in a slow kinetic phase. The chevron arms for both the fast and slow phases of refolding are shown to have upward curvatures, suggesting that at least two pathways each defined by a different intermediate are operational during these kinetic phases of structure formation. Refolding switches from one pathway to the other as the GdnHCl concentration is increased.     

Temperature-dependent folding pathways of Pin1 WW domain: an all-atom molecular dynamics simulation of a Gō model

We study the folding thermodynamics and kinetics of the Pin1 WW domain, a three-stranded beta-sheet protein, by using all-atom (except nonpolar hydrogens) discontinuous molecular dynamics simulations at various temperatures with a Gō model. The protein exhibits a two-state folding kinetics near the folding transition temperature. A good agreement between our simulations and the experimental measurements by the Gruebele group has been found, and the simulation sheds new insights into the structure of transition state, which is hard to be straightforwardly captured in experiments. The simulation also reveals that the folding pathways at approximately the transition temperature and at low temperatures are much different, and an intermediate state at a low temperature is predicted. The transition state of this small beta-protein at its folding transition temperature has a well-established hairpin 1 made of beta1 and beta2 strands while its low-temperature kinetic intermediate has a formed hairpin 2 composed of beta2 and beta3 strands. Theoretical results are compared with other simulation results as well as available experimental data. This study confirms that specific side-chain packing in an all-atom Gō model can yield a reasonable prediction of specific folding kinetics for a given protein. Different folding behaviors at different temperatures are interpreted in terms of the interplay of entropy and enthalpy in folding process.     

Thermodynamics and Kinetics of Single-Chain Monellin Folding with Structural Insights into Specific Collapse in the Denatured State Ensemble

Proteins, which behave as random coils in high denaturant concentrations undergo collapse transition similar to polymers on denaturant dilution. We study collapse in the denatured ensemble of single-chain monellin (MNEI) using a coarse-grained protein model and molecular dynamics simulations. The model is validated by quantitatively comparing the computed guanidinium chloride and pH-dependent thermodynamic properties of MNEI folding with the experiments. The computed properties such as the fraction of the protein in the folded state and radius of gyration (R_g) as function of [GuHCl] are in good agreement with the experiments. The folded state of MNEI is destabilized with an increase in pH due to the deprotonation of the residues Glu24 and Cys42. On decreasing [GuHCl], the protein in the unfolded ensemble showed specific compaction. The R_g of the protein decreased steadily with [GuHCl] dilution due to increase in the number of native contacts in all the secondary structural elements present in the protein. MNEI folding kinetics is complex with multiple folding pathways and transiently stable intermediates are populated in these pathways. In strong stabilizing conditions, the protein in the unfolded ensemble showed transition to a more compact unfolded state where R_g decreased by ≈ 17% due to the formation of specific native contacts in the protein. The intermediate populated in the dominant MNEI folding pathway satisfies the structural features of the dry molten globule inferred from experiments.