An on-pathway hidden intermediate and the early rate-limiting transition state of Rd-apocytochrome b562 characterized by protein engineering

The folding pathway of Rd-apocytochrome b562, a four-helix bundle protein, was characterized using Trp and Ala/Gly pair mutations. We found that the Trp mutants (F65W) of both the fully folded Rd-apocytochrome b562 and a partially unfolded intermediate with the N-terminal helix (helix I) unfolded, fold with identical folding rates, providing direct evidence for the conclusion that the rate-limiting transition state folds before the partially unfolded intermediate; and that this hidden intermediate is an on-pathway intermediate. We further characterized the helical structures formed in the rate-limiting transition state by measuring the folding/unfolding rates for Ala/Gly pair mutations at solvent-exposed positions. Little change in folding rates occurred for the Ala/Gly pair mutations at positions in helix I and the C-terminal regions of helix II and IV. In contrast, a significant difference in folding rates was observed for the Ala/Gly pair mutations in helix III and the N-terminal regions of helix II and IV, suggesting that helix III and the N-terminal regions of helix II and IV are formed in the rate-limiting transition state. These results complement those obtained from earlier studies and help to define the folding pathway of Rd-apocytochrome b562 in more detail.     

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.     

Structure and function in bacteriorhodopsin: the effect of the interhelical loops on the protein folding kinetics

The loops connecting the seven transmembrane helices of bacteriorhodopsin have each been replaced in turn by structureless linkers of Gly-Gly-Ser repeat sequences, and the effect on the protein folding kinetics has been determined. An SDS-denatured state of each loop mutant bacterio-opsin was folded in l-alpha-1,2-dihexanoylphosphatidylcholine/l-alpha-1,2-dimyristoylphosphatidylcholine micelles, containing retinal, to give functional bacteriorhodopsin. Stopped-flow mixing was used to initiate the folding reaction, giving a time resolution of milliseconds, and changes in protein fluorescence were used to monitor folding. All loop mutant proteins folded according to the same reaction scheme as wild-type protein. The folding kinetics of the AB, BC and DE loop mutants were the same as wild-type protein, despite the blue-shifted chromophore band of the BC loop mutant bR state. A partially folded apoprotein intermediate state of the AB loop mutant did however appear to decay in the absence of retinal. The most significant effects on the folding kinetics were seen for mutant protein with structureless linkers in place of the CD, EF and FG loops. The rate-limiting apoprotein folding step of the CD loop mutant was about ten times slower than wild-type, whilst that of the EF loop mutant was almost four times slower than wild-type. Wild-type behaviour was observed for the other folding and retinal binding events of the CD and EF loop mutant proteins. These effects of the CD and EF loop mutations on apoprotein folding correlate with the fact that these two loop mutants also have the least stable, partially folded apoprotein intermediate of all the loop mutants, and are the most affected by a decrease in lipid lateral pressure. In contrast, the FG loop mutant exhibited wild-type apoprotein folding, but altered covalent binding of retinal and final folding to bacteriorhodopsin. This correlates with the fact that the FG loop mutant bacteriorhodopsin is the most susceptible to denaturation by SDS of all the loop mutants, but its partially folded apoprotein intermediate is more stable than that of the CD and EF mutants. Thus the CD and EF loops may contribute to the transition state for the rate-limiting apoprotein folding step and the FG loop to that for final folding and covalent binding of retinal.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.     

Sequential barriers and an obligatory metastable intermediate define the apparent two-state folding pathway of the ubiquitin-like PB1 domain of NBR1

