Thermodynamic characterisation of two transition states along parallel protein folding pathways

Recent work has shown that a β-sandwich domain from the human muscle protein titin (TI I27) unfolds via more than one pathway, providing experimental evidence for a long-standing theoretical prediction in protein folding. Here we present a thermodynamic analysis of two transition states along different folding pathways for this protein. The unusual upwards curvature previously observed in the denaturant-dependent unfolding kinetics is increased at both high and low temperatures, indicating that the high denaturant pathway is becoming more accessible. The transition states in each pathway are structurally distinct and have very different heat capacities. Interestingly the nucleation-condensation pathway is dominant at all physiologically relevant temperatures, supporting the suggestion that pathways with diffuse rather than localised transition states have been selected for by evolution to prevent misfolding.     

Folding pathway of FKBP12 and characterisation of the transition state.

The folding pathway of human FKBP12, a 12 kDa FK506-binding protein (immunophilin), has been characterised. Unfolding and refolding rate constants have been determined over a wide range of denaturant concentrations and data are shown to fit to a two-state model of folding in which only the denatured and native states are significantly populated, even in the absence of denaturant. This simple model for folding, in which no intermediate states are significantly populated, is further supported from stopped-flow circular dichroism experiments in which no fast "burst" phases are observed. FKBP12, with 107 residues, is the largest protein to date which folds with simple two-state kinetics in water (kF=4 s(-1)at 25 degrees C). The topological crossing of two loops in FKBP12, a structural element suggested to cause kinetic traps during folding, seems to have little effect on the folding pathway.The transition state for folding has been characterised by a series of experiments on wild-type FKBP12. Information on the thermodynamic nature of, the solvent accessibility of, and secondary structure in, the transition state was obtained from experiments measuring the unfolding and refolding rate constants as a function of temperature, denaturant concentration and trifluoroethanol concentration. In addition, unfolding and refolding studies in the presence of ligand provided information on the structure of the ligand-binding pocket in the transition state. The data suggest a compact transition state relative to the unfolded state with some 70 % of the surface area buried. The ligand-binding site, which is formed mainly by two loops, is largely unstructured in the transition state. The trifluoroethanol experiments suggest that the alpha-helix may be formed in the transition state. These results are compared with results from protein engineering studies and molecular dynamics simulations (see the accompanying paper).     

Shape of the free energy barriers for protein folding probed by multiple perturbation analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.     

Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.     

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.     

The folding pathway of barnase: the rate-limiting transition state and a hidden intermediate under native conditions

The nature of the rate-limiting transition state at zero denaturant (TS(1)) and whether there are hidden intermediates are the two major unsolved problems in defining the folding pathway of barnase. In earlier studies, it was shown that TS(1) has small phi values throughout the structure of the protein, suggesting that the transition state has either a defined partially folded secondary structure with all side chains significantly exposed or numerous different partially unfolded structures with similar stability. To distinguish the two possibilities, we studied the effect of Gly mutations on the folding rate of barnase to investigate the secondary structure formation in the transition state. Two mutations in the same region of a beta-strand decreased the folding rate by 20- and 50-fold, respectively, suggesting that the secondary structures in this region are dominantly formed in the rate-limiting transition state. We also performed native-state hydrogen exchange experiments on barnase at pD 5.0 and 25 degrees C and identified a partially unfolded state. The structure of the intermediate was investigated using protein engineering and NMR. The results suggest that the intermediate has an omega loop unfolded. This intermediate is more folded than the rate-limiting transition state previously characterized at high denaturant concentrations (TS(2)). Therefore, it exists after TS(2) in folding. Consistent with this conclusion, the intermediate folds with the same rate and denaturant dependence as the wild-type protein, but unfolds faster with less dependence on the denaturant concentration. These and other results in the literature suggest that barnase folds through partially unfolded intermediates that exist after the rate-limiting step. Such folding behavior is similar to those of cytochrome c and Rd-apocyt b(562). Together, we suggest that other small apparently two-state proteins may also fold through hidden intermediates.     

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.     

