Abp1p and Fyn SH3 domains fold through similar low-populated intermediate states

Src homology 3 (SH3) domains are small modules that are thought to fold via a two-state mechanism, without the accumulation of significant populations of intermediate states. Relaxation dispersion NMR studies of the folding of G48V and G48M mutants of the Fyn SH3 domain have established that, at least for these modules, folding proceeds through the formation of a transient on-pathway intermediate with an equilibrium population of 1-2% that can be readily detected [Korzhnev, D. M., et al. (2004) Nature 430, 586-590]. To investigate the generality of this result, we present an (15)N relaxation dispersion NMR study of a pair of additional SH3 domains, including a G48V mutant of a stabilized Abp1p SH3 domain that shares 36% sequence identity with the Fyn SH3 module, and a A39V/N53P/V55L mutant Fyn SH3 domain. A transient folding intermediate is detected for both of the proteins studied here, and the dispersion data are well fit to a folding model of the form F <--> I <--> U, where F, I, and U correspond to folded, intermediate, and unfolded states, respectively. The temperature dependencies of the folding/unfolding rate constants were obtained so that the thermodynamic properties of each of F, I, and U could be established. The detection of I states in folding pathways of all SH3 domains examined to date via relaxation dispersion NMR spectroscopy indicates that such intermediates may well be a conserved feature in the folding of such domains in general but that their transient nature along with their low population makes detection difficult using more well-established approaches to the study of folding.     

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.     

Population shuffling between ground and high energy excited states

Stochastic processes powered by thermal energy lead to protein motions traversing time-scales from picoseconds to seconds. Fundamental to protein functionality is the utilization of these dynamics for tasks such as catalysis, folding, and allostery. A hierarchy of motion is hypothesized to connect and synergize fast and slow dynamics toward performing these essential activities. Population shuffling predicts a “top-down” temporal hierarchy, where slow time-scale conformational interconversion leads to a shuffling of the free energy landscape for fast time-scale events. Until now, population shuffling was only applied to interconverting ground states. Here, we extend the framework of population shuffling to be applicable for a system interconverting between low energy ground and high energy excited states, such as the SH3 domain mutants G48M and A39V/N53P/V55L from the Fyn tyrosine kinase, providing another tool for accessing the structural dynamics of high energy excited states. Our results indicate that the higher energy gauche− rotameric state for the leucine χ₂ dihedral angle contributes significantly to the distribution of rotameric states in both the major and minor forms of the SH3 domain. These findings are corroborated with unrestrained molecular dynamics (MD) simulations on both the major and minor states of the SH3 domain demonstrating high correlations between experimental and back-calculated leucine χ₂ rotameric populations. Taken together, we demonstrate how fast time-scale rotameric side-chain population distributions can be extracted from slow time-scale conformational exchange data further extending the scope and the applicability of the population shuffling model.     

Revealing a Concealed Intermediate that Forms after the Rate-limiting Step of Refolding of the SH3 Domain of PI3 Kinase

Kinetic and equilibrium studies of the folding and unfolding of the SH3 domain of the PI3 kinase, have been used to identify a folding intermediate that forms after the rate-limiting step on the folding pathway. Folding and unfolding, in urea as well as in guanidine hydrochloride (GdnHCl), were studied by monitoring changes in the intrinsic fluorescence or in the far-UV circular dichroism (CD) of the protein. The two probes yield non-coincident equilibrium transitions for unfolding in urea, indicating that an intermediate, I, exists in equilibrium with native (N) and unfolded (U) protein, during unfolding. Hence, the equilibrium unfolding data were analyzed according to a three-state N ↔ I ↔ U mechanism. An intermediate is observed also in kinetic unfolding studies, and its presence leads to the unfolding reaction in urea as well as in GdnHCl, occurring in two steps. The fast step is complete within the initial 11 ms of unfolding and manifests itself in a burst phase change in fluorescence. At high concentrations of GdnHCl, the entire change in fluorescence during unfolding occurs during the 11 ms burst phase. CD measurements indicate, however, that I retains N-like secondary structure. An analysis of the kinetic and thermodynamic data, according to a minimal three-state N ↔ I ↔ U mechanism, positions I after the rate-limiting transition state, TS1, of folding, on the reaction coordinate of folding in GdnHCl. Hence, I is not revealed when folding is commenced from U, regardless of the nature of the probe used to follow the folding reaction. Interrupted unfolding experiments, in which the protein is unfolded transiently in GdnHCl for various lengths of time before being refolded, showed that I refolds to N much faster than does U, confirms the analysis of the direct folding and unfolding experiments, that I is formed after the rate-limiting step of refolding in GdnHCl.     

