Sequence- and interactome-based prediction of viral protein hotspots targeting host proteins: a case study for HIV Nef

Virus proteins alter protein pathways of the host toward the synthesis of viral particles by breaking and making edges via binding to host proteins. In this study, we developed a computational approach to predict viral sequence hotspots for binding to host proteins based on sequences of viral and host proteins and literature-curated virus-host protein interactome data. We use a motif discovery algorithm repeatedly on collections of sequences of viral proteins and immediate binding partners of their host targets and choose only those motifs that are conserved on viral sequences and highly statistically enriched among binding partners of virus protein targeted host proteins. Our results match experimental data on binding sites of Nef to host proteins such as MAPK1, VAV1, LCK, HCK, HLA-A, CD4, FYN, and GNB2L1 with high statistical significance but is a poor predictor of Nef binding sites on highly flexible, hoop-like regions. Predicted hotspots recapture CD8 cell epitopes of HIV Nef highlighting their importance in modulating virus-host interactions. Host proteins potentially targeted or outcompeted by Nef appear crowding the T cell receptor, natural killer cell mediated cytotoxicity, and neurotrophin signaling pathways. Scanning of HIV Nef motifs on multiple alignments of hepatitis C protein NS5A produces results consistent with literature, indicating the potential value of the hotspot discovery in advancing our understanding of virus-host crosstalk.     

Searching for protein signatures using a multilevel alphabet

Short motifs are known to play diverse roles in proteins, such as in mediating the interactions with other molecules, binding to membranes, or conducting a specific biological function. Standard approaches currently employed to detect short motifs in proteins search for enrichment of amino acid motifs considering mostly the sequence information. Here, we presented a new approach to search for common motifs (protein signatures) which share both physicochemical and structural properties, looking simultaneously at different features. Our method takes as an input an amino acid sequence and translates it to a new alphabet that reflects its intrinsic structural and chemical properties. Using the MEME search algorithm, we identified the proteins signatures within subsets of protein which encompass common sequence and structural information. We demonstrated that we can detect enriched structural motifs, such as the amphipathic helix, from large datasets of linear sequences, as well as predicting common structural properties (such as disorder, surface accessibility, or secondary structures) of known functional‐motifs. Finally, we applied the method to the yeast protein interactome and identified novel putative interacting motifs. We propose that our approach can be applied for de novo protein function prediction given either sequence or structural information. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Sequence Alignment Reveals Possible MAPK Docking Motifs on HIV Proteins

Over the course of HIV infection, virus replication is facilitated by the phosphorylation of HIV proteins by human ERK1 and ERK2 mitogen-activated protein kinases (MAPKs). MAPKs are known to phosphorylate their substrates by first binding with them at a docking site. Docking site interactions could be viable drug targets because the sequences guiding them are more specific than phosphorylation consensus sites. In this study we use multiple bioinformatics tools to discover candidate MAPK docking site motifs on HIV proteins known to be phosphorylated by MAPKs, and we discuss the possibility of targeting docking sites with drugs. Using sequence alignments of HIV proteins of different subtypes, we show that MAPK docking patterns previously described for human proteins appear on the HIV matrix, Tat, and Vif proteins in a strain dependent manner, but are absent from HIV Rev and appear on all HIV Nef strains. We revise the regular expressions of previously annotated MAPK docking patterns in order to provide a subtype independent motif that annotates all HIV proteins. One revision is based on a documented human variant of one of the substrate docking motifs, and the other reduces the number of required basic amino acids in the standard docking motifs from two to one. The proposed patterns are shown to be consistent with in silico docking between ERK1 and the HIV matrix protein. The motif usage on HIV proteins is sufficiently different from human proteins in amino acid sequence similarity to allow for HIV specific targeting using small-molecule drugs.     

Characterization of protein hubs by inferring interacting motifs from protein interactions

