Functional Diversification after Gene Duplication: Paralog Specific Regions of Structural Disorder and Phosphorylation in p53, p63, and p73

Conformational and functional flexibility promote protein evolvability. High evolvability allows related proteins to functionally diverge and perhaps to neostructuralize. p53 is a multifunctional protein frequently referred to as the Guardian of the Genome–a hub for e.g. incoming and outgoing signals in apoptosis and DNA repair. p53 has been found to be structurally disordered, an extreme form of conformational flexibility. Here, p53, and its paralogs p63 and p73, were studied for further insights into the evolutionary dynamics of structural disorder, secondary structure, and phosphorylation. This study is focused on the post gene duplication phase for the p53 family in vertebrates, but also visits the origin of the protein family and the early domain loss and gain events. Functional divergence, measured by rapid evolutionary dynamics of protein domains, structural properties, and phosphorylation propensity, is inferred across vertebrate p53 proteins, in p63 and p73 from fish, and between the three paralogs. In particular, structurally disordered regions are redistributed among paralogs, but within clades redistribution of structural disorder also appears to be an ongoing process. Despite its deemed importance as the Guardian of the Genome, p53 is indeed a protein with high evolvability as seen not only in rearranged structural disorder, but also in fluctuating domain sequence signatures among lineages.     

Protein intrinsic disorder and human papillomaviruses: increased amount of disorder in E6 and E7 oncoproteins from high risk HPVs

It is recognized now that many functional proteins or their long segments are devoid of stable secondary and/or tertiary structure and exist instead as very dynamic ensembles of conformations. They are known by different names including natively unfolded, intrinsically disordered, intrinsically unstructured, rheomorphic, pliable, and different combinations thereof. Many important functions and activities have been associated with these intrinsically disordered proteins (IDPs), including molecular recognition, signaling, and regulation. It is also believed that disorder of these proteins allows function to be readily modified through phosphorylation, acetylation, ubiquitination, hydroxylation, and proteolysis. Bioinformatics analysis revealed that IDPs comprise a large fraction of different proteomes. Furthermore, it is established that the intrinsic disorder is relatively abundant among cancer-related and other disease-related proteins and IDPs play a number of key roles in oncogenesis. There are more than 100 different types of human papillomaviruses (HPVs), which are the causative agents of benign papillomas/warts, and cofactors in the development of carcinomas of the genital tract, head and neck, and epidermis. With respect to their association with cancer, HPVs are grouped into two classes, known as low (e.g., HPV-6 and HPV-11) and high-risk (e.g., HPV-16 and HPV-18) types. The entire proteome of HPV includes six nonstructural proteins [E1, E2, E4, E5, E6, and E7 (the latter two are known to function as oncoproteins in the high-risk HPVs)] and two structural proteins (L1 and L2). To understand whether intrinsic disorder plays a role in the oncogenic potential of different HPV types, we have performed a detailed bioinformatics analysis of proteomes of high-risk and low-risk HPVs with the major focus on E6 and E7 oncoproteins. The results of this analysis are consistent with the conclusion that high-risk HPVs are characterized by the increased amount of intrinsic disorder in transforming proteins E6 and E7.     

Flexible nets. The roles of intrinsic disorder in protein interaction networks

Proteins participate in complex sets of interactions that represent the mechanistic foundation for much of the physiology and function of the cell. These protein-protein interactions are organized into exquisitely complex networks. The architecture of protein-protein interaction networks was recently proposed to be scale-free, with most of the proteins having only one or two connections but with relatively fewer 'hubs' possessing tens, hundreds or more links. The high level of hub connectivity must somehow be reflected in protein structure. What structural quality of hub proteins enables them to interact with large numbers of diverse targets? One possibility would be to employ binding regions that have the ability to bind multiple, structurally diverse partners. This trait can be imparted by the incorporation of intrinsic disorder in one or both partners. To illustrate the value of such contributions, this review examines the roles of intrinsic disorder in protein network architecture. We show that there are three general ways that intrinsic disorder can contribute: First, intrinsic disorder can serve as the structural basis for hub protein promiscuity; secondly, intrinsically disordered proteins can bind to structured hub proteins; and thirdly, intrinsic disorder can provide flexible linkers between functional domains with the linkers enabling mechanisms that facilitate binding diversity. An important research direction will be to determine what fraction of protein-protein interaction in regulatory networks relies on intrinsic disorder.     

