Search

Intrinsically disordered proteins and unfolded proteins have fluctuating conformational ensembles that are fundamental to their biological function and impact protein folding, stability and misfolding. Despite the importance of protein dynamics and conformational sampling, time-dependent data types are not fully exploited when defining and refining disordered protein ensembles. Here we introduce a computational framework using an elastic network model and normal mode displacements to generate a dynamic disordered ensemble consistent with NMR-derived dynamics parameters, including transverse [Formula: see text] relaxation rates and Lipari-Szabo order parameters ([Formula: see text] values). We illustrate our approach using the unfolded state of the drkN SH3 domain to show that the dynamical ensembles give better agreement than a static ensemble for a wide range of experimental validation data including NMR chemical shifts, J-couplings, nuclear Overhauser effects, paramagnetic relaxation enhancements, residual dipolar couplings, hydrodynamic radii, single-molecule fluorescence Förster resonance energy transfer, and small-angle X-ray scattering.

Many pairwise additive force fields are in active use for intrinsically disordered proteins (IDPs) and regions (IDRs), some of which modify energetic terms to improve the description of IDPs/IDRs but are largely in disagreement with solution experiments for the disordered states. This work considers a new direction—the connection to configurational entropy—and how it might change the nature of our understanding of protein force field development to equally well encompass globular proteins, IDRs/IDPs, and disorder-to-order transitions. We have evaluated representative pairwise and many-body protein and water force fields against experimental data on representative IDPs and IDRs, a peptide that undergoes a disorder-to-order transition, for seven globular proteins ranging in size from 130 to 266 amino acids. We find that force fields with the largest statistical fluctuations consistent with the radius of gyration and universal Lindemann values for folded states simultaneously better describe IDPs and IDRs and disorder-to-order transitions. Hence, the crux of what a force field should exhibit to well describe IDRs/IDPs is not just the balance between protein and water energetics but the balance between energetic effects and configurational entropy of folded states of globular proteins.

Understanding signaling and other complex biological processes requires elucidating the critical roles of intrinsically disordered proteins (IDPs) and regions (IDRs), which represent ∼30% of the proteome and enable unique regulatory mechanisms. In this review, we describe the structural heterogeneity of disordered proteins that underpins these mechanisms and the latest progress in obtaining structural descriptions of conformational ensembles of disordered proteins that are needed for linking structure and dynamics to function. We describe the diverse interactions of IDPs that can have unusual characteristics such as "ultrasensitivity" and "regulated folding and unfolding". We also summarize the mounting data showing that large-scale assembly and protein phase separation occurs within a variety of signaling complexes and cellular structures. In addition, we discuss efforts to therapeutically target disordered proteins with small molecules. Overall, we interpret the remodeling of disordered state ensembles due to binding and post-translational modifications within an expanded framework for allostery that provides significant insights into how disordered proteins transmit biological information.

Increased understanding of intrinsically disordered proteins (IDPs) and protein regions has revolutionized our view of the relationship between protein structure and function. Data now support that IDPs can be functional in the absence of a single, fixed, three-dimensional structure. Due to their dynamic morphology, IDPs have the ability to display a range of kinetics and affinity depending on what the system requires, as well as the potential for large-scale association. Although several studies have shed light on the functional properties of IDPs, the class of intrinsically disordered transcription factors (TFs) is still poorly characterized biophysically due to their combination of ordered and disordered sequences. In addition, TF modulation by small molecules has long been considered a difficult or even impossible task, limiting functional probe development. However, with evolving technology, it is becoming possible to characterize TF structure-function relationships in unprecedented detail and explore avenues not available or not considered in the past. Here we provide an introduction to the biophysical properties of intrinsically disordered TFs and we discuss recent computational and experimental efforts toward understanding the role of intrinsically disordered TFs in biology and disease. We describe a series of successful TF targeting strategies that have overcome the perception of the "undruggability" of TFs, providing new leads on drug development methodologies. Lastly, we discuss future challenges and opportunities to enhance our understanding of the structure-function relationship of intrinsically disordered TFs.

Post-translational modifications (PTMs), which are found largely in intrinsically disordered protein regions (IDRs), regulate protein activity, stability and interactions with partners. They are therefore critical for controlling essentially all cellular processes. A single modification event can have dramatic effects; however, proteins are often modified on multiple sites to collectively modulate the biological outcome. Multiple PTMs can mediate the same, complementary or opposing effects and the result of their interplay is determined by a complex combination of the number, positioning and type of modifications. Multiple PTMs can also synergize to shift the conformational or binding equilibria of the modified protein to modulate its interaction with partners or formation of higher order assembly. Recognition of such PTM crosstalk is crucial for understanding the underlying mechanisms of complex regulatory processes.

