Simple topological features reflect dynamics and modularity in protein interaction networks

A partially shared joint clustering framework for detecting protein complexes from multiple state-specific signed interaction networks

Detecting protein complexes is critical for studying cellular organizations and functions. The accumulation of protein-protein interaction (PPI) data enables the identification of protein complexes computationally. Although a great number of computational methods have been proposed to identify protein complexes from PPI networks, most of them ignore the signs of PPIs that reflect the ways proteins interact (activation or inhibition). As not all PPIs imply co-complex relationships, taking into account the signs of PPIs can benefit the identification of protein complexes. Moreover, PPI networks are not static, but vary with the change of cell states or environments. However, existing methods are primarily designed for single-network clustering, and rarely consider joint clustering of multiple PPI networks. In this study, we propose a novel partially shared signed network clustering (PS-SNC) model for identifying protein complexes from multiple state-specific signed PPI networks jointly. PS-SNC can not only consider the signs of PPIs, but also identify the common and unique protein complexes in different states. Experimental results on synthetic and real datasets show that our PS-SNC model can achieve better performance than other state-of-the-art protein complex detection methods. Extensive analysis on real datasets demonstrate the effectiveness of PS-SNC in revealing novel insights about the underlying patterns of different cell lines.

Protein conformational switch discerned via network centrality properties

Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensemble-informed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.

Identification of Crucial Candidate Genes and Pathways in Glioblastoma Multiform by Bioinformatics Analysis

The present study aimed to investigate the molecular mechanisms underlying glioblastoma multiform (GBM) and its biomarkers. The differentially expressed genes (DEGs) were diagnosed using the limma software package. The ToppGene (ToppFun) was used to perform pathway and Gene Ontology (GO) enrichment analysis of the DEGs. Protein-protein interaction (PPI) networks, extracted modules, miRNA-target genes regulatory network and TF-target genes regulatory network were used to obtain insight into the actions of DEGs. Survival analysis for DEGs was carried out. A total of 590 DEGs, including 243 up regulated and 347 down regulated genes, were diagnosed between scrambled shRNA expression and Lin7A knock down. The up-regulated genes were enriched in ribosome, mitochondrial translation termination, translation, and peptide biosynthetic process. The down-regulated genes were enriched in focal adhesion, VEGFR3 signaling in lymphatic endothelium, extracellular matrix organization, and extracellular matrix. The current study screened the genes in the PPI network, extracted modules, miRNA-target genes regulatory network, and TF-target genes regulatory network with higher degrees as hub genes, which included NPM1, CUL4A, YIPF1, SHC1, AKT1, VLDLR, RPL14, P3H2, DTNA, FAM126B, RPL34, and MYL5. Survival analysis indicated that the high expression of RPL36A and MRPL35 were predicting longer survival of GBM, while high expression of AP1S1 and AKAP12 were predicting shorter survival of GBM. High expression of RPL36A and AP1S1 were associated with pathogenesis of GBM, while low expression of ALPL was associated with pathogenesis of GBM. In conclusion, the current study diagnosed DEGs between scrambled shRNA expression and Lin7A knock down samples, which could improve our understanding of the molecular mechanisms in the progression of GBM, and these crucial as well as new diagnostic markers might be used as therapeutic targets for GBM.

Self-organizing scale-free patterns in a phase-modulated periodic connecting system

OBJECTIVE

The regularity of scale-free patterns in rank-size relations has been observed in word frequency, city size distribution, firm size distribution, and gene expression. Because of the common emergence of this regularity, understanding its mechanisms has been of great interest. For obtaining the scale-free pattern regularity, various models based on the rich-get-richer mechanism have been proposed; however, the overarching procedure of searching for the "rich" is in disagreement with the locally interacting behaviors seen in the aforementioned natural and social phenomena.

RESULTS

We implemented a computational model of a resource distribution system inspired by observations of word connectivity, which is created by local constraints with periodic and phase modulatory features. Here, we empirically demonstrated that a phase-modulated periodic connecting system can reach a dynamic equilibrium state as the most probable case, with the self-organizing scale-free patterns. The regularity could be a result of the configurational balance in spatiotemporal inequity during the resource distribution process with an adaptive constrained connectivity. Our results suggest that investigations of interferences of oscillating fluctuations in the system will elucidate the autoregulatory dynamic behavior.

Exploiting locational and topological overlap model to identify modules in protein interaction networks

Background