How do chemical denaturants affect the mechanical folding and unfolding of proteins?

We present the first single-molecule atomic force microscopy study on the effect of chemical denaturants on the mechanical folding/unfolding kinetics of a small protein GB1 (the B1 immunoglobulin-binding domain of protein G from Streptococcus). Upon increasing the concentration of the chemical denaturant guanidinium chloride (GdmCl), we observed a systematic decrease in the mechanical stability of GB1, indicating the softening effect of the chemical denaturant on the mechanical stability of proteins. This mechanical softening effect originates from the reduced free-energy barrier between the folded state and the unfolding transition state, which decreases linearly as a function of the denaturant concentration. Chemical denaturants, however, do not alter the mechanical unfolding pathway or shift the position of the transition state for mechanical unfolding. We also found that the folding rate constant of GB1 is slowed down by GdmCl in mechanical folding experiments. By combining the mechanical folding/unfolding kinetics of GB1 in GdmCl solution, we developed the “mechanical chevron plot” as a general tool to understand how chemical denaturants influence the mechanical folding/unfolding kinetics and free-energy diagram in a quantitative fashion. This study demonstrates great potential in combining chemical denaturation with single-molecule atomic force microscopy techniques to reveal invaluable information on the energy landscape underlying protein folding/unfolding reactions.     

Effects of heme on the structure of the denatured state and folding kinetics of cytochrome b562

Heme-linked proteins, such as cytochromes, are popular subjects for protein folding studies. There is the underlying question of whether the heme affects the structure of the denatured state by cross-linking it and forming other interactions, which would perturb the folding pathway. We have studied wild-type and mutant cytochrome b562 from Escherichia coli, a 106 residue four-alpha-helical bundle. The holo protein apparently refolds with a half-life of 4 micros in its ferrous state. We have analysed the folding of the apo protein using continuous-flow fluorescence as well as stopped-flow fluorescence and CD. The apo protein folded much more slowly with a half-life of 270 micros that was unaffected by the presence of exogenous heme. We examined the nature of the denatured states of both holo and apo proteins by NMR methods over a range of concentrations of guanidine hydrochloride. The starting point for folding of the holo protein in concentrations of denaturant around the denaturation transition was a highly ordered native-like species with heme bound. Fully denatured holo protein at higher concentrations of denaturant consisted of denatured apo protein and free heme. Our results suggest that the very fast folding species of denatured holo protein is in a compact state, whereas the normal folding pathway from fully denatured holo protein consists of the slower folding of the apo protein followed by the binding of heme. These data should be considered in the analysis of folding of heme proteins.     

Structural dynamics, stability and folding of proteins

The present concepts of protein folding in vitro are reviewed. According to these concepts, amino acid sequence of protein, which has appeared a result of evolutionary selection, determines the native structure of protein, the pathway of protein folding, and the existence of free energy barrier between native and denatured states of protein. The latter means that protein macromolecule can exist in either native or denatured state. And all macromolecules in the native state are identical but for structural fluctuations due to Brownian motion of their atoms. Identity of all molecules in native state is of primary importance for their correct functioning. The dependence of protein stability, which is measured as the difference between free energy of protein in native and denatured states, on temperature and denaturant concentration is discussed. The modern approaches characterizing transition state and nucleation are regarded. The role of intermediate and misfolded states in amorphous aggregate and amyloid fibril formation is discussed.     

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.     

Effect of solvent conditions upon refolding pathways and intermediates for a simple lattice protein

