Information Flow Analysis of Interactome Networks

Recent studies of cellular networks have revealed modular organizations of genes and proteins. For example, in interactome networks, a module refers to a group of interacting proteins that form molecular complexes and/or biochemical pathways and together mediate a biological process. However, it is still poorly understood how biological information is transmitted between different modules. We have developed information flow analysis, a new computational approach that identifies proteins central to the transmission of biological information throughout the network. In the information flow analysis, we represent an interactome network as an electrical circuit, where interactions are modeled as resistors and proteins as interconnecting junctions. Construing the propagation of biological signals as flow of electrical current, our method calculates an information flow score for every protein. Unlike previous metrics of network centrality such as degree or betweenness that only consider topological features, our approach incorporates confidence scores of protein–protein interactions and automatically considers all possible paths in a network when evaluating the importance of each protein. We apply our method to the interactome networks of Saccharomyces cerevisiae and Caenorhabditis elegans. We find that the likelihood of observing lethality and pleiotropy when a protein is eliminated is positively correlated with the protein's information flow score. Even among proteins of low degree or low betweenness, high information scores serve as a strong predictor of loss-of-function lethality or pleiotropy. The correlation between information flow scores and phenotypes supports our hypothesis that the proteins of high information flow reside in central positions in interactome networks. We also show that the ranks of information flow scores are more consistent than that of betweenness when a large amount of noisy data is added to an interactome. Finally, we combine gene expression data with interaction data in C. elegans and construct an interactome network for muscle-specific genes. We find that genes that rank high in terms of information flow in the muscle interactome network but not in the entire network tend to play important roles in muscle function. This framework for studying tissue-specific networks by the information flow model can be applied to other tissues and other organisms as well.     

Identification of potential host proteins for influenza A virus based on topological and biological characteristics by proteome-wide network approach

A proteome-wide network approach was performed to characterize significant patterns of influenza A virus (IAV)-human interactions, and to further identify potentially valuable targets for prophylactic and therapeutic interventions. Topological analysis demonstrated a strong tendency for IAV to interplay with highly connected and central proteins located in sparsely connected sub-networks. Additionally, functional analysis based on biological process revealed a number of functional groups overrepresented for IAV interactions, in which regulation of cell death and apoptosis, and phosphorus metabolic process is the most highly enriched. In order to investigate whether these topological and biological features are significant enough to distinguish IAV targets from human proteome, a discrimination model was constructed based on these features using support vector machine coupled with genetic algorithm. The average result of overall prediction accuracy is 71.04% by leave-one-out across validation test. The optimized classifier was then applied to 9706 human proteins. As a result, 1418 novel genes were identified from human interactome, some of which were experimentally validated by others' works to be important for IAV infection. The findings presented in this study might be important in discovering new drug targets for therapeutic treatments as well as revealing topological features and functional properties specific for viral infection.     

Protein Interaction Network of Arabidopsis thaliana Female Gametophyte Development Identifies Novel Proteins and Relations

Although the female gametophyte in angiosperms consists of just seven cells, it has a complex biological network. In this study, female gametophyte microarray data from Arabidopsis thaliana were integrated into the Arabidopsis interactome database to generate a putative interaction map of the female gametophyte development including proteome map based on biological processes and molecular functions of proteins. Biological and functional groups as well as topological characteristics of the network were investigated by analyzing phytohormones, plant defense, cell death, transporters, regulatory factors, and hydrolases. This approach led to the prediction of critical members and bottlenecks of the network. Seventy-four and 24 upregulated genes as well as 171 and 3 downregulated genes were identified in subtracted networks based on biological processes and molecular function respectively, including novel genes such as the pathogenesis-related protein 4, ER type Ca²⁺ ATPase 3, dihydroflavonol reductase, and ATP disulfate isomerase. Biologically important relationships between genes, critical nodes, and new essential proteins such as AT1G26830, AT5G20850, CYP74A, AT1G42396, PR4 and MEA were found in the interactome''s network. The positions of novel genes, both upregulated and downregulated, and their relationships with biological pathways, in particular phytohormones, were highlighted in this study.     

Probabilistic model of the human protein-protein interaction network

A catalog of all human protein-protein interactions would provide scientists with a framework to study protein deregulation in complex diseases such as cancer. Here we demonstrate that a probabilistic analysis integrating model organism interactome data, protein domain data, genome-wide gene expression data and functional annotation data predicts nearly 40,000 protein-protein interactions in humans-a result comparable to those obtained with experimental and computational approaches in model organisms. We validated the accuracy of the predictive model on an independent test set of known interactions and also experimentally confirmed two predicted interactions relevant to human cancer, implicating uncharacterized proteins into definitive pathways. We also applied the human interactome network to cancer genomics data and identified several interaction subnetworks activated in cancer. This integrative analysis provides a comprehensive framework for exploring the human protein interaction network.     