Clustering molecular network is a typical method in system biology, which is effective in predicting protein complexes or functional modules. However, few studies have realized that biological molecules are spatial-temporally regulated to form a dynamic cellular network and only a subset of interactions take place at the same location in cells.

Results

In this study, considering the subcellular localization of proteins, we first construct a co-localization human protein interaction network (PIN) and systematically investigate the relationship between subcellular localization and biological functions. After that, we propose a Locational and Topological Overlap Model (LTOM) to preprocess the co-localization PIN to identify functional modules. LTOM requires the topological overlaps, the common partners shared by two proteins, to be annotated in the same localization as the two proteins. We observed the model has better correspondence with the reference protein complexes and shows more relevance to cancers based on both human and yeast datasets and two clustering algorithms, ClusterONE and MCL.

Conclusion

Taking into consideration of protein localization and topological overlap can improve the performance of module detection from protein interaction networks.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2598-7) contains supplementary material, which is available to authorized users.

Biomarker panels for characterizing microbial community biofilm formation as composite molecular process

Microbial consortia execute collaborative molecular processes with contributions from individual species, on such basis enabling optimized molecular function. Such collaboration and synergies benefit metabolic flux specifically in extreme environmental conditions as seen in acid mine drainage, with biofilms as relevant microenvironment. However, knowledge about community species composition is not sufficient for deducing presence and efficiency of composite molecular function. For this task molecular resolution of the consortium interactome is to be retrieved, with molecular biomarkers particularly suited for characterizing composite molecular processes involved in biofilm formation and maintenance. A microbial species set identified in 18 copper environmental sites provides a data matrix for deriving a cross-species molecular process model of biofilm formation composed of 191 protein coding genes contributed from 25 microbial species. Computing degree and stress centrality of biofilm molecular process nodes allows selection of network hubs and central connectors, with the top ranking molecular features proposed as biomarker candidates for characterizing biofilm homeostasis. Functional classes represented in the biomarker panel include quorum sensing, chemotaxis, motility and extracellular polysaccharide biosynthesis, complemented by chaperones. Abundance of biomarker candidates identified in experimental data sets monitoring different biofilm conditions provides evidence for the selected biomarkers as sensitive and specific molecular process proxies for capturing biofilm microenvironments. Topological criteria of process networks covering an aggregate function of interest support the selection of biomarker candidates independent of specific community species composition. Such panels promise efficient screening of environmental samples for presence of microbial community composite molecular function.

Protein complex prediction by date hub removal

Proteins physically interact with each other and form protein complexes to perform their biological functions. The prediction of protein complexes from protein-protein interaction (PPI) network is usually difficult when the complexes are overlapping with each other in a dense region of the network. To address the problem of predicting overlapping complexes, a previously proposed network-decomposition approach is promising. It decomposes a PPI network by e.g. removing proteins with high degree (hubs) which may participate in different complexes. This motivates us to examine a list of proteins, which bind their different partners at different time or at different location (viz. date hubs), manually collected from literature, for network decomposition. Results show that the CMC complex discovery algorithm after removing date hubs recalls more overlapping complexes that were missed earlier. Further improvement in performance is achieved when we predict date hub proteins based on simple network features and remove them from PPI networks.

Unified Alignment of Protein-Protein Interaction Networks

Paralleling the increasing availability of protein-protein interaction (PPI) network data, several network alignment methods have been proposed. Network alignments have been used to uncover functionally conserved network parts and to transfer annotations. However, due to the computational intractability of the network alignment problem, aligners are heuristics providing divergent solutions and no consensus exists on a gold standard, or which scoring scheme should be used to evaluate them. We comprehensively evaluate the alignment scoring schemes and global network aligners on large scale PPI data and observe that three methods, HUBALIGN, L-GRAAL and NATALIE, regularly produce the most topologically and biologically coherent alignments. We study the collective behaviour of network aligners and observe that PPI networks are almost entirely aligned with a handful of aligners that we unify into a new tool, Ulign. Ulign enables complete alignment of two networks, which traditional global and local aligners fail to do. Also, multiple mappings of Ulign define biologically relevant soft clusterings of proteins in PPI networks, which may be used for refining the transfer of annotations across networks. Hence, PPI networks are already well investigated by current aligners, so to gain additional biological insights, a paradigm shift is needed. We propose such a shift come from aligning all available data types collectively rather than any particular data type in isolation from others.

Interaction and localization diversities of global and local hubs in human protein–protein interaction networks

Hubs, the highly connected nodes in protein–protein interaction networks (PPINs), are associated with several characteristic properties and are known to perform vital roles in cells.