The characterization of protein interactions is essential for understanding biological systems. While genome-scale methods are available for identifying interacting proteins, they do not pinpoint the interacting motifs (e.g., a domain, sequence segments, a binding site, or a set of residues). Here, we develop and apply a method for delineating the interacting motifs of hub proteins (i.e., highly connected proteins). The method relies on the observation that proteins with common interaction partners tend to interact with these partners through a common interacting motif. The sole input for the method are binary protein interactions; neither sequence nor structure information is needed. The approach is evaluated by comparing the inferred interacting motifs with domain families defined for 368 proteins in the Structural Classification of Proteins (SCOP). The positive predictive value of the method for detecting proteins with common SCOP families is 75% at sensitivity of 10%. Most of the inferred interacting motifs were significantly associated with sequence patterns, which could be responsible for the common interactions. We find that yeast hubs with multiple interacting motifs are more likely to be essential than hubs with one or two interacting motifs, thus rationalizing the previously observed correlation between essentiality and the number of interacting partners of a protein. We also find that yeast hubs with multiple interacting motifs evolve slower than the average protein, contrary to the hubs with one or two interacting motifs. The proposed method will help us discover unknown interacting motifs and provide biological insights about protein hubs and their roles in interaction networks.     

Interactome-Wide Prediction of Protein-Protein Binding Sites Reveals Effects of Protein Sequence Variation in Arabidopsis thaliana

The specificity of protein-protein interactions is encoded in those parts of the sequence that compose the binding interface. Therefore, understanding how changes in protein sequence influence interaction specificity, and possibly the phenotype, requires knowing the location of binding sites in those sequences. However, large-scale detection of protein interfaces remains a challenge. Here, we present a sequence- and interactome-based approach to mine interaction motifs from the recently published Arabidopsis thaliana interactome. The resultant proteome-wide predictions are available via www.ab.wur.nl/sliderbio and set the stage for further investigations of protein-protein binding sites. To assess our method, we first show that, by using a priori information calculated from protein sequences, such as evolutionary conservation and residue surface accessibility, we improve the performance of interface prediction compared to using only interactome data. Next, we present evidence for the functional importance of the predicted sites, which are under stronger selective pressure than the rest of protein sequence. We also observe a tendency for compensatory mutations in the binding sites of interacting proteins. Subsequently, we interrogated the interactome data to formulate testable hypotheses for the molecular mechanisms underlying effects of protein sequence mutations. Examples include proteins relevant for various developmental processes. Finally, we observed, by analysing pairs of paralogs, a correlation between functional divergence and sequence divergence in interaction sites. This analysis suggests that large-scale prediction of binding sites can cast light on evolutionary processes that shape protein-protein interaction networks.     

Prediction of short linear protein binding regions

Short linear motifs in proteins (typically 3-12 residues in length) play key roles in protein-protein interactions by frequently binding specifically to peptide binding domains within interacting proteins. Their tendency to be found in disordered segments of proteins has meant that they have often been overlooked. Here we present SLiMPred (short linear motif predictor), the first general de novo method designed to computationally predict such regions in protein primary sequences independent of experimentally defined homologs and interactors. The method applies machine learning techniques to predict new motifs based on annotated instances from the Eukaryotic Linear Motif database, as well as structural, biophysical, and biochemical features derived from the protein primary sequence. We have integrated these data sources and benchmarked the predictive accuracy of the method, and found that it performs equivalently to a predictor of protein binding regions in disordered regions, in addition to having predictive power for other classes of motif sites such as polyproline II helix motifs and short linear motifs lying in ordered regions. It will be useful in predicting peptides involved in potential protein associations and will aid in the functional characterization of proteins, especially of proteins lacking experimental information on structures and interactions. We conclude that, despite the diversity of motif sequences and structures, SLiMPred is a valuable tool for prioritizing potential interaction motifs in proteins.     

The mechanism of action of the SSB interactome reveals it is the first OB‐fold family of genome guardians in prokaryotes

The single‐stranded DNA binding protein (SSB) is essential to all aspects of DNA metabolism in bacteria. This protein performs two distinct, but closely intertwined and indispensable functions in the cell. SSB binds to single‐stranded DNA (ssDNA) and at least 20 partner proteins resulting in their regulation. These partners comprise a family of genome guardians known as the SSB interactome. Essential to interactome regulation is the linker/OB‐fold network of interactions. This network of interactions forms when one or more PXXP motifs in the linker of SSB bind to an OB‐fold in a partner, with interactome members involved in competitive binding between the linker and ssDNA to their OB‐fold. Consequently, when linker‐binding occurs to an OB‐fold in an interactome partner, proteins are loaded onto the DNA. When linker/OB‐fold interactions occur between SSB tetramers, cooperative ssDNA‐binding results, producing a multi‐tetrameric complex that rapidly protects the ssDNA. Within this SSB‐ssDNA complex, there is an extensive and dynamic network of linker/OB‐fold interactions that involves multiple tetramers bound contiguously along the ssDNA lattice. The dynamic behavior of these tetramers which includes binding mode changes, sliding as well as DNA wrapping/unwrapping events, are likely coupled to the formation and disruption of linker/OB‐fold interactions. This behavior is essential to facilitating downstream DNA processing events. As OB‐folds are critical to the essence of the linker/OB‐fold network of interactions, and they are found in multiple interactome partners, the SSB interactome is classified as the first family of prokaryotic, oligosaccharide/oligonucleotide binding fold (OB‐fold) genome guardians.     

