Eukaryotic Linear Motifs Research Articles

Fast and scalable querying of eukaryotic linear motifs with gget elm.

Eukaryotic linear motifs (ELMs), or Short Linear Motifs, are protein interaction modules that play an essential role in cellular processes and signaling networks and are often involved in diseases like cancer. The ELM database is a collection of manually curated motif knowledge from scientific papers. It has become a crucial resource for investigating motif biology and recognizing candidate ELMs in novel amino acid sequences. Users can search amino acid sequences or UniProt Accessions on the ELM resource web interface. However, as with many web services, there are limitations in the swift processing of large-scale queries through the ELM web interface or API calls, and, therefore, integration into protein function analysis pipelines is limited. To allow swift, large-scale motif analyses on protein sequences using ELMs curated in the ELM database, we have extended the gget suite of Python and command line tools with a new module, gget elm, which does not rely on the ELM server for efficiently finding candidate ELMs in user-submitted amino acid sequences and UniProt Accessions. gget elm increases accessibility to the information stored in the ELM database and allows scalable searches for motif-mediated interaction sites in the amino acid sequences. The manual and source code are available at https://github.com/pachterlab/gget.

Novel Insights into Phytoplasma Effectors

Effectoromics has become integral to the identification of pathogen targets and/or host-resistant proteins for the genetic improvement of plants in agriculture and horticulture. Phytoplasmas are the causal agents of more than 100 plant diseases in economically important crops such as vegetables, spices, medicinal plants, ornamentals, palms, fruit trees, etc. To date, around 20 effectors in phytoplasmas have been experimentally validated but the list of putative effectors comprises hundreds of different proteins. Very few families (tribes) have been identified based on homology, such as the SAP05-like, SAP11-like, SAP54-like and TENGU-like families. The lack of conservation in amino acid sequences slows the progress of effectoromics in phytoplasmas since many effectors must be studied individually. Here, 717 phytoplasma effector candidates and 21 validated effectors were characterized in silico to identify common features. We identified functional domains in 153 effectors, while 585 had no known domains. The most frequently identified domain was the sequence-variable mosaic domain (SVM domain), widely distributed in 87 phytoplasma effectors. Searching for de novo amino acid motifs, 50 were found in the phytoplasma effector dataset; 696 amino acid sequences of effectors had at least 1 motif while 42 had no motif at all. These data allowed us to organize effectors into 15 tribes, uncovering, for the first time, evolutionary relationships largely masked by lack of sequence conservation among effectors. We also identified 42 eukaryotic linear motifs (ELMs) in phytoplasma effector sequences. Since the motifs are related to common functions, this novel organization of phytoplasma effectors may help further advance effectoromics research to combat phytoplasma infection in agriculture and horticulture.

Implications of intrinsic disorder and functional proteomics in the merkel cell polyomavirus life cycle.

Infection with merkel cell polyomavirus (MCPyV) is implicated in the development of merkel cell carcinoma (MCC), a rare but aggressive skin cancer. MCC has a mortality rate near 50%, and incidence has been rapidly increasing in recent decades, making development of improved treatment strategies critical to addressing its growing social burden. The parallel increasing necessity for novel research to better understand MCPyV pathogenesis has prompted numerous studies in recent years, yet the role of intrinsic disorder in MCPyV proteins remains unexplored. This study carries out computational characterization of intrinsic disorder within the MCPyV proteome and suggests mechanisms that may contribute to the oncogenicity of the virus to invade and hijack host immune systems. Our analysis finds that significant levels of intrinsic disorder are present in proteins LT, ALTO, 57kT, and VP1, and suggests that regions of sT may also contain large, disordered regions. The investigation further shows correlation of disorder propensity with the outputs for functional predictors of eukaryotic linear motifs (ELMs), molecular recognition features (MoRFs), and propensity for liquid-liquid phase separation (LLPS). Our findings indicate that MCPyV may use disorder and phase condensation to alter viral function that may accentuate or provide the basis for oncogenic activities. It is intended that this study will inform future experimental validation efforts around the phase separation capacity of MCPyV and its host protein-protein interactions. Furthermore, we hope to inform other investigators on the potential role of disorder in the MCPyV life cycle toward ultimately progressing the development of novel therapeutic agents.