The 90-residue N-terminal Phox and Bem1p (PB1) domain of NBR1 forms an α/β ubiquitin-like fold. Kinetic analysis using stopped-flow fluorescence reveals two-state kinetics; however, nonlinear effects in the denaturant dependence of the unfolding data demonstrate changes in the position of the rate-limiting barrier along the folding coordinate as the folding conditions change. The kinetics of wt-PB1 and several mutants show that this curvature is consistent with a single-pathway mechanism involving sequential transition states (TS1 and TS2) separated by a transiently populated high-energy intermediate, rather than movement of the transition state on a broad energy plateau. We show that the two transition states within the sequential model represent structurally and thermodynamically distinct species. TS1 is a collapsed state (α_TS1 = 0.71) with a large enthalpic barrier to formation that is rate-limiting under conditions that strongly favour folding. TS2 is highly native-like (α_TS2 = 0.93) and represents a late entropic barrier to formation of the native state. In support of the sequential transition state mechanism, we show that the G62A helix 2 substitution stabilises TS1 and the intermediate to such an extent that the latter becomes significantly populated, leading to the observation of a fast kinetic phase representing the initial U → I transition, with TS2 (α_TS2 = 0.87) becoming rate-limiting. The folding rate is not retarded by populating an intermediate, which would be expected for a misfold state, but is accelerated, suggesting that the I state is productive and on-pathway. The results show that the apparent two-state folding of the wt-PB1 domain occurs along a well-defined pathway involving structurally and thermodynamically distinct sequential transition states and an obligatory metastable intermediate that represents a productive local minimum in the energy landscape that increases the efficiency of barrier crossing through favourable effects on the entropy of activation.     

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 α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). 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.     

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.     

Energetic basis of structural stability in the molten globule state: alpha-lactalbumin

The denatured states of alpha-lactalbumin, which have features of a molten globule state, have been studied to elucidate the energetics of the molten globule state and its contribution to the stability of the native conformation. Analysis of calorimetric and CD data shows that the heat capacity increment of alpha-lactalbumin denaturation highly correlates with the degree of disorder of the residual structure of the state. As a result, the denaturational transition of alpha-lactalbumin from the native to a highly ordered compact denatured state, and from the native to the disordered unfolded state are described by different thermodynamic functions. The enthalpy and entropy of the denaturation of alpha-lactalbumin to compact denatured state are always greater than the enthalpy and entropy of its unfolding. This difference represents the unfolding of the molten globule state. Calorimetric measurements of the heat effect associated with the unfolding of the molten globule state reveal that it is negative in sign over the temperature range of molten globule stability. This observation demonstrates the energetic specificity of the molten globule state, which, in contrast to a protein with unique tertiary structure, is stabilized by the dominance of negative entropy and enthalpy of hydration over the positive conformational entropy and enthalpy of internal interactions. It is concluded that at physiological temperatures the entropy of dehydration is the dominant factor providing stability for the compact intermediate state on the folding pathway, while for the stability of the native state, the conformational enthalpy is the dominant factor.     

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.     

Progressive stabilization of intermediate and transition states in protein folding reactions by introducing surface hydrophobic residues

It can be argued from the principle of solvent exclusion that the introduction of hydrophobic residues onto the surface of a protein will not destabilize the folded state because the nonpolar side chain will be at least as exposed in the unfolded state as it is when the protein chain is folded. A comparison of the folding pathway of wild type and 11 site-directed mutants of CD2.d1 shows this to be true. In fact, owing to partial burial of nonpolar groups as folding proceeds, we find that the rapidly formed intermediate state and, to a greater extent, the transition state are generally stabilized by hydrophobic surface mutations. This effect is slightly moderated in the folded state presumably by the perturbation of van der Waals' contacts and/or local electrostatic interactions that have a greater influence in this fully compact structure. The fact that in all but one case we find that stabilization of the rapidly collapsed intermediate is accompanied by a faster acquisition of the folded state refutes the argument that I states are generally "off pathway" conformations or ensembles that lead to the inhibition of otherwise more rapid folding trajectories.     

Free energy landscape and folding mechanism of a beta-hairpin in explicit water: a replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.     

Exploring the temperature-pressure phase diagram of staphylococcal nuclease

