Targeting Intrinsically Disordered Transcription Factors: Changing the Paradigm.

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.

Characterization of the PilN, PilO and PilP type IVa pilus subcomplex.

Type IVa pili are bacterial nanomachines required for colonization of surfaces, but little is known about the organization of proteins in this system. The Pseudomonas aeruginosa pilMNOPQ operon encodes five key members of the transenvelope complex facilitating pilus function. While PilQ forms the outer membrane secretin pore, the functions of the inner membrane-associated proteins PilM/N/O/P are less well defined. Structural characterization of a stable C-terminal fragment of PilP (PilP(Δ71)) by NMR revealed a modified ß-sandwich fold, similar to that of Neisseria meningitidis PilP, although complementation experiments showed that the two proteins are not interchangeable likely due to divergent surface properties. PilP is an inner membrane putative lipoprotein, but mutagenesis of the putative lipobox had no effect on the localization and function of PilP. A larger fragment, PilP(Δ18-6His), co-purified with a PilN(Δ44)/PilO(Δ51) heterodimer as a stable complex that eluted from a size exclusion chromatography column as a single peak with a molecular weight equivalent to two heterotrimers with 1:1:1 stoichiometry. Although PilO forms both homodimers and PilN-PilO heterodimers, PilP(Δ18-6His) did not interact stably with PilO(Δ51) alone. Together these data demonstrate that PilN/PilO/PilP interact directly to form a stable heterotrimeric complex, explaining the dispensability of PilP's lipid anchor for localization and function.

Slow dynamics in folded and unfolded states of an SH3 domain.

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

Structural characterization of proteins with an attached ATCUN motif by paramagnetic relaxation enhancement NMR spectroscopy.

The use of a short, three-residue Cu(2+)-binding sequence, the ATCUN motif, is presented as an approach for extracting long-range distance restraints from relaxation enhancement NMR spectroscopy. The ATCUN motif is prepended to the N-termini of proteins and binds Cu(2+) with a very high affinity. Relaxation rates of amide protons in ATCUN-tagged protein in the presence and absence of Cu(2+) can be converted into distance restraints and used for structure refinement by using a new routine, PMAG, that has been written for the structure calculation program CNS. The utility of the approach is demonstrated with an application to ATCUN-tagged ubiquitin. Excellent agreement between measured relaxation rates and those calculated on the basis of the X-ray structure of the protein have been obtained.

Calculation of ensembles of structures representing the unfolded state of an SH3 domain.

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 (<60) dominate the final ensembles compared to the total number in the starting ensembles, suggesting that the drkN SH3 domain U(exch) state can be described by a limited number of lower-energy conformations. The calculated U(exch) state ensembles are much more compact than a "random coil" chain, with significant native-like residual structure observed. In particular, a sizable population of conformers having the n-src loop and distal beta-hairpin structures exist in the calculated U(exch) state ensembles, and Trp36 is involved in a large number of interactions, both native and non-native.

Solution structure of a Nedd4 WW domain-ENaC peptide complex.

Nedd4 is a ubiquitin protein ligase composed of a C2 domain, three (or four) WW domains and a ubiquitin ligase Hect domain. Nedd4 was demonstrated to bind the epithelial sodium channel (alphabetagammaENaC), by association of its WW domains with PY motifs (XPPXY) present in each ENaC subunit, and to regulate the cell surface stability of the channel. The PY motif of betaENaC is deleted or mutated in Liddle syndrome, a hereditary form of hypertension caused by elevated ENaC activity. Here we report the solution structure of the third WW domain of Nedd4 complexed to the PY motif-containing region of betaENaC (TLPIPGTPPPNYDSL, referred to as betaP2). A polyproline type II helical conformation is adopted by the PPPN sequence. Unexpectedly, the C-terminal sequence YDSL forms a helical turn and both the tyrosine and the C-terminal leucine contact the WW domain. This is unlike other proline-rich peptides complexed to WW domains, which bind in an extended conformation and lack molecular interactions with residues C-terminal to the tyrosine or the structurally equivalent residue in non-PY motif WW domain targets. The Nedd4 WW domain-ENaC betaP2 peptide structure expands our understanding of the mechanisms involved in WW domain-ligand recognition and the molecular basis of Liddle syndrome.

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states.

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.

Multidimensional NMR methods for protein structure determination.

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.

Global folds of proteins with low densities of NOEs using residual dipolar couplings: application to the 370-residue maltodextrin-binding protein.

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.

Sequential assignment of proline-rich regions in proteins: application to modular binding domain complexes.