Uncovering potential host proteins and pathways that may interact with eukaryotic short linear motifs in viral proteins of MERS, SARS and SARS2 coronaviruses that infect humans.

A coronavirus pandemic caused by a novel coronavirus (SARS-CoV-2) has spread rapidly worldwide since December 2019. Improved understanding and new strategies to cope with novel coronaviruses are urgently needed. Viruses (especially RNA viruses) encode a limited number and size (length of polypeptide chain) of viral proteins and must interact with the host cell components to control (hijack) the host cell machinery. To achieve this goal, the extensive mimicry of SLiMs in host proteins provides an effective strategy. However, little is known regarding SLiMs in coronavirus proteins and their potential targets in host cells. The objective of this study is to uncover SLiMs in coronavirus proteins that are present within host cells. These SLiMs have a high possibility of interacting with host intracellular proteins and hijacking the host cell machinery for virus replication and dissemination. In total, 1,479 SLiM hits were identified in the 16 proteins of 590 coronaviruses infecting humans. Overall, 106 host proteins were identified that may interact with SLiMs in 16 coronavirus proteins. These SLiM-interacting proteins are composed of many intracellular key regulators, such as receptors, transcription factors and kinases, and may have important contributions to virus replication, immune evasion and viral pathogenesis. A total of 209 pathways containing proteins that may interact with SLiMs in coronavirus proteins were identified. This study uncovers potential mechanisms by which coronaviruses hijack the host cell machinery. These results provide potential therapeutic targets for viral infections.

Intrinsically Disordered Regions/Proteins Compensate For Genomic Economization In Mycobacterium Tuberculosis

Central to the emergence of Mycobacterium tuberculosis (M. tb) as an extremely successful pathogen is its predilection to persist within the hostile host environment. In this context, the role of disordered proteins and linear motif rich intrinsically disordered proteins (IDPs) to perturb and hijack host cell networks for a productive infection holds special importance. M.tb has a relatively high percentage of IDPs in its proteome among prokaryotes, the significance of which is not known. Moreover, M.tb has evolved through reductive evolution, with the exception of PE‐PPE/PGRS protein family, which has been expanding in pathogenic mycobacteria. Here, we analyzed the M.tb proteome and its entire secretome. PE‐PGRS subfamily and the secretome were enriched in disorder and disordered binding sites, implying that structural disorder might be important for host‐pathogen interactions. We studied the PPE37 protein in more detail that carries an N‐terminal PPE domain capable of iron binding followed by two transmembrane domains and a disordered C‐terminal segment harboring number of eukaryotic linear motifs (ELMs). N‐terminal segment caused proliferation and differentiation of monocytic THP‐1 cells into CD11c, DC‐SIGN positive semi‐mature dendritic cells. The C‐terminal segment localizes to macrophage nucleus and induced caspase‐3 dependent apoptosis of host cells. We also studied the structural shift in PE35‐PPE68 and PE32‐PPE65 protein pairs of pathogenesis‐related RD1 and RD8 regions by analyzing their secondary structure through Fourier transformed infrared spectroscopy. These disordered proteins displayed a remarkable structural shift from disorder to order while engaging in the formation of complexes. These proteins are immunogenic individually and enhance pro‐pathogen response; however, their corresponding complexes enhanced the responses manifold in immunized mice. These results support the role of IDPs in performing contrasting activities to modulate the host processes possibly through molecular mimicry and cross talk. The results from our study further reinforce structural dynamics as a strategy to subvert host defence and to compensate genomic content loss due to reductional evolution.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Analysis of coevolution in nonstructural proteins of chikungunya virus.

