SLiMPrints: conservation-based discovery of functional motif fingerprints in intrinsically disordered protein regions

Large portions of higher eukaryotic proteomes are intrinsically disordered, and abundant evidence suggests that these unstructured regions of proteins are rich in regulatory interaction interfaces. A major class of disordered interaction interfaces are the compact and degenerate modules known as short linear motifs (SLiMs). As a result of the difficulties associated with the experimental identification and validation of SLiMs, our understanding of these modules is limited, advocating the use of computational methods to focus experimental discovery. This article evaluates the use of evolutionary conservation as a discriminatory technique for motif discovery. A statistical framework is introduced to assess the significance of relatively conserved residues, quantifying the likelihood a residue will have a particular level of conservation given the conservation of the surrounding residues. The framework is expanded to assess the significance of groupings of conserved residues, a metric that forms the basis of SLiMPrints (short linear motif fingerprints), a de novo motif discovery tool. SLiMPrints identifies relatively overconstrained proximal groupings of residues within intrinsically disordered regions, indicative of putatively functional motifs. Finally, the human proteome is analysed to create a set of highly conserved putative motif instances, including a novel site on translation initiation factor eIF2A that may regulate translation through binding of eIF4E.     

Interactome-wide prediction of short, disordered protein interaction motifs in humans

Many of the specific functions of intrinsically disordered protein segments are mediated by Short Linear Motifs (SLiMs) interacting with other proteins. Well known examples include SLiMs that interact with 14-3-3, PDZ, SH2, SH3, and WW domains but the true extent and diversity of SLiM-mediated interactions is largely unknown. Here, we attempt to expand our knowledge of human SLiMs by applying in silico SLiM prediction to the human interactome. Combining data from seven different interaction databases, we analysed approximately 6000 protein-centred and 1600 domain-centred human interaction datasets of 3+ unrelated proteins that interact with a common partner. Results were placed in context through comparison to randomised datasets of similar size and composition. The search returned thousands of evolutionarily conserved, intrinsically disordered occurrences of hundreds of significantly enriched recurring motifs, including many that have never been previously identified (). In addition to True Positive results for at least 25 different known SLiMs, a striking number of "off-target" proteins/domains also returned significantly enriched known motifs. Often, this was due to the non-independence of the datasets, with many proteins sharing interaction partners or contributing interactions to multiple domain datasets. The majority of these motif classes, however, were also found to be significantly enriched in one or more randomised datasets. This highlights the need for care when interpreting motif predictions of this nature but also raises the possibility that SLiM occurrences may be successfully identified independently of interaction data. Although not as compositionally biased as previous studies, patterns matching known SLiMs tended to cluster into a few large groups of similar sequence, while novel predictions tended to be more distinctive and less abundant. Whether this is due to ascertainment bias or a true functional composition bias of SLiMs is not clear and warrants further investigation.     

A structure filter for the Eukaryotic Linear Motif Resource

Background  

Many proteins are highly modular, being assembled from globular domains and segments of natively disordered polypeptides. Linear motifs, short sequence modules functioning independently of protein tertiary structure, are most abundant in natively disordered polypeptides but are also found in accessible parts of globular domains, such as exposed loops. The prediction of novel occurrences of known linear motifs attempts the difficult task of distinguishing functional matches from stochastically occurring non-functional matches. Although functionality can only be confirmed experimentally, confidence in a putative motif is increased if a motif exhibits attributes associated with functional instances such as occurrence in the correct taxonomic range, cellular compartment, conservation in homologues and accessibility to interacting partners. Several tools now use these attributes to classify putative motifs based on confidence of functionality.     

A comparative study of short linear motif compositions of the influenza A virus ribonucleoproteins

