Linear motif-mediated interactions have contributed to the evolution of modularity in complex protein interaction networks

Short Linear Motifs in Colorectal Cancer Interactome and Tumorigenesis

Colorectal tumorigenesis is driven by alterations in genes and proteins responsible for cancer initiation, progression, and invasion. This multistage process is based on a dense network of protein-protein interactions (PPIs) that become dysregulated as a result of changes in various cell signaling effectors. PPIs in signaling and regulatory networks are known to be mediated by short linear motifs (SLiMs), which are conserved contiguous regions of 3-10 amino acids within interacting protein domains. SLiMs are the minimum sequences required for modulating cellular PPI networks. Thus, several in silico approaches have been developed to predict and analyze SLiM-mediated PPIs. In this review, we focus on emerging evidence supporting a crucial role for SLiMs in driver pathways that are disrupted in colorectal cancer (CRC) tumorigenesis and related PPI network alterations. As a result, SLiMs, along with short peptides, are attracting the interest of researchers to devise small molecules amenable to be used as novel anti-CRC targeted therapies. Overall, the characterization of SLiMs mediating crucial PPIs in CRC may foster the development of more specific combined pharmacological approaches.

Thousands of protein linear motif classes may still be undiscovered

Linear motifs are short protein subsequences that mediate protein interactions. Hundreds of motif classes including thousands of motif instances are known. Our theory estimates how many motif classes remain undiscovered. As commonly done, we describe motif classes as regular expressions specifying motif length and the allowed amino acids at each motif position. We measure motif specificity for a pair of motif classes by quantifying how many motif-discriminating positions prevent a protein subsequence from matching the two classes at once. We derive theorems for the maximal number of motif classes that can simultaneously maintain a certain number of motif-discriminating positions between all pairs of classes in the motif universe, for a given amino acid alphabet. We also calculate the fraction of all protein subsequences that would belong to a motif class if all potential motif classes came into existence. Naturally occurring pairs of motif classes present most often a single motif-discriminating position. This mild specificity maximizes the potential number of coexisting motif classes, the expansion of the motif universe due to amino acid modifications and the fraction of amino acid sequences that code for a motif instance. As a result, thousands of linear motif classes may remain undiscovered.

IMCC: A Novel Quantitative Approach Revealing Variation of Global Modular Map and Local Inter-Module Coordination Among Differential Drug's Targeted Cerebral Ischemic Networks

Stroke is a common disease characterized by multiple genetic dysfunctions. In this complex disease, detecting the strength of inter-module coordination (genetic community interaction) and subsequent modular rewiring is essential to characterize the reactive biosystematic variation (biosystematic perturbation) brought by multiple-target drugs, whose effects are achieved by hitting on a series of targets (target profile) jointly. Here, a quantitative approach for inter-module coordination and its transition, named as IMCC, was developed. Applying IMCC to mouse cerebral ischemia–related gene microarray, we investigated a holistic view of modular map and its rewiring from ischemic stroke to drugs (baicalin, BA; ursodeoxycholic acid, UA; and jasminoidin, JA) perturbation states and locally identified the cooperative pathological module pair and its dissection. Our result suggested the global modular map in cerebral ischemia exhibited a characteristic “core–periphery” architecture, and this architecture was rewired by the effective drugs heterogeneously: BA and UA converged modules into an intensively connected integrity, whereas JA diverged partial modules and widened the remaining inter-module paths. Locally, the PMP dissociation brought by drugs contributed to the reversion of the pathological condition: the focus of the cellular function shift from survival after nervous system injury into development and repair, including neurotrophin regulation, hormone releasing, and chemokine signaling activation. The core targets and mechanisms were validated by in vivo experiments. Overall, our result highlights the holistic inter-module coordination rearrangement rather than a target or a single module that brings phenotype alteration. This strategy may lead to systematically explore detailed variation of inter-module pharmacological action mode of multiple-target drugs, which is the principal problem of module pharmacology for network-based drug discovery.

Domain-mediated interactions for protein subfamily identification

Within a protein family, proteins with the same domain often exhibit different cellular functions, despite the shared evolutionary history and molecular function of the domain. We hypothesized that domain-mediated interactions (DMIs) may categorize a protein family into subfamilies because the diversified functions of a single domain often depend on interacting partners of domains. Here we systematically identified DMI subfamilies, in which proteins share domains with DMI partners, as well as with various functional and physical interaction networks in individual species. In humans, DMI subfamily members are associated with similar diseases, including cancers, and are frequently co-associated with the same diseases. DMI information relates to the functional and evolutionary subdivisions of human kinases. In yeast, DMI subfamilies contain proteins with similar phenotypic outcomes from specific chemical treatments. Therefore, the systematic investigation here provides insights into the diverse functions of subfamilies derived from a protein family with a link-centric approach and suggests a useful resource for annotating the functions and phenotypic outcomes of proteins.