BackgroundRNA viruses are characterized by high rate of mutations mainly due to the lack of proofreading repair activities associated with its RNA-dependent RNA-polymerase (RdRp). In case of arboviruses, this phenomenon has lead to the existence of mixed population of genomic variants within the host called quasi-species. The stability of strains within the quasi-species lies on mutations that are positively selected which in turn depend on whether these mutations are beneficial in either or both hosts. Coevolution of amino acids (aa) is one phenomenon that leads to establishment of favorable traits in viruses and leading to their fitness.ResultsFourteen CHIKV clinical samples collected over three years were subjected to RT-PCR, the four non-structural genes amplified and subjected to various genetic analyses. Coevolution analysis showed 30 aa pairs coevolving in nsP1, 23 aa pairs coevolving in nsP2, 239 in nsP3 and 46 aa coevolving pairs in nsP4 when each non-structural protein was considered independently. Further analysis showed that 705 amino acids pairs of the non-structural polyproteins coevolved together with a correlation coefficient of ≥0.5. Functional relevance of these coevolving amino acids in all the nonstructural proteins of CHIKV were predicted using Eukaryotic Linear Motifs (ELMs) of human.ConclusionsThe present study was undertaken to study co-evolving amino acids in the non-structural proteins of chikungunya virus (CHIKV), an important arbovirus. It was observed that several amino acids residues were coevolving and shared common functions.

Evolution of bopA Gene in Burkholderia: A Case of Convergent Evolution as a Mechanism for Bacterial Autophagy Evasion.

Autophagy is an important defense mechanism targeting intracellular bacteria to restrict their survival and growth. On the other hand, several intracellular pathogens have developed an antiautophagy mechanism to facilitate their own replication or intracellular survival. Up to now, no information about the origin or evolution of the antiautophagic genes in bacteria is available. BopA is an effector protein secreted by Burkholderia pseudomallei via the type three secretion system, and it has been shown to play a pivotal role in their escape from autophagy. The evolutionary origin of bopA was examined in this work. Sequence similarity searches for BopA showed that no homolog of BopA was detected in eukaryotes. However, eukaryotic linear motifs were detected in BopA. The phylogenetic tree of the BopA proteins in our analysis is congruent with the species phylogeny derived from housekeeping genes. Moreover, there was no obvious difference in GC content values of bopA gene and their respective genomes. Integrated information on the taxonomic distribution, phylogenetic relationships, and GC content of the bopA gene of Burkholderia revealed that this gene was acquired via convergent evolution, not from eukaryotic host through horizontal gene transfer (HGT) event. This work has, for the first time, characterized the evolutionary mechanism of bacterial evasion of autophagy. The results of this study clearly demonstrated the role of convergent evolution in the evolution of how bacteria evade autophagy.

Molecular principles of human virus protein-protein interactions.

Viruses, from the human protein-protein interaction network perspective, target hubs, bottlenecks and interconnected nodes enriched in certain biological pathways. However, not much is known about the general characteristic features of the human proteins interacting with viral proteins (referred to as hVIPs) as well as the motifs and domains utilized by human-virus protein-protein interactions (referred to as Hu-Vir PPIs). Our study has revealed that hVIPs are mostly disordered proteins, whereas viral proteins are mostly ordered proteins. Protein disorder in viral proteins and hVIPs varies from one subcellular location to another. In any given viral-human PPI pair, at least one of the two proteins is structurally disordered suggesting that disorder associated conformational flexibility as one of the characteristic features of virus-host interaction. Further analyses reveal that hVIPs are (i) slowly evolving proteins, (ii) associated with high centrality scores in human-PPI network, (iii) involved in multiple pathways, (iv) enriched in eukaryotic linear motifs (ELMs) associated with protein modification, degradation and regulatory processes, (v) associated with high number of splice variants and (vi) expressed abundantly across multiple tissues. These aforementioned findings suggest that conformational flexibility, spatial diversity, abundance and slow evolution are the characteristic features of the human proteins targeted by viral proteins. Hu-Vir PPIs are mostly mediated via domain-motif interactions (DMIs) where viral proteins employ motifs that mimic host ELMs to bind to domains in human proteins. DMIs are shared among viruses belonging to different families indicating a possible convergent evolution of these motifs to help viruses to adopt common strategies to subvert host cellular pathways. Hu-Vir PPI data, DDI and DMI data for human-virus PPI can be downloaded from http://cdfd.org.in/labpages/computational_biology_datasets.html. Supplementary data are available at Bioinformatics online.