Protein-protein interactions through short linear motifs (SLiMs) are an emerging concept that is different from interactions between globular domains. The SLiMs encode a functional interaction interface in a short (three to ten residues) poorly conserved sequence. This characteristic makes them much more likely to arise/disappear spontaneously via mutations, and they may be more evolutionarily labile than globular domains. The diversity of SLiM composition may provide functional diversity for a viral protein from different viral strains. This study is designed to determine the different SLiM compositions of ribonucleoproteins (RNPs) from influenza A viruses (IAVs) from different hosts and with different levels of virulence. The 96 consensus sequences (regular expressions) of SLiMs from the ELM server were used to conduct a comprehensive analysis of the 52,513 IAV RNP sequences. The SLiM compositions of RNPs from IAVs from different hosts and with different levels of virulence were compared. The SLiM compositions of 845 RNPs from highly virulent/pandemic IAVs were also analyzed. In total, 292 highly conserved SLiMs were found in RNPs regardless of the IAV host range. These SLiMs may be basic motifs that are essential for the normal functions of RNPs. Moreover, several SLiMs that are rare in seasonal IAV RNPs but are present in RNPs from highly virulent/pandemic IAVs were identified.The SLiMs identified in this study provide a useful resource for experimental virologists to study the interactions between IAV RNPs and host intracellular proteins. Moreover, the SLiM compositions of IAV RNPs also provide insights into signal transduction pathways and protein interaction networks with which IAV RNPs might be involved. Information about SLiMs might be useful for the development of anti-IAV drugs.     

Evolution of Linear Motifs within the Papillomavirus E7 Oncoprotein

Many protein functions can be traced to linear sequence motifs of less than five residues, which are often found within intrinsically disordered domains. In spite of their prevalence, their role in protein evolution is only beginning to be understood. The study of papillomaviruses has provided many insights on the evolution of protein structure and function. We have chosen the papillomavirus E7 oncoprotein as a model system for the evolution of functional linear motifs. The multiple functions of E7 proteins from paradigmatic papillomavirus types can be explained to a large extent in terms of five linear motifs within the intrinsically disordered N-terminal domain and two linear motifs within the globular homodimeric C-terminal domain. We examined the motif inventory of E7 proteins from over 200 known papillomavirus types and found that the motifs reported for paradigmatic papillomavirus types are absent from many uncharacterized E7 proteins. Several motif pairs occur more often than expected, suggesting that linear motifs may evolve and function in a cooperative manner. The E7 linear motifs have appeared or disappeared multiple times during papillomavirus evolution, confirming the evolutionary plasticity of short functional sequences. Four of the motifs appeared several times during papillomavirus evolution, providing direct evidence for convergent evolution. Interestingly, the evolution pattern of a motif is independent of its location in a globular or disordered domain. The correlation between the presence of some motifs and virus host specificity and tissue tropism suggests that linear motifs play a role in the adaptive evolution of papillomaviruses.     

DisCons: a novel tool to quantify and classify evolutionary conservation of intrinsic protein disorder

Background

Analyzing the amino acid sequence of an intrinsically disordered protein (IDP) in an evolutionary context can yield novel insights on the functional role of disordered regions and sequence element(s). However, in the case of many IDPs, the lack of evolutionary conservation of the primary sequence can hamper the study of functionality, because the conservation of their disorder profile and ensuing function(s) may not appear in a traditional analysis of the evolutionary history of the protein.

Results

Here we present DisCons (Disorder Conservation), a novel pipelined tool that combines the quantification of sequence- and disorder conservation to classify disordered residue positions. According to this scheme, the most interesting categories (for functional purposes) are constrained disordered residues and flexible disordered residues. The former residues show conservation of both the sequence and the property of disorder and are associated mainly with specific binding functionalities (e.g., short, linear motifs, SLiMs), whereas the latter class correspond to segments where disorder as a feature is important for function as opposed to the identity of the underlying sequence (e.g., entropic chains and linkers). DisCons therefore helps with elucidating the function(s) arising from the disordered state by analyzing individual proteins as well as large-scale proteomics datasets.

Conclusions

DisCons is an openly accessible sequence analysis tool that identifies and highlights structurally disordered segments of proteins where the conformational flexibility is conserved across homologs, and therefore potentially functional. The tool is freely available both as a web application and as stand-alone source code hosted at http://pedb.vib.be/discons.     

Equivalent binding sites reveal convergently evolved interaction motifs