Systematic prediction of FFAT motifs across eukaryote proteomes identifies nucleolar and eisosome proteins with the predicted capacity to form bridges to the endoplasmic reticulum

The endoplasmic reticulum (ER), the most pervasive organelle, exchanges information and material with many other organelles, but the extent of its inter-organelle connections and the proteins that form bridges are not well known. The integral ER membrane protein VAMP-associated protein (VAP) is found in multiple bridges, interacting with many proteins that contain a short linear motif consisting of "two phenylalanines in an acidic tract" (FFAT). The VAP-FFAT interaction is the most common mechanism by which cytoplasmic proteins, particularly inter-organelle bridges, target the ER. Therefore, predicting new FFAT motifs may both find new individual peripheral ER proteins and identify new routes of communication involving the ER. Here we searched for FFAT motifs across whole proteomes. The excess of eukaryotic proteins with FFAT motifs over background was ≥0.8%, suggesting this is the minimum number of peripheral ER proteins. In yeast, where VAP was previously known to bind 4 proteins with FFAT motifs, a detailed analysis of a subset of proteins predicted 20 FFAT motifs. Extrapolating these findings to the whole proteome estimated the number of FFAT motifs in yeast at approximately 50-55 (0.9% of proteome). Among these previously unstudied FFAT motifs, most have known functions outside the ER, so could be involved in inter-organelle communication. Many of these can target well-characterised membrane contact sites, however some are in nucleoli and eisosomes, organelles previously unknown to have molecular bridges to the ER. We speculate that the nucleolar and eisosomal proteins with predicted motifs may function while bridging to the ER, indicating novel ER-nucleolus and ER-eisosome routes of inter-organelle communication.

Link clustering explains non-central and contextually essential genes in protein interaction networks

Recent studies have shown that many essential genes (EGs) change their essentiality across various contexts. Finding contextual EGs in pathogenic conditions may facilitate the identification of therapeutic targets. We propose link clustering as an indicator of contextual EGs that are non-central in protein-protein interaction (PPI) networks. In various human and yeast PPI networks, we found that 29–47% of EGs were better characterized by link clustering than by centrality. Importantly, non-central EGs were prone to change their essentiality across different human cell lines and between species. Compared with central EGs and non-EGs, non-central EGs had intermediate levels of expression and evolutionary conservation. In addition, non-central EGs exhibited a significant impact on communities at lower hierarchical levels, suggesting that link clustering is associated with contextual essentiality, as it depicts locally important nodes in network structures.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.

Controllability of protein-protein interaction phosphorylation-based networks: Participation of the hub 14-3-3 protein family

Posttranslational regulation of protein function is an ubiquitous mechanism in eukaryotic cells. Here, we analyzed biological properties of nodes and edges of a human protein-protein interaction phosphorylation-based network, especially of those nodes critical for the network controllability. We found that the minimal number of critical nodes needed to control the whole network is 29%, which is considerably lower compared to other real networks. These critical nodes are more regulated by posttranslational modifications and contain more binding domains to these modifications than other kinds of nodes in the network, suggesting an intra-group fast regulation. Also, when we analyzed the edges characteristics that connect critical and non-critical nodes, we found that the former are enriched in domain-to-eukaryotic linear motif interactions, whereas the later are enriched in domain-domain interactions. Our findings suggest a possible structure for protein-protein interaction networks with a densely interconnected and self-regulated central core, composed of critical nodes with a high participation in the controllability of the full network, and less regulated peripheral nodes. Our study offers a deeper understanding of complex network control and bridges the controllability theorems for complex networks and biological protein-protein interaction phosphorylation-based networked systems.

K-shell Analysis Reveals Distinct Functional Parts in an Electron Transfer Network and Its Implications for Extracellular Electron Transfer