Misofolding of mammalian prion proteins (PrP) is believed to be the cause of a group of rare and fatal neurodegenerative diseases. Despite intense scrutiny however, the mechanism of the misfolding reaction remains unclear. We perform nuclear Magnetic Resonance and thermodynamic stability measurements on the C-terminal domains (residues 90–231) of two PrP variants exhibiting different pH-induced susceptibilities to aggregation: the susceptible hamster prion (GHaPrP) and its less susceptible rabbit homolog (RaPrP). The pKa of histidines in these domains are determined from titration experiments, and proton-exchange rates are measured at pH 5 and pH 7. A single buried highly conserved histidine, H187/H186 in GHaPrP/RaPrP, exhibited a markedly down shifted pKa ~5 for both proteins. However, noticeably larger pH-induced shifts in exchange rates occur for GHaPrP versus RaPrP. Analysis of the data indicates that protonation of the buried histidine destabilizes both PrP variants, but produces a more drastic effect in the less stable GHaPrP. This interpretation is supported by urea denaturation experiments performed on both PrP variants at neutral and low pH, and correlates with the difference in disease susceptibility of the two species, as expected from the documented linkage between destabilization of the folded state and formation of misfolded and aggregated species.

Chloride channel gating and trafficking of the cystic fibrosis transmembrane conductance regulator (CFTR) are regulated by phosphorylation. Intrinsically disordered segments of the protein are responsible for phospho-regulation, particularly the regulatory (R) region that is a target for several kinases and phosphatases. The R region remains disordered following phosphorylation, with different phosphorylation states sampling various conformations. Recent studies have demonstrated the crucial role that intramolecular and intermolecular interactions of the R region play in CFTR regulation. Different partners compete for the same binding segment, with the R region containing multiple overlapping binding elements. The non-phosphorylated R region interacts with the nucleotide binding domains and inhibits channel activity by blocking heterodimerization. Phosphorylation shifts the equilibrium such that the R region is excluded from the dimer interface, facilitating gating and processing by stimulating R region interactions with other domains and proteins. The dynamic conformational sampling and transient binding of the R region to multiple partners enables complex control of CFTR channel activity and trafficking.

© 2016 American Chemical Society. The traditional structure-function paradigm has provided significant insights for well-folded proteins in which structures can be easily and rapidly revealed by X-ray crystallography beamlines. However, approximately one-third of the human proteome is comprised of intrinsically disordered proteins and regions (IDPs/IDRs) that do not adopt a dominant well-folded structure, and therefore remain "unseen" by traditional structural biology methods. This Perspective considers the challenges raised by the "Dark Proteome", in which determining the diverse conformational substates of IDPs in their free states, in encounter complexes of bound states, and in complexes retaining significant disorder requires an unprecedented level of integration of multiple and complementary solution-based experiments that are analyzed with state-of-the art molecular simulation, Bayesian probabilistic models, and high-throughput computation. We envision how these diverse experimental and computational tools can work together through formation of a "computational beamline" that will allow key functional features to be identified in IDP structural ensembles.

Cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel in the apical surface of epithelial cells in the airway and gastrointestinal tract, and mutation of CFTR is the underlying cause of cystic fibrosis. However, the precise molecular details of the structure and function of CFTR in native and disease states remains elusive and cystic fibrosis researchers are hindered by a lack of high specificity, high affinity binding reagents for use in structural and biological studies. Here, we describe a panel of synthetic antigen-binding fragments (Fabs) isolated from a phage-displayed library that are specific for intracellular domains of CFTR that include the nucleotide-binding domains (NBD1 and NBD2), the R-region, and the regulatory insertion loop of NBD1. Binding assays performed under conditions that promote the native fold of the protein demonstrated that all Fabs recognized full-length CFTR. However, only the NBD1-specific Fab recognized denatured CFTR by western blot, suggesting a conformational epitope requirement for the other Fabs. Surface plasmon resonance experiments showed that the R-region Fab binds with high affinity to both the phosphorylated and unphosphorylated R-region. In addition, NMR analysis of bound versus unbound R-region revealed a distinct conformational effect upon Fab binding. We further defined residues involved with antibody recognition using an overlapping peptide array. In summary, we describe methodology complementary to previous hybridoma-based efforts to develop antibody reagents to CFTR, and introduce a synthetic antibody panel to aid structural and biological studies.