Schizophrenia Gene Networks and Pathways and Their Applications for Novel Candidate Gene Selection

Background

Schizophrenia (SZ) is a heritable, complex mental disorder. We have seen limited success in finding causal genes for schizophrenia from numerous conventional studies. Protein interaction network and pathway-based analysis may provide us an alternative and effective approach to investigating the molecular mechanisms of schizophrenia.

Methodology/Principal Findings

We selected a list of schizophrenia candidate genes (SZGenes) using a multi-dimensional evidence-based approach. The global network properties of proteins encoded by these SZGenes were explored in the context of the human protein interactome while local network properties were investigated by comparing SZ-specific and cancer-specific networks that were extracted from the human interactome. Relative to cancer genes, we observed that SZGenes tend to have an intermediate degree and an intermediate efficiency on a perturbation spreading throughout the human interactome. This suggested that schizophrenia might have different pathological mechanisms from cancer even though both are complex diseases. We conducted pathway analysis using Ingenuity System and constructed the first schizophrenia molecular network (SMN) based on protein interaction networks, pathways and literature survey. We identified 24 pathways overrepresented in SZGenes and examined their interactions and crosstalk. We observed that these pathways were related to neurodevelopment, immune system, and retinoic X receptor (RXR). Our examination of SMN revealed that schizophrenia is a dynamic process caused by dysregulation of the multiple pathways. Finally, we applied the network/pathway approach to identify novel candidate genes, some of which could be verified by experiments.

Conclusions/Significance

This study provides the first comprehensive review of the network and pathway characteristics of schizophrenia candidate genes. Our preliminary results suggest that this systems biology approach might prove promising for selection of candidate genes for complex diseases. Our findings have important implications for the molecular mechanisms for schizophrenia and, potentially, other psychiatric disorders.     

The protein–protein interaction network of the human Sirtuin family

Protein–protein interaction networks are useful for studying human diseases and to look for possible health care through a holistic approach. Networks are playing an increasing and important role in the understanding of physiological processes such as homeostasis, signaling, spatial and temporal organizations, and pathological conditions. In this article we show the complex system of interactions determined by human Sirtuins (Sirt) largely involved in many metabolic processes as well as in different diseases. The Sirtuin family consists of seven homologous Sirt-s having structurally similar cores but different terminal segments, being rather variable in length and/or intrinsically disordered. Many studies have determined their cellular location as well as biological functions although molecular mechanisms through which they act are actually little known therefore, the aim of this work was to define, explore and understand the Sirtuin-related human interactome. As a first step, we have integrated the experimentally determined protein–protein interactions of the Sirtuin-family as well as their first and second neighbors to a Sirtuin-related sub-interactome. Our data showed that the second-neighbor network of Sirtuins encompasses 25% of the entire human interactome, and exhibits a scale-free degree distribution and interconnectedness among top degree nodes. Moreover, the Sirtuin sub interactome showed a modular structure around the core comprising mixed functions. Finally, we extracted from the Sirtuin sub-interactome subnets related to cancer, aging and post-translational modifications for information on key nodes and topological space of the subnets in the Sirt family network.     

Examination of the relationship between essential genes in PPI network and hub proteins in reverse nearest neighbor topology

Background  

In many protein-protein interaction (PPI) networks, densely connected hub proteins are more likely to be essential proteins. This is referred to as the "centrality-lethality rule", which indicates that the topological placement of a protein in PPI network is connected with its biological essentiality. Though such connections are observed in many PPI networks, the underlying topological properties for these connections are not yet clearly understood. Some suggested putative connections are the involvement of essential proteins in the maintenance of overall network connections, or that they play a role in essential protein clusters. In this work, we have attempted to examine the placement of essential proteins and the network topology from a different perspective by determining the correlation of protein essentiality and reverse nearest neighbor topology (RNN).     

Unraveling protein interaction networks with near-optimal efficiency

The functional characterization of genes and their gene products is the main challenge of the genomic era. Examining interaction information for every gene product is a direct way to assemble the jigsaw puzzle of proteins into a functional map. Here we demonstrate a method in which the information gained from pull-down experiments, in which single proteins act as baits to detect interactions with other proteins, is maximized by using a network-based strategy to select the baits. Because of the scale-free distribution of protein interaction networks, we were able to obtain fast coverage by focusing on highly connected nodes (hubs) first. Unfortunately, locating hubs requires prior global information about the network one is trying to unravel. Here, we present an optimized 'pay-as-you-go' strategy that identifies highly connected nodes using only local information that is collected as successive pull-down experiments are performed. Using this strategy, we estimate that 90% of the human interactome can be covered by 10,000 pull-down experiments, with 50% of the interactions confirmed by reciprocal pull-down experiments.     