Identification and sequence analysis of cDNAs encoding a 110-kilodalton actin filament-associated pp60src substrate.

Transformation of chicken embryo cells by oncogenic forms of pp60src (e.g., pp60v-src or pp60527F) is linked with a concomitant increase in the steady-state levels of tyrosine-phosphorylated cellular proteins. Activated forms of the Src protein-tyrosine kinase stably associate with tyrosine-phosphorylated proteins, including a protein of 110 kDa, pp110. Previous reports have established that stable complex formation between pp110 and pp60src requires the structural integrity of the Src SH2 and SH3 domains, whereas tyrosine phosphorylation of pp110 requires only the structural integrity of the SH3 domain. In normal chicken embryo cells, pp110 colocalizes with actin stress filaments, and in Src-transformed cells, pp110 is found associated with podosomes (rosettes). Here, we report the identification and characterization of cDNAs encoding pp110. The predicted open reading frame encodes a polypeptide of 635 amino acids which exhibits little sequence similarity with other protein sequences present in the available sequence data bases. Thus, pp110 is a distinctive cytoskeleton-associated protein. On the basis of its association with actin stress filaments, we propose the term AFAP-110, for actin filament-associated protein of 110 kDa. In vitro analysis of AFAP-110 binding to bacterium-encoded glutathione S-transferase (GST) fusion proteins revealed that AFAP-110 present in normal cell extracts binds efficiently to Src SH3/SH2-containing fusion proteins, less efficiently to Src SH3-containing proteins, and poorly to SH2-containing fusion proteins. In contrast, AFAP-110 in Src-transformed cell extracts bound to GST-SH3/SH2 and GST-SH2 fusion proteins. Analysis of AFAP-110 cDNA sequences revealed the presence of sequence motifs predicted to bind to SH2 and SH3 domains, respectively. We suggest that AFAP-110 may represent a cellular protein capable of interacting with SH3-containing proteins and, upon tyrosine phosphorylation, binds tightly to SH2-containing proteins, such as pp60src or pp59fyn. The potential roles of AFAP-110 as an SH3/SH2 cytoskeletal binding protein are discussed.     

Structural disorder in eukaryotes

Based on early bioinformatic studies on a handful of species, the frequency of structural disorder of proteins is generally thought to be much higher in eukaryotes than in prokaryotes. To refine this view, we present here a comparative prediction study and analysis of 194 fully described eukaryotic proteomes and 87 reference prokaryotes for structural disorder. We found that structural disorder does distinguish eukaryotes from prokaryotes, but its frequency spans a very wide range in the two superkingdoms that largely overlap. The number of disordered binding regions and different Pfam domain types also contribute to distinguish eukaryotes from prokaryotes. Unexpectedly, the highest levels--and highest variability--of predicted disorder is found in protists, i.e. single-celled eukaryotes, often surpassing more complex eukaryote organisms, plants and animals. This trend contrasts with that of the number of domain types, which increases rather monotonously toward more complex organisms. The level of structural disorder appears to be strongly correlated with lifestyle, because some obligate intracellular parasites and endosymbionts have the lowest levels, whereas host-changing parasites have the highest level of predicted disorder. We conclude that protists have been the evolutionary hot-bed of experimentation with structural disorder, in a period when structural disorder was actively invented and the major functional classes of disordered proteins established.     

Advantages of proteins being disordered