Assessment of the Effects of Increased Relaxation Dispersion Data on the Extraction of 3-site Exchange Parameters Characterizing the Unfolding of an SH3 Domain

Recently a suite of six CPMG relaxation dispersion experiments has been described for quantifying millisecond time-scale exchange processes in proteins. The methodology has been applied to study the folding reaction of a G48M Fyn SH3 domain mutant that exchanges between the native state, and low populated unfolded and intermediate states. A complex non-linear global optimization protocol allows extraction of the kinetics and thermodynamics of the 3-site exchange process from the experimental data, as well as reconstruction of the amide group chemical shifts of the excited states. We show here, through a series of Monte-Carlo simulations on various synthetic data sets, that the 3-site exchange parameters extracted for this system on the basis of ¹⁵N single-quantum (SQ) dispersion profiles exclusively, recorded at a single temperature, are significantly in error. While a temperature dependent ¹⁵N study improves the robustness of extracted parameters, as does a combined analysis of ¹⁵N and ¹H SQ data sets measured at a single temperature, the best agreement is observed in cases where the full complement of six dispersion profiles per residue is analyzed.     

Side-chain interactions in the folding pathway of a Fyn SH3 domain mutant studied by relaxation dispersion NMR spectroscopy

A major challenge to the study of protein folding is the fact that intermediate states along the reaction pathway are generally unstable and thus difficult to observe. Recently developed NMR relaxation dispersion experiments present an avenue to accessing such states, providing kinetic, thermodynamic, and structural information for intermediates with small (greater than or equal to approximately 1%) populations at equilibrium. We have employed these techniques to study the three-state folding reaction of the G48M Fyn SH3 domain. Using (13)C-, (1)H-, and (15)N-based methods, we have characterized backbone and side-chain interactions in the folded, unfolded, intermediate, and transition states, thereby mapping the energy landscape of the protein. We find that the intermediate, populated to approximately 1%, contains nativelike structure in a central beta-sheet, and is disordered at the amino and carboxy termini. The intermediate is stabilized by side-chain van der Waals contacts, yet (13)C chemical shifts indicate that methyl-containing residues remain disordered. This state has a partially structured backbone and a collapsed yet mobile hydrophobic core and thus closely resembles a molten globule. Nonpolar side-chain contacts are formed in the unfolded-intermediate transition state; these interactions are disrupted in the intermediate-folded transition state, possibly allowing side chains to rearrange as they adopt the native packing configuration. This work illustrates the power of novel relaxation dispersion experiments in characterizing excited states that are "invisible" in even the most sensitive of NMR experiments.     

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N¹H/N²H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. ¹⁵N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.     

Folding mechanism of an ankyrin repeat protein: scaffold and active site formation of human CDK inhibitor p19(INK4d)

The p19(INK4d) protein consists of five ankyrin repeats (ANK) and controls the human cell cycle by inhibiting the cyclin D-dependent kinases (CDK) 4 and 6. We investigated the folding of p19(INK4d) by urea-induced unfolding transitions, kinetic analyses of unfolding and refolding, including double-mixing experiments and a special assay for folding intermediates. Folding is a sequential two-step reaction via a hyperfluorescent on-pathway intermediate. This intermediate is present under all conditions, during unfolding, refolding and at equilibrium. The folding mechanism was confirmed by a quantitative global fit of a consistent set of equilibrium and kinetic data revealing the thermodynamics and intrinsic folding rates of the different states. Surprisingly, the N<-->I transition is much faster compared to the I<-->U transition. The urea-dependence of the intrinsic folding rates causes population of the intermediate at equilibrium close to the transition midpoint. NMR detected hydrogen/deuterium exchange and the analysis of truncated variants showed that the C-terminal repeats ANK3-5 are already folded in the on-pathway intermediate, whereas the N-terminal repeats 1 and 2 are not folded. We suggest that during refolding, repeats ANK3-ANK5 first form the scaffold for the subsequent assembly of repeats ANK1 and ANK2. The binding function of p19(INK4d) resides in the latter repeats. We propose that the graded stability and the facile unfolding of repeats 1 and 2 is a prerequisite for the down-regulation of the inhibitory activity of p19(INK4d) during the cell-cycle.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Destabilization of the Escherichia coli RNase H kinetic intermediate: switching between a two-state and three-state folding mechanism