Shewanella oneidensis MR-1 is capable of extracellular electron transfer (EET) and hence has attracted considerable attention. The EET pathways mainly consist of c-type cytochromes, along with some other proteins involved in electron transfer processes. By whole genome study and protein interactions inquisition, we constructed a large-scale electron transfer network containing 2276 interactions among 454 electron transfer related proteins in S. oneidensis MR-1. Using the k-shell decomposition method, we identified and analyzed distinct parts of the electron transfer network. We found that there was a negative correlation between the k s (k-shell values) and the average DR_100 (disordered regions per 100 amino acids) in every shell, which suggested that disordered regions of proteins played an important role during the formation and extension of the electron transfer network. Furthermore, proteins in the top three shells of the network are mainly located in the cytoplasm and inner membrane; these proteins can be responsible for transfer of electrons into the quinone pool in a wide variety of environmental conditions. In most of the other shells, proteins are broadly located throughout the five cellular compartments (cytoplasm, inner membrane, periplasm, outer membrane, and extracellular), which ensures the important EET ability of S. oneidensis MR-1. Specifically, the fourth shell was responsible for EET and the c-type cytochromes in the remaining shells of the electron transfer network were involved in aiding EET. Taken together, these results show that there are distinct functional parts in the electron transfer network of S. oneidensis MR-1, and the EET processes could achieve high efficiency through cooperation through such an electron transfer network.

Hammock: a hidden Markov model-based peptide clustering algorithm to identify protein-interaction consensus motifs in large datasets

Motivation: Proteins often recognize their interaction partners on the basis of short linear motifs located in disordered regions on proteins’ surface. Experimental techniques that study such motifs use short peptides to mimic the structural properties of interacting proteins. Continued development of these methods allows for large-scale screening, resulting in vast amounts of peptide sequences, potentially containing information on multiple protein-protein interactions. Processing of such datasets is a complex but essential task for large-scale studies investigating protein-protein interactions.

Results: The software tool presented in this article is able to rapidly identify multiple clusters of sequences carrying shared specificity motifs in massive datasets from various sources and generate multiple sequence alignments of identified clusters. The method was applied on a previously published smaller dataset containing distinct classes of ligands for SH3 domains, as well as on a new, an order of magnitude larger dataset containing epitopes for several monoclonal antibodies. The software successfully identified clusters of sequences mimicking epitopes of antibody targets, as well as secondary clusters revealing that the antibodies accept some deviations from original epitope sequences. Another test indicates that processing of even much larger datasets is computationally feasible.

Availability and implementation: Hammock is published under GNU GPL v. 3 license and is freely available as a standalone program (from http://www.recamo.cz/en/software/hammock-cluster-peptides/) or as a tool for the Galaxy toolbox (from https://toolshed.g2.bx.psu.edu/view/hammock/hammock). The source code can be downloaded from https://github.com/hammock-dev/hammock/releases.

Contact:muller@mou.cz

Supplementaryinformation:Supplementary data are available at Bioinformatics online.

How do disordered regions achieve comparable functions to structured domains?

The traditional structure to function paradigm conceives of a protein's function as emerging from its structure. In recent years, it has been established that unstructured, intrinsically disordered regions (IDRs) in proteins are equally crucial elements for protein function, regulation and homeostasis. In this review, we provide a brief overview of how IDRs can perform similar functions to structured proteins, focusing especially on the formation of protein complexes and assemblies and the mediation of regulated conformational changes. In addition to highlighting instances of such functional equivalence, we explain how differences in the biological and physicochemical properties of IDRs allow them to expand the functional and regulatory repertoire of proteins. We also discuss studies that provide insights into how mutations within functional regions of IDRs can lead to human diseases.

LMPID: a manually curated database of linear motifs mediating protein-protein interactions

Linear motifs (LMs), used by a subset of all protein-protein interactions (PPIs), bind to globular receptors or domains and play an important role in signaling networks. LMPID (Linear Motif mediated Protein Interaction Database) is a manually curated database which provides comprehensive experimentally validated information about the LMs mediating PPIs from all organisms on a single platform. About 2200 entries have been compiled by detailed manual curation of PubMed abstracts, of which about 1000 LM entries were being annotated for the first time, as compared with the Eukaryotic LM resource. The users can submit their query through a user-friendly search page and browse the data in the alphabetical order of the bait gene names and according to the domains interacting with the LM. LMPID is freely accessible at http://bicresources.jcbose. ac.in/ssaha4/lmpid and contains 1750 unique LM instances found within 1181 baits interacting with 552 prey proteins. In summary, LMPID is an attempt to enrich the existing repertoire of resources available for studying the LMs implicated in PPIs and may help in understanding the patterns of LMs binding to a specific domain and develop prediction model to identify novel LMs specific to a domain and further able to predict inhibitors/modulators of PPI of interest.