The temperature dependence of the pressure-induced equilibrium unfolding of staphylococcal nuclease (Snase) was determined by fluorescence of the single tryptophan residue, FTIR absorption for the amide I' and tyrosine O-H bands, and small-angle X-ray scattering (SAXS). The results from these three techniques were similar, although the stability as measured by fluorescence was slightly lower than that measured by FTIR and SAXS. The resulting phase diagram exhibits the well-known curvature for heat and cold denaturation of proteins, due to the large decrease in heat capacity upon folding. The volume change for unfolding became less negative with increasing temperatures, consistent with a larger thermal expansivity for the unfolded state than for the folded state. Fluorescence-detected pressure-jump kinetics measurements revealed that the curvature in the phase diagram is due primarily to the rate constant for folding, indicating a loss in heat capacity for the transition state relative to the unfolded state. The similar temperature dependence of the equilibrium and activation volume changes for folding indicates that the thermal expansivities of the folded and transition states are similar. This, along with the fact that the activation volume for folding is positive over the temperature range examined, the nonlinear dependence of the folding rate constant upon temperature implicates significant dehydration in the rate-limiting step for folding of Snase.     

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.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.     

Kinetic analysis of folding and unfolding the 56 amino acid IgG-binding domain of streptococcal protein G.

The 56 amino acid B domain of protein G (GB) is a stable globular folding unit with no disulfide cross-links. The physical properties of GB offer extraordinary flexibility for evaluating the energetics of the folding reaction. The protein is monomeric and very soluble in both folded and unfolded forms. The folding reaction has been previously examined by differential scanning calorimetry (Alexander et al., 1992) and found to exhibit two-state unfolding behavior over a wide pH range with an unfolding transition near 90 degrees C (GB1) at neutral pH. Here, the kinetics of folding and unfolding two naturally occurring versions of GB have been measured using stopped-flow mixing methods and analyzed according to transition-state theory. GB contains no prolines, and the kinetics of folding and unfolding can be fit to a single, first-order rate constant over the temperature range of 5-35 degrees C. The major thermodynamic changes going from the unfolded state to the transition state are (1) a large decrease in heat capacity (delta Cp), indicating that the transition state is compact and solvent inaccessible relative to the unfolded state; (2) a large loss of entropy; and (3) a small increase in enthalpy. The most surprising feature of the folding of GB compared to that of previously studied proteins is that its folding approximates a rapid diffusion controlled process with little increase in enthalpy going from the unfolded to the transition state.     

COSMIC analysis of the major alpha-helix of barnase during folding

The structures of transition states and intermediates in protein folding may be analysed by protein engineering methods that remove simple interactions that stabilize the folded state. We have now extended the range and reliability of the procedure by using the COSMIC (Combination of Sequential Mutant Interaction Cycles) technique, in which a series of double-mutant cycles is constructed. In each cycle, the side-chains of two amino acid residues that interact in the folded state are mutated separately and together. Kinetic and equilibrium measurements on folding for each cycle show unambiguously whether or not two residues interact during protein folding. A series of such cycles has been constructed to leapfrog along the major alpha-helix of barnase, comprising residues 6 to 18. The helix is found to be intact from its C terminus to residue 12 but begins to unwind towards the N terminus in both the transition state for unfolding and in a folding intermediate.     

Mechanism of folding and assembly of a small tetrameric protein domain from tumor suppressor p53

We have analyzed the folding pathway of the tetramerization domain of the tumor suppressor protein p53. Structures of transition states were determined from phi-values for 25 mutations, including leucine to norvaline, and the analysis encompassed nearly every residue in the domain. Denatured monomers fold and dimerize, through a transition state with little native structure, to form a transient, highly structured dimeric intermediate. The intermediate dimerizes, through a native-like transition state with the primary dimers fully folded but with interdimer interactions only partially formed, to form the native tetramer as a 'dimer of dimers'.