MOTIVATION: Much research has been devoted to the characterization of interaction interfaces found in complexes with known structure. In this context, the interactions of non-homologous domains at equivalent binding sites are of particular interest, as they can reveal convergently evolved interface motifs. Such motifs are an important source of information to formulate rules for interaction specificity and to design ligands based on the common features shared among diverse partners. RESULTS: We develop a novel method to identify non-homologous structural domains which bind at equivalent sites when interacting with a common partner. We systematically apply this method to all pairs of interactions with known structure and derive a comprehensive database for these interactions. Of all non-homologous domains, which bind with a common interaction partner, 4.2% use the same interface of the common interaction partner (excluding immunoglobulins and proteases). This rises to 16% if immunoglobulin and proteases are included. We demonstrate two applications of our database: first, the systematic screening for viral protein interfaces, which can mimic native interfaces and thus interfere; and second, structural motifs in enzymes and its inhibitors. We highlight several cases of virus protein mimicry: viral M3 protein interferes with a chemokine dimer interface. The virus has evolved the motif SVSPLP, which mimics the native SSDTTP motif. A second example is the regulatory factor Nef in HIV which can mimic a kinase when interacting with SH3. Among others the virus has evolved the kinase's PxxP motif. Further, we elucidate motif resemblances in Baculovirus p35 and HIV capsid proteins. Finally, chymotrypsin is subject to scrutiny wrt. its structural similarity to subtilisin and wrt. its inhibitor's similar recognition sites. SUPPLEMENTARY INFORMATION: A database is online at scoppi.biotec.tu-dresden.de/abac/.     

The MDMX Acidic Domain Uses Allovalency to Bind Both p53 and MDMX

Autoinhibition of p53 binding to MDMX requires two short-linear motifs (SLiMs) containing adjacent tryptophan (WW) and tryptophan-phenylalanine (WF) residues. NMR spectroscopy was used to show the WW and WF motifs directly compete for the p53 binding site on MDMX and circular dichroism spectroscopy was used to show the WW motif becomes helical when it is bound to the p53 binding domain (p53BD) of MDMX. Binding studies using isothermal titration calorimetry showed the WW motif is a stronger inhibitor of p53 binding than the WF motif when they are both tethered to p53BD by the natural disordered linker. We also investigated how the WW and WF motifs interact with the DNA binding domain (DBD) of p53. Both motifs bind independently to similar sites on DBD that overlap the DNA binding site. Taken together our work defines a model for complex formation between MDMX and p53 where a pair of disordered SLiMs bind overlapping sites on both proteins.     

Prediction of Protein Binding Regions in Disordered Proteins

Many disordered proteins function via binding to a structured partner and undergo a disorder-to-order transition. The coupled folding and binding can confer several functional advantages such as the precise control of binding specificity without increased affinity. Additionally, the inherent flexibility allows the binding site to adopt various conformations and to bind to multiple partners. These features explain the prevalence of such binding elements in signaling and regulatory processes. In this work, we report ANCHOR, a method for the prediction of disordered binding regions. ANCHOR relies on the pairwise energy estimation approach that is the basis of IUPred, a previous general disorder prediction method. In order to predict disordered binding regions, we seek to identify segments that are in disordered regions, cannot form enough favorable intrachain interactions to fold on their own, and are likely to gain stabilizing energy by interacting with a globular protein partner. The performance of ANCHOR was found to be largely independent from the amino acid composition and adopted secondary structure. Longer binding sites generally were predicted to be segmented, in agreement with available experimentally characterized examples. Scanning several hundred proteomes showed that the occurrence of disordered binding sites increased with the complexity of the organisms even compared to disordered regions in general. Furthermore, the length distribution of binding sites was different from disordered protein regions in general and was dominated by shorter segments. These results underline the importance of disordered proteins and protein segments in establishing new binding regions. Due to their specific biophysical properties, disordered binding sites generally carry a robust sequence signal, and this signal is efficiently captured by our method. Through its generality, ANCHOR opens new ways to study the essential functional sites of disordered proteins.     

Marked Variability in the Extent of Protein Disorder within and between Viral Families