Molecular Architecture of Spinal Cord Injury Protein Interaction Network

Spinal cord injury (SCI) is associated with complex pathophysiological processes that follow the primary traumatic event and determine the extent of secondary damage and functional recovery. Numerous reports have used global and hypothesis-driven approaches to identify protein changes that contribute to the overall pathology of SCI in an effort to identify potential therapeutic interventions. In this study, we use a semi-automatic annotation approach to detect terms referring to genes or proteins dysregulated in the SCI literature and develop a curated SCI interactome. Network analysis of the SCI interactome revealed the presence of a rich-club organization corresponding to a “powerhouse” of highly interacting hub-proteins. Studying the modular organization of the network have shown that rich-club proteins cluster into modules that are specifically enriched for biological processes that fall under the categories of cell death, inflammation, injury recognition and systems development. Pathway analysis of the interactome and the rich-club revealed high similarity indicating the role of the rich-club proteins as hubs of the most prominent pathways in disease pathophysiology and illustrating the centrality of pro-and anti-survival signal competition in the pathology of SCI. In addition, evaluation of centrality measures of single nodes within the rich-club have revealed that neuronal growth factor (NGF), caspase 3, and H-Ras are the most central nodes and potentially an interesting targets for therapy. Our integrative approach uncovers the molecular architecture of SCI interactome, and provide an essential resource for evaluating significant therapeutic candidates.     

Roles for the two-hybrid system in exploration of the yeast protein interactome

Comprehensive analysis of protein-protein interactions is a challenging endeavor of functional proteomics and has been best explored in the budding yeast. The yeast protein interactome analysis was achieved first by using the yeast two-hybrid system in a proteome-wide scale and next by large-scale mass spectrometric analysis of affinity-purified protein complexes. While these interaction data have led to a number of novel findings and the emergence of a single huge network containing thousands of proteins, they suffer many false signals and fall short of grasping the entire interactome. Thus, continuous efforts are necessary in both bioinformatics and experimentation to fully exploit these data and to proceed another step forward to the goal. Computational tools to integrate existing biological knowledge buried in literature and various functional genomic data with the interactome data are required for biological interpretation of the huge protein interaction network. Novel experimental methods have to be developed to detect weak, transient interactions involving low abundance proteins as well as to obtain clues to the biological role for each interaction. Since the yeast two-hybrid system can be used for the mapping of the interaction domains and the isolation of interaction-defective mutants, it would serve as a technical basis for the latter purpose, thereby playing another important role in the next phase of protein interactome research.     

Peeling the yeast protein network

A set of highly connected proteins (or hubs) plays an important role for the integrity of the protein interaction network of Saccharomyces cerevisae by connecting the network's intrinsic modules. The importance of the hubs' central placement is further confirmed by their propensity to be lethal. However, although highly emphasized, little is known about the topological coherence among the hubs. Applying a core decomposition method which allows us to identify the inherent layer structure of the protein interaction network, we find that the probability of nodes both being essential and evolutionary conserved successively increases toward the innermost cores. While connectivity alone is often not a sufficient criterion to assess a protein's functional, evolutionary and topological relevance, we classify nodes as globally and locally central depending on their appearance in the inner or outer cores. The observation that globally central proteins participate in a substantial number of protein complexes which display an elevated degree of evolutionary conservation allows us to hypothesize that globally central proteins serve as the evolutionary backbone of the proteome. Even though protein interaction data are extensively flawed, we find that our results are very robust against inaccurately determined protein interactions.     

Why do hubs in the yeast protein interaction network tend to be essential: reexamining the connection between the network topology and essentiality

The centrality-lethality rule, which notes that high-degree nodes in a protein interaction network tend to correspond to proteins that are essential, suggests that the topological prominence of a protein in a protein interaction network may be a good predictor of its biological importance. Even though the correlation between degree and essentiality was confirmed by many independent studies, the reason for this correlation remains illusive. Several hypotheses about putative connections between essentiality of hubs and the topology of protein-protein interaction networks have been proposed, but as we demonstrate, these explanations are not supported by the properties of protein interaction networks. To identify the main topological determinant of essentiality and to provide a biological explanation for the connection between the network topology and essentiality, we performed a rigorous analysis of six variants of the genomewide protein interaction network for Saccharomyces cerevisiae obtained using different techniques. We demonstrated that the majority of hubs are essential due to their involvement in Essential Complex Biological Modules, a group of densely connected proteins with shared biological function that are enriched in essential proteins. Moreover, we rejected two previously proposed explanations for the centrality-lethality rule, one relating the essentiality of hubs to their role in the overall network connectivity and another relying on the recently published essential protein interactions model.     

Computational analysis of human protein interaction networks

