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.     

Hydration and packing along the folding pathway of SH3 domains by pressure-dependent NMR

The volumetric properties associated with protein folding transitions reflect changes in protein packing and hydration of the states that participate in the folding reaction. Here, NMR spin relaxation techniques are employed to probe the folding-unfolding kinetics of two SH3 domains as a function of pressure so that the changes in partial molar volumes along the folding pathway can be measured. The two domains fold with rates that differ by approximately 3 orders of magnitude, so their folding dynamics must be probed using different NMR relaxation experiments. In the case of the drkN SH3 domain that folds via a two-state mechanism on a time scale of seconds, nitrogen magnetization exchange spectroscopy is employed, while for the G48M mutant of the Fyn SH3 domain where the folding occurs on the millisecond time scale (three-step reaction), relaxation dispersion experiments are utilized. The NMR methodology is extremely sensitive to even small changes in equilibrium and rate constants, so reliable estimates of partial molar volumes can be obtained using low pressures (1-120 bar), thus minimizing perturbations to any of the states along the folding reaction coordinate. The volumetric data that were obtained are consistent with a similar folding mechanism for both SH3 domains, involving early chain compaction to states that are at least partially hydrated. This work emphasizes the role of NMR spin relaxation in studying dynamic processes over a wide range of time scales.     

Multiple folding pathways of the SH3 domain

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

Parallel folding pathways in the SH3 domain protein

The transition-state ensemble (TSE) is the set of protein conformations with an equal probability to fold or unfold. Its characterization is crucial for an understanding of the folding process. We determined the TSE of the src-SH3 domain protein by using extensive molecular dynamics simulations of the Go model and computing the folding probability of a generated set of TSE candidate conformations. We found that the TSE possesses a well-defined hydrophobic core with variable enveloping structures resulting from the superposition of three parallel folding pathways. The most preferred pathway agrees with the experimentally determined TSE, while the two least preferred pathways differ significantly. The knowledge of the different pathways allows us to design the interactions between amino acids that guide the protein to fold through the least preferred pathway. This particular design is akin to a circular permutation of the protein. The finding motivates the hypothesis that the different experimentally observed TSEs in homologous proteins and circular permutants may represent potentially available pathways to the wild-type protein.     

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.     

The apomyoglobin folding pathway revisited: structural heterogeneity in the kinetic burst phase intermediate

Extensive analysis of accurate quench-flow hydrogen exchange results indicates that the burst phase kinetic intermediate in the folding of apomyoglobin (apoMb) from urea is structurally heterogeneous. The structural variability is associated with the partial folding of the E helix during the burst phase (<6.4ms) of the folding process. Analysis of the effects of exchange-out of amide proton labels during the labeling pulse ( approximately pH 10) of the quench-flow process indicates that three of the amide protons in the E helix are in fact largely protected in the burst phase of folding, while the remainder of the E helix has a substantial complement of amide protons that show biphasic kinetics, i.e. are protected partly during the burst phase and partly during the slow phase of folding. The locations of these amide protons can be used to map the sites of structural heterogeneity in the kinetic molten globule. These sites include, besides the E helix, the ends of the A and B helices and part of the C helix. Our results give significant support to the hypothesis that the kinetic molten globule intermediate of apoMb is native-like.     

Hydrophobic core packing in the SH3 domain folding transition state.

How tightly packed is the hydrophobic core of a folding transition state structure? We have addressed this question by characterizing the effects on folding kinetics of > 40 substitutions of both large and small amino acids in the hydrophobic core of the Fyn SH3 domain. Our results show that residues at three positions, which we designate as the 'core folding nucleus', are tightly packed in the transition state, and substitutions at these positions cause the largest changes in the folding rate. The other six positions examined appear to be loosely packed; thus, substitutions at these positions with larger hydrophobic residues generally accelerate folding, presumably by increasing the rate of nonspecific hydrophobic collapse. Surprisingly, the folding rate can be greatly accelerated by residues that also significantly destabilize the native state structure. Furthermore, mutants with identical thermodynamic stability can differ by up to 55-fold in their folding rates. These results highlight the importance of hydrophobic core composition, as opposed to only topology, in determining the folding rate of a protein. They also provide a new explanation for the 'abnormal' phi-values observed in many protein folding kinetics studies.     

