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.     

Sequence optimization for native state stability determines the evolution and folding kinetics of a small protein

Investigating the relative importance of protein stability, function, and folding kinetics in driving protein evolution has long been hindered by the fact that we can only compare modern natural proteins, the products of the very process we seek to understand, to each other, with no external references or baselines. Through a large-scale all-atom simulation of protein evolution, we have created a large diverse alignment of SH3 domain sequences which have been selected only for native state stability, with no other influencing factors. Although the average pairwise identity between computationally evolved and natural sequences is only 17%, the residue frequency distributions of the computationally evolved sequences are similar to natural SH3 sequences at 86% of the positions in the domain, suggesting that optimization for the native state structure has dominated the evolution of natural SH3 domains. Additionally, the positions which play a consistent role in the transition state of three well-characterized SH3 domains (by phi-value analysis) are structurally optimized for the native state, and vice versa. Indeed, we see a specific and significant correlation between sequence optimization for native state stability and conservation of transition state structure.     

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.     

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.     

Study of the structure and dynamics of a chimeric variant of the SH3 domain (SHA-Bergerac) by NMR spectroscopy

A structural-dynamic study of one of the chimeric proteins (SHA) belonging to the SH3-Bergerac family and containing the KATANGKTYE sequence instead of the N47D48 β-turn in the spectrin SH3-domain was carried out by high resolution NMR spectroscopy. The spatial structure of the protein was determined and its dynamics in solution was investigated on the basis of the NMR data. The elongation of the SHA polypeptide chain in comparison with the WT-SH3 original protein (by ～17%) exerts practically no effect on the general topology of the molecule. The presence of a stable β-hairpin in the region of insertion was confirmed. This hairpin was shown to have a higher mobility in comparison with other regions of the protein.     

Cooperative interactions and a non-native buried Trp in the unfolded state of an SH3 domain

The presence of residual structure in the unfolded state of the N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) has been investigated using far- and near-UV circular dichroism (CD), fluorescence, and NMR spectroscopy. The unfolded (U(exch)) state of the drkN SH3 domain is significantly populated and exists in equilibrium with the folded (F(exch)) state under non-denaturing conditions near physiological pH. Denaturation experiments have been performed on the drkN SH3 domain in order to monitor the change in ellipticity, fluorescence intensity, and chemical shift between the U(exch) state and chemically or thermally denatured states. Differences between the unfolded and chemically or thermally denatured states highlight specific areas of residual structure in the unfolded state that are cooperatively disrupted upon denaturation. Results provide evidence for cooperative interactions in the unfolded state involving residues of the central beta-sheet, particularly the beta4 strand. Denaturation as well as hydrogen-exchange experiments demonstrate a non-native burial of the Trp ring within this "cooperative" core of the unfolded state. These findings support the presence of non-native hydrophobic clusters, organised by Trp rings, within disordered states.     

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.     

Src homology domains of v-Src stabilize an active conformation of the tyrosine kinase catalytic domain

To examine the interactions between Src homology,domains and the tyrosine kinase catalytic domain of v-Src, various combinations of domains have been expressed in bacteria as fusion proteins. Constructs containing the isolated catalytic domain, SH2 + catalytic domain, and SH3 + SH2 + catalytic domains were active in autophosphorylation assays. For the catalytic domain of v-Src, but not for v-Abl, addition of exogenous Src SH3-SH2 domains stimulated the autophosphorylation activity. In contrast to results for autophosphorylation, constructs containing Src homology domains were more active towards a synthetic peptide substrate than the isolated catalytic domain. The ability of the SH2 and SH3 domains of v-Src to stabilize an active enzyme conformation was also confirmed by refolding after denaturation in guanidinium hydrochloride. Collectively the data suggest that, in addition to their roles in intermolecular protein-protein interactions, the Src homology regions of v-Src exert a positive influence on tyrosine kinase function, potentially by maintaining an active conformation of the catalytic domain.     