Detecting Protein Complexes from Signed Protein-Protein Interaction Networks

Identification of protein complexes is fundamental for understanding the cellular functional organization. With the accumulation of physical protein-protein interaction (PPI) data, computational detection of protein complexes from available PPI networks has drawn a lot of attentions. While most of the existing protein complex detection algorithms focus on analyzing the physical protein-protein interaction network, none of them take into account the "signs"  (i.e., activation-inhibition relationships) of physical interactions. As the "signs"  of interactions reflect the way proteins communicate, considering the "signs"  of interactions can not only increase the accuracy of protein complex identification, but also deepen our understanding of the mechanisms of cell functions. In this study, we proposed a novel Signed Graph regularized Nonnegative Matrix Factorization (SGNMF) model to identify protein complexes from signed PPI networks. In our experiments, we compared the results collected by our model on signed PPI networks with those predicted by the state-of-the-art complex detection techniques on the original unsigned PPI networks. We observed that considering the "signs"  of interactions significantly benefits the detection of protein complexes. Furthermore, based on the predicted complexes, we predicted a set of signed complex-complex interactions for each dataset, which provides a novel insight of the higher level organization of the cell. All the experimental results and codes can be downloaded from http://mail.sysu.edu.cn/home/stsddq@mail.sysu.edu.cn/dai/others/SGNMF.zip.

Genome-Wide Detection and Analysis of Multifunctional Genes

Many genes can play a role in multiple biological processes or molecular functions. Identifying multifunctional genes at the genome-wide level and studying their properties can shed light upon the complexity of molecular events that underpin cellular functioning, thereby leading to a better understanding of the functional landscape of the cell. However, to date, genome-wide analysis of multifunctional genes (and the proteins they encode) has been limited. Here we introduce a computational approach that uses known functional annotations to extract genes playing a role in at least two distinct biological processes. We leverage functional genomics data sets for three organisms—H. sapiens, D. melanogaster, and S. cerevisiae—and show that, as compared to other annotated genes, genes involved in multiple biological processes possess distinct physicochemical properties, are more broadly expressed, tend to be more central in protein interaction networks, tend to be more evolutionarily conserved, and are more likely to be essential. We also find that multifunctional genes are significantly more likely to be involved in human disorders. These same features also hold when multifunctionality is defined with respect to molecular functions instead of biological processes. Our analysis uncovers key features about multifunctional genes, and is a step towards a better genome-wide understanding of gene multifunctionality.

Almost every aspect of cellular function depends on protein activity. In spite of being fine-tuned to carry out highly specific functions, proteins can also multitask. Experimental studies have identified genes and proteins endowed with more than one molecular function, or participating in very different biological processes. These studies suggest that the degree of functional plasticity exhibited by proteins might go well beyond a simple “one protein—one function” relationship. However, systematic studies of the properties of multifunctional genes (and their encoded proteins) have been limited. Here we present a computational framework to identify putative multifunctional genes, and compare their properties with those of other genes. We find that multifunctional genes are significantly different from other genes with respect to their physicochemical properties, expression profiles, and interaction properties. We also observe that multifunctional genes tend to be more conserved, and that a greater fraction of them are associated with human disorders. Taken together, these results represent a step towards a more complete understanding of the role multifunctional genes play in the functional organization of the cell.

Biological Relevance and Therapeutic Potential of the Hypusine Modification System

Hypusine modification of the eukaryotic initiation factor 5A (eIF-5A) is emerging as a crucial regulator in cancer, infections, and inflammation. Although its contribution in translational regulation of proline repeat-rich proteins has been sufficiently demonstrated, its biological role in higher eukaryotes remains poorly understood. To establish the hypusine modification system as a novel platform for therapeutic strategies, we aimed to investigate its functional relevance in mammals by generating and using a range of new knock-out mouse models for the hypusine-modifying enzymes deoxyhypusine synthase and deoxyhypusine hydroxylase as well as for the cancer-related isoform eIF-5A2. We discovered that homozygous depletion of deoxyhypusine synthase and/or deoxyhypusine hydroxylase causes lethality in adult mice with different penetrance compared with haploinsufficiency. Network-based bioinformatic analysis of proline repeat-rich proteins, which are putative eIF-5A targets, revealed that these proteins are organized in highly connected protein-protein interaction networks. Hypusine-dependent translational control of essential proteins (hubs) and protein complexes inside these networks might explain the lethal phenotype observed after deletion of hypusine-modifying enzymes. Remarkably, our results also demonstrate that the cancer-associated isoform eIF-5A2 is dispensable for normal development and viability. Together, our results provide the first genetic evidence that the hypusine modification in eIF-5A is crucial for homeostasis in mammals. Moreover, these findings highlight functional diversity of the hypusine system compared with lower eukaryotes and indicate eIF-5A2 as a valuable and safe target for therapeutic intervention in cancer.