Many protein-protein interactions involve amino acid sequences containing proline-rich motifs and even polyproline stretches. The lack of amide protons in such regions complicates assignment, since 1HN-based triple-resonance assignment strategies cannot be employed. Two such systems that we are currently studying include an SH2 domain from the protein Crk with a region containing 9 prolines in a 14 amino acid sequence, as well as a WW domain that interacts with a proline-rich target. A modified version of the HACAN pulse scheme, originally described by Bax and co-workers [Wang et al. (1995) J. Biomol. NMR, 5, 376-382], and an experiment which correlates the intra-residue 1Halpha, 13Calpha/13Cbeta chemical shifts with the 15N shift of the subsequent residue are presented and applied to the two systems listed above, allowing sequential assignment of the molecules.

Multiple modes of peptide recognition by the PTB domain of the cell fate determinant Numb.

The phosphotyrosine-binding (PTB) domain of the cell fate determinant Numb is involved in the formation of multiple protein complexes in vivo and can bind a diverse array of peptide sequences in vitro. To investigate the structural basis for the promiscuous nature of this protein module, we have determined its solution structure by NMR in a complex with a peptide containing an NMSF sequence derived from the Numb-associated kinase (Nak). The Nak peptide was found to adopt a significantly different structure from that of a GPpY sequence-containing peptide previously determined. In contrast to the helical turn adopted by the GPpY peptide, the Nak peptide forms a beta-turn at the NMSF site followed by another turn near the C-terminus. The Numb PTB domain appears to recognize peptides that differ in both primary and secondary structures by engaging various amounts of the binding surface of the protein. Our results suggest a mechanism through which a single PTB domain might interact with multiple distinct target proteins to control a complex biological process such as asymmetric cell division.

Diversity in protein recognition by PTB domains.

Phosphotyrosine-binding (PTB) domains were originally identified as modular domains that recognize phosphorylated Asn-Pro-Xxx-p Tyr-containing proteins. Recent binding and structural studies of PTB domain complexes with target peptides have revealed a number of deviations from the previously described mode of interaction, with respect to both the sequences of possible targets and their structures within the complexes. This diversity of recognition by PTB domains extends and strengthens our general understanding of modular binding domain recognition.

The 'dynamics' in the thermodynamics of binding.

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.

Novel mode of ligand binding by the SH2 domain of the human XLP disease gene product SAP/SH2D1A.

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.

A simple in vivo assay for increased protein solubility.

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.

Analysis of deuterium relaxation-derived methyl axis order parameters and correlation with local structure.

Methyl axis (S2axis) and backbone NH (S2NH) order parameters derived from eight proteins have been analyzed. Similar distribution profiles for Ala S2axis and S2NH order parameters were observed. A good correlation between the two S2axis values of Val and Leu methyl groups is noted, although differences between order parameters can arise. The relation of S2axis or S2NH to solvent accessibility and packing density has also been investigated. Correlations are weak, likely reflecting the importance of collective, non-local motions in proteins. The lack of correlation between these simple structural parameters and dynamics emphasizes the importance of motional studies to fully characterize proteins.

High populations of non-native structures in the denatured state are compatible with the formation of the native folded state.

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.

Structure of a Numb PTB domain-peptide complex suggests a basis for diverse binding specificity.

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.

Correlation between binding and dynamics at SH2 domain interfaces.

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.

NMR studies of tandem WW domains of Nedd4 in complex with a PY motif-containing region of the epithelial sodium channel.

Nedd4 (neuronal precursor cell-expressed developmentally down-regulated 4) is a ubiquitin-protein ligase containing multiple WW domains. We have previously demonstrated the association between the WW domains of Nedd4 and PPxY (PY) motifs of the epithelial sodium channel (ENaC). In this paper, we report the assignment of backbone 1H alpha, 1HN, 15N, 13C', 13C alpha, and aliphatic 13C resonances of a fragment of rat Nedd4 (rNedd4) containing the two C-terminal WW domains, WW(II+III), complexed to a PY motif-containing peptide derived from the beta subunit of rat ENaC, the betaP2 peptide. The secondary structures of these two WW domains, determined from chemical shifts of 13C alpha and 13C beta resonances, are virtually identical to those of the WW domains of the Yes-associated protein YAP65 and the peptidyl-prolyl isomerase Pin1. Triple resonance experiments that detect the 1H alpha chemical shift were necessary to complete the chemical shift assignment, owing to the large number of proline residues in this fragment of rNedd4. A new experiment, which correlates sequential residues via their 15N nuclei and also detects 1H alpha chemical shifts, is introduced and its utility for the chemical shift assignment of sequential proline residues is discussed. Data collected on the WW(II+III)-betaP2 complex indicate that these WW domains have different affinities for the betaP2 peptide.