The past decade has witnessed great advances in our understanding of protein structure‐function relationships in terms of the ubiquitous existence of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). The structural disorder of IDPs/IDRs enables them to play essential functions that are complementary to those of ordered proteins. In addition, IDPs/IDRs are persistent in evolution. Therefore, they are expected to possess some advantages over ordered proteins. In this review, we summarize and survey nine possible advantages of IDPs/IDRs: economizing genome/protein resources, overcoming steric restrictions in binding, achieving high specificity with low affinity, increasing binding rate, facilitating posttranslational modifications, enabling flexible linkers, preventing aggregation, providing resistance to non‐native conditions, and allowing compatibility with more available sequences. Some potential advantages of IDPs/IDRs are not well understood and require both experimental and theoretical approaches to decipher. The connection with protein design is also briefly discussed.     

Conformational role for the C-terminal tail of the intrinsically disordered high mobility group A (HMGA) chromatin factors

The architectural factors HMGA are highly connected hubs in the chromatin network and affect key cellular functions. HMGA have a causal involvement in cancer development; in fact, truncated or chimeric HMGA forms, resulting from chromosomal rearrangements, lack the constitutively phosphorylated acidic C-terminal tail and display increased oncogenic potential, suggesting a functional role for this domain. HMGA belong to the intrinsically disordered protein category, and this prevents the use of classical approaches to obtain structural data. Therefore, we combined limited proteolysis, ion mobility separation-mass spectrometry (IMS-MS), and electrospray ionization-mass spectrometry (ESI-MS) to obtain structural information regarding full length and C-terminal truncated HMGA forms. Limited proteolysis indicates that HMGA acidic tail shields the inner portions of the protein. IMS-MS and ESI-MS show that HMGA proteins can assume a compact form and that the degree of compactness is dependent upon the presence of the acidic tail and its constitutive phosphorylations. Moreover, we demonstrate that C-terminal truncated forms and wild type proteins are post-translationally modified in a different manner. Therefore, we propose that the acidic tail and its phosphorylation could affect HMGA post-translational modification status and likely their activity. Finally, the mass spectrometry-based approach adopted here proves to be a valuable new tool to obtain structural data regarding intrinsically disordered proteins.     

Computational studies reveal phosphorylation-dependent changes in the unstructured R domain of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is a cAMP-dependent chloride channel that is mutated in cystic fibrosis, an inherited disease of high morbidity and mortality. The phosphorylation of its ∼ 200 amino acid R domain by protein kinase A is obligatory for channel gating under normal conditions. The R domain contains more than ten PKA phosphorylation sites. No individual site is essential but phosphorylation of increasing numbers of sites enables progressively greater channel activity. In spite of numerous studies of the role of the R domain in CFTR regulation, its mechanism of action remains largely unknown. This is because neither its structure nor its interactions with other parts of CFTR have been completely elucidated. Studies have shown that the R domain lacks well-defined secondary structural elements and is an intrinsically disordered region of the channel protein. Here, we have analyzed the disorder pattern and employed computational methods to explore low-energy conformations of the R domain. The specific disorder and secondary structure patterns detected suggest the presence of molecular recognition elements (MoREs) that may mediate phosphorylation-regulated intra- and inter-domain interactions. Simulations were performed to generate an ensemble of accessible R domain conformations. Although the calculated structures may represent more compact conformers than occur in vivo, their secondary structure propensities are consistent with predictions and published experimental data. Equilibrium simulations of a mimic of a phosphorylated R domain showed that it exhibited an increased radius of gyration. In one possible interpretation of these findings, by changing its size, the globally unstructured R domain may act as an entropic spring to perturb the packing of membrane-spanning sequences that constitute the ion permeability pathway and thereby activate channel gating.     

Predicted and measured disorder in peripherin/rds, a retinal tetraspanin