Comparative analysis of Leishmania exoproteomes: Implication for host–pathogen interactions

Leishmaniasis is a vector-borne disease caused by the protozoa Leishmania. We have analyzed and compared the sequences of three experimental exoproteomes of Leishmania promastigotes from different species to determine their specific features and to identify new candidate proteins involved in interactions of Leishmania with the host. The exoproteomes differ from the proteomes by a decrease in the average molecular weight per protein, in disordered amino acid residues and in basic proteins. The exoproteome of the visceral species is significantly enriched in sites predicted to be phosphorylated as well as in features frequently associated with molecular interactions (intrinsic disorder, number of disordered binding regions per protein, interaction and/or trafficking motifs) compared to the other species. The visceral species might thus have a larger interaction repertoire with the host than the other species. Less than 10% of the exoproteomes contain heparin-binding and RGD sequences, and ~ 30% the host targeting signal RXLXE/D/Q. These latter proteins might thus be exported inside the host cell during the intracellular stage of the infection. Furthermore we have identified nine protein families conserved in the three exoproteomes with specific combinations of Pfam domains and selected eleven proteins containing at least three interaction and/or trafficking motifs including two splicing factors, phosphomannomutase, 2,3-bisphosphoglycerate-independent phosphoglycerate mutase, the paraflagellar rod protein-1D and a putative helicase. Their role in host–Leishmania interactions warrants further investigation but the putative ATP-dependent DEAD/H RNA helicase, which contains numerous interaction motifs, a host targeting signal and two disordered regions, is a very promising candidate.     

Disorder and sequence repeats in hub proteins and their implications for network evolution

Protein interaction networks display approximate scale-free topology, in which hub proteins that interact with a large number of other proteins determine the overall organization of the network. In this study, we aim to determine whether hubs are distinguishable from other networked proteins by specific sequence features. Proteins of different connectednesses were compared in the interaction networks of Saccharomyces cerevisiae, Drosophila melanogaster, Caenorhabditis elegans, and Homo sapienswith respect to the distribution of predicted structural disorder, sequence repeats, low complexity regions, and chain length. Highly connected proteins ("hub proteins") contained significantly more of, and greater proportion of, these sequence features and tended to be longer overall as compared to less connected proteins. These sequence features provide two different functional means for realizing multiple interactions: (1) extended interaction surface and (2) flexibility and adaptability, providing a mechanism for the same region to bind distinct partners. Our view contradicts the prevailing view that scaling in protein interactomes arose from gene duplication and preferential attachment of equivalent proteins. We propose an alternative evolutionary network specialization process, in which certain components of the protein interactome improved their fitness for binding by becoming longer or accruing regions of disorder and/or internal repeats and have therefore become specialized in network organization.     

DYNLL/LC8: a light chain subunit of the dynein motor complex and beyond

The LC8 family members of dynein light chains (DYNLL1 and DYNLL2 in vertebrates) are highly conserved ubiquitous eukaryotic homodimer proteins that interact, besides dynein and myosin 5a motor proteins, with a large (and still incomplete) number of proteins involved in diverse biological functions. Despite an earlier suggestion that LC8 light chains function as cargo adapters of the above molecular motors, they are now recognized as regulatory hub proteins that interact with short linear motifs located in intrinsically disordered protein segments. The most prominent LC8 function is to promote dimerization of their binding partners that are often scaffold proteins of various complexes, including the intermediate chains of the dynein motor complex. Structural and functional aspects of this intriguing hub protein will be highlighted in this minireview.     

Insilco analysis of functionally important residues in folate receptors

Lack of crystal structure data of folate binding proteins has left so many questions unanswered (for example, important residues in active site, binding domain, important amino acid residues involved in interactions between ligand and receptor). With sequence alignment and PROSITE motif identification, we attempted to answer evolutionarily significant residues that are of functional importance for ligand binding and that form catalytic sites. We have analyzed 46 different FRs and FBP sequences of various organisms obtained from Genbank. Multiple sequence alignment identified 44 highly conserved identical amino acid residues with 10 cysteine residues and 12 motifs including ECSPNLGPW (which might help in the structural stability of FR).     

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.     

