Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include ¹H-¹⁵N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O₂-induced ¹³C paramagnetic shifts, hydrogen-exchange protection factors, and ¹⁵N R₂ data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu²⁺-Ni²⁺ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.     

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.     

Ensemble modeling of protein disordered states: Experimental restraint contributions and validation

Disordered states of proteins include the biologically functional intrinsically disordered proteins and the unfolded states of normally folded proteins. In recent years, ensemble‐modeling strategies using various experimental measurements as restraints have emerged as powerful means for structurally characterizing disordered states. However, these methods are still in their infancy compared with the structural determination of folded proteins. Here, we have addressed several issues important to ensemble modeling using our ENSEMBLE methodology. First, we assessed how calculating ensembles containing different numbers of conformers affects their structural properties. We find that larger ensembles have very similar properties to smaller ensembles fit to the same experimental restraints, thus allowing a considerable speed improvement in our calculations. In addition, we analyzed the contributions of different experimental restraints to the structural properties of calculated ensembles, enabling us to make recommendations about the experimental measurements that should be made for optimal ensemble modeling. The effects of different restraints, most significantly from chemical shifts, paramagnetic relaxation enhancements and small‐angle X‐ray scattering, but also from other data, underscore the importance of utilizing multiple sources of experimental data. Finally, we validate our ENSEMBLE methodology using both cross‐validation and synthetic experimental restraints calculated from simulated ensembles. Our results suggest that secondary structure and molecular size distribution can generally be modeled very accurately, whereas the accuracy of calculated tertiary structure is dependent on the number of distance restraints used. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Distribution of molecular size within an unfolded state ensemble using small-angle X-ray scattering and pulse field gradient NMR techniques

The size distribution of molecules within an unfolded state of the N-terminal SH3 domain of drk (drkN SH3) has been studied by small-angle X-ray scattering (SAXS) and pulsed-field-gradient NMR (PFG-NMR) methods. An empirical model to describe this distribution in the unfolded state ensemble has been proposed based on (i) the ensemble-averaged radius of gyration and hydrodynamic radius derived from the SAXS and PFG-NMR data, respectively, and (ii) a histogram of the size distribution of structures obtained from preliminary analyses of structural parameters recorded on the unfolded state. Results show that this unfolded state, U(exch), which exists in equilibrium with the folded state, F(exch), under non-denaturing conditions, is relatively compact, with the average size of conformers within the unfolded state ensemble only 30-40% larger than the folded state structure. In addition, the model predicts a significant overlap in the size range of structures comprising the U(exch) state with those in a denatured state obtained by addition of 2 M guanidinium chloride.     

NOE data demonstrating a compact unfolded state for an SH3 domain under non-denaturing conditions.

The N-terminal SH3 domain of drk (drkN SH3) is unstable, existing in equilibrium between a folded state (Fexch) and an unfolded state (Uexch) under non-denaturing buffer conditions. Using a15N/2H-labeled sample, long range amide NOEs can be observed in the Uexchstate as a result of reduced relaxation, in some cases correlating protons over 40 residues apart. These long range NOEs disappear upon addition of 2 M guanidinium chloride, demonstrating that there are substantial differences between the Uexchand the guanidine denatured states. Calculations using the long range NOEs of the Uexchstate yield highly compact structures having non-native turns and a non-native buried tryptophan residue. These structures agree with experimental stopped-flow fluorescence data and analytical ultracentrifugation results. Since protein stability depends on the structural and dynamic properties of both the folded and unfolded states, this study provides insights into the stability of the drkN SH3 domain. These results provide the first strong NOE-based evidence for compact unfolded states of proteins and suggest that some unfolded states under physiological conditions have specific interactions leading to compact structures.     

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.     

Aromatic and methyl NOEs highlight hydrophobic clustering in the unfolded state of an SH3 domain

