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.     

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.     

Structural and kinetic characterization of the simplified SH3 domain FP1

The simplified SH3 domain sequence, FP1, obtained in phage display selection experiments has an amino acid composition that is 95% Ile, Lys, Glu, Ala, Gly. Here we use NMR to investigate the tertiary structure of FP1. We find that the overall topology of FP1 resembles that of the src SH3 domain, the hydrogen-deuterium exchange and chemical shift perturbation profiles are similar to those of naturally occurring SH3 domains, and the (15)N relaxation rates are in the range of naturally occurring small proteins. Guided by the structure, we further simplify the FP1 sequence and compare the effects on folding kinetics of point mutations in FP1 and the wild-type src SH3 domain. The results suggest that the folding transition state of FP1 is similar to but somewhat less polarized than that of the wild-type src SH3 domain.     

Backbone dynamics in an intramolecular prolylpeptide-SH3 complex from the diphtheria toxin repressor, DtxR

The diphtheria toxin repressor contains an SH3-like domain that forms an intramolecular complex with a proline-rich (Pr) peptide segment and stabilizes the inactive state of the repressor. Upon activation of diphtheria toxin repressor (DtxR) by transition metals, this intramolecular complex must dissociate as the SH3 domain and Pr segment form different interactions in the active repressor. Here we investigate the dynamics of this intramolecular complex using backbone amide nuclear spin relaxation rates determined using NMR spectroscopy and molecular dynamics trajectories. The SH3 domain in the unbound and bound states showed typical dynamics in that the secondary structures were fairly ordered with high generalized order parameters and low effective correlation times, while residues in the loops connecting β-strands exhibited reduced generalized order parameters and required additional motional terms to adequately model the relaxation rates. Residues forming the Pr segment exhibited low-order parameters with internal rotational correlation times on the order of 0.6 ns-1 ns. Further analysis showed that the SH3 domain was rich in millisecond time scale motions while the Pr segment exhibited motions on the 100 μs time scale. Molecular dynamics simulations indicated structural rearrangements that may contribute to the observed relaxation rates and, together with the observed relaxation rate data, suggested that the Pr segment exhibits a binding ↔ unbinding equilibrium. The results here provide new insights into the nature of the intramolecular complex and provide a better understanding of the biological role of the SH3 domain in regulating DtxR activity.     

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.     

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.     

Identification of a collapsed intermediate with non-native long-range interactions on the folding pathway of a pair of Fyn SH3 domain mutants by NMR relaxation dispersion spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central beta-sheet is at least partially formed, but without interactions between the NH2- and COOH-terminal beta-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F<-->I<-->U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (< 25 degrees C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F<-->I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.     

Structural insight into the binding diversity between the Tyr-phosphorylated human ephrinBs and Nck2 SH2 domain

The binding interaction between the Nck2 SH2 domain and the phosphorylated ephrinB initiates a critical pathway for the reverse signaling network mediated by Eph receptor-ephrinB. Previously, the NMR structure and Tyr phosphorylations of the human ephrinB cytoplasmic domain have been studied. To obtain a complete story, it would be of significant interest to determine the structure of the Nck2 SH2 domain that shows a low sequence identity to other SH2 domains with known structures. Here, we report the determination of the solution structure of the human Nck2 SH2 domain and investigate its interactions with three phosphorylated ephrinB fragments by NMR spectroscopy. The results indicate that: 1) although the human Nck2 SH2 domain adopts a core tertiary fold common to all SH2 domains, it owns some unique properties such as a shorter C-terminal helix and unusual electrostatic potential surface. However, the most striking finding is that the C-terminal tail of the human Nck2 SH2 domain adopts a short antiparallel beta-sheet that, to the best of our knowledge, has never been identified in other SH2 domains. The truncation study suggests that one function of the C-terminal tail is to control the folding/solubility of the SH2 domain. 2) In addition to [Tyr(P)304]ephrinB2(301-322) and [Tyr(P)316]ephrinB2(301-322), here we identified [Tyr(P)330]ephrinB2(324-333) also capable of binding to the SH2 domain. The detailed NMR study indicated that the binding mechanisms for the three ephrinB fragments might be different. The binding with [Tyr(P)304]-ephrinB2(301-322) and [Tyr(P)316]ephrinB2(301-322) might be mostly involved in the residues over the N-half of the SH2 domain and provoked a significant increase in the backbone and side chain dynamics of the SH2 domain on the microsecond-millisecond time scale. In contrast, binding with [Tyr(P)330]ephrinB2(324-333) might have most residues over both halves engaged but induced less profound conformational dynamics on the mus-ms time scale.     

Characterization of early and late transition states of the folding pathway of a SH2 domain