Vertebrate photoreceptor outer segment (OS) morphogenesis requires peripherin/rds (P/rds). We have characterized this protein's C-terminus and present evidence that suggests it is intrinsically disordered. We propose that structural flexibility may underlie the multifunctionality proposed for this domain previously. The extremely short C-termini present in other tetraspanin family members suggest that intrinsic disorder may also play a role for those proteins.     

Protein disorder in the centrosome correlates with complexity in cell types number

Here we study the properties and the evolution of proteins that constitute the Centrosome, the complex molecular assembly that regulates the division and differentiation of animal cells. We found that centrosomal proteins are predicted to be significantly enriched in disordered and coiled-coil regions, more phosphorylated and longer than control proteins of the same organism. Interestingly, the ratio of these properties in centrosomal and control proteins tends to increase with the number of cell-types. We reconstructed indels evolution, finding that indels significantly increase disorder in both centrosomal and control proteins, at a rate that is typically larger along branches associated with a large growth in cell-types number, and larger for centrosomal than for control proteins. Substitutions show a similar trend for coiled-coil, but they contribute less to the evolution of disorder. Our results suggest that the increase in cell-types number in animal evolution is correlated with the gain of disordered and coiled-coil regions in centrosomal proteins, establishing a connection between organism and molecular complexity. We argue that the structural plasticity conferred to the Centrosome by disordered regions and phosphorylation plays an important role in its mechanical properties and its regulation in space and time.     

Expanding the Proteome of an RNA Virus by Phosphorylation of an Intrinsically Disordered Viral Protein

The human proteome contains myriad intrinsically disordered proteins. Within intrinsically disordered proteins, polyproline-II motifs are often located near sites of phosphorylation. We have used an unconventional experimental paradigm to discover that phosphorylation by protein kinase A (PKA) occurs in the intrinsically disordered domain of hepatitis C virus non-structural protein 5A (NS5A) on Thr-2332 near one of its polyproline-II motifs. Phosphorylation shifts the conformational ensemble of the NS5A intrinsically disordered domain to a state that permits detection of the polyproline motif by using ¹⁵N-, ¹³C-based multidimensional NMR spectroscopy. PKA-dependent proline resonances were lost in the presence of the Src homology 3 domain of c-Src, consistent with formation of a complex. Changing Thr-2332 to alanine in hepatitis C virus genotype 1b reduced the steady-state level of RNA by 10-fold; this change was lethal for genotype 2a. The lethal phenotype could be rescued by changing Thr-2332 to glutamic acid, a phosphomimetic substitution. Immunofluorescence and transmission electron microscopy showed that the inability to produce Thr(P)-2332-NS5A caused loss of integrity of the virus-induced membranous web/replication organelle. An even more extreme phenotype was observed in the presence of small molecule inhibitors of PKA. We conclude that the PKA-phosphorylated form of NS5A exhibits unique structure and function relative to the unphosphorylated protein. We suggest that post-translational modification of viral proteins containing intrinsic disorder may be a general mechanism to expand the viral proteome without a corresponding expansion of the genome.     

Structural characterization of the N-terminal kinase-interacting domain of an Hsp90-cochaperone Cdc37 by CD and solution NMR spectroscopy

Cdc37 is a protein kinase-targeting molecular chaperone, which cooperates with Hsp90 to assist the folding, assembly and maturation of various signaling kinases. It consists of three distinct domains: the N-terminal, middle, and C-terminal domain. While the middle domain is an Hsp90-binding domain, the N-terminal domain is recognized as a kinase-interacting domain. The N-terminal domain contains a well-conserved Ser residue at position 13, and the phosphorylation at this site has been shown to be a prerequisite for the interaction between Cdc37 and signaling kinases. Although the phosphorylation of Ser13 might induce some conformational change in Cdc37 molecule, little is known about the structure of the N-terminal domain of Cdc37. We examined the structural and dynamic properties of several fragment proteins corresponding to the N-terminal region of Cdc37 by circular dichroism and solution NMR spectroscopy. We found that the N-terminal domain of Cdc37 exhibits highly dynamic structure, and it exists in the equilibrium between α-helical and more disordered structures. We also found that phosphorylation at Ser13 did not significantly change the overall structure of N-terminal fragment protein of Cdc37. The results suggested that more complicated mechanisms might be necessary to explain the phosphorylation-activated interaction of Cdc37 with various kinases.     