Complete folding is not a prerequisite for protein function, as disordered and partially folded states of proteins frequently perform essential biological functions. In order to understand their functions at the molecular level, we utilized diverse experimental measurements to calculate ensemble models of three non-homologous, intrinsically disordered proteins: I-2, spinophilin and DARPP-32, which bind to and regulate protein phosphatase 1 (PP1). The models demonstrate that these proteins have dissimilar propensities for secondary and tertiary structure in their unbound forms. Direct comparison of these ensemble models with recently determined PP1 complex structures suggests a significant role for transient, pre-formed structure in the interactions of these proteins with PP1. Finally, we generated an ensemble model of partially disordered I-2 bound to PP1 that provides insight into the relationship between flexibility and biological function in this dynamic complex.

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 (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) 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 alpha-helical structure and both native-like and nonnative tertiary structure.

Post-translational modifications (PTMs) produce significant changes in the structural properties of intrinsically disordered proteins (IDPs) by affecting their energy landscapes. PTMs can induce a range of effects, from local stabilization or destabilization of transient secondary structure to global disorder-to-order transitions, potentially driving complete state changes between intrinsically disordered and folded states or dispersed monomeric and phase-separated states. Here, we discuss diverse biological processes that are dependent on PTM regulation of IDPs. We also present recent tools for generating homogenously modified IDPs for studies of PTM-mediated IDP regulatory mechanisms.

Synuclein-γ (SNCG) plays oncogenic roles in breast carcinogenesis. Although the expression of SNCG is abnormally high in advanced and metastatic breast carcinomas, SNCG is not expressed in normal or benign breast tissues. SNCG is an intrinsically disordered protein known to interact with BubR1, a mitotic checkpoint kinase. The SNCG-BubR1 interaction inhibits mitotic checkpoint control upon spindle damage caused by anticancer drugs, such as nocodazole and taxol. Antimicrotubule drugs that cause mitotic arrest and subsequent apoptosis of cancer cells are frequently used to treat breast cancer patients with advanced or metastatic diseases. However, patient response rates to this class of chemotherapeutic agents vary significantly. In this study, we have designed a novel peptide (ANK) and shown its interaction with SNCG using fluorometry, surface plasmon resonance, and isothermal titration calorimetry. Binding of the ANK peptide did not induce folding of SNCG, suggesting that SNCG can function biologically in its intrinsically disordered state. Microinjection of the ANK peptide in breast cancer cell line overexpressing SNCG (MCF7-SNCG) exhibited a similar cell killing response by nocodazole as in the SNCG-negative MCF7 cells. Overexpression of enhanced green fluorescent protein–tagged ANK reduces SNCG-mediated resistance to paclitaxel treatment by ∼3.5-fold. Our coimmunoprecipitation and colocalization results confirmed the intracellular association of the ANK peptide with SNCG. This is likely due to the disruption of the interaction of SNCG with BubR1 interaction. Our findings shed light on the molecular mechanism of the ANK peptide in releasing SNCG-mediated drug resistance. [Cancer Res 2007;67(2):626–33]

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.

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

SummaryInteractions between the WW domains of Drosophila Nedd4 (dNedd4) and Commissureless (Comm) PY motifs promote axon crossing at the CNS midline and muscle synaptogenesis. Here we report the solution structure of the dNedd4 WW3∗ domain complexed to the second PY motif (227′TGLPSYDEALH237′) of Comm. Unexpectedly, there are interactions between WW3∗ and ligand residues both N- and C-terminal to the PY motif. Residues Y232′–L236′ form a helical turn, following the PPII helical PY motif. Mutagenesis and binding studies confirm the importance of these extensive contacts, not simultaneously observed in other WW domain complexes, and identify a variable loop in WW3∗ responsible for its high-affinity interaction. These studies expand our general understanding of the molecular determinants involved in WW domain-ligand recognition. In addition, they provide insights into the specific regulation of dNedd4-mediated ubiquitination of Comm and subsequent internalization of Comm or the Comm/Roundabout complex, critical for CNS and muscle development.

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.