Mining the modular structure of protein interaction networks

BACKGROUND

Cluster-based descriptions of biological networks have received much attention in recent years fostered by accumulated evidence of the existence of meaningful correlations between topological network clusters and biological functional modules. Several well-performing clustering algorithms exist to infer topological network partitions. However, due to respective technical idiosyncrasies they might produce dissimilar modular decompositions of a given network. In this contribution, we aimed to analyze how alternative modular descriptions could condition the outcome of follow-up network biology analysis.

METHODOLOGY

We considered a human protein interaction network and two paradigmatic cluster recognition algorithms, namely: the Clauset-Newman-Moore and the infomap procedures. We analyzed to what extent both methodologies yielded different results in terms of granularity and biological congruency. In addition, taking into account Guimera's cartographic role characterization of network nodes, we explored how the adoption of a given clustering methodology impinged on the ability to highlight relevant network meso-scale connectivity patterns.

RESULTS

As a case study we considered a set of aging related proteins and showed that only the high-resolution modular description provided by infomap, could unveil statistically significant associations between them and inter/intra modular cartographic features. Besides reporting novel biological insights that could be gained from the discovered associations, our contribution warns against possible technical concerns that might affect the tools used to mine for interaction patterns in network biology studies. In particular our results suggested that sub-optimal partitions from the strict point of view of their modularity levels might still be worth being analyzed when meso-scale features were to be explored in connection with external source of biological knowledge.

Merging and scoring molecular interactions utilising existing community standards: tools, use-cases and a case study

The evidence that two molecules interact in a living cell is often inferred from multiple different experiments. Experimental data is captured in multiple repositories, but there is no simple way to assess the evidence of an interaction occurring in a cellular environment. Merging and scoring of data are commonly required operations after querying for the details of specific molecular interactions, to remove redundancy and assess the strength of accompanying experimental evidence. We have developed both a merging algorithm and a scoring system for molecular interactions based on the proteomics standard initiative-molecular interaction standards. In this manuscript, we introduce these two algorithms and provide community access to the tool suite, describe examples of how these tools are useful to selectively present molecular interaction data and demonstrate a case where the algorithms were successfully used to identify a systematic error in an existing dataset.

Dynamic modularity of host protein interaction networks in Salmonella Typhi infection

Background

Salmonella Typhi is a human-restricted pathogen, which causes typhoid fever and remains a global health problem in the developing countries. Although previously reported host expression datasets had identified putative biomarkers and therapeutic targets of typhoid fever, the underlying molecular mechanism of pathogenesis remains incompletely understood.

Methods

We used five gene expression datasets of human peripheral blood from patients suffering from S. Typhi or other bacteremic infections or non-infectious disease like leukemia. The expression datasets were merged into human protein interaction network (PIN) and the expression correlation between the hubs and their interacting proteins was measured by calculating Pearson Correlation Coefficient (PCC) values. The differences in the average PCC for each hub between the disease states and their respective controls were calculated for studied datasets. The individual hubs and their interactors with expression, PCC and average PCC values were treated as dynamic subnetworks. The hubs that showed unique trends of alterations specific to S. Typhi infection were identified.

Results

We identified S. Typhi infection-specific dynamic subnetworks of the host, which involve 81 hubs and 1343 interactions. The major enriched GO biological process terms in the identified subnetworks were regulation of apoptosis and biological adhesions, while the enriched pathways include cytokine signalling in the immune system and downstream TCR signalling. The dynamic nature of the hubs CCR1, IRS2 and PRKCA with their interactors was studied in detail. The difference in the dynamics of the subnetworks specific to S. Typhi infection suggests a potential molecular model of typhoid fever.

Conclusions

Hubs and their interactors of the S. Typhi infection-specific dynamic subnetworks carrying distinct PCC values compared with the non-typhoid and other disease conditions reveal new insight into the pathogenesis of S. Typhi.

Presence and utility of intrinsically disordered regions in kinases

We reveal presence of intrinsically disordered regions in human kinome and build a kinase–kinase interaction network identifying a novel SRC–SMAD relationship.