The N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) exists in equilibrium between a folded (F(exch)) state and a relatively compact unfolded (U(exch)) state under nondenaturing conditions. Selectively labeled samples of the domain have been analyzed by NOESY NMR experiments to probe residual hydrophobic clustering in the U(exch) state. The labeling strategy included selective protonation of aromatic rings or delta-methyl groups on Ile and Leu residues in a highly deuterated background. Combined with long mixing times, the methods permitted observation of significant numbers of long-range interactions between hydrophobic side chains, providing evidence for multiple conformers involving non-native hydrophobic clusters around the Trp 36 indole. Comparison of these data with previously reported HN-HN NOEs yields structural insight into the diversity of structures within the U(exch) ensemble in the drkN SH3 domain. Many of the HN-HN NOEs are consistent with models containing compact residual nativelike secondary structure and greater exposure of the Trp 36 indole to solvent, similar to kinetic intermediates formed in the hierarchic condensation model of folding. However, the methyl and aromatic NOE data better fit conformations with non-native burial of the Trp indole surrounded by hydrophobic groups and more loosely formed beta-structure; these structural characteristics are more consistent with those of kinetic intermediates formed during the hydrophobic collapse mechanism of folding. This suite of NOE data provides a more complete picture of the structures that span the U(exch) state ensemble, from conformers with non-native structure but long-range contacts to those that are highly nativelike. Together, the results are also consistent with the folding funnel view involving multiple folding pathways for this molecule.     

Formation of a copper specific binding site in non-native states of beta-2-microglobulin

A debilitating complication of long-term hemodialysis is the deposition of beta-2-microglobulin (beta2m) as amyloid plaques in the joint space. We have recently shown that Cu(2+) can be a contributing, if not causal, factor at concentrations encountered during dialysis therapy. The basis for this effect is destabilization and incorporation of beta2m into amyloid fibers upon binding of Cu(2+). In this work, we demonstrate that while beta2m binds Cu(2+) specifically in the native state, it is binding of Cu(2+) by non-native states of beta2m which is responsible for destabilization. Mutagenesis of potential coordinating groups for Cu(2+) shows that native state binding of Cu(2+) is mediated by residues and structures that are different than those which bind in non-native states. An increased affinity for copper by non-native states compared to that of the native state gives rise to overall destabilization. Using mass spectrometry, NMR, and fluorescence techniques, we show that native state binding is localized to H31 and W60 and is highly specific for Cu(2+) over Zn(2+) and Ni(2+). Binding of Cu(2+) in non-native states of beta2m is mediated by residues H13, H51, and H84, but not H31. Although denatured beta2m has characteristics of a globally unfolded state, it nevertheless demonstrates the following strong specificity of binding: Cu(2+) > Zn(2+) > Ni(2+). This requires the existence of a well-defined structure in the unfolded state of this protein. As Cu(2+) effects are reported in many other amyloidoses, e.g., PrP, alpha-synuclein, and Abeta, our results may be extended to the emerging field of divalent ion-associated amyloidosis.     

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.     

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by F?rster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.     

Backbone 1H and 15N resonance assignments of the N-terminal SH3 domain of drk in folded and unfolded states using enhanced-sensitivity pulsed field gradient NMR techniques

Summary The backbone ¹H and ¹⁵N resonances of the N-terminal SH3 domain of the Drosophila signaling adapter protein, drk, have been assigned. This domain is in slow exchange on the NMR timescale between folded and predominantly unfolded states. Data were collected on both states simultaneously, on samples of the SH3 in near physiological buffer exhibiting an approximately 1:1 ratio of the two states. NMR methods which exploit the chemical shift dispersion of the ¹⁵N resonances of unfolded states and pulsed field gradient water suppression approaches for avoiding saturation and dephasing of amide protons which rapidly exchange with solvent were utilized for the assignment.Abbreviations 2D, 3D two-, three-dimensional - drkN SH3 N-terminal SH3 domain of Drosophila drk - HSQC heteronuclear single-quantum spectroscopy - NOE nuclear Overhauser enhancement - SH3 Src homology domain 3 - TOCSY total correlation spectroscopy     

Structural comparison of the unstable drkN SH3 domain and a stable mutant