The bacterial outer membrane enzyme PagP transfers a palmitate chain from a phospholipid to lipid A. In a number of pathogenic Gram-negative bacteria, PagP confers resistance to certain cationic antimicrobial peptides produced during the host innate immune response. The global fold of Escherichia coli PagP was determined in both dodecylphosphocholine and n -octyl-β- d -glucoside detergent micelles using solution NMR spectroscopy. PagP consists of an eight-stranded anti-parallel β-barrel preceded by an amphipathic α helix. The β-barrel is well defined, whereas NMR relaxation measurements reveal considerable mobility in the loops connecting individual β-strands. Three amino acid residues critical for enzymatic activity localize to extracellular loops near the membrane interface, positioning them optimally to interact with the polar headgroups of lipid A. Hence, the active site of PagP is situated on the outer surface of the outer membrane. Because the phospholipids that donate palmitate in the enzymatic reaction are normally found only in the inner leaflet of the outer membrane, PagP activity may depend on the aberrant migration of phospholipids into the outer leaflet. This finding is consistent with an emerging paradigm for outer membrane enzymes in providing an adaptive response toward disturbances in the outer membrane.

A triple-resonance NMR pulse scheme is presented for measuring aspartic and glutamic acid side-chain pK(a) values in unfolded protein states where chemical shift overlap is limiting. The experiment correlates side-chain carboxyl carbon chemical shifts of these residues with the backbone amide proton chemical shift of the following residue. The methodology is applied to an (15)N, (13)C labeled sample of the N-terminal SH3 domain of the Drosophila protein drk, which exists in equilibrium between folded (F(exch)) and unfolded (U(exch)) states under nondenaturing conditions. Residue-specific pK(a) values of side-chain carboxyl groups are presented for the first time for an unfolded protein (drk U(exch) state), determined from a pH titration. Results indicate that deviations from pK(a) values measured for model compounds are likely due to local effects, while long-range electrostatic interactions appear to be of minor importance for this protein.

Intrinsically disordered proteins play important roles in cell signalling, transcription, translation and cell cycle regulation. Although they lack stable tertiary structure, many intrinsically disordered proteins undergo disorder-to-order transitions upon binding to partners. Similarly, several folded proteins use regulated order-to-disorder transitions to mediate biological function. In principle, the function of intrinsically disordered proteins may be controlled by post-translational modifications that lead to structural changes such as folding, although this has not been observed. Here we show that multisite phosphorylation induces folding of the intrinsically disordered 4E-BP2, the major neural isoform of the family of three mammalian proteins that bind eIF4E and suppress cap-dependent translation initiation. In its non-phosphorylated state, 4E-BP2 interacts tightly with eIF4E using both a canonical YXXXXLΦ motif (starting at Y54) that undergoes a disorder-to-helix transition upon binding and a dynamic secondary binding site. We demonstrate that phosphorylation at T37 and T46 induces folding of residues P18-R62 of 4E-BP2 into a four-stranded β-domain that sequesters the helical YXXXXLΦ motif into a partly buried β-strand, blocking its accessibility to eIF4E. The folded state of pT37pT46 4E-BP2 is weakly stable, decreasing affinity by 100-fold and leading to an order-to-disorder transition upon binding to eIF4E, whereas fully phosphorylated 4E-BP2 is more stable, decreasing affinity by a factor of approximately 4,000. These results highlight stabilization of a phosphorylation-induced fold as the essential mechanism for phospho-regulation of the 4E-BP:eIF4E interaction and exemplify a new mode of biological regulation mediated by intrinsically disordered proteins.

The global fold of maltose-binding protein in complex with the substrate beta-cyclodextrin was determined by solution NMR methods. The two-domain protein is comprised of a single polypeptide chain of 370 residues, with a molecular mass of 42 kDa. Distance information in the form of H(N)-H(N), H(N)-CH(3) and CH(3)-CH(3) NOEs was recorded on (15)N, (2)H and (15)N, (13)C, (2)H-labeled proteins with methyl protonation in Val, Leu, and Ile (C(delta1) only) residues. Distances to methyl protons, critical for the structure determination, comprised 77 % of the long-range restraints. Initial structures were calculated on the basis of 1943 NOEs, 48 hydrogen bond and 555 dihedral angle restraints. A global pair-wise backbone rmsd of 5.5 A was obtained for these initial structures with rmsd values for the N and C domains of 2.4 and 3.8 A, respectively. Direct refinement against one-bond (1)H(N)-(15)N, (13)C(alpha)-(13)CO, (15)N-(13)CO, two-bond (1)H(N)-(13)CO and three-bond (1)H(N)-(13)C(alpha) dipolar couplings resulted in structures with large numbers of dipolar restraint violations. As an alternative to direct refinement against measured dipolar couplings we have developed an approach where discrete orientations are calculated for each peptide plane on the basis of the dipolar couplings described above. The orientation which best matches that in initial NMR structures calculated from NOE and dihedral angle restraints exclusively is used to refine further the structures using a new module written for CNS. Modeling studies from four different proteins with diverse structural motifs establishes the utility of the methodology. When applied to experimental data recorded on MBP the precision of the family of structures generated improves from 5.5 to 2.2 A, while the rmsd with respect to the X-ray structure (1dmb) is reduced from 5.1 to 3.3 A.