Use of Host-like Peptide Motifs in Viral Proteins Is a Prevalent Strategy in Host-Virus Interactions

Viruses interact extensively with host proteins, but the mechanisms controlling these interactions are not well understood. We present a comprehensive analysis of eukaryotic linear motifs (ELMs) in 2,208 viral genomes and reveal that viruses exploit molecular mimicry of host-like ELMs to possibly assist in host-virus interactions. Using a statistical genomics approach, we identify a large number of potentially functional ELMs and observe that the occurrence of ELMs is often evolutionarily conserved but not uniform across virus families. Some viral proteins contain multiple types of ELMs, in striking similarity to complex regulatory modules in host proteins, suggesting that ELMs may act combinatorially to assist viral replication. Furthermore, a simple evolutionary model suggests that the inherent structural simplicity of ELMs often enables them to tolerate mutations and evolve quickly. Our findings suggest that ELMs may allow fast rewiring of host-virus interactions, which likely assists rapid viral evolution and adaptation to diverse environments.

Ordered Disorder of the Astrocytic Dystrophin-Associated Protein Complex in the Norm and Pathology

The abundance and potential functional roles of intrinsically disordered regions in aquaporin-4, Kir4.1, a dystrophin isoforms Dp71, α-1 syntrophin, and α-dystrobrevin; i.e., proteins constituting the functional core of the astrocytic dystrophin-associated protein complex (DAPC), are analyzed by a wealth of computational tools. The correlation between protein intrinsic disorder, single nucleotide polymorphisms (SNPs) and protein function is also studied together with the peculiarities of structural and functional conservation of these proteins. Our study revealed that the DAPC members are typical hybrid proteins that contain both ordered and intrinsically disordered regions. Both ordered and disordered regions are important for the stabilization of this complex. Many disordered binding regions of these five proteins are highly conserved among vertebrates. Conserved eukaryotic linear motifs and molecular recognition features found in the disordered regions of five protein constituting DAPC likely enhance protein-protein interactions that are required for the cellular functions of this complex. Curiously, the disorder-based binding regions are rarely affected by SNPs suggesting that these regions are crucial for the biological functions of their corresponding proteins.

Intrinsic Disorder in PTEN and its Interactome Confers Structural Plasticity and Functional Versatility

IDPs, while structurally poor, are functionally rich by virtue of their flexibility and modularity. However, how mutations in IDPs elicit diseases, remain elusive. Herein, we have identified tumor suppressor PTEN as an intrinsically disordered protein (IDP) and elucidated the molecular principles by which its intrinsically disordered region (IDR) at the carboxyl-terminus (C-tail) executes its functions. Post-translational modifications, conserved eukaryotic linear motifs and molecular recognition features present in the C-tail IDR enhance PTEN's protein-protein interactions that are required for its myriad cellular functions. PTEN primary and secondary interactomes are also enriched in IDPs, most being cancer related, revealing that PTEN functions emanate from and are nucleated by the C-tail IDR, which form pliable network-hubs. Together, PTEN higher order functional networks operate via multiple IDP-IDP interactions facilitated by its C-tail IDR. Targeting PTEN IDR and its interaction hubs emerges as a new paradigm for treatment of PTEN related pathologies.

SiteBinder: An Improved Approach for Comparing Multiple Protein Structural Motifs