Structure–function relationships in HIV-1 Nef

The accessory Nef protein of HIV and SIV is essential for viral pathogenesis, yet it is perplexing in its multitude of molecular functions. In this review we analyse the structure–function relationships of motifs recently proposed to play roles in aspects of Nef modification, signalling and trafficking, and thereby to impinge on the ability of the virus to survive in, and to manipulate, its cellular host. Based on the full-length structure assembly of HIV Nef, we correlate surface accessibility with secondary structure elements and sequence conservation. Motifs involved in Nef-mediated CD4 and MHC I downregulation are located in flexible regions of Nef, suggesting that the formation of the transient trafficking complexes involved in these processes depends on the recognition of primary sequences. In contrast, the interaction sites for signalling molecules that contain SH3 domains or the p21-activated kinases are associated with the well folded core domain, suggesting the recognition of highly structured protein surfaces.     

Human phospho‐signaling networks of SARS‐CoV‐2 infection are rewired by population genetic variants

SARS‐CoV‐2 infection hijacks signaling pathways and induces protein–protein interactions between human and viral proteins. Human genetic variation may impact SARS‐CoV‐2 infection and COVID‐19 pathology; however, the genetic variation in these signaling networks remains uncharacterized. Here, we studied human missense single nucleotide variants (SNVs) altering phosphorylation sites modulated by SARS‐CoV‐2 infection, using machine learning to identify amino acid substitutions altering kinase‐bound sequence motifs. We found 2,033 infrequent phosphorylation‐associated SNVs (pSNVs) that are enriched in sequence motif alterations, potentially reflecting the evolution of signaling networks regulating host defenses. Proteins with pSNVs are involved in viral life cycle and host responses, including RNA splicing, interferon response (TRIM28), and glucose homeostasis (TBC1D4) with potential associations with COVID‐19 comorbidities. pSNVs disrupt CDK and MAPK substrate motifs and replace these with motifs of Tank Binding Kinase 1 (TBK1) involved in innate immune responses, indicating consistent rewiring of signaling networks. Several pSNVs associate with severe COVID‐19 and hospitalization (STARD13, ARFGEF2). Our analysis highlights potential genetic factors contributing to inter‐individual variation of SARS‐CoV‐2 infection and COVID‐19 and suggests leads for mechanistic and translational studies.Subject Categories: Computational Biology, Microbiology, Virology & Host Pathogen Interaction, Post-translational Modifications & Proteolysis

An integrative proteogenomic study reveals that human phospho‐signaling networks responding to SARS‐CoV‐2 infection are enriched in genetic variants that modify kinase binding motifs, suggesting that genetic variation may impact infection and COVID‐19 pathology.     

Composite structural motifs of binding sites for delineating biological functions of proteins

Most biological processes are described as a series of interactions between proteins and other molecules, and interactions are in turn described in terms of atomic structures. To annotate protein functions as sets of interaction states at atomic resolution, and thereby to better understand the relation between protein interactions and biological functions, we conducted exhaustive all-against-all atomic structure comparisons of all known binding sites for ligands including small molecules, proteins and nucleic acids, and identified recurring elementary motifs. By integrating the elementary motifs associated with each subunit, we defined composite motifs that represent context-dependent combinations of elementary motifs. It is demonstrated that function similarity can be better inferred from composite motif similarity compared to the similarity of protein sequences or of individual binding sites. By integrating the composite motifs associated with each protein function, we define meta-composite motifs each of which is regarded as a time-independent diagrammatic representation of a biological process. It is shown that meta-composite motifs provide richer annotations of biological processes than sequence clusters. The present results serve as a basis for bridging atomic structures to higher-order biological phenomena by classification and integration of binding site structures.     

Towards functional repertoire of the earliest proteins