Understanding the folding mechanism of an alpha-helical hairpin

The alpha-helical hairpin is the fundamental building block of the widespread helix-turn-helix DNA binding motif. With two antiparallel helices connected by a reverse turn, the alpha-helical hairpin structure may be regarded as a "supersecondary structural element" and, therefore, could exhibit rather unique folding properties. So far, the folding mechanism of alpha-helical hairpins has not been studied in detail and remains elusive. Herein, we examine the effects of the turn, the hydrophobic cluster, and a disulfide cross-linker on the folding kinetics of a designed alpha-helical hairpin, Z34C, using an infrared temperature-jump (T-jump) method in conjunction with site-specific mutagenesis. Our results show that Z34C folds with an ultrafast rate ( approximately 4.0 x 10(5) s(-1)) and support a folding mechanism in which the rate-limiting step corresponds to the formation of the reverse turn. On the other hand, the hydrophobic cluster and the disulfide cross-linker appear to largely stabilize the native state but not the folding transition state.     

Stability and folding of the SH3 domain of bruton's tyrosine kinase

Bruton's tyrosine kinase (BTK) plays an important role in B cell development. Deletion of C-terminal 14 amino acids of the SH3 domain of BTK results in X-linked agammaglobulinemia (XLA), an inherited disease. We report here on the stability and folding of SH3 domain of BTK. Peptides corresponding to residues 216–273 (58 residues) and 216–259 (44 residues) of BTK SH3 domain were synthesized by solid phase methods; the first peptide constitutes the entire SH3 domain of BTK while the latter peptide lacks 14 amino acid residues of the C-terminal. The 58 amino acid peptide forms mainly a β-barrel type folding unit. Although small and lacking disulfide bonds, this peptide is extremely stable to thermal denaturation. Based on circular dichroism measurements, its melting temperature was found to be high, 82°C at pH 6.0. However, the Gibbs free energy (ΔG) of the intrinsic stability and thermodynamic spontaneity of unfolding were found to be low, 2.6 kcal/mol by Gdn·HCl denaturation experiments, as compared to 12 kcal/mol obtained for larger single domain proteins, indicating poor stability of SH3 domain. Addition of 500 mM of Na₂SO₄ increased the free energy change ΔG to 4.0 kcal/mol, suggesting an ionic strength effect. The truncated peptide fails to fold correctly and adopts random coil conformation in contrast to 58 amino acid β-barrel peptide, which exhibits high thermal stability but normal or low stability at ambient temperature. These results, to our knowledge the first to delineate the importance of C-terminal in structural integrity of SH3 domains, indicate also that improper folding and/or poor stability of mutant SH3 domain in BTK likely causes XLA. Proteins 28:465–471 © 1996 Wiley-Liss, Inc.     

Hierarchy of structure loss in MD simulations of src SH3 domain unfolding.

To complement experimental studies of the src SH3 domain folding, we studied 30 independent, high-temperature, molecular dynamics simulations of src SH3 domain unfolding. These trajectories were observed to differ widely from each other. Thus, rather than analyzing individual trajectories, we sought to identify the recurrent features of the high-temperature unfolding process. The conformations from all simulations were combined and then divided into groups based on the number of native contacts. Average occupancies of each side-chain hydrophobic contact and hydrogen bond in the protein were then determined. In the symmetric funnel limit, the occupancies of all contacts should decrease in concert with the loss in total number of native contacts. If there is a lack of symmetry or hierarchy to the unfolding process, the occupancies of some contacts should decrease more slowly, and others more rapidly. Despite the heterogeneity of the individual trajectories, the ensemble averaging revealed an order to the unfolding process: contacts between the N and C-terminal strands are the first to disappear, whereas contacts within the distal beta-hairpin and a hydrogen-bonding network involving the distal loop beta-turn and the diverging turn persist well after the majority of the native contacts are lost. This hierarchy of events resembles but is somewhat less pronounced than that observed in our experimental studies of the folding of src SH3 domain.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Thermal unfolding of small proteins with SH3 domain folding pattern