There is a paramount need to develop new techniques and tools that will extract as much information as possible from the ever growing repository of protein 3D structures. We report here on the development of a software tool for the multiple superimposition of large sets of protein structural motifs. Our superimposition methodology performs a systematic search for the atom pairing that provides the best fit. During this search, the RMSD values for all chemically relevant pairings are calculated by quaternion algebra. The number of evaluated pairings is markedly decreased by using PDB annotations for atoms. This approach guarantees that the best fit will be found and can be applied even when sequence similarity is low or does not exist at all. We have implemented this methodology in the Web application SiteBinder, which is able to process up to thousands of protein structural motifs in a very short time, and which provides an intuitive and user-friendly interface. Our benchmarking analysis has shown the robustness, efficiency, and versatility of our methodology and its implementation by the successful superimposition of 1000 experimentally determined structures for each of 32 eukaryotic linear motifs. We also demonstrate the applicability of SiteBinder using three case studies. We first compared the structures of 61 PA-IIL sugar binding sites containing nine different sugars, and we found that the sugar binding sites of PA-IIL and its mutants have a conserved structure despite their binding different sugars. We then superimposed over 300 zinc finger central motifs and revealed that the molecular structure in the vicinity of the Zn atom is highly conserved. Finally, we superimposed 12 BH3 domains from pro-apoptotic proteins. Our findings come to support the hypothesis that there is a structural basis for the functional segregation of BH3-only proteins into activators and enablers.

Attributes of short linear motifs

Traditionally, protein-protein interactions were thought to be mediated by large, structured domains. However, it has become clear that the interactome comprises a wide range of binding interfaces with varying degrees of flexibility, ranging from rigid globular domains to disordered regions that natively lack structure. Enrichment for disorder in highly connected hub proteins and its correlation with organism complexity hint at the functional importance of disordered regions. Nevertheless, they have not yet been extensively characterised. Shifting the attention from globular domains to disordered regions of the proteome might bring us closer to elucidating the dense and complex connectivity of the interactome. An important class of disordered interfaces are the compact mono-partite, short linear motifs (SLiMs, or eukaryotic linear motifs (ELMs)). They are evolutionarily plastic and interact with relatively low affinity due to the limited number of residues that make direct contact with the binding partner. These features confer to SLiMs the ability to evolve convergently and mediate transient interactions, which is imperative to network evolution and to maintain robust cell signalling, respectively. The ability to discriminate biologically relevant SLiMs by means of different attributes will improve our understanding of the complexity of the interactome and aid development of bioinformatics tools for motif discovery. In this paper, the curated instances currently available in the Eukaryotic Linear Motif (ELM) database are analysed to provide a clear overview of the defining attributes of SLiMs. These analyses suggest that functional SLiMs have higher levels of conservation than their surrounding residues, frequently evolve convergently, preferentially occur in disordered regions and often form a secondary structure when bound to their interaction partner. These results advocate searching for small groupings of residues in disordered regions with higher relative conservation and a propensity to form the secondary structure. Finally, the most interesting conclusions are examined in regard to their functional consequences.

Molecular Recognition Features Facilitate the Binding Diversity of Hub Proteins

Molecular recognition features (MoRFs) are short, unstructured regions that perform molecular binding and recognition functions via disorder-to-order transitions. These short disordered peptides have various names in the literature, such as ELMs (eukaryotic linear motifs), SLiMs (short linear motifs) and ANCHORS (predicted disordered binding sites based on energy calculation). The MoRF mechanisms are hypothesized to involve in the binding diversity of hub proteins in two ways: one-to-many signaling in disordered hubs (e.g. p53) and many-to-one signaling in ordered hubs (e.g. 14-3-3). The former hub uses the same disordered region to form different secondary structures to bind four structurally different partners. The latter hub binds to five distinct disordered regions via one structured binding groove. Our work explored hundreds of hub protein examples by analyzing the complexes deposited in Protein Data Bank that have same binding MoRFs or similar binding grooves. Further computational experiments and analyses provided us more detailed and specific explanations regarding how disordered regions facilitate the binding diversity in different complex structures. Exploring more examples absolutely gave us much clearer picture to uncover the conformational changes that occur accompanying with binding, and showing that, in general, flexibility allows both subtle and complex structural variation thereby enabling the same sequence to morph into differently shaped binding sites and different sequences to fit into the same binding site. These results not only verify the two MoRF mechanisms we observed in the single examples of one-to-many and many-to-one signaling previously, but also generalize the specific examples to many hundreds of other well-characterized protein-protein interactions.

Exploring the Binding Diversity of Intrinsically Disordered Proteins