The conserved protein sequence motifs present in all prokaryotic proteomes, “omnipresent motifs,” presumably, correspond to the earliest proteins of the Last Universal Cellular Ancestor, from which all the proteomes have descended. Fifteen proteomes, each representing one of the total 15 diverse phyla of 131 Eubacteria and Archea, from which the omnipresent elements have been originally derived, are exhaustively screened. All those proteins which harbor the omnipresent motifs are identified. Six “omnipresent” protein types are revealed which are located in all 15 proteomes: ABC cassettes, FtsH proteases, translation initiation factors, translation elongation factors, isoleucyl-tRNA synthases, and RNA polymerases β’. In addition to the omnipresent motifs, these proteins also contain other highly conserved motifs, standing for additional modules of the proteins. Remarkably, the identified tentative earliest proteins are responsible for only three basic functions: supply of monomers (ABC transporters and proteases), protein synthesis (initiation and elongation factors, aminoacyl-tRNA synthases), and RNA synthesis (polymerases). No enzymes involved in metabolic activities are present in the list of the earliest proteins derived by this approach. Some of the omnipresent sequence motifs are found, indeed, in the metabolic enzymes (e.g. NTP binding motifs), but these enzymes do not make a sequence matching collection of 15 sequences, i.e. they are not omnipresent. Future analysis of less conserved sequence motifs may reveal at what degree of conservation (stage of evolution) the metabolic enzymes could have entered the scene.     

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.     

Leucine-specific, functional interactions between human immunodeficiency virus type 1 Nef and adaptor protein complexes

The human immunodeficiency virus type 1 virulence protein Nef interacts with the endosomal sorting machinery via a leucine-based motif. Similar sequences within the cytoplasmic domains of cellular transmembrane proteins bind to the adaptor protein (AP) complexes of coated vesicles to modulate protein traffic, but the molecular basis of the interactions between these motifs and the heterotetrameric complexes is controversial. To identify the target of the Nef leucine motif, the native sequence was replaced with either leucine- or tyrosine-based AP-binding sequences from cellular proteins, and the interactions with AP subunits were correlated with function. Tyrosine motifs predictably modulated the interactions between Nef and the mu subunits of AP-1, AP-2, and AP-3; heterologous leucine motifs caused little change in these interactions. Conversely, leucine motifs mediated a ternary interaction between Nef and hemicomplexes containing the sigma1 plus gamma subunits of AP-1 or the sigma3 plus delta subunits of AP-3, whereas tyrosine motifs did not. Similarly, only leucine motifs supported the Nef-mediated association of AP-1 and AP-3 with endosomal membranes in cells treated with brefeldin A. Functionally, Nef proteins containing leucine motifs down-regulated CD4 from the cell surface and enhanced viral replication, whereas those containing tyrosine motifs were inactive. Apparently, the interaction of Nef with the mu subunits of AP complexes is insufficient for function. A leucine-specific mode of interaction that likely involves AP hemicomplexes is further required for Nef activity. The mu and hemicomplex interactions may cooperate to yield high avidity binding of AP complexes to Nef. This binding likely underlies the unusual ability of Nef to induce the stabilization of these complexes on endosomal membranes, an activity that correlates with enhancement of viral replication.     

Evaluating Caveolin Interactions: Do Proteins Interact with the Caveolin Scaffolding Domain through a Widespread Aromatic Residue-Rich Motif?

Caveolins are coat proteins of caveolae, small flask-shaped pits of the plasma membranes of most cells. Aside from roles in caveolae formation, caveolins recruit, retain and regulate many caveolae-associated signalling molecules. Caveolin-protein interactions are commonly considered to occur between a ∼20 amino acid region within caveolin, the caveolin scaffolding domain (CSD), and an aromatic-rich caveolin binding motif (CBM) on the binding partner (фXфXXXXф, фXXXXфXXф or фXфXXXXфXXф, where ф is an aromatic and X an unspecified amino acid). The CBM resembles a typical linear motif - a short, simple sequence independently evolved many times in different proteins for a specific function. Here we exploit recent improvements in bioinformatics tools and in our understanding of linear motifs to critically examine the role of CBMs in caveolin interactions. We find that sequences conforming to the CBM occur in 30% of human proteins, but find no evidence for their statistical enrichment in the caveolin interactome. Furthermore, sequence- and structure-based considerations suggest that CBMs do not have characteristics commonly associated with true interaction motifs. Analysis of the relative solvent accessible area of putative CBMs shows that the majority of their aromatic residues are buried within the protein and are thus unlikely to interact directly with caveolin, but may instead be important for protein structural stability. Together, these findings suggest that the canonical CBM may not be a common characteristic of caveolin-target interactions and that interfaces between caveolin and targets may be more structurally diverse than presently appreciated.