Search

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.

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.

Eph receptor tyrosine kinases (RTKs) mediate numerous developmental processes. Their activity is regulated by auto-phosphorylation on two tyrosines within the juxtamembrane segment (JMS) immediately N-terminal to the kinase domain (KD). Here, we probe the molecular details of Eph kinase activation through mutational analysis, X-ray crystallography and NMR spectroscopy on auto-inhibited and active EphB2 and EphA4 fragments. We show that a Tyr750Ala gain-of-function mutation in the KD and JMS phosphorylation independently induce disorder of the JMS and its dissociation from the KD. Our X-ray analyses demonstrate that this occurs without major conformational changes to the KD and with only partial ordering of the KD activation segment. However, conformational exchange for helix alphaC in the N-terminal KD lobe and for the activation segment, coupled with increased inter-lobe dynamics, is observed upon kinase activation in our NMR analyses. Overall, our results suggest that a change in inter-lobe dynamics and the sampling of catalytically competent conformations for helix alphaC and the activation segment rather than a transition to a static active conformation underlies Eph RTK activation.

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.

SH3 is a common binding domain that mediates protein-protein interactions. The N-terminal SH3 domain of the Drosophila adaptor protein Drk is marginally stable under physiological conditions, showing nearly 50/50 equilibrium between the folded and unfolded states. Due to the dynamic ensemble nature of this unfolded state under non-denaturing conditions, the structural characterization of SH3 is difficult using ensemble measurements. Experimental data have been used to define sets of heterogeneous conformations but ensemble averaging minimizes the information content [1]. Importantly, little is known about interconversion rates within the unfolded ensemble.We performed fluorescence correlation spectroscopy (FCS) experiments on the wild-type and on a single-site mutant stabilizing the folded state. Based on this data, we estimated the size (hydrodynamic radii), of the folded (Rf), unfolded (Ru) and denatured (Rd) states. The size of unfolded state measured by FCS is in excellent agreement with the NMR-measured value, Ru/Rf = 1.30 ± 0.01 [1]. Quite surprisingly, our results indicate that the unfolded state under non-denaturing conditions is less compact than the chemically-denatured state, as Rd/Rf = 1.21 ± 0.03. This agrees with previously reported NMR and SAXS data when adequate shape factors for different conformations are considered. Using an environment sensitive dye and fitting the full correlation curve reveals a fast (∼200 ns) process within the unfolded ensemble that is absent in the folded state.Single-molecule FRET experiments were performed to quantify end-to-end distance distributions. We probed directly the conformational heterogeneity of the wild-type protein as well as the response to denaturants and osmolytes. Quite surprisingly, N- and the C-termini remained in close proximity at high denaturant concentrations for ca. 40% of the conformations, suggesting that DrkN-SH3behaves at least partially like a disordered circular chain.Reference:1. J.A. Marsh and J.D. Forman-Kay, Proteins 80(2):556-572 (2012)

The Smad2 Mad homology 2 (MH2) domain binds to a diverse group of proteins which do not share a common sequence motif. We have used NMR to investigate the structure of one of these interacting proteins, the Smad binding domain (SBD) of Smad anchor for receptor activation (SARA). Our results indicate that the unbound SBD is highly disordered and forms no stable secondary or tertiary structures. Additionally we have used fluorescence binding studies to study the interaction between the MH2 domain and SBD and find that no region of the SBD dominates the interaction between the MH2 and the SBD. Our results are consistent with a series of hydrophobic patches on the MH2 that are able to recognize disordered regions of proteins. These findings elucidate a mechanism by which a single domain (MH2) can specifically recognize a diverse set of proteins which are unrelated by sequence, lead to a clearer picture of how MH2 domains function in the transforming growth factor-beta-signaling pathway and suggest possible mechanisms for controlling interactions with MH2 domains.

Understanding protein stability is a significant challenge requiring characterization of interactions within both folded and unfolded states. Of these, electrostatic interactions influence ionization equilibria of acidic and basic groups and diversify their pK a values. The pH dependence of the thermodynamic stability (ΔG FU ) of a protein arises as a consequence of differential pK a values between folded and unfolded states. Previous attempts to calculate pH-dependent contributions to stability have been limited by the lack of experimental unfolded state pK a values. Using recently developed NMR spectroscopic methods, we have determined residue-specific pK a values for a thermodynamically unstable Src homology 3 domain in both states, enabling the calculation of the pH dependence of stability based on simple analytical expressions. The calculated pH stability profile obtained agrees very well with experiment, unlike profiles derived from two current models of electrostatic interactions within unfolded states. Most importantly, per-residue contributions to the pH dependence of ΔG FU derived from the data provide insights into specific electrostatic interactions in both the folded and unfolded states and their roles in protein stability. These interactions include a hydrogen bond between the Asp-8 side-chain and the Lys-21 backbone amide group in the folded state, which represents a highly conserved interaction in Src homology 3 domains.