Intrinsically disordered regions in eukaryotic proteomes contain key signaling and regulatory modules and mediate interactions with many proteins. Many viral proteomes encode disordered proteins and modulate host factors through the use of short linear motifs (SLiMs) embedded within disordered regions. However, the degree of viral protein disorder across different viruses is not well understood, so we set out to establish the constraints acting on viruses, in terms of their use of disordered protein regions. We surveyed predicted disorder across 2,278 available viral genomes in 41 families, and correlated the extent of disorder with genome size and other factors. Protein disorder varies strikingly between viral families (from 2.9% to 23.1% of residues), and also within families. However, this substantial variation did not follow the established trend among their hosts, with increasing disorder seen across eubacterial, archaebacterial, protists, and multicellular eukaryotes. For example, among large mammalian viruses, poxviruses and herpesviruses showed markedly differing disorder (5.6% and 17.9%, respectively). Viral families with smaller genome sizes have more disorder within each of five main viral types (ssDNA, dsDNA, ssRNA+, dsRNA, retroviruses), except for negative single-stranded RNA viruses, where disorder increased with genome size. However, surveying over all viruses, which compares tiny and enormous viruses over a much bigger range of genome sizes, there is no strong association of genome size with protein disorder. We conclude that there is extensive variation in the disorder content of viral proteomes. While a proportion of this may relate to base composition, to extent of gene overlap, and to genome size within viral types, there remain important additional family and virus-specific effects. Differing disorder strategies are likely to impact on how different viruses modulate host factors, and on how rapidly viruses can evolve novel instances of SLiMs subverting host functions, such as innate and acquired immunity.     

The next wave of interactomics: Mapping the SLiM-based interactions of the intrinsically disordered proteome

Short linear motifs (SLiMs) are a unique and ubiquitous class of protein interaction modules that perform key regulatory functions and drive dynamic complex formation. For decades, interactions mediated by SLiMs have accumulated through detailed low-throughput experiments. Recent methodological advances have opened this previously underexplored area of the human interactome to high-throughput protein–protein interaction discovery. In this article, we discuss that SLiM-based interactions represent a significant blind spot in the current interactomics data, introduce the key methods that are illuminating the elusive SLiM-mediated interactome of the human cell on a large scale, and discuss the implications for the field.     

Close encounters of the third kind: disordered domains and the interactions of proteins

Protein–protein interactions are thought to be mediated by domains, which are autonomous folding units of proteins. Recently, a second type of interaction has been suggested, mediated by short segments termed linear motifs, which are related to recognition elements of intrinsically disordered regions. Here, we propose a third kind of protein–protein recognition mechanism, mediated by disordered regions longer than 20–30 residues. Bioinformatics predictions and well‐characterized examples, such as the kinase‐inhibitory domain of Cdk inhibitors and the Wiskott–Aldrich syndrome protein (WASP)‐homology domain 2 of actin‐binding proteins, show that these disordered regions conform to the definition of domains rather than motifs, i.e., they represent functional, evolutionary, and structural units. Their functions are distinct from those of short motifs and ordered domains, and establish a third kind of interaction principle. With these points, we argue that these long disordered regions should be recognized as a distinct class of biologically functional protein domains.     

Deciphering the cause of evolutionary variance within intrinsically disordered regions in human proteins

Why the intrinsically disordered regions evolve within human proteome has became an interesting question for a decade. Till date, it remains an unsolved yet an intriguing issue to investigate why some of the disordered regions evolve rapidly while the rest are highly conserved across mammalian species. Identifying the key biological factors, responsible for the variation in the conservation rate of different disordered regions within the human proteome, may revisit the above issue. We emphasized that among the other biological features (multifunctionality, gene essentiality, protein connectivity, number of unique domains, gene expression level and expression breadth) considered in our study, the number of unique protein domains acts as a strong determinant that negatively influences the conservation of disordered regions. In this context, we justified that proteins having a fewer types of domains preferably need to conserve their disordered regions to enhance their structural flexibility which in turn will facilitate their molecular interactions. In contrast, the selection pressure acting on the stretches of disordered regions is not so strong in the case of multi-domains proteins. Therefore, we reasoned that the presence of conserved disordered stretches may compensate the functions of multiple domains within a single domain protein. Interestingly, we noticed that the influence of the unique domain number and expression level acts differently on the evolution of disordered regions from that of well-structured ones.     