Large amounts of human protein interaction data have been produced by experiments and prediction methods. However, the experimental coverage of the human interactome is still low in contrast to predicted data. To gain insight into the value of publicly available human protein network data, we compared predicted datasets, high-throughput results from yeast two-hybrid screens, and literature-curated protein-protein interactions. This evaluation is not only important for further methodological improvements, but also for increasing the confidence in functional hypotheses derived from predictions. Therefore, we assessed the quality and the potential bias of the different datasets using functional similarity based on the Gene Ontology, structural iPfam domain-domain interactions, likelihood ratios, and topological network parameters. This analysis revealed major differences between predicted datasets, but some of them also scored at least as high as the experimental ones regarding multiple quality measures. Therefore, since only small pair wise overlap between most datasets is observed, they may be combined to enlarge the available human interactome data. For this purpose, we additionally studied the influence of protein length on data quality and the number of disease proteins covered by each dataset. We could further demonstrate that protein interactions predicted by more than one method achieve an elevated reliability.     

Systems properties of proteins encoded by imprinted genes

Genomically imprinted genes show parentally fixed mono-allelic expression and are important for the mammalian development. Dysregulation of genomic imprinting leads to several complex pathological conditions. Though the genetic and epigenetic regulation of imprinted genes has been well studied, their protein aspects are largely ignored. Here, we systematically studied a sub-network centered on proteins encoded by imprinted genes within human interactome. Using concepts of network biology, we uncover a highly connected, transitive and central network module of imprinted gene-products and their interacting partners (IGPN). The network is enriched in development, metabolism and cell cycle related functions and its malfunctioning ascribes error intolerance to human interactome network. Further, detailed analysis revealed that its higher centrality is determined by ‘date’ interactions among the proteins belonging to different functional classes than the ‘party’ interactions within the same functional class. Interestingly, a significant proportion of this network genetically associates with disease phenotypes. Moreover, the network comprises of gene-sets that are upregulated in leukemia, psychosis, obesity/diabetes and downregulated in autism. We conclude that imprinted gene-products are part of a functionally and topologically important module of human interactome and errors in this sub-network are intolerant to, otherwise robust, human interactome. The findings might also shed light on how imprinted genes, which are rather very few, coordinate at protein level to pleiotropically regulate growth and metabolism during embryonic and post-natal development.     

Analysis and identification of essential genes in humans using topological properties and biological information

Genes that are indispensable for survival are termed essential genes. The analysis and identification of essential genes are very important for understanding the minimal requirements of cellular survival and for practical purposes. Proteins do not exert their function in isolation of one another but rather interact together in PPI networks. A global analysis of protein interaction networks provides an effective way to elucidate the relationships between proteins. With the recent large-scale identifications of essential genes and the production of large amounts of PPIs in humans, we are able to investigate the topological properties and biological properties of essential genes. However, until recently, no one has ever investigated human essential genes using topological and biological properties. In this study, for the first time, 28 topological properties and 22 biological properties were used to investigate the characteristics of essential and non-essential genes in humans. Most of the properties were statistically discriminative between essential and non-essential genes. The F-score was used to estimate the essentiality of each property. The GO-enrichment analysis was performed to investigate the functions of the essential and non-essential genes. Finally, based on the topological features and the biological characteristics, a machine-learning classifier was constructed to predict the essential genes. The results of the jackknife test and 10-fold cross validation test are encouraging, indicating that our classifier is an effective human essential gene discovery method.     

A systems biology approach reveals the role of a novel methyltransferase in response to chemical stress and lipid homeostasis

Using small molecule probes to understand gene function is an attractive approach that allows functional characterization of genes that are dispensable in standard laboratory conditions and provides insight into the mode of action of these compounds. Using chemogenomic assays we previously identified yeast Crg1, an uncharacterized SAM-dependent methyltransferase, as a novel interactor of the protein phosphatase inhibitor cantharidin. In this study we used a combinatorial approach that exploits contemporary high-throughput techniques available in Saccharomyces cerevisiae combined with rigorous biological follow-up to characterize the interaction of Crg1 with cantharidin. Biochemical analysis of this enzyme followed by a systematic analysis of the interactome and lipidome of CRG1 mutants revealed that Crg1, a stress-responsive SAM-dependent methyltransferase, methylates cantharidin in vitro. Chemogenomic assays uncovered that lipid-related processes are essential for cantharidin resistance in cells sensitized by deletion of the CRG1 gene. Lipidome-wide analysis of mutants further showed that cantharidin induces alterations in glycerophospholipid and sphingolipid abundance in a Crg1-dependent manner. We propose that Crg1 is a small molecule methyltransferase important for maintaining lipid homeostasis in response to drug perturbation. This approach demonstrates the value of combining chemical genomics with other systems-based methods for characterizing proteins and elucidating previously unknown mechanisms of action of small molecule inhibitors.