The N-terminal SH3 domain of the Drosophila adapter protein Drk (drkN SH3 domain) is marginally stable (DeltaG(U) = 1 kcal/mol) and exists in equilibrium between folded and highly populated unfolded states. The single substitution T22G, however, completely stabilizes the protein (DeltaG(U) = 4.0 kcal/mol). To probe the causes of instability of the wild-type (WT) protein and the dramatic stabilization of the mutant, we determined and compared nuclear magnetic resonance structures of the folded WT and mutant drkN SH3 domains. Residual dipolar coupling (RDC) and carbonyl chemical-shift anisotropy (C'-CSA) restraints measured for the WT and T22G domains were used for calculating the structures. The structures for the WT and mutant are highly similar. Thr22 of the WT and Gly22 of the mutant are at the i + 2 position of the diverging, type-II beta-turn. Interestingly, not only Gly22 but also Thr22 successfully adopt an alpha(L) conformation, required at this position of the turn, despite the fact that positive phi values are energetically unfavorable and normally disallowed for threonine residues. Forcing the Thr22 residue into this unnatural conformation increases the free energy of the folded state of the WT domain relative to its T22G mutant. Evidence for residual helix formation in the diverging turn region has been previously reported for the unfolded state of the WT drkN SH3 domain, and this, in addition to other residual structure, has been proposed to play a role in decreasing the free energy of the unfolded state of the protein. Together these data provide evidence that both increasing the free energy of the folded state and decreasing the free energy of the unfolded state of the protein contribute to instability of the WT drkN SH3 domain.     

Exploring the Minimally Frustrated Energy Landscape of Unfolded ACBP

The unfolded state of globular proteins is not well described by a simple statistical coil due to residual structural features, such as secondary structure or transiently formed long-range contacts. The principle of minimal frustration predicts that the unfolded ensemble is biased toward productive regions in the conformational space determined by the native structure. Transient long-range contacts, both native-like and non-native-like, have previously been shown to be present in the unfolded state of the four-helix-bundle protein acyl co-enzyme binding protein (ACBP) as seen from both perturbations in nuclear magnetic resonance (NMR) chemical shifts and structural ensembles generated from NMR paramagnetic relaxation data. To study the nature of the contacts in detail, we used paramagnetic NMR relaxation enhancements, in combination with single-point mutations, to obtain distance constraints for the acid-unfolded ensemble of ACBP. We show that, even in the acid-unfolded state, long-range contacts are specific in nature and single-point mutations affect the free-energy landscape of the unfolded protein. Using this approach, we were able to map out concerted, interconnected, and productive long-range contacts. The correlation between the native-state stability and compactness of the denatured state provides further evidence for native-like contact formation in the denatured state. Overall, these results imply that, even in the earliest stages of folding, ACBP dynamics are governed by native-like contacts on a minimally frustrated energy landscape.     

Tryptophan solvent exposure in folded and unfolded states of an SH3 domain by 19F and 1H NMR

The isolated N-terminal SH3 domain of the Drosophila signal transduction protein Drk (drkN SH3) is a useful model for the study of residual structure and fluctuating structure in disordered proteins since it exists in slow exchange between a folded (Fexch) and compact unfolded (Uexch) state in roughly equal proportions under nondenaturing conditions. The single tryptophan residue, Trp36, is believed to play a key role in forming a non-native hydrophobic cluster in the Uexch state, with a number of long-range nuclear Overhauser contacts (NOEs) observed primarily to the indole proton. Substitution of Trp36 for 5-fluoro-Trp36 resulted in a substantial shift in the equilibrium to favor the Fexch state. A variety of 19F NMR measurements were performed to investigate the degree of solvent exposure and hydrophobicity associated with the 5-fluoro position in both the Fexch and Uexch states. Ambient T1 measurements and H2O/D2O solvent isotope effects indicated extensive protein contacts to the 5-fluoro position in the Fexch state and greater solvent exposure in the Uexch state. This was corroborated by the measurements of paramagnetic effects (chemical shift perturbations and T1 relaxation enhancement) from dissolved oxygen at a partial pressure of 20 atm. In contrast, paramagnetic effects from dissolved oxygen revealed less solvent exposure to the indole proton of Trp36 in the Uexch state than that observed for the Fexch state, consistent with the model in which Trp36 indole belongs to a non-native cluster. Thus, although the Uexch state may be described as a dynamically interconverting ensemble of conformers, there appears to be significant asymmetry in the environment of the indole group and the six-membered ring or backbone of Trp36. This implied lack of averaging of a side chain position is in contrast to the general view of fluctuating side chains within disordered states.     

Stiffness of the distal loop restricts the structural heterogeneity of the transition state ensemble in SH3 domains

Protein engineering experiments and Phi(F)-value analysis of SH3 domains reveal that their transition state ensemble (TSE) is conformationally restricted, i.e. the fluctuations in the transition state (TS) structures are small. In the TS of src SH3 and alpha-spectrin SH3 the distal loop and the associated hairpin are fully structured, while the rest of the protein is relatively disordered. If native structure predominantly determines the folding mechanism, the findings for SH3 folds raise the question: What are the features of the native topology that determine the nature of the TSE? We propose that the presence of stiff loops in the native state that connect local structural elements (such as the distal hairpin in SH3 domains) conformationally restricts TSE. We validate this hypothesis using the simulations of a "control" system (16 residue beta-hairpin forming C-terminal fragment of the GBl protein) and its variants. In these fragments the role of bending rigidity in determining the nature of the TSE can be directly examined without complications arising from interactions with the rest of the protein. The TSE structures in the beta-hairpins are determined computationally using cluster analysis and limited Phi(F)-value analysis. Both techniques prove that the conformational heterogeneity decreases as the bending rigidity of the loop increases. To extend this finding to SH3 domains a measure of bending rigidity based on loop curvature, which utilizes native structures in the Protein Data Bank (PDB), is introduced. Using this measure we show that, with few exceptions, the ordering of stiffness of the distal, n-src, and RT loops in the 29 PDB structures of SH3 domains is conserved. Combining the simulation results for beta-hairpins and the analysis of PDB structures for SH3 domains, we propose that the stiff distal loop restricts the conformational fluctuations in the TSE. We also predict that constraining the distal loop to be preformed in the denatured ensemble should not alter the nature of TSE. On the other hand, if the amino and carboxy terminals are cross-linked to form a circular polypeptide chain, the pathways and TSs are altered. These contrasting scenarios are illustrated using simulations of cross-linked WT beta-hairpin fragments. Computations of bending rigidities for immunoglobulin-like domain proteins reveal no clear separation in the stiffness of their loops. In the beta-sandwich proteins, which have large fractions of non-local native contacts, the nature of the TSE cannot be apparently determined using purely local structural characteristics. Nevertheless, the measure of loop stiffness still provides qualitative predictions of the ordered regions in the TSE of Ig27 and TenFn3.     

The effect of increasing the stability of non-native interactions on the folding landscape of the bacterial immunity protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.     

Native-like mean structure in the unfolded ensemble of small proteins

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.     

Slo1 tail domains,but not the Ca2+ bowl,are required for the beta 1 subunit to increase the apparent Ca2+ sensitivity of BK channels

Functional large-conductance Ca(2+)- and voltage-activated K(+) (BK) channels can be assembled from four alpha subunits (Slo1) alone, or together with four auxiliary beta1 subunits to greatly increase the apparent Ca(2+) sensitivity of the channel. We examined the structural features involved in this modulation with two types of experiments. In the first, the tail domain of the alpha subunit, which includes the RCK2 (regulator of K(+) conductance) domain and Ca(2+) bowl, was replaced with the tail domain of Slo3, a BK-related channel that lacks both a Ca(2+) bowl and high affinity Ca(2+) sensitivity. In the second, the Ca(2+) bowl was disrupted by mutations that greatly reduce the apparent Ca(2+) sensitivity. We found that the beta1 subunit increased the apparent Ca(2+) sensitivity of Slo1 channels, independently of whether the alpha subunits were expressed as separate cores (S0-S8) and tails (S9-S10) or full length, and this increase was still observed after the Ca(2+) bowl was mutated. In contrast, beta1 subunits no longer increased Ca(2+) sensitivity when Slo1 tails were replaced by Slo3 tails. The beta1 subunits were still functionally coupled to channels with Slo3 tails, as DHS-I and 17 beta-estradiol activated these channels in the presence of beta1 subunits, but not in their absence. These findings indicate that the increase in apparent Ca(2+) sensitivity induced by the beta1 subunit does not require either the Ca(2+) bowl or the linker between the RCK1 and RCK2 domains, and that Slo3 tails cannot substitute for Slo1 tails. The beta1 subunit also induced a decrease in voltage sensitivity that occurred with either Slo1 or Slo3 tails. In contrast, the beta1 subunit-induced increase in apparent Ca(2+) sensitivity required Slo1 tails. This suggests that the allosteric activation pathways for these two types of actions of the beta1 subunit may be different.     

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.