We have expanded our description of the energy landscape for folding of the SH3 domain of chicken alpha-spectrin by a detailed structural characterization of its denatured state ensemble (DSE). This DSE is significantly populated under mildly acidic conditions in equilibrium with the folded state. Evidence from heteronuclear nuclear magnetic resonance (NMR) experiments on (2)H, (15)N-labeled protein suggests the presence of conformers whose residual structure bears some resemblence to the structure of the folding transition state of this protein. NMR analysis in a mutant with an engineered, non-native alpha-helical tendency shows a significant amount of local non-native structure in the mutant, while the overall characteristics of the DSE are unchanged. Comparison with recent theoretical predictions of SH3 domain folding reactions reveals an interesting correlation with the predicted early events. Based on these results and recent data from other systems, we propose that the DSE of a protein will resemble the intermediate or transition state of its nearest rate-limiting step, as a consequence of simple energetic and kinetic principles.

The structures of the denatured states of the spectrin SH3 domain and a mutant designed to have a non-native helical tendency at the N terminus have been analyzed under mild acidic denaturing conditions by nuclear magnetic resonance methods with improved resolution. The wild-type denatured state has little residual structure. However, the denatured state of the mutant has an approximately 50% populated helical structure from residues 2 to 14, a region that forms part of the beta-sheet structure in the folded state. Comparison with a peptide corresponding to the same sequence shows that the helix is stabilized in the whole domain, likely by non-local interactions with other parts of the protein as suggested by changes in a region far from the mutated sequence. These results demonstrate that high populations of non-native secondary structure elements in the denatured state are compatible with the formation of the native folded structure.

The isolated N-terminal SH3 domain of the Drosophila adapter protein drk (drkN SH3 domain) exists in a dynamic equilibrium between a folded (F(exch)) and an unfolded (U(exch)) state under native-like buffer conditions [Zhang, O.,Forman-Kay, J. D. (1995) Biochemistry 34, 6784-6794]. The effect of binding a proline-rich peptide derived from the protein Sos, a biological target of the drkN SH3 domain, on this equilibrium has been investigated. The stabilization of the F(exch) folded state upon binding provides an example of the link between binding and protein folding or stabilization. We have compared NMR parameters of the U(exch) state with those of a denatured state in 2 M guanidine hydrochloride (U(Gdn)). Variable-temperature experiments demonstrate that interactions in a region encompassing residues Gln 23-Leu 28 in the U(exch) state are destabilized upon addition of guanidine hydrochloride. Data from an 15N HSQC-NOESY-HSQC experiment as well as recently developed methods provide more unambiguous structural information than described previously, showing the presence of preferential structure in both unfolded states. Backbone NOEs observed in both unfolded states as well as chemical shifts and coupling constants suggest a rapid equilibrium between extended structure and turn-like structures which may play a role in initiation of protein folding. However, differences in detailed structural features between the two unfolded states argue that caution is needed in interpretation of results from structural characterization of protein conformational states generated using denaturing conditions.

NMR-based structural studies of macromolecules focus to a large extent on the establishmentof interproton distances within the molecule based on the nuclear Overhauser effect (NOE).Despite the improvements in resolution resulting from multidimensional NMR experiments,the detailed characterization of disordered states of proteins or highly overlapped regions offolded molecules using current NMR methods remains challenging. A suite of triple-resonanceNOESY-type pulse schemes is presented which require uniform 15N and 13C labeling andmake use of the chemical shift dispersion of backbone 15N and 13C′ (carbonyl)resonances to increase the spectral resolution. In particular, for the case of partially folded andunfolded proteins, the experiments exploit the fact that the dispersion of 15N and 13C′resonances is comparable to that observed in folded states. Ambiguities that arise in theassignment of NOEs as a result of the severe chemical shift degeneracy in 1H and aliphatic13C nuclei are resolved, therefore, by recording the chemical shifts of 15N or 13C′either before or after the NOE mixing period. Applications of these methods to the study ofthe unfolded state of the N-terminal SH3 domain of drk (drkN SH3) and a partially foldedlarge fragment of staphylococcal nuclease (SNase), Δ131Δ, are presented. Inaddition, an application to folded SNase in complex with the ligands thymidine3′,5′-bisphosphate (pdTp) and Ca2+ is illustrated which allows the assignmentof NOEs between degenerate Hα protons or protons resonating close to water.