Escherichia coli RNase H folds through a partially folded kinetic intermediate that mirrors a rarely populated, partially unfolded form detectable by native-state hydrogen exchange under equilibrium conditions. Residue 53 is at the interface of two helices known to be structured in this intermediate. Kinetic refolding studies on mutant proteins varying in size and hydrophobicity at residue 53 support a contribution of hydrophobicity to the stabilities of the kinetic intermediate and the transition state. Packing interactions also play a significant role in the stability of these two states, though they play a much larger role in the native-state stability. One dramatic mutation, I53D, results in the conversion from a three-state to a two-state folding mechanism, which is explained most easily through a simple destabilization of the kinetic intermediate such that it is no longer stable with respect to the unfolded state. These results demonstrate that interactions that stabilize an intermediate can accelerate folding if these same interactions are present in the transition state. Our results are consistent with a hierarchical model of folding, where the intermediate consists of native-like interactions, is on-pathway, and is productive for folding.     

Characterization of the folding kinetics of a three-helix bundle protein via a minimalist Langevin model.

We use a simple off-lattice Langevin model of protein folding to characterize the folding and unfolding of a fast-folding, 46 residue three-helix bundle. Under conditions at which the C-terminal helix is 30 % stable, we observe a clear three-state folding mechanism. In the on-pathway intermediate state, the middle and C-terminal helices are folded and in contact with each other, while the N-terminal region remains disordered. Nevertheless, under these conditions this intermediate is thermodynamically unstable relative to its unfolded state. The first and highest folding barrier corresponds to the organization of the hinge between the middle and C-terminal helices. A subsequent major barrier corresponds to the organization of the hinge between the middle and N-terminal helices. Hyperstabilizing the hinge regions leads to twice the folding rate that is obtained from hyperstabilizing the helices, even though much fewer contacts are involved in hinge hyperstabilization than in helix hyperstabilization. Unfolding follows single-exponential kinetics, even at temperatures only slightly above the folding transition temperature.     

An NMR view of the folding process of a CheY mutant at the residue level

The folding of CheY mutant F14N/V83T was studied at 75 residues by NMR. Fluorescence, NMR, and sedimentation equilibrium studies at different urea and protein concentrations reveal that the urea-induced unfolding of this CheY mutant includes an on-pathway molten globule-like intermediate that can associate off-pathway. The populations of native and denatured forms have been quantified from a series of 15N-1H HSQC spectra recorded under increasing concentrations of urea. A thermodynamic analysis of these data provides a detailed picture of the mutant's unfolding at the residue level: (1) the transition from the native state to the molten globule-like intermediate is highly cooperative, and (2) the unfolding of this state is sequential and yields another intermediate showing a collapsed N-terminal domain and an unfolded C-terminal tail. This state presents a striking similarity to the kinetic transition state of the CheY folding pathway.     

The kinetic pathway of folding of barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.     

Identification of a Conserved Aggregation-Prone Intermediate State in the Folding Pathways of Spc-SH3 Amyloidogenic Variants

We compared the folding pathways of selected mutational variants of the α-spectrin SH3 domain (Spc-SH3) by using a continuum model that combines a full atomistic protein representation with the Gō potential. Experimental data show that the N47G mutant shows very little tendency to aggregate while the N47A and triple mutant D48G(2Y) are both amyloidogenic, with the latter being clearly more aggregation prone. We identified a strikingly similar native-like folding intermediate across the three mutants, in which strand β(1) is totally unstructured and more than half of the major hydrophobic core residues are highly solvent exposed. Results from extensive docking simulations show that the ability of the intermediates to dimerize is largely driven by strand β(1) and is consistent with the in vitro aggregation behavior reported for the corresponding mutants. They further suggest that residues 44 and 53, which are key players in the nucleation-condensation mechanism of folding, are also important triggers of the aggregation process.     

Equilibrium amide hydrogen exchange and protein folding kinetics