CSSI-PRO: a method for secondary structure type editing,assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone ¹H^α and ¹³C′ chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to α-helical/β-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Local secondary structure content predicts folding rates for simple,two-state proteins

Many single-domain proteins exhibit two-state folding kinetics, with folding rates that span more than six orders of magnitude. A quantity of much recent interest for such proteins is their contact order, the average separation in sequence between contacting residue pairs. Numerous studies have reached the surprising conclusion that contact order is well-correlated with the logarithm of the folding rate for these small, well-characterized molecules. Here, we investigate the physico-chemical basis for this finding by asking whether contact order is actually a composite number that measures the fraction of local secondary structure in the protein; viz. turns, helices, and hairpins. To pursue this question, we calculated the secondary structure content for 24 two-state proteins and obtained coefficients that predict their folding rates. The predicted rates correlate strongly with experimentally determined rates, comparable to the correlation with contact order. Further, these predicted folding rates are correlated strongly with contact order. Our results suggest that the folding rate of two-state proteins is a function of their local secondary structure content, consistent with the hierarchic model of protein folding. Accordingly, it should be possible to utilize secondary structure prediction methods to predict folding rates from sequence alone.     

The high resolution NMR structure of the third SH3 domain of CD2AP

CD2 associated protein (CD2AP) is an adaptor protein that plays an important role in cell to cell union needed for the kidney function. CD2AP interacts, as an adaptor protein, with different natural targets, such as CD2, nefrin, c-Cbl and podocin. These proteins are believed to interact to one of the three SH3 domains that are positioned in the N-terminal region of CD2AP. To understand the network of interactions between the natural targets and the three SH3 domains (SH3-A, B and C), we have started to determine the structures of the individual SH3 domains. Here we present the high-resolution structure of the SH3-C domain derived from NMR data. Full backbone and side-chain assignments were obtained from triple-resonance spectra. The structure was determined from distance restraints derived from high-resolution 600 and 800 MHz NOESY spectra, together with phi and psi torsion angle restraints based on the analysis of ¹HN, ¹⁵N, ¹Hα, ¹³Cα, ¹³CO and ¹³Cβ chemical shifts. Structures were calculated using CYANA and refined in water using RECOORD. The three-dimensional structure of CD2AP SH3-C contains all the features that are typically found in other SH3 domains, including the general binding site for the recognition of polyproline sequences.     

A lack of peptide binding and decreased thermostability suggests that the CASKIN2 scaffolding protein SH3 domain may be vestigial

Background

CASKIN2 is a neuronal signaling scaffolding protein comprised of multiple ankyrin repeats, two SAM domains, and one SH3 domain. The CASKIN2 SH3 domain for an NMR structural determination because its peptide-binding cleft appeared to deviate from the repertoire of aromatic enriched amino acids that typically bind polyproline-rich sequences.

Results

The structure demonstrated that two non-canonical basic amino acids (K290/R319) in the binding cleft were accommodated well in the SH3 fold. An K290Y/R319W double mutant restoring the typical aromatic amino acids found in the binding cleft resulted in a 20 °C relative increase in the thermal stability. Considering the reduced stability, we speculated that the CASKIN2 SH3 could be a nonfunctional remnant in this scaffolding protein.

Conclusions

While the NMR structure demonstrates that the CASKIN2 SH3 domain is folded, its cleft has suffered two substitutions that prevent it from binding typical polyproline ligands. This observation led us to additionally survey and describe other SH3 domains in the Protein Data Bank that may have similarly lost their ability to promote protein-protein interactions.

     

ALASKA: A Mathematica package for two-state kinetic analysis of protein folding reactions

A Mathematica package (ALASKA) has been developed to simplify the measurement of protein folding kinetics by analysis of ¹H NMR lineshape analysis. This package reads NMR data in ASCII format and can simulate an aromatic ¹ NMR spectrum with or without lineshape broadening from chemical exchange. We describe the analysis of a urea denaturation series of a fast-folding protein, the G46A/G48A variant of monomeric repressor.     

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.