IAPs contain an evolutionarily conserved ubiquitin-binding domain that regulates NF-kappaB as well as cell survival and oncogenesis

The covalent attachment of ubiquitin to target proteins influences various cellular processes, including DNA repair, NF-kappaB signalling and cell survival. The most common mode of regulation by ubiquitin-conjugation involves specialized ubiquitin-binding proteins that bind to ubiquitylated proteins and link them to downstream biochemical processes. Unravelling how the ubiquitin-message is recognized is essential because aberrant ubiquitin-mediated signalling contributes to tumour formation. Recent evidence indicates that inhibitor of apoptosis (IAP) proteins are frequently overexpressed in cancer and their expression level is implicated in contributing to tumorigenesis, chemoresistance, disease progression and poor patient-survival. Here, we have identified an evolutionarily conserved ubiquitin-associated (UBA) domain in IAPs, which enables them to bind to Lys 63-linked polyubiquitin. We found that the UBA domain is essential for the oncogenic potential of cIAP1, to maintain endothelial cell survival and to protect cells from TNF-alpha-induced apoptosis. Moreover, the UBA domain is required for XIAP and cIAP2-MALT1 to activate NF-kappaB. Our data suggest that the UBA domain of cIAP2-MALT1 stimulates NF-kappaB signalling by binding to polyubiquitylated NEMO. Significantly, 98% of all cIAP2-MALT1 fusion proteins retain the UBA domain, suggesting that ubiquitin-binding contributes to the oncogenic potential of cIAP2-MALT1 in MALT lymphoma. Our data identify IAPs as ubiquitin-binding proteins that contribute to ubiquitin-mediated cell survival, NF-kappaB signalling and oncogenesis.     

Partner choice during meiosis is regulated by Hop1-promoted dimerization of Mek1

Meiotic recombination differs from mitotic recombination in that DSBs are repaired using homologous chromosomes, rather than sister chromatids. This change in partner choice is due in part to a barrier to sister chromatid repair (BSCR) created by the meiosis-specific kinase, Mek1, in a complex with two other meiosis-specific proteins, Hop1 and Red1. HOP1 contains two functional domains, called the N and C domains. Analysis of a point mutation that specifically inactivates the C domain (hop1-K593A) reveals that the N domain is sufficient for Hop1 localization to chromosomes and for Red1 and Hop1 interactions. The C domain is needed for spore viability, for chromosome synapsis, and for preventing DMC1-independent DSB repair, indicating it plays a role in the BSCR. All of the hop1-K593A phenotypes can be bypassed by fusion of ectopic dimerization domains to Mek1, suggesting that the function of the C domain is to promote Mek1 dimerization. Hop1 is a DSB-dependent phosphoprotein, whose phosphorylation requires the presence of the C domain, but is independent of MEK1. These results suggest a model in which Hop1 phosphorylation in response to DSBs triggers dimerization of Mek1 via the Hop1 C domain, thereby enabling Mek1 to phosphorylate target proteins that prevent repair of DSBs by sister chromatids.     

Phosphoproteomic analysis of the mouse brain cytosol reveals a predominance of protein phosphorylation in regions of intrinsic sequence disorder