The classical Linderstrøm-Lang hydrogen exchange (HX) model is extended to describe the relationship between the HX behaviors (EX1 and EX2) and protein folding kinetics for the amide protons that can only exchange by global unfolding in a three-state system including native (N), intermediate (I), and unfolded (U) states. For these slowly exchanging amide protons, it is shown that the existence of an intermediate (I) has no effect on the HX behavior in an off-pathway three-state system (IUN). On the other hand, in an on-pathway three-state system (UIN), the existence of a stable folding intermediate has profound effect on the HX behavior. It is shown that fast refolding from the unfolded state to the stable intermediate state alone does not guarantee EX2 behavior. The rate of refolding from the intermediate state to the native state also plays a crucial role in determining whether EX1 or EX2 behavior should occur. This is mainly due to the fact that only amide protons in the native state are observed in the hydrogen exchange experiment. These new concepts suggest that caution needs to be taken if one tries to derive the kinetic events of protein folding from equilibrium hydrogen exchange experiments.     

The effect of glycosylation on the folding kinetics of erythropoietin

Glycosylation is a common posttranslational modification that generally increases protein solubility and thermodynamic stability. Less is known about how this modification influences protein folding, particularly folding processes involving intermediate species. In the present report, folding comparisons of a nonglycosylated erythropoietin (EPO) mutant are made with the fully glycosylated EPO, which was recently shown to fold by a three-state on-pathway mechanism. The absence of glycosylation did not alter the folding mechanism of EPO but did greatly decrease the stability of the intermediate species, change the rate-limiting step of the folding reaction, and accelerate the folding kinetics to both the intermediate state and the native state. Surprisingly, glycosylation stabilized the intermediate species to a greater extent than it increased the EPO equilibrium stability. These results suggest that glycosylation impedes the latter EPO folding steps rather than accelerating them by biasing particular folding pathways, as previously proposed for folding reactions initiated from unfolded ensembles with minimal residual structure. Due to the specific biological processes modulated by EPO glycosylation, however, there may be little evolutionary pressure to fold on a faster, more direct pathway at the expense of biological function, particularly given the protective role glycosylation has at preventing EPO aggregation. Lastly, evidence that is consistent with glycosylation destabilizing the unfolded state to some degree and contributing to the greater equilibrium stability of the glycosylated EPO is presented.     

Some NMR Experiments and a Structure Determination Employing a {15N,2H} Enriched Protein

We present the results of studies of an aqueous sample of a highly {15N,2H} enriched protein, the SH3 domain from Fyn. Measurements of 1H relaxation and interactions between H2O solvent and exchangeable protons are given, as well as a method for increasing the effective longitudinal relaxation of solvent exchangeable proton resonances. The long-range isotope shifts are measured, for 1H and 15N, which arise due to perdeuteration. Simulations, which employed a 7 or 8 spin relaxation matrix analysis, were compared to the experimental data from a time series of 2D NOESY datasets for some resonances. The agreement between experiment and simulation suggest that, with this 1H dilute sample, relatively long mixing times (up to 1.2 s) can be used to detect specific dipolar interactions between amide protons up to about 7Å apart. A set of 155 inter-amide NOEs and 7 side chain NOEs were thus identified in a series of 3D HSQC-NOESY-HSQC experiments. These data, alone and in combination with previously collected restraints, were used to calculate sets of structures using X-PLOR. These results are compared to the available X-ray and NMR structures of the Fyn SH3 domain.     

Probing the influence on folding behavior of structurally conserved core residues in P. aeruginosa apo-azurin

The effects on folding kinetics and equilibrium stability of core mutations in the apo-mutant C112S of azurin from Pseudomonas aeruginosa were studied. A number of conserved residues within the cupredoxin family were recognized by sequential alignment as constituting a common hydrophobic core: I7, F15, L33, W48, F110, L50, V95, and V31. Of these, I7, V31, L33, and L50 were mutated for the purpose of obtaining information on the transition state and a potential folding nucleus. In addition, residue V5 in the immediate vicinity of the common core, as well as T52, separate from the core, were mutated as controls. All mutants exhibited a nonlinear dependence of activation free energy of folding on denaturant concentration, although the refolding kinetics of the V31A/C112S mutant indicated that the V31A mutation destabilizes the transition state enough to allow folding via a parallel transition state ensemble. Phi-values could be calculated for three of the six mutants, V31A/C112S, L33A/C112S, and L50A/C112S, and the fractional values of 0.63, 0.33, and 0.50 (respectively) obtained at 0.5 M GdmCl suggest that these residues are important for stabilizing the transition state. Furthermore, a linear dependence of ln k(obs)(H2O) on DeltaG(U-N)(H2O) of the core mutations and the putative involvement of ground-state effects suggest the presence of native-like residual interactions in the denatured state that bias this ensemble toward a folding-competent state.