Disordered Binding Regions and Linear Motifs—Bridging the Gap between Two Models of Molecular Recognition

Intrinsically disordered proteins (IDPs) exist without the presence of a stable tertiary structure in isolation. These proteins are often involved in molecular recognition processes via their disordered binding regions that can recognize partner molecules by undergoing a coupled folding and binding process. The specific properties of disordered binding regions give way to specific, yet transient interactions that enable IDPs to play central roles in signaling pathways and act as hubs of protein interaction networks. An alternative model of protein-protein interactions with largely overlapping functional properties is offered by the concept of linear interaction motifs. This approach focuses on distilling a short consensus sequence pattern from proteins with a common interaction partner. These motifs often reside in disordered regions and are considered to mediate the interaction roughly independent from the rest of the protein. Although a connection between linear motifs and disordered binding regions has been established through common examples, the complementary nature of the two concepts has yet to be fully explored. In many cases the sequence based definition of linear motifs and the structural context based definition of disordered binding regions describe two aspects of the same phenomenon. To gain insight into the connection between the two models, prediction methods were utilized. We combined the regular expression based prediction of linear motifs with the disordered binding region prediction method ANCHOR, each specialized for either model to get the best of both worlds. The thorough analysis of the overlap of the two methods offers a bioinformatics tool for more efficient binding site prediction that can serve a wide range of practical implications. At the same time it can also shed light on the theoretical connection between the two co-existing interaction models.     

Structured proteins and proteins with internal disorder

The point of view that a uniquely folded protein tertiary structure is required for the protein functioning has been prevailing in the literature quite recently. However of lately it has been found that many proteins in a cell have no such structure in an isolated state, though they have a well-defined function in physiological conditions. These proteins were named as proteins with natural or internal disorder. The portion of disordered regions in such proteins may vary from a sequence of several amino acids to a completely disordered sequence containing from tens to hundreds of amino acids. The main difference of these proteins from the structured (globular) ones is that they have no unique tertiary structure in an isolated state and acquire it after interaction with their partners. Their conformation in such a complex depends on the interacting partner and not only on their own amino acid sequence, which is specific for structured (globular) proteins. The problem of structural and functional relations in the structured proteins and proteins with internal disorder is discussed in this review. The complexity of the problem and its potential solutions are illustrated by the example of elongation factors EFlA.     

Linear motifs: lost in (pre)translation

Pretranslational modification by alternative splicing, alternative promoter usage and RNA editing enables the production of multiple protein isoforms from a single gene. A large quantity of data now supports the notion that short linear motifs (SLiMs), which are protein interaction modules enriched within intrinsically disordered regions, are key for the functional diversification of these isoforms. The inclusion or removal of these SLiMs can switch the subcellular localisation of an isoform, promote cooperative associations, refine the affinity of an interaction, coordinate phase transitions within the cell, and even create isoforms of opposing function. This article discusses the novel functionality enabled by the addition or removal of SLiM-containing exons by pretranslational modifications, such as alternative splicing and alternative promoter usage, and how these alterations enable the creation and modulation of complex regulatory and signalling pathways.     

Structured proteins and proteins with intrinsic disorder

Until recently, the point of view that the unique tertiary structure is necessary for protein function has prevailed. However, recent data have demonstrated that many cell proteins do not possess such structure in isolation, although displaying a distinct function under physiological conditions. These proteins were named the naturally, or intrinsically, disordered proteins. The fraction of intrinsically disordered regions in such proteins may vary from several amino acid residues to a completely unordered sequence of several tens or even several hundreds of residues. The main distinction of these proteins from structured (globular) proteins is that they have no unique tertiary structure in isolation and acquire it only upon interaction with their partners. The conformation of these proteins in a complex is determined not only by their own amino acid sequence (as is typical of structured, or globular, proteins) but also by the interacting partner. This review discusses the structure-function relationships in structured and intrinsically disordered proteins. The intricateness of this problem and the possible ways to solve it are illustrated by the example of the EF1A elongation factor family.