We analyzed the mouse forebrain cytosolic phosphoproteome using sequential (protein and peptide) IMAC purifications, enzymatic dephosphorylation, and targeted tandem mass spectrometry analysis strategies. In total, using complementary phosphoenrichment and LC-MS/MS strategies, 512 phosphorylation sites on 540 non-redundant phosphopeptides from 162 cytosolic phosphoproteins were characterized. Analysis of protein domains and amino acid sequence composition of this data set of cytosolic phosphoproteins revealed that it is significantly enriched in intrinsic sequence disorder, and this enrichment is associated with both cellular location and phosphorylation status. The majority of phosphorylation sites found by MS were located outside of structural protein domains (97%) but were mostly located in regions of intrinsic sequence disorder (86%). 368 phosphorylation sites were located in long regions of disorder (over 40 amino acids long), and 94% of proteins contained at least one such long region of disorder. In addition, we found that 58 phosphorylation sites in this data set occur in 14-3-3 binding consensus motifs, linear motifs that are associated with unstructured regions in proteins. These results demonstrate that in this data set protein phosphorylation is significantly depleted in protein domains and significantly enriched in disordered protein sequences and that enrichment of intrinsic sequence disorder may be a common feature of phosphoproteomes. This supports the hypothesis that disordered regions in proteins allow kinases, phosphatases, and phosphorylation-dependent binding proteins to gain access to target sequences to regulate local protein conformation and activity.     

GlobPlot: Exploring protein sequences for globularity and disorder

A major challenge in the proteomics and structural genomics era is to predict protein structure and function, including identification of those proteins that are partially or wholly unstructured. Non-globular sequence segments often contain short linear peptide motifs (e.g. SH3-binding sites) which are important for protein function. We present here a new tool for discovery of such unstructured, or disordered regions within proteins. GlobPlot (http://globplot.embl.de) is a web service that allows the user to plot the tendency within the query protein for order/globularity and disorder. We show examples with known proteins where it successfully identifies inter-domain segments containing linear motifs, and also apparently ordered regions that do not contain any recognised domain. GlobPlot may be useful in domain hunting efforts. The plots indicate that instances of known domains may often contain additional N- or C-terminal segments that appear ordered. Thus GlobPlot may be of use in the design of constructs corresponding to globular proteins, as needed for many biochemical studies, particularly structural biology. GlobPlot has a pipeline interface--GlobPipe--for the advanced user to do whole proteome analysis. GlobPlot can also be used as a generic infrastructure package for graphical displaying of any possible propensity.     

Intrinsic disorder in cell-signaling and cancer-associated proteins

The number of intrinsically disordered proteins known to be involved in cell-signaling and regulation is growing rapidly. To test for a generalized involvement of intrinsic disorder in signaling and cancer, we applied a neural network predictor of natural disordered regions (PONDR VL-XT) to four protein datasets: human cancer-associated proteins (HCAP), signaling proteins (AfCS), eukaryotic proteins from SWISS-PROT (EU_SW) and non-homologous protein segments with well-defined (ordered) 3D structure (O_PDB_S25). PONDR VL-XT predicts >or=30 consecutive disordered residues for 79(+/-5)%, 66(+/-6)%, 47(+/-4)% and 13(+/-4)% of the proteins from HCAP, AfCS, EU_SW, and O_PDB_S25, respectively, indicating significantly more intrinsic disorder in cancer-associated and signaling proteins as compared to the two control sets. The disorder analysis was extended to 11 additional functionally diverse categories of human proteins from SWISS-PROT. The proteins involved in metabolism, biosynthesis, and degradation together with kinases, inhibitors, transport, G-protein coupled receptors, and membrane proteins are predicted to have at least twofold less disorder than regulatory, cancer-associated and cytoskeletal proteins. In contrast to 44.5% of the proteins from representative non-membrane categories, just 17.3% of the cancer-associated proteins had sequence alignments with structures in the Protein Data Bank covering at least 75% of their lengths. This relative lack of structural information correlated with the greater amount of predicted disorder in the HCAP dataset. A comparison of disorder predictions with the experimental structural data for a subset of the HCAP proteins indicated good agreement between prediction and observation. Our data suggest that intrinsically unstructured proteins play key roles in cell-signaling, regulation and cancer, where coupled folding and binding is a common mechanism.