Folding pathways and intermediates for a two-dimensional lattice protein have been investigated via computer simulation at various denaturant concentrations. The protein is represented as a chain of 8 hydrophobic (H) and 12 polar (P) beads on a square lattice sequenced in such a way that the native state is a compact hydrophobic core surrounded by a shell of polar beads. Two nonbonded H beads are said to attract each other with a potential of mean force of strength ϵ. Increasing |ϵ/kT| mimics decreasing the denaturant concentration in the solution. Dynamic Monte Carlo simulations have been performed in order to investigate the folding transition and the folding pathways. Sharp folding—unfolding transitions are observed and the folding process proceeds along well-defined pathways that are populated by partially folded intermediates. The folding pathways as well as the populations of the intermediates are strongly dependent upon the denaturant concentration. Generally, intermediates containing long open stretches of H beads are more populated at high denaturant concentration, whereas compact intermediates containing a substantial number of hydrophobic contacts are more populated at low denaturant concentrations. The folding process is also observed to be cooperative in nature in that the chain does not start folding until a key fold in the middle section of the chain is formed correctly. © 1997 John Wiley & Sons, Inc. Biopoly 42: 399–409, 1997     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.     

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.     

Trehalose favors a cutinase compact intermediate off-folding pathway

The folding of cutinase, an enzyme displaying lipolytic activity, has been studied in the presence of trehalose. Equilibrium unfolding data show that trehalose increases the free energy change between folded and unfolded states. Unfolding kinetics reveal the presence of an intermediate which is ca. 60% folded in terms of solvent exposure. Trehalose stabilizes this intermediate relative to the folded state. In contrast, the intermediate revealed by folding kinetics is more compact than the transition state, as shown by the positive slope observed at low denaturant concentration in the chevron plot, as well as the decrease in the observable rate constant for folding with the increase in trehalose concentration. This intermediate displays more than 50% of area buried from the solvent (relative to the native state) compared to around 40% for the transition state for folding and therefore appears to be off the folding pathway. Trehalose stabilizes and guanidine hydrochloride destabilizes this compact intermediate. Both unfolding and folding kinetics show that compact conformational states are stabilized by trehalose, in agreement with current models on the effect of compatible solutes. This effect occurs even for compact states that decelerate the folding as in the case of the intermediate revealed by folding kinetics.     

Wild type, mutant protein unfolding and phase transition detected by single-nanopore recording

Understanding protein folding remains a challenge. A difficulty is to investigate experimentally all the conformations in the energy landscape. Only single molecule methods, fluorescence and force spectroscopy, allow observing individual molecules along their folding pathway. Here we observe that single-nanopore recording can be used as a new single molecule method to explore the unfolding transition and to examine the conformational space of native or variant proteins. We show that we can distinguish unfolded states from partially folded ones with the aerolysin pore. The unfolding transition curves of the destabilized variant are shifted toward the lower values of the denaturant agent compared to the wild type protein. The dynamics of the partially unfolded wild type protein follows a first-order transition. The denaturation curve obtained with the aerolysin pore is similar to that obtained with the α-hemolysin pore. The nanopore geometry or net charge does not influence the folding transition but changes the dynamics.     

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.     

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

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

Role of a solvent-exposed aromatic cluster in the folding of Escherichia coli CspA

Escherichia coli CspA is a member of the cold shock protein family. All cold shock proteins studied to date fold rapidly by an apparent two-state mechanism. CspA contains an unusual cluster of aromatic amino acids on its surface that is necessary for nucleic acid binding and also provides stability to CspA (Hillier et al., 1998). To elucidate the role this aromatic cluster plays in the determining the folding rate and pathway of CspA, we have studied the folding kinetics of mutants containing either leucine or serine substituted for Phe 18, Phe20, and/or Phe31. The leucine substitutions are found to accelerate folding and the serine substitutions to decelerate folding. Because these residues exert effects on the free energy of the folding transition state, they may be necessary for nucleating folding. They are not responsible, however, for the very compact, native-like transition state ensemble seen in the cold shock proteins, as the refolding rates of the mutants all show a similar, weak dependence of unfolding rate on denaturant concentration. Using mutant cycle analysis, we show that there is energetic coupling among the three residues between the unfolded and transition states, suggesting that the cooperative nature of these interactions helps to determine the unfolding rate. Overall, our results suggest that separate evolutionary pressures can act simultaneously on the same group of residues to maintain function, stability, and folding rate.     

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.