Intrinsically disorder proteins often perform their functions through the binding of a short, loosely structured region from one protein onto the binding site of a protein partner. These short binding regions go by various names including eukaryotic linear motifs (ELMs), short linear motifs (SLiMs) and molecular recognition features (MoRFs). All of these represent a class of disordered protein that executes molecular recognition and binding functions typically via a disorder-to-order structural transition. Previous studies from our group showed 2 distinct examples of hub proteins, a disordered hub example (p53) and an ordered hub example (14-3-3), performing one-to-many signaling and many-to-one signaling, respectively. In the former, a single disordered region binds to multiple structured partners. In the latter, many distinct disordered regions bind to a single structured partner. Both alternatives use the MoRF mechanism described above. In this study, we tried to expand our previous work to find more examples and to determine whether the prevalence of the MoRF mechanism is common or rare through analyzing protein complexes deposited in PDB that consist of short nonglobular fragments bound to large globular partners. After examining and clustering the various MoRFs, we found 298 disordered hub examples and 246 ordered hub examples. Further experiments are providing detailed information about how intrinsic disorder facilitates binding to diverse partners. Exploring these examples is yielding a much clearer picture of the conformational changes that occur upon binding and showing that, in general, flexibility allows both subtle and complex structural variation thereby enabling different sequences to fit into the same binding site and the same sequence to morph into differently shaped binding sites. These results indicate that the observations made previously on the single examples of many-to-one and one-to-many signaling generalize to many hundreds of other well characterized protein-protein interactions.

Host sequence motifs shared by HIV predict response to antiretroviral therapy

BackgroundThe HIV viral genome mutates at a high rate and poses a significant long term health risk even in the presence of combination antiretroviral therapy. Current methods for predicting a patient's response to therapy rely on site-directed mutagenesis experiments and in vitro resistance assays. In this bioinformatics study we treat response to antiretroviral therapy as a two-body problem: response to therapy is considered to be a function of both the host and pathogen proteomes. We set out to identify potential responders based on the presence or absence of host protein and DNA motifs on the HIV proteome.ResultsAn alignment of thousands of HIV-1 sequences attested to extensive variation in nucleotide sequence but also showed conservation of eukaryotic short linear motifs on the protein coding regions. The reduction in viral load of patients in the Stanford HIV Drug Resistance Database exhibited a bimodal distribution after 24 weeks of antiretroviral therapy, with 2,000 copies/ml cutoff. Similarly, patients allocated into responder/non-responder categories based on consistent viral load reduction during a 24 week period showed clear separation. In both cases of phenotype identification, a set of features composed of short linear motifs in the reverse transcriptase region of HIV sequence accurately predicted a patient's response to therapy. Motifs that overlap resistance sites were highly predictive of responder identification in single drug regimens but these features lost importance in defining responders in multi-drug therapies.ConclusionHIV sequence mutates in a way that preferentially preserves peptide sequence motifs that are also found in the human proteome. The presence and absence of such motifs at specific regions of the HIV sequence is highly predictive of response to therapy. Some of these predictive motifs overlap with known HIV-1 resistance sites. These motifs are well established in bioinformatics databases and hence do not require identification via in vitro mutation experiments.

Understanding eukaryotic linear motifs and their role in cell signaling and regulation

It is now clear that a detailed picture of cell regulation requires a comprehensive understanding of the abundant short protein motifs through which signaling is channeled. The current body of knowledge has slowly accumulated through piecemeal experimental investigation of individual motifs in signaling. Computational methods contributed little to this process. A new generation of bioinformatics tools will aid the future investigation of motifs in regulatory proteins, and the disordered polypeptide regions in which they frequently reside. Allied to high throughput methods such as phosphoproteomics, signaling networks are becoming amenable to experimental deconstruction. In this review, we summarise the current state of linear motif biology, which uses low affinity interactions to create cooperative, combinatorial and highly dynamic regulatory protein complexes. The discrete deterministic properties implicit to these assemblies suggest that models for cell regulatory networks in systems biology should neither be overly dependent on stochastic nor on smooth deterministic approximations.