The epithelial Na+ channel (alphabetagammaENaC) regulates salt and fluid homeostasis and blood pressure. Each ENaC subunit contains a PY motif (PPXY) that binds to the WW domains of Nedd4, a Hect family ubiquitin ligase containing 3-4 WW domains and usually a C2 domain. It has been proposed that Nedd4-2, but not Nedd4-1, isoforms can bind to and suppress ENaC activity. Here we challenge this notion and show that, instead, the presence of a unique WW domain (WW3*) in either Nedd4-2 or Nedd4-1 determines high affinity interactions and the ability to suppress ENaC. WW3* from either Nedd4-2 or Nedd4-1 binds ENaC-PY motifs equally well (e.g. Kd approximately 10 microm for alpha- or betaENaC, 3-6-fold higher affinity than WW4), as determined by intrinsic tryptophan fluorescence. Moreover, dNedd4-1, which naturally contains a WW3* instead of WW2, is able to suppress ENaC function equally well as Nedd4-2. Homology models of the WW3*.betaENaC-PY complex revealed that a Pro and Ala conserved in all WW3*, but not other Nedd4-WW domains, help form the binding pocket for PY motif prolines. Extensive contacts are formed between the betaENaC-PY motif and the Pro in WW3*, and the small Ala creates a large pocket to accommodate the peptide. Indeed, mutating the conserved Pro and Ala in WW3* reduces binding affinity 2-3-fold. Additionally, we demonstrate that mutations in PY motif residues that form contacts with the WW domain based on our previously solved structure either abolish or severely reduce binding affinity to the WW domain and that the extent of binding correlates with the level of ENaC suppression. Independently, we show that a peptide encompassing the PY motif of sgk1, previously proposed to bind to Nedd4-2 and alter its ability to regulate ENaC, does not bind (or binds poorly) the WW domains of Nedd4-2. Collectively, these results suggest that high affinity of WW domain-PY-motif interactions rather than affiliation with Nedd4-1/Nedd-2 is critical for ENaC suppression by Nedd4 proteins.

Summary: ENSEMBLE is a computational approach for determining a set of conformations that represents the structural ensemble of a disordered protein based on input experimental data. The disordered protein can be an unfolded or intrinsically disordered state. Here, we introduce the latest version of the program, which has been enhanced to facilitate its general release and includes an intuitive user interface, as well as new approaches to treat data and analyse results. Availability and implementation: ENSEMBLE is a program implemented in C and embedded in a Perl wrapper. It is supported on main Linux distributions. Source codes and installation files, including a detailed example, can be freely downloaded at http://abragam.med.utoronto.ca/∼JFKlab. Contact: forman@sickkids.ca Supplementary information: Supplementary Material are available at Bioinformatics online.

X-linked lymphoproliferative disease is caused by mutations in the protein SAP, which consists almost entirely of a single SH2 domain. SAP interacts with the Tyr281 site of the T↔B cell signaling protein SLAM via its SH2 domain. Interestingly, binding is not dependent on phosphorylation but does involve interactions with residues N-terminal to the Tyr. We have used 15 N and 2 H NMR relaxation experiments to investigate the motional properties of the SAP SH2 domain backbone amides and side-chain methyl groups in the free protein and complexes with phosphorylated and non-phosphorylated peptides derived from the Tyr281 site of SLAM. The most mobile methyl groups are in side-chains with large RMSD values between the three crystal structures of SAP, suggesting that fast time-scale dynamics in side-chains is associated with conformational plasticity. The backbone amides of two residues which interact with the C-terminal part of the peptides experience fast time-scale motions in the free SH2 domain that are quenched upon binding of either the phosphorylated or non-phosphorylated peptide. Of most importance, the mobility of methyl groups in and around the binding site for residues in the N-terminus of the peptide is significantly restricted in the complexes, underscoring the dominance of this interaction with SAP and demonstrating a correlation between changes in rapid side-chain motion upon binding with local binding energy.

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 β-sheet, particularly the β4 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.

Structural studies of proteins are critical for understanding biological processes at the molecular level. Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for obtaining structural and dynamic information on proteins and protein-ligand complexes. In the present review, methodologies for NMR structure determination of proteins and macromolecular complexes are described. In addition, a number of recent advances that reduce the molecular weight limitations previously imposed on NMR studies of biomolecules are discussed, highlighting applications of these technologies to protein systems studied in our laboratories.

The SH2 domain protein SAP/SH2D1A, encoded by the X-linked lymphoproliferative (XLP) syndrome gene, associates with the hematopoietic cell surface receptor SLAM in a phosphorylation-independent manner. By screening a repertoire of synthetic peptides, the specificity of SAP/SH2D1A has been mapped and a consensus sequence motif for binding identified, T/S-x-x-x-x-V/I, where x represents any amino acid. Remarkably, this motif contains neither a Tyr nor a pTyr residue, a hallmark of conventional SH2 domain–ligand interactions. The structures of the protein, determined by NMR, in complex with two distinct peptides provide direct evidence in support of a ‘three-pronged’ binding mechanism for the SAP/SH2D1A SH2 domain in contrast to the ‘two-pronged’ binding for conventional SH2 domains. Differences in the structures of the two complexes suggest considerable flexibility in the SH2 domain, as further confirmed and characterized by hydrogen exchange studies. The structures also explain binding defects observed in disease-causing SAP/SH2D1A mutants and suggest that phosphorylation-independent interactions mediated by SAP/SH2D1A likely play an important role in the pathogenesis of XLP.