The thermal unfolding of three SH3 domains of the Tec family of tyrosine kinases was studied by differential scanning calorimetry and CD spectroscopy. The unfolding transition of the three protein domains in the acidic pH region can be described as a reversible two-state process. For all three SH3 domains maximum stability was observed in the pH region 4.5 < pH < 7.0 where these domains unfold at temperatures of 353K (Btk), 342K (Itk), and 344K (Tec). At these temperatures an enthalpy change of 196 kJ/mol, 178 kJ/mol, and 169 kJ/mol was measured for Btk-, Itk-, and Tec-SH3 domains, respectively. The determined changes in heat capacity between the native and the denatured state are in an usual range expected for small proteins. Our analysis revealed that all SH3 domains studied are only weakly stabilized and have free energies of unfolding which do not exceed 12–16 kJ/mol but show quite high melting temperatures. Comparing unfolding free energies measured for eukaryotic SH3 domains with those of the topologically identical Sso7d protein from the hyperthermophile Sulfolobus solfataricus, the increased melting temperature of the thermostable protein is due to a broadening as well as a significant lifting of its stability curve. However, at their physiological temperatures, 310K for mesophilic SH3 domains and 350K for Sso7d, eukaryotic SH3 domains and Sso7d show very similar stabilities. Proteins 31:309–319, 1998. © 1998 Wiley-Liss, Inc.     

Solution structure and folding characteristics of the C-terminal SH3 domain of c-Crk-II

Crk-II is a signaling adaptor protein that is involved in many cellular processes including apoptosis, proliferation, and differentiation. It has a modular domain architecture consisting of an Src homology 2 domain (SH2) followed by two Src homology 3 (SH3) domains. The structures and ligand-binding properties of the SH2 and the middle SH3 domains are well-characterized. Several studies suggest that the C-terminal SH3 domain plays an important regulatory role in the protein; however, no structural information is available on this domain, and relatively little is known about its binding partners. In the current work, we have solved the solution NMR structure of the C-terminal SH3 domain. The domain adopts the standard SH3 fold comprising a five-stranded beta barrel. In agreement with alignment and modeling studies, the structure indicates that the canonical-binding surface of the SH3 domain is unusually polar and suggests that this domain may not bind typical PXXP ligands or that it may bind them with reduced affinity. Thermodynamic and kinetic studies show that the domain folds in a reversible two-state manner and that the stability of the fold is similar to that observed for other SH3 domains. These studies offer some insight into the likely structural and thermodynamic consequences of point mutations in the cSH3 domain that are known to deregulate Crk-II function. Our results set the stage for a better understanding the role of the cSH3 domain in the context of the full-length protein.     

Topography for independent binding of alpha-helical and PPII-helical ligands to a peroxisomal SH3 domain

While the function of most small signaling domains is confined to binary ligand interactions, the peroxisomal Pex13p SH3 domain has the unique capacity of binding to two different ligands, Pex5p and Pex14p. We have used this domain as a model to decipher its structurally independent ligand binding sites. By the combined use of X-ray crystallography, NMR spectroscopy, and circular dichroism, we show that the two ligands bind in unrelated conformations to patches located at opposite surfaces of this SH3 domain. Mutations in the Pex13p SH3 domain that abolish interactions within the Pex13p-Pex5p interface specifically impair PTS1-dependent protein import into yeast peroxisomes.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

Experiment and theory highlight role of native state topology in SH3 folding.

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.     

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

To investigate the relationships between protein topology, amino acid sequence and folding mechanisms, the folding transition state of the Sso7d protein has been characterised both experimentally and theoretically. Although Sso7d protein has a similar topology to that of the SH3 domains, the structure of its transition state is different from that of alpha-spectrin and src SH3 domains previously studied. The folding algorithm, Fold-X, including an energy function with specific sequence features, accounts for these differences and reproduces with a good agreement the set of experimental phi(double dagger-U) values obtained for the three proteins. Our analysis shows that taking into account sequence features underlying protein topology is critical for an accurate prediction of the folding process.     

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.