Intrinsically disordered proteins play crucial roles in regulatory processes and often function as protein interaction hubs. Here, we present a detailed characterization of a full-length disordered hub protein region involved in multiple dynamic complexes. We performed NMR, CD, and fluorescence binding studies on the nonphosphorylated and highly PKA-phosphorylated human cystic fibrosis transmembrane conductance regulator (CFTR) regulatory region, a ∼200-residue disordered segment involved in phosphorylation-dependent regulation of channel trafficking and gating. Our data provide evidence for dynamic, phosphorylation-dependent, multisite interactions of various segments of the regulatory region for its intra- and intermolecular partners, including the CFTR nucleotide binding domains 1 and 2, a 42-residue peptide from the C terminus of CFTR, the SLC26A3 sulphate transporter and antisigma factor antagonist (STAS) domain, and 14-3-3β. Because of its large number of binding partners, multivalent binding of individually weak sites facilitates rapid exchange between free and bound states to allow the regulatory region to engage with different partners and generate a graded or rheostat-like response to phosphorylation. Our results enrich the understanding of how disordered binding segments interact with multiple targets. We present structural models consistent with our data that illustrate this dynamic aspect of phospho-regulation of CFTR by the disordered regulatory region.

Protein flexibility spans a broad spectrum, from highly stable folded to intrinsically disordered states. In this review, we discuss how various techniques, including X-ray crystallography, nuclear magnetic resonance spectroscopy and ensemble-modeling strategies employing various experimental measurements, have enabled detailed structural and dynamic characterizations of proteins in their free and bound states. This has revealed a variety of possible binding scenarios in which flexibility can either decrease or increase upon binding. Furthermore, dynamic free-state ensembles have repeatedly been observed to contain transiently formed conformations that partially or completely resemble bound states. These results demonstrate an intimate connection between protein flexibility and protein interactions and illustrate the huge diversity of structure and dynamics in both free proteins and protein complexes.

A heteronuclear correlation experiment is described which permits simultaneous characterization of both 15N longitudinal decay rates and slow conformational exchange rates. Data pertaining to the exchange between folded and unfolded forms of an SH3 domain is used to illustrate the technique. Because the unfolded form of the molecule, on average, shows significantly higher NH exchange rates than the folded form, and approach which minimizes the degree of water saturation is employed, enabling the extraction of accurate rate constants.

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.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a multi-domain membrane chloride channel whose activity is regulated by ATP at two nucleotide-binding domains (NBD1 and NBD2) and by phosphorylation of the regulatory (R) region. The NBDs and the R region have functionally relevant motions that are critical for channel gating. Nuclear magnetic resonance (NMR) spectroscopy is a highly useful technique for obtaining information on the structure and interactions of CFTR and is extremely powerful for probing dynamics. NMR approaches for studying CFTR are reviewed, using our previous NBD1 and the R region results to provide examples. These NMR data are yielding insights into the dynamic properties and interactions that facilitate normal CFTR regulation as well as pathological effects of mutations, including the most common disease mutant, deletion of F508 in NBD1.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a multi-domain membrane chloride channel whose activity is regulated by ATP at two nucleotide-binding domains (NBD1 and NBD2) and by phosphorylation of the regulatory (R) region. The NBDs and the R region have functionally relevant motions that are critical for channel gating. Nuclear magnetic resonance (NMR) spectroscopy is a highly useful technique for obtaining information on the structure and interactions of CFTR and is extremely powerful for probing dynamics. NMR approaches for studying CFTR are reviewed, using our previous NBD1 and the R region results to provide examples. These NMR data are yielding insights into the dynamic properties and interactions that facilitate normal CFTR regulation as well as pathological effects of mutations, including the most common disease mutant, deletion of F508 in NBD1.