Intrinsically disordered proteins (IDPs) have been implicated in the regulation of many important cellular processes, such as regulation of cell cycle progression, of transcription and of translation. IDPs are thought to function primarily in the mediation of protein-protein interactions, but detailed analysis of the molecular mechanisms that render disorder beneficial in protein interactions is required. According to a common view, IDPs are usually thought to undergo disorder-to-order transitions upon interaction with their binding partners. It is becoming increasingly clear that predominantly disordered protein complexes are functionally relevant. We have shown previously that the intrinsically disordered cyclin dependent kinase (CDK) inhibitor Sic1 interacts with the substrate adapter, Cdc4, of its ubiquitin ligase via multiple phosphorylated binding motifs in a dynamic complex. Cdc4 is the substrate recognition subunit of a culin ubiquitin ligase and targets Sic1 for degradation at the G1/S phase transition of the yeast cell cycle. Individual binding motifs in Sic1 are ordered transiently without a global disorder-to-order transition upon binding Cdc4. The dynamic complex allows for engagement of several phosphorylation sites in a dynamic interface resulting in an affinity that depends on the number of phosphorylation sites in an ultrasensitive manner. The dynamic nature of the complex allows for ‘counting’ of phosphorylation sites via largely electrostatic interactions. The Sic1-Cdc4 interaction therefore acts as a sensor of the concentration of active kinase and the cell cycle status. We continue to use NMR spectroscopy and other biophysical methods to study dynamic interactions in the ubiquitin proteasome pathway with the objective of unraveling nature's repertoire of disorder in protein function. The combination of disorder with multi-site phosphorylation may serve as a general means to set thresholds in regulated protein-protein interactions.

(15)N relaxation dispersion experiments were applied to the isolated N-terminal SH3 domain of the Drosophila protein drk (drkN SH3) to study microsecond to second time scale exchange processes. The drkN SH3 domain exists in equilibrium between folded (F(exch)) and unfolded (U(exch)) states under nondenaturing conditions in a ratio of 2:1 at 20 degrees C, with an average exchange rate constant, k(ex), of 2.2 s(-1) (slow exchange on the NMR chemical shift time scale). Consequently a discrete set of resonances is observed for each state in NMR spectra. Within the U(exch) ensemble there is a contiguous stretch of residues undergoing conformational exchange on a micros/ms time scale, likely due to local, non-native hydrophobic collapse. For these residues both the F(exch)--U(exch) conformational exchange process and the micros/ms exchange event within the U(exch) state contribute to the (15)N line width and can be analyzed using CPMG-based (15)N relaxation dispersion measurements. The contribution of both processes to the apparent relaxation rate can be deconvoluted numerically by combining the experimental (15)N relaxation dispersion data with results from an (15)N longitudinal relaxation experiment that accurately quantifies exchange rates in slow exchanging systems (Farrow, N. A.; Zhang, O.; Forman-Kay, J. D.; Kay, L. E. J. Biomol. NMR 1994, 4, 727-734). A simple, generally applicable analytical expression for the dependence of the effective transverse relaxation rate constant on the pulse spacing in CPMG experiments has been derived for a two-state exchange process in the slow exchange limit, which can be used to fit the experimental data on the global folding/unfolding transition. The results illustrate that relaxation dispersion experiments provide an extremely sensitive tool to probe conformational exchange processes in unfolded states and to obtain information on the free energy landscape of such systems.