Albeit SH2 domains are abundant protein–protein interaction modules with fundamental roles in the regulation of several physiological and molecular pathways in the cell, the available information about the determinants of their thermodynamic stability and folding properties are still very limited. In this work, we provide a quantitative characterization of the folding pathway of the C‐terminal SH2 domain of SHP2, conducted through a combination of site‐directed mutagenesis and kinetic (un)folding experiments (Φ‐value analysis). The energetic profile of the folding reaction of the C‐SH2 domain is described by a three‐state mechanism characterized by the presence of two transition states and a high‐energy intermediate. The production of 29 site‐directed variants allowed us to calculate the degree of native‐like interactions occurring in the early and late events of the folding reaction. Data analysis highlights the presence of a hydrophobic folding nucleus surrounded by a lower degree of structure in the early events of folding, further consolidated as the reaction proceeds towards the native state. Interestingly, residues physically located in the functional region of the domain reported unusual Φ‐values, a hallmark of the presence of transient misfolding. We compared our results with previous ones obtained for the N‐terminal SH2 domain of SHP2. Notably, a conserved complex folding mechanism implying the presence of a folding intermediate arise from comparison, and the relative stability of such intermediate appears to be highly sequence dependent. Data are discussed under the light of previous works on SH2 domains.     

Molecular dynamics simulations from putative transition states of alpha-spectrin SH3 domain

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff-Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental Phi - values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 ns. Complete folding was observed in one case, while in the other simulations partial folding or unfolding events were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states.     

Folding of the alphaII-spectrin SH3 domain under physiological salt conditions

The SH3 domain has often been used as a model for protein folding due to its typical two-state behaviour. However, recent experimental data at low pH as well as molecular dynamic simulations have indicated that the folding process of SH3 probably is more complicated, and may involve intermediate states. Using both kinetic and equilibrium measurements we have obtained evidence that under native-like conditions the folding of the spectrin SH3 domain does not follow a classic two-state behaviour. The curvature we observed in the Chevron plots is a strong indication of a non-linear activation energy relationship due to the presence of high-energy intermediates. In addition, circular dichroism measurements indicated that refolding after thermal denaturation did not follow the same pattern as thermal unfolding but rather implied less cooperativity and that the refolding transition increased with increasing protein concentration. Further, NMR experiments indicated that upon refolding the SH3 domain gave rise to more than one conformation. Therefore, our results suggest that the folding of the SH3 domain of αII-spectrin does not follow a classical two-state process under high-salt conditions and neutral pH. Heterogeneous folding pathways, which can include folding intermediates as well as misfolded intermediates, might give a more reasonable insight into the folding behaviour of the αII-spectrin SH3 domain.     

The folding transition state between SH3 domains is conformationally restricted and evolutionarily conserved.

The protein engineering analysis of the alpha-spectrin SH3 domain at three different stability conditions (pH 7.0, 3.5 and 2.5) reveals a folding transition state structured around the distal loop beta-hairpin and the 310-helix. This region is impervious to overall changes in protein stability, suggesting a transition state ensemble with little conformational variability. Comparison with the Src SH3 domain (36% sequence homology) indicates that the transition state in this protein family may be conserved. Discrepancies at some positions can be rationalized in terms of the different interactions made by the different side chains in both domains. Br?nsted plot analysis confirms the straight phi(doubledagger-U) results and shows two folding subdomains for this small protein. These results, together with previous data on circular permutants of the alpha-spectrin SH3 domain, indicate that polypeptide topology and chain connectivity play a major role in the folding reaction of this protein family.     

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.     

NMR dynamics-derived insights into the binding properties of a peptide interacting with an SH2 domain

The signal transduction protein phospholipase C-gamma1 (PLC-gamma1) is activated when its C-terminal SH2 domain (PLCC) binds the phosphorylated Tyr-1021 site (pTyr-1021) in the beta-platelet-derived growth factor receptor (PDGFR). To better understand the contributions that dynamics make to binding, we have used NMR relaxation experiments to investigate the motional properties of backbone amide and side chain methyl groups in a peptide derived from the pTyr-1021 site of PDGFR, both free and in complex with the PLCC SH2 domain. The free peptide has relaxation properties that are typical for a small, unstructured polymer, while the backbone of the bound peptide is least flexible for residues in the central portion of the binding site with the amplitude of pico- to nanosecond time scale motions increasing toward the C-terminus of the peptide. The increase in large amplitude motion toward the end of the pY1021 peptide is consistent with the bound peptide existing as an ensemble of states with C-terminal residues having the broadest distribution of backbone conformations, while residues in the central binding site are the most restricted. Deuterium spin relaxation experiments establish that the protein-peptide interface is highly dynamic, and this mobility may play an important role in modulating the affinity of the interaction.     

Altered dynamics in Lck SH3 upon binding to the LBD1 domain of Herpesvirus saimiri Tip