Muscle relaxation is triggered by the dephosphorylation of Ser19 in the myosin regulatory light chain. This reaction is catalyzed by the holoenzyme myosin phosphatase (MP), which includes the catalytic subunit protein phosphatase 1 (PP1) and the regulatory targeting subunit (MYPT). MYPT1 (myosin phosphatase targeting subunit 1) is responsible for both targeting the holoenzyme to subcellular compartments in the muscle and directing PP1 specificity toward myosin. To understand the molecular events leading to the MYPT1-PP1 holoenzyme formation, we used NMR spectroscopy to determine the structural and dynamic characteristics of unbound MYPT1. This allowed the conformations of MYPT1 in the free, unbound state to be directly compared to the PP1-bound state. Our results show that MYPT1(1-98) behaves like a two-domain protein in solution. The first 40 residues of MYPT1(1-98), the disordered region, are intrinsically disordered and highly dynamic, whereas residues 41-98, the folded ankyrin-repeat region, are well-structured and rigid. Furthermore, the integrated use of NMR and biophysical data enabled us to calculate an ensemble model for MYPT1(1-98). The most prominent structural feature of the MYPT1(1-98) ensemble is a 25% populated transient α-helix in the disordered region of MYPT1(1-98). This α-helix becomes fully populated when bound to PP1 and, as we show, likely plays a central role in the formation of the MYPT1-PP1 holoenzyme complex. Finally, this combined analysis shows that the structural and dynamic behaviors exhibited by MYPT1 for PP1 are distinct from those of any other previously analyzed PP1 regulatory protein. Collectively, these data enable us to present a new model of the molecular events that drive MYPT1-PP1 holoenzyme formation and demonstrate that there are structural differences in unbound PP1 regulators that have not been previously observed. Thus, this work adds significant insights to the currently limited data for molecular structures and dynamics of PP1 regulators.

Intrinsically disordered proteins (IDPs), which lack folded structure and are disordered under nondenaturing conditions, have been shown to perform important functions in a large number of cellular processes. These proteins have interesting structural properties that deviate from the random-coil-like behavior exhibited by chemically denatured proteins. In particular, IDPs are often observed to exhibit significant compaction. In this study, we have analyzed the hydrodynamic radii of a number of IDPs to investigate the sequence determinants of this compaction. Net charge and proline content are observed to be strongly correlated with increased hydrodynamic radii, suggesting that these are the dominant contributors to compaction. Hydrophobicity and secondary structure, on the other hand, appear to have negligible effects on compaction, which implies that the determinants of structure in folded and intrinsically disordered proteins are profoundly different. Finally, we observe that polyhistidine tags seem to increase IDP compaction, which suggests that these tags have significant perturbing effects and thus should be removed before any structural characterizations of IDPs. Using the relationships observed in this analysis, we have developed a sequence-based predictor of hydrodynamic radius for IDPs that shows substantial improvement over a simple model based upon chain length alone.

Amide solvent exchange rates are regarded as a valuable source of information on structure/dynamics of unfolded (disordered) proteins. Proton-based saturation transfer experiments, normally used to measure solvent exchange, are known to meet some serious difficulties. The problems mainly arise from the need to (1) manipulate water magnetization and (2) discriminate between multiple magnetization transfer pathways that occur within the proton pool. Some of these issues are specific to unfolded proteins. For example, the compensation scheme used to cancel the Overhauser effect in the popular CLEANEX experiment is not designed for use with unfolded proteins. In this report we describe an alternative experimental strategy, where amide (15)N is used as a probe of solvent exchange. The experiment is performed in 50% H(2)O-50% D(2)O solvent and is based on the (HACACO)NH pulse sequence. The resulting spectral map is fully equivalent to the conventional HSQC. To fulfill its purpose, the experiment monitors the conversion of deuterated species, (15)N(D), into protonated species, (15)N(H), as effected by the solvent exchange. Conceptually, this experiment is similar to EXSY which prompted the name of (15)N(H/D)-SOLEXSY (SOLvent EXchange SpectroscopY). Of note, our experimental scheme, which relies on nitrogen rather than proton to monitor solvent exchange, is free of the complications described above. The developed pulse sequence was used to measure solvent exchange rates in the chemically denatured state of the drkN SH3 domain. The results were found to correlate well with the CLEANEX-PM data, r = 0.97, thus providing a measure of validation for both techniques. When the experimentally measured exchange rates are converted into protection factors, most of the values fall in the range 0.5-2, consistent with random-coil behavior. However, elevated values, ca. 5, are obtained for residues R38 and A39, as well as the side-chain indole of W36. This is surprising, given that high protection factors imply hydrogen bonding or hydrophobic burial not expected to occur in a chemically denatured state of a protein. We, therefore, hypothesized that elevated protection factors are an artefact arising from the calculation of the reference (random-coil) exchange rates. To confirm this hypothesis, we prepared samples of several short peptides derived from the sequence of the drkN SH3 domain; these samples were used to directly measure the reference exchange rates. The revised protection factors obtained in this manner proved to be close to 1.0. These results also have implications for the more compact unfolded state of drkN SH3, which appears to be fully permeable to water as well, with no manifestations of hydrophobic burial.