The N-terminal SH3 domain of drk (drkN SH3 domain) exists in equilibrium between a folded (F(exch)) and an unfolded (U(exch)) form under non-denaturing conditions. In order to further our previous descriptions of the U(exch) state, we have developed a protocol for calculating ensembles of structures, based on experimental spectroscopic data, which broadly represent the unfolded state. A large number of unfolding trajectories were generated, starting from the folded state structure of the protein, in order to provide a reasonable sampling of the conformational space accessible to this sequence. Unfolded state ensembles have been "calculated" using a newly developed program ENSEMBLE, which optimizes the population weights assigned to each structure based on experimental properties of the U(exch) state. Pseudo-energy terms for nuclear Overhauser effects, J-coupling constants, (13)C chemical shifts, translational diffusion coefficients and tryptophan ring burial based on NMR and fluorescence data have been implemented. The population weight assignment procedure was performed for different starting ensembles. Small numbers of structures (

Low solubility is a major stumbling block in the detailed structural and functional characterization of many proteins and isolated protein domains. The production of some proteins in a soluble form may only be possible through alteration of their sequences by mutagenesis. The feasibility of this approach has been demonstrated in a number of cases where amino acid substitutions were shown to increase protein solubility without altering structure or function. However, identifying residues to mutagenize to increase solubility is difficult, especially in the absence of structural knowledge. For this reason, we have developed a method by which soluble mutants of an insoluble protein can be easily distinguished in vivo in Escherichia coli. This method is based on our observation that cells expressing fusions of an insoluble protein to chloramphenicol acetyltransferase (CAT) exhibit decreased resistance to chloramphenicol compared to fusions with soluble proteins. We found that a soluble mutant of an insoluble protein fused to CAT could be selected by plating on high levels of chloramphenicol.

The N-terminal SH3 domain of drk (drkN SH3) is unstable, existing in equilibrium between a folded state (F exch ) and an unfolded state (U exch ) under non-denaturing buffer conditions. Using a 15 N/ 2 H-labeled sample, long range amide NOEs can be observed in the U exch state 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 U exch and the guanidine denatured states. Calculations using the long range NOEs of the U exch state 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.

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.

The phosphotyrosine-binding (PTB) domain of Numb, a protein involved in asymmetric cell division, has recently been shown to bind to the adapter protein Lnx through an LDNPAY sequence, to the Numb-associated kinase (Nak) through a sequence that does not contain an NPXY motif and to GP(p)Y-containing peptides obtained from library screening. We show here that these diverse peptide sequences bind with comparable affinities to the Numb PTB domain at a common binding site on the surface of the protein. The NMR structure of the Numb PTB domain in complex with a GPpY-containing peptide reveals a novel mechanism of binding with the peptide in a helical turn that does not hydrogen bond to the PTB domain beta-sheet. These results suggest that PTB domains can potentially have multiple modes of peptide recognition and provide a structural basis from which the multiple functions of the Numb PTB domain during asymmetric cell division could arise.

Protein recognition is a key determinant in regulating biological processes. Structures of complexes of interacting proteins provide significant insights into the mechanism of specific recognition. However, studies performed by modifying residues within a protein interface demonstrate that binding is not fully explained by these static pictures. Thus, structural data alone was not predictive of affinities in binding studies of phospholipase Cgamma1 and Syp phosphatase SH2 domains with phosphopeptides. NMR relaxation experiments probing dynamics of methyl groups of these complexes indicate a correlation between binding energy and restriction of motion at the interfacial region responsible for specific binding.

Comprehensive NOE results from detailed structural characterization of a 131 residue partially folded fragment of staphylococcal nuclease (Delta131Delta) made possible by NMR methods with improved resolution are presented. The resulting NOE patterns reflect sampling of both alpha and beta regions of phi, phi conformational space, yet demonstrate significant preferences for both native-like and non-native-like turn and potentially helical conformations. Together with data from studies of the unfolded state of the drkN SH3 domain, NOE patterns observed for partially folded or unfolded proteins are summarized. It is surprising that few long-range NOEs were observed in Delta131Delta. The two longest-range NOEs are both native-like; one of these, an (i,i+5) NOE, provides evidence for a Schellman capping motif for helix termination. Many aliphatic-aliphatic and aliphatic-amide NOEs, which are not normally observed in folded proteins, were detected. We have ruled out significant contributions from spin-diffusion for a number of these NOEs and suggest that one source may be sampling of non-prolyl cis peptide bond configurations in the disordered state of Delta131Delta.

Measurements of 15N NMR relaxation parameters have been used to characterize the backbone dynamics of folded and denatured states of the N-terminal SH3 domain from the adapter protein drk, in high salt or guanidinium chloride, respectively. Values of the spectral density function evaluated at a number of frequencies are compared. The levels of backbone dynamics in the folded protein show little variation across the molecule and are of similar magnitude to those determined previously for the folded state of the protein in exchange with an unfolded state at low salt concentrations [Farrow et al. (1995) Biochemistry 34, 868-878]. The denatured state of the domain exhibits both more extensive and more heterogeneous dynamics than the folded state. In particular the profile of the spectral density function evaluated at zero-frequency for the unfolded state of the domain indicates that residues in the middle of the protein sequence are considerably less mobile than those at the termini. These data suggest that the molecule is not behaving as an extended polymer and that concerted motions of the central portions of the molecule are occurring, consistent with a reasonably compact conformation in this region. The backbone dynamics of the folded and unfolded states were studied at two temperatures. The level of high-frequency motions in the folded molecule is largely unaffected by changes in temperature, whereas an increase in temperature results in increased high-frequency motion in the unfolded state.

Protein folding kinetic data have been obtained for the marginally stable N-terminal SH3 domain of the Drosophila protein drk as a function of pH in order to investigate the electrostatic properties of Asp8 in the folding transition state ensemble. The slow exchange between folded and unfolded forms of the protein gives rise to separate NMR resonances for both folded and unfolded states at equilibrium. As a result, kinetic data can be derived from magnetization transfer between these two states without the need for denaturants. Using the fact that ionization of Asp8 dominates the electrostatic behavior of the protein between pH 2 and 3, along with pKa values for titrating groups in both folded and unfolded states that have been determined in a previous study, values of 2.9 +/- 0.1 and 3.3 +/- 0.2 are obtained for the pKa of Asp8 in the transition state for the wild-type protein and for a His7Ala mutant, respectively. The data are consistent with the partial formation in the transition state ensemble of an Asp8 side chain carboxylate-a Lys21 backbone amide interaction that represents a highly conserved contact in folded SH3 domains.

Cystic fibrosis is caused by mutations in CFTR (cystic fibrosis transmembrane conductance regulator), leading to folding and processing defects and to chloride channel gating misfunction. CFTR is regulated by ATP binding to its cytoplasmic nucleotide-binding domains, NBD1 and NBD2, and by phosphorylation of the NBD1 regulatory insert (RI) and the regulatory extension (RE)/R region. These regulatory effects are transmitted to the rest of the channel via NBD interactions with intracellular domain coupling helices (CL), particularly CL4. Using a sensitive method for detecting inter-residue correlations between chemical shift changes in NMR spectra, an allosteric network was revealed within NBD1, with a construct lacking RI. The CL4-binding site couples to the RI-deletion site and the C-terminal residues of NBD1 that precede the R region in full-length CFTR. Titration of CL4 peptide into NBD1 perturbs the conformational ensemble in these sites with similar titration patterns observed in F508del, the major CF-causing mutant, and in suppressor mutants F494N, V510D and Q637R NBD1, as well as in a CL4-NBD1 fusion construct. Reciprocally, the C-terminal mutation, Q637R, perturbs dynamics in these three sites. This allosteric network suggests a mechanism synthesizing diverse regulatory NBD1 interactions and provides biophysical evidence for the allosteric coupling required for CFTR function.

The phosphotyrosine-binding (PTB) domain of Drosophila Shc (dShc) binds in vitro to phosphopeptides containing the sequence motif NPXpY, and physically associates with the activated Drosophila epidermal growth factor receptor homologue (DER) in vivo. The structural elements, specificity and binding kinetics of the dShc PTB domain have now been characterized. The dShc PTB domain appeared similar to the insulin-like receptor substrate-1 PTB domain in secondary structure as suggested by Fourier transform infrared spectroscopy. Surface plasmon resonance measurements indicated that the dShc PTB domain bound with high affinity to phosphopeptides (Der) derived from the Tyr1228 site of the DER receptor. The kinetics of the dShc PTB domain-Der phosphopeptide interaction differed from those of a typical SH2 domain-ligand interaction, in that the PTB domain displayed slower on/off rates. Competition binding assays using truncated versions of the Der peptides revealed that high affinity binding to the dShc PTB domain requires, in addition to the NPXpY motif, the presence of hydrophobic residues at both positions −5 and −7 relative to phosphotyrosine. The dShc PTB domain showed a similar binding specificity to the human Shc (hShc) PTB domain, but subtle differences were noted; such that the hShc PTB domain bound preferentially to a phosphopeptide from the mammalian nerve growth factor receptor, whereas the dShc PTB domain bound preferentially to phosphopeptides from the Drosophila DER receptor. The invertebrate dShc PTB domain therefore possesses a binding specificity for tyrosine-phosphorylated peptides that is optimally suited for recognition of the activated DER receptor.

Protein-protein interfaces can consist of interactions between large numbers of residues of each molecule; some of these interactions are critical in determining binding affinity and conferring specificity, while others appear to play only a marginal role. Src-homology-2 (SH2) domains bind to proteins containing phosphorylated tyrosines, with additional specificity provided by interactions with residues C-terminal to the phosphotyrosine (pTyr) residue. While the C-terminal SH2 domain of phospholipase C-gamma 1 (PLCC SH2) interacts with eight residues of a pTyr-containing peptide from its high affinity binding site on the beta-platelet-derived growth factor receptor, it can still bind tightly to a phosphopeptide containing only three residues. Novel deuterium (2H) based nuclear magnetic resonance (NMR) spin relaxation experiments which probe the nanosecond-picosecond time scale dynamics of methyl containing side chain residues have established that certain regions of the PLCC SH2 domain contacting the residues C-terminal to the pTyr have a high degree of mobility in both the free and peptide complexed states. In contrast, there is significant restriction of motion in the pTyr binding site. These results suggest a correlation between the dynamic behavior of certain groups in the PLCC SH2 complex and their contribution to high affinity binding and binding specificity.

Recent advances in NMR methodology have made it a powerful approach for the study of biomolecular structure and dynamics (Bax and Grzesiek, 1993; Bax, 1994; Farrow et al., 1994a). In addition to NMR being a structural tool, as a solution spectroscopy it is exquisitely sensitive to dynamic processes - not only fast, low amplitude motions which can often be described by analysis of X-ray crystallographic B factors (Ringe and Petsko, 1985), but also slower, larger amplitude motions. These may include conformational exchange on millisecond time-scales or longer between states as dissimilar as folded and unfolded states of proteins and reflecting motions of tens of angstroms. We have exploited this distinguishing capability of NMR spectroscopy to describe the dynamic processes observed during structural studies of two isolated domains of signal transduction proteins, a Src Homology 2 (SH2) domain of phospholipase Cγ in complex with a phosphopeptide from the platelet-derived growth factor receptor (PDGFR) and an isolated Src Homology 3 (SH3) domain from the drosophila protein Drk1.

The isolated N-terminal Src homology 3 (SH3) domain of Drosophila drk exists in equilibrium between folded and unfolded states in aqueous buffer near neutral pH. Nuclear magnetic resonance spectra recorded on both states simultaneously exhibit an approximate 1:1 ratio of protein conformations. The folded form is similar to other known SH3 structures, especially the N-terminal SH3 domain of the mammalian homologue GRB2. A stretch of sequential amide-amide nuclear Overhauser effect cross-peaks for resonances of the unfolded state is observed in a region corresponding to beta-strands in the folded state. The results suggest that turn-like conformations may be preferentially sampled in the folding pathway for this predominantly beta-structured SH3 domain. In addition, a stable turn at Leu-28 is observed in the unfolded but not in the folded state. Comparison of this unfolded form with a denatured state in 2 M guanidine hydrochloride shows that, while both are highly disordered, these states are not identical and more residual structure is present under nondenaturing conditions.

Two-dimensional NMR 15N relaxation studies have been used to characterize the backbone dynamics and folding transition of the N-terminal SH3 domain of the protein drk (drkN SH3). The isolated drkN SH3 domain exists in equilibrium between a folded and an unfolded state in aqueous buffer and near neutral pH [Zhang et al. (1994) J. Biomol. NMR 4, 845]. The backbone dynamics of both the folded and unfolded states in this exchanging system have been determined and the rates of the folding transition estimated at 14 degrees C. For 12 residues, the values of the spectral density functions of the backbone amide bond vectors at a number of frequencies have been established. Results show that while the unfolded state has considerably greater dynamic behavior, the overall motional properties are consistent with it having a reasonably compact structure in solution. In contrast to the folded state, considerable variations are seen in the values of the spectral densities of the unfolded state as a function of residue number; these variations do not appear to be strongly correlated with structural elements in the folded state. The mean value of the exchange rate associated with the folding transition was determined to be 0.89 s-1, similar to previous measurements of the rate of formation of beta-structure.

The SH3 domain is a commonly found modular binding domain that functions to mediate protein-protein interactions and has been a frequent target of protein folding and other biophysical studies. The N-terminal SH3 domain (6.8 kDa) of the Drosophila adapter protein Drk is marginally stable under physiological conditions and exists in nearly 50/50 equilibrium between the folded and unfolded states, with a slow conversion rate of 2.2 s−1(1). Due to the dynamic ensemble nature of this unfolded state under non-denaturing conditions it is difficult to characterize structurally. Experimental data have been used to define sets of heterogeneous conformations but ensemble averaging minimizes the information content (2). Importantly, little is known about inter-conversion rates within the unfolded ensemble.Fluorescence correlation spectroscopy (FCS) experiments performed on the wild-type and on a mutant that stabilizes the folded state (T22G) provided the hydrodynamic radii of the folded (Rh,f), unfolded (Rh,u) and denatured (Rh,d) states at 20°C. Quite surprisingly, the results indicate that the unfolded state under non-denaturing conditions is less compact than the chemically-denatured state (Rh,u/Rh,f=1.31±0.05 vs. Rh,d/Rh,f=1.21±0.03). However, this can be explained by introducing adequate shape factors in previously reported NMR/SAXS data. In addition the FCS measured compactness of unfolded state is in an excellent agreement with PFG-NMR value of Rh,u/Rh,f=1.30±0.01 (2). Furthermore, FCS analysis reveals a fast 200 ns process within the unfolded ensemble. To resolve slow (milli)second conformational dynamics, dual cysteine mutant (E2C-60C) was labelled with the Bodipy Fl/Alexa 647 FRET pair and single proteins were encapsulated in surface-tethered liposomes for long observations. This study provides further insights into the unfolded state heterogeneity and its importance in protein folding in general.1.M Tollinger et al., J. Am. Chem. Soc.;123(46):11341-11352 (2001).2.WY Choy et al., J. Mol. Biol.;316(1):101-112 (2002).

The backbone dynamics of the C-terminal SH2 domain of phospholipase C gamma 1 have been investigated. Two forms of the domain were studied, one in complex with a high-affinity binding peptide derived from the platelet-derived growth factor receptor and the other in the absence of this peptide. 2-D 1H-15N NMR methods, employing pulsed field gradients, were used to determine steady-state 1H-15N NOE values and T1 and T2 15N relaxation times. Backbone dynamics were characterized by the overall correlation time (tau m), order parameters (S2), effective correlation times for internal motions (tau e), and, if required, terms to account for motions on a microsecond-to-millisecond-time scale. An extended two-time-scale formalism was used for residues having relaxation data and that could not be fit adequately using a single-time-scale formalism. The overall correlation times of the uncomplexed and complexed forms of SH2 were found to be 9.2 and 6.5 ns, respectively, suggesting that the uncomplexed form is in a monomer-dimer equilibrium. This was subsequently confirmed by hydrodynamic measurements. Analysis of order parameters reveals that residues in the so-called phosphotyrosine-binding loop exhibited higher than average disorder in both forms of SH2. Although localized differences in order parameters were observed between the uncomplexed and complexed forms of SH2, overall, higher order parameters were not found in the peptide-bound form, indicating that on average, picosecond-time-scale disorder is not reduced upon binding peptide. The relaxation data of the SH2-phosphopeptide complex were fit with fewer exchange terms than the uncomplexed form. This may reflect the monomer-dimer equilibrium that exists in the uncomplexed form or may indicate that the complexed form has lower conformational flexibility on a microsecond-to-millisecond-time scale.

The solution structure of the C-terminal SH2 domain of phospholipase C-gamma 1 (PLC-gamma 1), in complex with a phosphopeptide corresponding to its Tyr-1021 high affinity binding site on the platelet-derived growth factor receptor, has been determined by nuclear magnetic resonance spectroscopy. The topology of the SH2-phosphopeptide complex is similar to previously reported Src and Lck SH2 complexes. However, the binding site for residues C-terminal to the phosphotyrosine (pTyr) is an extended groove that contacts peptide residues at the +1 to +6 positions relative to the pTyr. This striking difference from Src and Lck reflects the fact that the PLC-gamma 1 complex involves binding of a phosphopeptide with predominantly hydrophobic residues C-terminal to the pTyr and therefore serves as a prototype for a second class of SH2-phosphopeptide interactions.

The complete assignment of 1H and 15N backbone resonances and near-complete 1H side-chain resonance assignments have been obtained for the reduced form of a mutant of human thioredoxin (105 residues) in which the three non-active site cysteines have been substituted by alanines: C62A, C69A, C73A. The assignments were made primarily on the basis of three-dimensional 15N-separated nuclear Overhauser and Hartmann-Hahn spectroscopy, in conjunction with two-dimensional homonuclear and heteronuclear correlation experiments. Based on comparisons of short-range and interstrand nuclear Overhauser effects, patterns of amide exchange, and chemical-shift differences, the structure appears essentially unchanged from that of the previously determined solution structure of the native protein [Forman-Kay, J.D. et al. (1991) Biochemistry, 30, 2685-2698]. An assay for thioredoxin shows that the C62A, C69A, C73A mutant retains activity. The assignment of the spectrum for this mutant of human thioredoxin constitutes the basis for future studies aimed at comparing the details of the active-site conformation in the reduced and oxidized forms of the protein.

Intrinsically disordered proteins play important roles in mediating regulatory interactions. While many of these proteins fold upon binding to targets, others appear to be only transiently and locally ordered, existing in a dynamic complex with protein partners. The interaction of the disordered Sic1 (a cyclin dependent kinase inhibitor) with Cdc4 (a component of an SCF ubiquitin ligase) involves the dynamic exchange of multiple linear binding motifs within Sic1 on and off of a single Cdc4 binding site. The disordered nature of the Sic1 chain enables Cdc4 to recognize multiply phosphorylated Sic1 by sensing a mean electrostatic field rather than only the field due to local binding of a linear motif. This sensitivity to multiple phosphorylations may facilitate a switch-like biological response of Sic1 ubiquitination and degradation to cyclin kinase activity. The regulatory (R) region of the cystic fibrosis transmembrane conductance regulator (CFTR) is another example, with transient and local ordering of multiple segments upon binding to multiple target proteins. The R region dynamic complex with the first nucleotide binding domain (NBD1) of CFTR appears to facilitate a rheostat-like function in activating CFTR channel conductance with an ability to integrate multiple regulatory binding inputs. We have measured NMR and SAXS structural data and used them as input to our program ENSEMBLE to generate sets of coordinates that represent significantly populated conformers within the isolated disordered states of Sic1 and the CFTR R region. The resulting ensembles can be docked to structures of the folded binding domains to further enable characterization of these dynamic complexes, leading to a better understanding of the role of disordered protein interactions in biological regulation.

The solution structure of recombinant human thioredoxin (105 residues) has been determined by nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Approximate interproton distance restraints were derived from nuclear Overhauser effect (NOE) measurements. In addition, a large number of stereospecific assignments for beta-methylene protons and torsion angle restraints for phi, psi, and chi 1 were obtained by using a conformational grid search on the basis of the intraresidue and sequential NOE data in conjunction with 3JHN alpha and 3J alpha beta coupling constants. The structure calculations were based on 1983 approximate interproton distance restraints, 52 hydrogen-bonding restraints for 26 hydrogen bonds, and 98 phi, 71 psi, and 72 chi 1 torsion angle restraints. The 33 final simulated annealing structures obtained had an average atomic rms distribution of the individual structures about the mean coordinate positions of 0.40 +/- 0.06 A for the backbone atoms and 0.78 +/- 0.05 A for all atoms. The solution structure of human thioredoxin consists of a five-stranded beta-sheet surrounded by four alpha-helices, with an active site protrusion containing the two redox-active cysteines. The overall structure is similar to the crystal and NMR structures of oxidized [Katti, S. K., LeMaster, D. M., & Eklund, H. (1990) J. Mol. Biol. 212, 167-184] and reduced [Dyson, J. H., Gippert, G. P., Case, D. A., Holmgren, A., & Wright, P. (1990) Biochemistry 29, 4129-4136] Escherichia coli thioredoxin, respectively, despite the moderate 25% amino acid sequence homology. Several differences, however, can be noted. The human alpha 1 helix is a full turn longer than the corresponding helix in E. coli thioredoxin and is characterized by a more regular helical geometry. The helix labeled alpha 3 in human thioredoxin has its counterpart in the 3(10) helix of the E. coli protein and is also longer in the human protein. In contrast to these structural differences, the conformation of the active site loop in both proteins is very similar, reflecting the perfect sequence identity for a stretch of eight amino acid residues around the redox-active cysteines.

Assumptions of restricted flexibility upon binding conflict with emerging data showing that motion can increase, decrease or stay the same within molecular complexes. Now, calculations of entropic contributions from dynamics at specific positions in a complex suggest that increases in motion can dominate the free energy of association in certain cases.

Two-dimensional {sup 1}H NMR spectroscopy has been applied to a structural analysis of the reduced form of a recombinant human thioredoxin, a ubiquitous dithiol oxidoreductase recently isolated from an immunocompetent lymphoblastoid cell line. The sequential assignment of the spectrum, including all proline residues, has been accomplished by using experiments to demonstrate through-bond and through-space connectivities. The secondary structure has been determined by a qualitative interpretation of nuclear Overhauser effects, NH exchange data, and {sup 3}J{sub HN{alpha}} coupling constants. The secondary structure was found to be similar to that of the X-ray structure of Escherichia coli thioredoxin, consisting of a mixed five-stranded {beta}-sheet surrounded by four {alpha}-helices. The assignment and structural characterization of human thioredoxin was facilitated by the increased resolution and sensitivity afforded by a magnetic field strength of 600 MHz and required the use of two temperatures and two pH conditions to resolve ambiguities caused by a duplication of resonances. This duplication, extending from Phe-41 to Val-59, and including Lys-3-Ile-5, Val-24, Val-25, Asn-39, and Ile-101-Glu-103, appears to be due to heterogeneity arising from the presence or absence of the N-terminal methionine.

Intrinsically disordered proteins (IDPs), which lack persistent structure, are a challenge to structural biology due to the inapplicability of standard methods for characterization of folded proteins as well as their deviation from the dominant structure/function paradigm. Their widespread presence and involvement in biological function, however, has spurred the growing acceptance of the importance of IDPs and the development of new tools for studying their structure, dynamics and function. The interplay of folded and disordered domains or regions for function and the existence of a continuum of protein states with respect to conformational energetics, motional timescales and compactness is shaping a unified understanding of structure-dynamics-disorder/function relationships. On the 20th anniversary of this journal, Structure, we provide a historical perspective on the investigation of IDPs and summarize the sequence features and physical forces that underlie their unique structural, functional and evolutionary properties.

SummaryIntrinsically disordered proteins can form highly dynamic complexes with partner proteins. One such dynamic complex involves the intrinsically disordered Sic1 with its partner Cdc4 in regulation of yeast cell cycle progression. Phosphorylation of six N-terminal Sic1 sites leads to equilibrium engagement of each phosphorylation site with the primary binding pocket in Cdc4, the substrate recognition subunit of a ubiquitin ligase. ENSEMBLE calculations using experimental nuclear magnetic resonance and small-angle X-ray scattering data reveal significant transient structure in both phosphorylation states of the isolated ensembles (Sic1 and pSic1) that modulates their electrostatic potential, suggesting a structural basis for the proposed strong contribution of electrostatics to binding. A structural model of the dynamic pSic1-Cdc4 complex demonstrates the spatial arrangements in the ubiquitin ligase complex. These results provide a physical picture of a protein that is predominantly disordered in both its free and bound states, enabling aspects of its structure/function relationship to be elucidated.

Background: The Src homology 2 (SH2) domains of cytoplasmic signaling proteins generally bind phosphotyrosine (pTyr) sites in the context of carboxy-terminal residues. SAP (also known as SH2D1A or DSHP), the product of the gene that is mutated in human X-linked lymphoproliferative (XLP) disease, comprises almost exclusively a single SH2 domain, which may modulate T-cell signaling by engaging T-cell co-activators such as SLAM, thereby blocking binding of other signaling proteins that contain SH2 domains. The SAP–SLAM interaction can occur in a phosphorylation-independent manner. Results: To characterize the interaction between SAP and SLAM, we synthesized peptides corresponding to the SAP-binding site at residue Y281 in SLAM. Both phosphorylated and non-phosphorylated versions of an 11-residue SLAM peptide bound SAP, with dissociation constants of 150 nM and 330 nM, respectively. SLAM phosphopeptides that were truncated either at the amino or carboxyl terminus bound with high affinity to SAP, suggesting that the SAP SH2 domain recognizes both amino-terminal and carboxy-terminal sequences relative to the pTyr residue. These results were confirmed by nuclear magnetic resonance (NMR) studies on 15 N- and 13 C-labeled SAP complexed with three SLAM peptides: an amino-terminally truncated phosphopeptide, a carboxy-terminally truncated phosphopeptide and a non-phosphorylated Tyr-containing full-length peptide. Conclusions: The SAP SH2 domain has a unique specificity. Not only does it bind peptides in a phosphorylation-independent manner, it also recognizes a pTyr residue either preceded by amino-terminal residues or followed by carboxy-terminal residues. We propose that the three ‘prongs' of a peptide ligand (the amino and carboxyl termini and the pTyr) can engage the SAP SH2 domain, accounting for its unusual properties. These data point to the flexibility of modular protein-interaction domains.

Intrinsically disordered proteins (IDPs) play critical yet often poorly understood roles in a variety of cellular processes. Despite an apparent antagonism between structural disorder and protein recognition, the large number of IDPs involved in protein regulation suggests that they could in fact provide advantages in recognition over well-folded proteins. Sic1 is an IDP inhibitor of a cyclin dependent kinase (CDK) in yeast, which interacts with a single site on its acceptor Cdc4 only upon phosphorylation of its multiple dispersed CDK sites. The multiple phosphorylation events can in principle be the basis for ultrasensitivity in protein-protein binding, however its molecular basis remains elusive1.We performed a systematic study of the Sic1's fluctuating conformations in different salt and denaturant concentrations using single-molecule Forster energy transfer (smFRET) and fluorescence correlation spectroscopy (FCS). smFRET data suggests that Sic1 protein is comprised of a continuum of conformers with varying end-end distances. Denaturant titration measurements suggest that these conformers are characterized by non-cooperative interactions. From FCS experiments, the exchange between Sic1 conformational states was found to occur on both ultrafast (10-100 ns timescale) and slow (10-100 ms) timescale. Burst smFRET experiments show that Sic1’s end-end distances do not vary significantly upon addition of salt, which suggests that charge-shielding by salt may only affect the structure locally, around charged amino acids. Our single-molecule data resolves conformational heterogeneity and dynamics in a model IDP protein and represent the first step towards the validation of the polyelectrostatic model of the Sic1-Cdc4 interaction2.(1) Nash, P., Tang, X., Orlicky, S., Chen, Q., Gertler, F. B., Mendenhall, M. D., Sicheri, F., Pawson, T., Tyers, M. Nature 2001, 414, 514.(2) Borg, M., Mittag, T., Pawson, T., Tyers, M., Forman-Kay, J. D., Chan, H. S. Proc. Natl. Acad. Sci. USA 2007, 104, 9650.