The Tip protein from Herpesvirus saimiri interacts with the SH3 domain from the Src-family kinase Lck via a proline-containing sequence termed LBD1. Src-family kinase SH3 domains related to Lck have been shown to be dynamic in solution and partially unfold under physiological conditions. The rate of such partial unfolding is reduced by viral protein binding. To determine if the Lck SH3 domain displayed similar behavior, the domain was investigated with hydrogen exchange and mass spectrometry. Lck SH3 was found to be highly dynamic in solution. While other SH3 domains require as much as 10,000 sec to become totally deuterated, Lck SH3 became almost completely labeled within 200 sec. A partial unfolding event involving 8-10 residues was observed with a half-life of approximately 10 sec. Tip LBD1 binding did not cause gross structural changes in Lck SH3 but globally stabilized the domain and reduced the rate of partial unfolding by a factor of five. The region of partial unfolding in Lck SH3 was found to be similar to that identified for other SH3 domains that partially unfold. Although the sequence conservation between Lck SH3 and other closely related SH3 domains is high, the dynamics do not appear to be conserved.     

Circularization changes the folding transition state of the src SH3 domain

Native state topology has been implicated as a major determinant of protein-folding mechanisms. Here, we test experimentally the robustness of the src SH3-domain folding transition state to changes in topology by covalently constraining regions of the protein with disulfide crosslinks and then performing kinetic analysis on point mutations in the context of these modified proteins. Circularization (crosslinking the N and C termini) of the src SH3 domain makes the protein topologically symmetric and causes delocalization of structure in the transition state ensemble suggesting a change in the folding mechanism. In contrast, crosslinking a single structural element (the distal beta-hairpin) which is an essential part of the transition state, results in a protein that folds 30 times faster, but does not change the distribution of structure in the transition state. As the transition states of distantly related SH3 domains were previously found to be very similar, we conclude that the free energy landscape of this protein family contains deep features which are relatively insensitive to sequence variations but can be altered by changes in topology.     

Cooperative propagation of local stability changes from low-stability and high-stability regions in a SH3 domain

Site-directed mutagenesis has been used to produce local stability changes at two regions of the binding site surface of the alpha-spectrin SH3 domain (Spc-SH3) differing in their intrinsic stability. Mutations were made at residue 56, located at the solvent-exposed side of the short 3(10) helix, and at residue 21 in the tip of the flexible RT-loop. NMR chemical-shift analysis and X-ray crystallography indicated negligible changes produced by the mutations in the native structure limited to subtle rearrangements near the mutated residue and at flexible loops. Additionally, mutations do not alter importantly the SH3 binding site structure, although produce significant changes in its affinity for a proline-rich decapeptide. The changes in global stability measured by differential scanning calorimetry are consistent the local energy changes predicted by theoretical models, with the most significant effects observed for the Ala-Gly mutations. Propagation of the local stability changes throughout the domain structure has been studied at a per-residue level of resolution by NMR-detected amide hydrogen-deuterium exchange (HX). Stability propagation is remarkably efficient in this small domain, apparently due to its intrinsically low stability. Nevertheless, the HX-core of the domain is not fully cooperative, indicating the existence of co-operative subunits within the core, which is markedly polarized. An equilibrium phi-analysis of the changes in the apparent Gibbs energies of HX per residue produced by the mutations has allowed us to characterize structurally the conformational states leading to HX. Some of these states resemble notably the folding transition state of the Spc-SH3 domain, suggesting a great potential of this approach to explore the folding energy landscape of proteins. An energy perturbation propagates more effectively from a flexible region to the core than in the opposite direction, because the former affects a broader region of the energy landscape than the latter. This might be of importance in understanding the special thermodynamic signature of the SH3-peptide interaction and the relevance of the dual character of SH3 binding sites.     

NMR characterizations of an amyloidogenic conformational ensemble of the PI3K SH3 domain

Amyloid formation is associated with structural changes of native polypeptides to monomeric intermediate states and their self-assembly into insoluble aggregates. Characterizations of the amyloidogenic intermediate state are, therefore, of great importance in understanding the early stage of amyloidogenesis. Here, we present NMR investigations of the structural and dynamic properties of the acid-unfolded amyloidogenic intermediate state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain--a model peptide. The monomeric amyloidogenic state of the SH3 domain studied at pH 2.0 (35 degrees C) was shown to be substantially disordered with no secondary structural preferences. (15)N NMR relaxation experiments indicated that the unfolded polypeptide is highly flexible on a subnanosecond timescale when observed under the amyloidogenic condition (pH 2.0, 35 degrees C). However, more restricted motions were detected in residues located primarily in the beta-strands as well as in a loop in the native fold. In addition, nonnative long-range interactions were observed between the residues with the reduced flexibility by paramagnetic relaxation enhancement (PRE) experiments. These indicate that the acid-unfolded state of the SH3 domain adopts a partly folded conformation through nonnative long-range contacts between the dynamically restricted residues at the amyloid-forming condition.