Recognition requires protein flexibility because it facilitates conformational rearrangements and induced-fit mechanisms upon target binding. Intrinsic disorder is an extreme on the continuous spectrum of possible protein dynamics and its role in recognition may seem counterintuitive. However, conformational disorder is widely found in many eukaryotic regulatory proteins involved in processes such as signal transduction and transcription. Disordered protein regions may in fact confer advantages over folded proteins in binding. Rapidly interconverting and diverse conformers may create mean electrostatic fields instead of presenting discrete charges. The resultant "polyelectrostatic" interactions allow for the utilization of post-translational modifications as a means to change the net charge and thereby modify the electrostatic interaction of a disordered region. Plasticity of disordered protein states enables steric advantages over folded proteins and allows for unique binding configurations. Disorder may also have evolutionary advantages, as it facilitates alternative splicing, domain shuffling and protein modularity. As proteins exist in a continuous spectrum of disorder, so do their complexes. Indeed, disordered regions in complexes may control the degree of motion between domains, mask binding sites, be targets of post-translational modifications, permit overlapping binding motifs, and enable transient binding of different binding partners, making them excellent candidates for signal integrators and explaining their prevalence in eukaryotic signaling pathways. "Dynamic" complexes arise if more than two transient protein interfaces are involved in complex formation of two binding partners in a dynamic equilibrium. "Disordered" complexes, in contrast, do not involve significant ordering of interacting protein segments but rely exclusively on transient contacts. The nature of these interactions is not well understood yet but advancements in the structural characterization of disordered states will help us gain insights into their function and their implications for health and disease.

The regulatory (R) region of the cystic fibrosis transmembrane conductance regulator (CFTR) is intrinsically disordered and must be phosphorylated at multiple sites for full CFTR channel activity, with no one specific phosphorylation site required. In addition, nucleotide binding and hydrolysis at the nucleotide-binding domains (NBDs) of CFTR are required for channel gating. We report NMR studies in the absence and presence of NBD1 that provide structural details for the isolated R region and its interaction with NBD1 at residue-level resolution. Several sites in the R region with measured fractional helical propensity mediate interactions with NBD1. Phosphorylation reduces the helicity of many R-region sites and reduces their NBD1 interactions. This evidence for a dynamic complex with NBD1 that transiently engages different sites of the R region suggests a structural explanation for the dependence of CFTR activity on multiple PKA phosphorylation sites.

The N-terminal SH3 domain of the Drosophila modular protein Drk undergoes slow exchange between a folded (Fexch) and highly populated unfolded (Uexch) state under nondenaturing buffer conditions, enabling both Fexch and Uexch states to be simultaneously monitored. The addition of dissolved oxygen, equilibrated to a partial pressure of either 30 atm or 60 atm, provides the means to study solvent exposure with atomic resolution via 13C NMR paramagnetic shifts in 1H,13C HSQC (heteronuclear single quantum coherence) spectra. Absolute differences in these paramagnetic shifts between the Fexch and Uexch states allow the discrimination of regions of the protein which undergo change in solvent exposure upon unfolding. Contact with dissolved oxygen for both the Fexch and Uexch states could also be assessed through 13C paramagnetic shifts which were normalized based on the corresponding paramagnetic shifts seen in the free amino acids. In the Fexch state, the 13C nuclei belonging to the hydrophobic core of the protein exhibited very weak normalized paramagnetic shifts while those with greater solvent accessible surface area exhibited significantly larger normalized shifts. The Uexch state displayed less varied 13C paramagnetic shifts although distinct regions of protection from solvent exposure could be identified by a lack of such shifts. These regions, which included Phe9, Thr12, Ala13, Lys21, Thr22, Ile24, Ile27, and Arg38, overlapped with those found to have residual nativelike and non-native structures in previous studies and in some cases provided novel information. Thus, the paramagnetic shifts from dissolved oxygen are highly useful in the study of a transient structure or clustering in disordered systems, where conventional NMR measurements (couplings, chemical shift deviations from random coil values, and NOEs) may give little information.

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.

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.

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.

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.

(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 (

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.

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.