Toward the dynamic interactome: it's about time

Beyond modules and hubs: the potential of gene coexpression networks for investigating molecular mechanisms of complex brain disorders

In a research environment dominated by reductionist approaches to brain disease mechanisms, gene network analysis provides a complementary framework in which to tackle the complex dysregulations that occur in neuropsychiatric and other neurological disorders. Gene-gene expression correlations are a common source of molecular networks because they can be extracted from high-dimensional disease data and encapsulate the activity of multiple regulatory systems. However, the analysis of gene coexpression patterns is often treated as a mechanistic black box, in which looming 'hub genes' direct cellular networks, and where other features are obscured. By examining the biophysical bases of coexpression and gene regulatory changes that occur in disease, recent studies suggest it is possible to use coexpression networks as a multi-omic screening procedure to generate novel hypotheses for disease mechanisms. Because technical processing steps can affect the outcome and interpretation of coexpression networks, we examine the assumptions and alternatives to common patterns of coexpression analysis and discuss additional topics such as acceptable datasets for coexpression analysis, the robust identification of modules, disease-related prioritization of genes and molecular systems and network meta-analysis. To accelerate coexpression research beyond modules and hubs, we highlight some emerging directions for coexpression network research that are especially relevant to complex brain disease, including the centrality-lethality relationship, integration with machine learning approaches and network pharmacology.

Dynamic protein interaction modules in human hepatocellular carcinoma progression

Background

Gene expression profiles have been frequently integrated with the human protein interactome to uncover functional modules under specific conditions like disease state. Beyond traditional differential expression analysis, differential co-expression analysis has emerged as a robust approach to reveal condition-specific network modules, with successful applications in a few human disease studies. Hepatocellular carcinoma (HCC), which is often interrelated with the Hepatitis C virus, typically develops through multiple stages. A comprehensive investigation of HCC progression-specific differential co-expression modules may advance our understanding of HCC's pathophysiological mechanisms.

Results

Compared with differentially expressed genes, differentially co-expressed genes were found more likely enriched with Hepatitis C virus binding proteins and cancer-mutated genes, and they were clustered more densely in the human reference protein interaction network. These observations indicated that a differential co-expression approach could outperform the standard differential expression network analysis in searching for disease-related modules. We then proposed a differential co-expression network approach to uncover network modules involved in HCC development. Specifically, we discovered subnetworks that enriched differentially co-expressed gene pairs in each HCC transition stage, and further resolved modules with coherent co-expression change patterns over all HCC developmental stages. Our identified network modules were enriched with HCC-related genes and implicated in cancer-related biological functions. In particular, APC and YWHAZ were highlighted as two most remarkable genes in the network modules, and their dynamic interaction partnership was resolved in HCC development.

Conclusions

We demonstrated that integration of differential co-expression with the protein interactome could outperform the traditional differential expression approach in discovering network modules of human diseases. In our application of this approach to HCC's gene expression data, we successfully identified subnetworks with marked differential co-expression in individual HCC stage transitions and network modules with coherent co-expression change patterns over all HCC developmental stages. Our results shed light on subtle HCC mechanisms, including temporal activation and dismissal of pivotal functions and dynamic interaction partnerships of key genes.

The role of shape complementarity in the protein-protein interactions

We use a dissipative particle dynamic simulation to investigate the effects of shape complementarity on the protein-protein interactions. By monitoring different kinds of protein shape-complementarity modes, we gave a clear mechanism to reveal the role of the shape complementarity in the protein-protein interactions, i.e., when the two proteins with shape complementarity approach each other, the conformation of lipid chains between two proteins would be restricted significantly. The lipid molecules tend to leave the gap formed by two proteins to maximize the configuration entropy, and therefore yield an effective entropy-induced protein-protein attraction, which enhances the protein aggregation. In short, this work provides an insight into understanding the importance of the shape complementarity in the protein-protein interactions especially for protein aggregation and antibody-antigen complexes. Definitely, the shape complementarity is the third key factor affecting protein aggregation and complex, besides the electrostatic-complementarity and hydrophobic complementarity.

Detecting protein complexes from active protein interaction networks constructed with dynamic gene expression profiles

Background

Protein interaction networks (PINs) are known to be useful to detect protein complexes. However, most available PINs are static, which cannot reflect the dynamic changes in real networks. At present, some researchers have tried to construct dynamic networks by incorporating time-course (dynamic) gene expression data with PINs. However, the inevitable background noise exists in the gene expression array, which could degrade the quality of dynamic networkds. Therefore, it is needed to filter out contaminated gene expression data before further data integration and analysis.

Results

Firstly, we adopt a dynamic model-based method to filter noisy data from dynamic expression profiles. Then a new method is proposed for identifying active proteins from dynamic gene expression profiles. An active protein at a time point is defined as the protein the expression level of whose corresponding gene at that time point is higher than a threshold determined by a standard variance involved threshold function. Furthermore, a noise-filtered active protein interaction network (NF-APIN) is constructed. To demonstrate the efficiency of our method, we detect protein complexes from the NF-APIN, compared with those from other dynamic PINs.

Conclusion

A dynamic model based method can effectively filter out noises in dynamic gene expression data. Our method to compute a threshold for determining the active time points of noise-filtered genes can make the dynamic construction more accuracy and provide a high quality framework for network analysis, such as protein complex prediction.

Protein signatures using electrostatic molecular surfaces in harmonic space

We developed a novel method based on the Fourier analysis of protein molecular surfaces to speed up the analysis of the vast structural data generated in the post-genomic era. This method computes the power spectrum of surfaces of the molecular electrostatic potential, whose three-dimensional coordinates have been either experimentally or theoretically determined. Thus we achieve a reduction of the initial three-dimensional information on the molecular surface to the one-dimensional information on pairs of points at a fixed scale apart. Consequently, the similarity search in our method is computationally less demanding and significantly faster than shape comparison methods. As proof of principle, we applied our method to a training set of viral proteins that are involved in major diseases such as Hepatitis C, Dengue fever, Yellow fever, Bovine viral diarrhea and West Nile fever. The training set contains proteins of four different protein families, as well as a mammalian representative enzyme. We found that the power spectrum successfully assigns a unique signature to each protein included in our training set, thus providing a direct probe of functional similarity among proteins. The results agree with established biological data from conventional structural biochemistry analyses.

Simple topological features reflect dynamics and modularity in protein interaction networks

The availability of large-scale protein-protein interaction networks for numerous organisms provides an opportunity to comprehensively analyze whether simple properties of proteins are predictive of the roles they play in the functional organization of the cell. We begin by re-examining an influential but controversial characterization of the dynamic modularity of the S. cerevisiae interactome that incorporated gene expression data into network analysis. We analyse the protein-protein interaction networks of five organisms, S. cerevisiae, H. sapiens, D. melanogaster, A. thaliana, and E. coli, and confirm significant and consistent functional and structural differences between hub proteins that are co-expressed with their interacting partners and those that are not, and support the view that the former tend to be intramodular whereas the latter tend to be intermodular. However, we also demonstrate that in each of these organisms, simple topological measures are significantly correlated with the average co-expression of a hub with its partners, independent of any classification, and therefore also reflect protein intra- and inter- modularity. Further, cross-interactomic analysis demonstrates that these simple topological characteristics of hub proteins tend to be conserved across organisms. Overall, we give evidence that purely topological features of static interaction networks reflect aspects of the dynamics and modularity of interactomes as well as previous measures incorporating expression data, and are a powerful means for understanding the dynamic roles of hubs in interactomes.

A better understanding of protein interaction networks would be a great aid in furthering our knowledge of the molecular biology of the cell. Towards this end, large-scale protein-protein physical interaction data have been determined for organisms across the evolutionary spectrum. However, the resulting networks give a static view of interactomes, and our knowledge about protein interactions is rarely time or context specific. A previous prominent but controversial attempt to characterize the dynamic modularity of the interactome was based on integrating physical interaction data with gene activity measurements from transcript expression data. This analysis distinguished between proteins that are co-expressed with their interacting partners and those that are not, and argued that the former are intramodular and the latter are intermodular. By analyzing the interactomes of five organisms, we largely confirm the biological significance of this characterization through a variety of statistical tests and computational experiments. Surprisingly, however, we find that similar results can be obtained using just network information without additionally integrating expression data, suggesting that purely topological characteristics of interaction networks strongly reflect certain aspects of the dynamics and modularity of interactomes.

A fluorescent reporter for mapping cellular protein-protein interactions in time and space

A method based on a combination of the Split-Ubiquitin system with two spectrally different fluorescent proteins (SPLIFF) is shown to enable measurement of protein interactions in vivo with high spatial and temporal resolution in yeast.

SPLIFF visualizes protein interactions with high spatial and temporal resolution.

Spc72p and Kar9p interact with the MAP Stu2p at opposite poles of microtubules.

Histone chaperone Nap1p and Kcc4 kinase interact preferentially at the bud site.

F-BAR protein Hof1p associates with the polarisome during cell fusion and cytokinesis.

We introduce a fluorescent reporter for monitoring protein–protein interactions in living cells. The method is based on the Split-Ubiquitin method and uses the ratio of two auto-fluorescent reporter proteins as signal for interaction (SPLIFF). The mating of two haploid yeast cells initiates the analysis and the interactions are followed online by two-channel time-lapse microscopy of the diploid cells during their first cell cycle. Using this approach we could with high spatio-temporal resolution visualize the differences between the interactions of the microtubule binding protein Stu2p with two of its binding partners, monitor the transient association of a Ran-GTPase with its receptors at the nuclear pore, and distinguish between protein interactions at the polar cortical domain at different phases of polar growth. These examples further demonstrate that protein–protein interactions identified from large-scale screens can be effectively followed up by high-resolution single-cell analysis.

Identification of highly synchronized subnetworks from gene expression data

Background

There has been a growing interest in identifying context-specific active protein-protein interaction (PPI) subnetworks through integration of PPI and time course gene expression data. However the interaction dynamics during the biological process under study has not been sufficiently considered previously.

Methods

Here we propose a topology-phase locking (TopoPL) based scoring metric for identifying active PPI subnetworks from time series expression data. First the temporal coordination in gene expression changes is evaluated through phase locking analysis; The results are subsequently integrated with PPI to define an activity score for each PPI subnetwork, based on individual member expression, as well topological characteristics of the PPI network and of the expression temporal coordination network; Lastly, the subnetworks with the top scores in the whole PPI network are identified through simulated annealing search.

Results

Application of TopoPL to simulated data and to the yeast cell cycle data showed that it can more sensitively identify biologically meaningful subnetworks than the method that only utilizes the static PPI topology, or the additive scoring method. Using TopoPL we identified a core subnetwork with 49 genes important to yeast cell cycle. Interestingly, this core contains a protein complex known to be related to arrangement of ribosome subunits that exhibit extremely high gene expression synchronization.

Conclusions

Inclusion of interaction dynamics is important to the identification of relevant gene networks.

Construction and application of dynamic protein interaction network based on time course gene expression data

In recent years, researchers have tried to inject dynamic information into static protein interaction networks (PINs). The paper first proposes a three-sigma method to identify active time points of each protein in a cellular cycle, where three-sigma principle is used to compute an active threshold for each gene according to the characteristics of its expression curve. Then a dynamic protein interaction network (DPIN) is constructed, which includes the dynamic changes of protein interactions. To validate the efficiency of DPIN, MCL, CPM, and core attachment algorithms are applied on two different DPINs, the static PIN and the time course PIN (TC-PIN) to detect protein complexes. The performance of each algorithm on DPINs outperforms those on other networks in terms of matching with known complexes, sensitivity, specificity, f-measure, and accuracy. Furthermore, the statistics of three-sigma principle show that 23-45% proteins are active at a time point and most proteins are active in about half of cellular cycle. In addition, we find 94% essential proteins are in the group of proteins that are active at equal or great than 12 timepoints of GSE4987, which indicates the potential existence of feedback mechanisms that can stabilize the expression level of essential proteins and might provide a new insight for predicting essential proteins from dynamic protein networks.

Dissecting the Characteristics and Dynamics of Human Protein Complexes at Transcriptome Cascade Using RNA-Seq Data

Human protein complexes play crucial roles in various biological processes as the functional module. However, the expression features of human protein complexes at the transcriptome cascade are poorly understood. Here, we used the RNA-Seq data from 16 disparate tissues and four types of human cancers to explore the characteristics and dynamics of human protein complexes. We observed that many individual components of human protein complexes can be generated by multiple distinct transcripts. Similar with yeast, the human protein complex constituents are inclined to co-express in diverse tissues. The dominant isoform of the genes involved in protein complexes tend to encode the complex constituents in each tissue. Our results indicate that the protein complex dynamics not only correlate with the presence or absence of complexes, but may also be related to the major isoform switching for complex subunits. Between any two cancers of breast, colon, lung and prostate, we found that only a few of the differentially expressed transcripts associated with complexes were identical, but 5-10 times more protein complexes involved in differentially expressed transcripts were common. Collectively, our study reveals novel properties and dynamics of human protein complexes at the transcriptome cascade in diverse normal tissues and different cancers.

Inference of dynamic networks using time-course data

Cells execute their functions through dynamic operations of biological networks. Dynamic networks delineate the operation of biological networks in terms of temporal changes of abundances or activities of nodes (proteins and RNAs), as well as formation of new edges and disappearance of existing edges over time. Global genomic and proteomic technologies can be used to decode dynamic networks. However, using these experimental methods, it is still challenging to identify temporal transition of nodes and edges. Thus, several computational methods for estimating dynamic topological and functional characteristics of networks have been introduced. In this review, we summarize concepts and applications of these computational methods for inferring dynamic networks and further summarize methods for estimating spatial transition of biological networks.

Visualizing time-related data in biology, a review

Time is of the essence in biology as in so much else. For example, monitoring disease progression or the timing of developmental defects is important for the processes of drug discovery and therapy trials. Furthermore, an understanding of the basic dynamics of biological phenomena that are often strictly time regulated (e.g. circadian rhythms) is needed to make accurate inferences about the evolution of biological processes. Recent advances in technologies have enabled us to measure timing effects more accurately and in more detail. This has driven related advances in visualization and analysis tools that try to effectively exploit this data. Beyond timeline plots, notable attempts at more involved temporal interpretation have been made in recent years, but awareness of the available resources is still limited within the scientific community. Here, we review some advances in biological visualization of time-driven processes and consider how they aid data analysis and interpretation.

Systems biology for molecular life sciences and its impact in biomedicine

Modern systems biology is already contributing to a radical transformation of molecular life sciences and biomedicine, and it is expected to have a real impact in the clinical setting in the next years. In this review, the emergence of systems biology is contextualized with a historic overview, and its present state is depicted. The present and expected future contribution of systems biology to the development of molecular medicine is underscored. Concerning the present situation, this review includes a reflection on the "inflation" of biological data and the urgent need for tools and procedures to make hidden information emerge. Descriptions of the impact of networks and models and the available resources and tools for applying them in systems biology approaches to molecular medicine are provided as well. The actual current impact of systems biology in molecular medicine is illustrated, reviewing two cases, namely, those of systems pharmacology and cancer systems biology. Finally, some of the expected contributions of systems biology to the immediate future of molecular medicine are commented.

Temporal dynamics of protein complexes in PPI networks: a case study using yeast cell cycle dynamics

Complexes of physically interacting proteins are one of the fundamental functional units responsible for driving key biological mechanisms within the cell. With the advent of high-throughput techniques, significant amount of protein interaction (PPI) data has been catalogued for organisms such as yeast, which has in turn fueled computational methods for systematic identification and study of protein complexes. However, many complexes are dynamic entities - their subunits are known to assemble at a particular cellular space and time to perform a particular function and disassemble after that - and while current computational analyses have concentrated on studying the dynamics of individual or pairs of proteins in PPI networks, a crucial aspect overlooked is the dynamics of whole complex formations. In this work, using yeast as our model, we incorporate 'time' in the form of cell-cycle phases into the prediction of complexes from PPI networks and study the temporal phenomena of complex assembly and disassembly across phases. We hypothesize that 'staticness' (constitutive expression) of proteins might be related to their temporal "reusability" across complexes, and test this hypothesis using complexes predicted from large-scale PPI networks across the yeast cell cycle phases. Our results hint towards a biological design principle underlying cellular mechanisms - cells maintain generic proteins as 'static' to enable their "reusability" across multiple temporal complexes. We also demonstrate that these findings provide additional support and alternative explanations to findings from existing works on the dynamics in PPI networks.

Systematic investigation of protein-small molecule interactions

Cell signaling is extensively wired between cellular components to sustain cell proliferation, differentiation, and adaptation. The interaction network is often manifested in how protein function is regulated through interacting with other cellular components including small molecule metabolites. While many biochemical interactions have been established as reactions between protein enzymes and their substrates and products, much less is known at the system level about how small metabolites regulate protein functions through allosteric binding. In the past decade, study of protein-small molecule interactions has been lagging behind other types of interactions. Recent technological advances have explored several high-throughput platforms to reveal many "unexpected" protein-small molecule interactions that could have profound impact on our understanding of cell signaling. These interactions will help bridge gaps in existing regulatory loops of cell signaling and serve as new targets for medical intervention. In this review, we summarize recent advances of systematic investigation of protein-metabolite/small molecule interactions, and discuss the impact of such studies and their potential impact on both biological researches and medicine.

A survey of computational methods for protein complex prediction from protein interaction networks

Complexes of physically interacting proteins are one of the fundamental functional units responsible for driving key biological mechanisms within the cell. Their identification is therefore necessary to understand not only complex formation but also the higher level organization of the cell. With the advent of "high-throughput" techniques in molecular biology, significant amount of physical interaction data has been cataloged from organisms such as yeast, which has in turn fueled computational approaches to systematically mine complexes from the network of physical interactions among proteins (PPI network). In this survey, we review, classify and evaluate some of the key computational methods developed till date for the identification of protein complexes from PPI networks. We present two insightful taxonomies that reflect how these methods have evolved over the years toward improving automated complex prediction. We also discuss some open challenges facing accurate reconstruction of complexes, the crucial ones being the presence of high proportion of errors and noise in current high-throughput datasets and some key aspects overlooked by current complex detection methods. We hope this review will not only help to condense the history of computational complex detection for easy reference but also provide valuable insights to drive further research in this area.

Biological function through network topology: a survey of the human diseasome

Molecular network data are increasingly becoming available, necessitating the development of well performing computational tools for their analyses. Such tools enabled conceptually different approaches for exploring human diseases to be undertaken, in particular, those that study the relationship between a multitude of biomolecules within a cell. Hence, a new field of network biology has emerged as part of systems biology, aiming to untangle the complexity of cellular network organization. We survey current network analysis methods that aim to give insight into human disease.

Dynamic hubs show competitive and static hubs non-competitive regulation of their interaction partners

Date hub proteins have 1 or 2 interaction interfaces but many interaction partners. This raises the question of whether all partner proteins compete for the interaction interface of the hub or if the cell carefully regulates aspects of this process? Here, we have used real-time rendering of protein interaction networks to analyse the interactions of all the 1 or 2 interface hubs of Saccharomyces cerevisiae during the cell cycle. By integrating previously determined structural and gene expression data, and visually hiding the nodes (proteins) and their edges (interactions) during their troughs of expression, we predict when interactions of hubs and their partners are likely to exist. This revealed that 20 out of all 36 one- or two- interface hubs in the yeast interactome fell within two main groups. The first was dynamic hubs with static partners, which can be considered as ‘competitive hubs’. Their interaction partners will compete for the interaction interface of the hub and the success of any interaction will be dictated by the kinetics of interaction (abundance and affinity) and subcellular localisation. The second was static hubs with dynamic partners, which we term ‘non-competitive hubs’. Regulatory mechanisms are finely tuned to lessen the presence and/or effects of competition between the interaction partners of the hub. It is possible that these regulatory processes may also be used by the cell for the regulation of other, non-cell cycle processes.

Identification of markers associated with global changes in DNA methylation regulation in cancers

DNA methylation exhibits different patterns in different cancers. DNA methylation rates at different genomic loci appear to be highly correlated in some samples but not in others. We call such phenomena conditional concordant relationships (CCRs). In this study, we explored DNA methylation patterns in 12 common cancers using data of 2434 patient samples collected by The Cancer Genome Atlas project. We developed an exploratory method to characterize CCRs in the methylation data and identified the 200 gene markers whose on-and-off statuses in DNA methylation are most significantly associated with drastic changes in CCRs throughout the genome. Clustering analysis of the methylation data of the 200 markers showed that they are tightly associated with cancer subtypes. We also generated a library of the significant CCRs that may be of interest to future studies of the regulation network of DNA methylation in cancer.

Static network structure can be used to model the phenotypic effects of perturbations in regulatory networks

MOTIVATION

Biological processes are dynamic, whereas the networks that depict them are typically static. Quantitative modeling using differential equations or logic-based functions can offer quantitative predictions of the behavior of biological systems, but they require detailed experimental characterization of interaction kinetics, which is typically unavailable. To determine to what extent complex biological processes can be modeled and analyzed using only the static structure of the network (i.e. the direction and sign of the edges), we attempt to predict the phenotypic effect of perturbations in biological networks from the static network structure.

RESULTS

We analyzed three networks from different sources: The EGFR/MAPK and PI3K/AKT network from a detailed experimental study, the TNF regulatory network from the STRING database and a large network of all NCI-curated pathways from the Protein Interaction Database. Altogether, we predicted the effect of 39 perturbations (e.g. by one or two drugs) on 433 target proteins/genes. In up to 82% of the cases, an algorithm that used only the static structure of the network correctly predicted whether any given protein/gene is upregulated or downregulated as a result of perturbations of other proteins/genes.

CONCLUSION

While quantitative modeling requires detailed experimental data and heavy computations, which limit its scalability for large networks, a wiring-based approach can use available data from pathway and interaction databases and may be scalable. These results lay the foundations for a large-scale approach of predicting phenotypes based on the schematic structure of networks.

Protein-protein interaction networks: unraveling the wiring of molecular machines within the cell

Mapping and understanding of the protein interaction networks with their key modules and hubs can provide deeper insights into the molecular machinery underlying complex phenotypes. In this article, we present the basic characteristics and definitions of protein networks, starting with a distinction of the different types of associations between proteins. We focus the review on protein-protein interactions (PPIs), a subset of associations defined as physical contacts between proteins that occur by selective molecular docking in a particular biological context. We present such definition as opposed to other types of protein associations derived from regulatory, genetic, structural or functional relations. To determine PPIs, a variety of binary and co-complex methods exist; however, not all the technologies provide the same information and data quality. A way of increasing confidence in a given protein interaction is to integrate orthogonal experimental evidences. The use of several complementary methods testing each single interaction assesses the accuracy of PPI data and tries to minimize the occurrence of false interactions. Following this approach there have been important efforts to unify primary databases of experimentally proven PPIs into integrated databases. These meta-databases provide a measure of the confidence of interactions based on the number of experimental proofs that report them. As a conclusion, we can state that integrated information allows the building of more reliable interaction networks. Identification of communities, cliques, modules and hubs by analysing the topological parameters and graph properties of the protein networks allows the discovery of central/critical nodes, which are candidates to regulate cellular flux and dynamics.

Towards the identification of protein complexes and functional modules by integrating PPI network and gene expression data

Background

Identification of protein complexes and functional modules from protein-protein interaction (PPI) networks is crucial to understanding the principles of cellular organization and predicting protein functions. In the past few years, many computational methods have been proposed. However, most of them considered the PPI networks as static graphs and overlooked the dynamics inherent within these networks. Moreover, few of them can distinguish between protein complexes and functional modules.

Results

In this paper, a new framework is proposed to distinguish between protein complexes and functional modules by integrating gene expression data into protein-protein interaction (PPI) data. A series of time-sequenced subnetworks (TSNs) is constructed according to the time that the interactions were activated. The algorithm TSN-PCD was then developed to identify protein complexes from these TSNs. As protein complexes are significantly related to functional modules, a new algorithm DFM-CIN is proposed to discover functional modules based on the identified complexes. The experimental results show that the combination of temporal gene expression data with PPI data contributes to identifying protein complexes more precisely. A quantitative comparison based on f-measure reveals that our algorithm TSN-PCD outperforms the other previous protein complex discovery algorithms. Furthermore, we evaluate the identified functional modules by using “Biological Process” annotated in GO (Gene Ontology). The validation shows that the identified functional modules are statistically significant in terms of “Biological Process”. More importantly, the relationship between protein complexes and functional modules are studied.

Conclusions

The proposed framework based on the integration of PPI data and gene expression data makes it possible to identify protein complexes and functional modules more effectively. Moveover, the proposed new framework and algorithms can distinguish between protein complexes and functional modules. Our findings suggest that functional modules are closely related to protein complexes and a functional module may consist of one or multiple protein complexes. The program is available at http://netlab.csu.edu.cn/bioinfomatics/limin/DFM-CIN/index.html.

Constructing structural networks of signaling pathways on the proteome scale

Proteins function through their interactions, and the availability of protein interaction networks could help in understanding cellular processes. However, the known structural data are limited and the classical network node-and-edge representation, where proteins are nodes and interactions are edges, shows only which proteins interact; not how they interact. Structural networks provide this information. Protein-protein interface structures can also indicate which binding partners can interact simultaneously and which are competitive, and can help forecasting potentially harmful drug side effects. Here, we use a powerful protein-protein interactions prediction tool which is able to carry out accurate predictions on the proteome scale to construct the structural network of the extracellular signal-regulated kinases (ERK) in the mitogen-activated protein kinase (MAPK) signaling pathway. This knowledge-based method, PRISM, is motif-based, and is combined with flexible refinement and energy scoring. PRISM predicts protein interactions based on structural and evolutionary similarity to known protein interfaces.

Topological analysis and interactive visualization of biological networks and protein structures

Computational analysis and interactive visualization of biological networks and protein structures are common tasks for gaining insight into biological processes. This protocol describes three workflows based on the NetworkAnalyzer and RINalyzer plug-ins for Cytoscape, a popular software platform for networks. NetworkAnalyzer has become a standard Cytoscape tool for comprehensive network topology analysis. In addition, RINalyzer provides methods for exploring residue interaction networks derived from protein structures. The first workflow uses NetworkAnalyzer to perform a topological analysis of biological networks. The second workflow applies RINalyzer to study protein structure and function and to compute network centrality measures. The third workflow combines NetworkAnalyzer and RINalyzer to compare residue networks. The full protocol can be completed in ∼2 h.

Protein interactions: mapping interactome networks to support drug target discovery and selection

Proteins are biomolecular structures that build the microscopic working machinery of any living system. Proteins within the cells and biological systems do not act alone, but rather team up into macromolecular structures enclosing intricate physicochemical dynamic connections to undertake biological functions. A critical step towards unraveling the complex molecular relationships in living systems is the mapping of protein-to-protein physical "interactions". The complete map of protein interactions that can occur in a living organism is called the "interactome". Achieving an adequate atlas of all the protein interactions within a living system should allow to build its interaction network and to identity the "central nodes" that can be critical for the function, the homeostasis, and the movement of such system. Focusing on human studies, the data about the human interactome are most relevant for current biomedical research, because it is clear that the location of the proteins in the interactome network will allow to evaluate their centrality and to redefine the potential value of each protein as a drug target. This chapter presents our current knowledge on the human protein-protein interactome and explains how such knowledge can help us to select adequate targets for drugs.

Comparative interaction networks: bridging genotype to phenotype

Over the past decade, biomedical research has witnessed an exponential increase in the throughput of the characterization of biological systems. Here we review the recent progress in large-scale methods to determine protein-protein, genetic and chemical-genetic interaction networks. We discuss some of the limitations and advantages of the different methods and give examples of how these networks are being used to study the evolutionary process. Comparative studies have revealed that different types of protein-protein interactions diverge at different rates with high conservation of co-complex membership but rapid divergence of more promiscuous interactions like those that mediate post-translational modifications. These evolutionary trends have consistent genetic consequences with highly conserved epistatic interactions within complex subunits but faster divergence of epistatic interactions across complexes or pathways. Finally, we discuss how these evolutionary observations are being used to interpret cross-species chemical-genetic studies and how they might shape therapeutic strategies. Together, these interaction networks offer us an unprecedented level of detail into how genotypes are translated to phenotypes, and we envision that they will be increasingly useful in the interpretation of genetic and phenotypic variation occurring within populations as well as the rational design of combinatorial therapeutics.

Network systems biology for targeted cancer therapies

The era of targeted cancer therapies has arrived. However, due to the complexity of biological systems, the current progress is far from enough. From biological network modeling to structural/dynamic network analysis, network systems biology provides unique insight into the potential mechanisms underlying the growth and progression of cancer cells. It has also introduced great changes into the research paradigm of cancer-associated drug discovery and drug resistance.

Heat shock partially dissociates the overlapping modules of the yeast protein-protein interaction network: a systems level model of adaptation

Network analysis became a powerful tool giving new insights to the understanding of cellular behavior. Heat shock, the archetype of stress responses, is a well-characterized and simple model of cellular dynamics. S. cerevisiae is an appropriate model organism, since both its protein-protein interaction network (interactome) and stress response at the gene expression level have been well characterized. However, the analysis of the reorganization of the yeast interactome during stress has not been investigated yet. We calculated the changes of the interaction-weights of the yeast interactome from the changes of mRNA expression levels upon heat shock. The major finding of our study is that heat shock induced a significant decrease in both the overlaps and connections of yeast interactome modules. In agreement with this the weighted diameter of the yeast interactome had a 4.9-fold increase in heat shock. Several key proteins of the heat shock response became centers of heat shock-induced local communities, as well as bridges providing a residual connection of modules after heat shock. The observed changes resemble to a 'stratus-cumulus' type transition of the interactome structure, since the unstressed yeast interactome had a globally connected organization, similar to that of stratus clouds, whereas the heat shocked interactome had a multifocal organization, similar to that of cumulus clouds. Our results showed that heat shock induces a partial disintegration of the global organization of the yeast interactome. This change may be rather general occurring in many types of stresses. Moreover, other complex systems, such as single proteins, social networks and ecosystems may also decrease their inter-modular links, thus develop more compact modules, and display a partial disintegration of their global structure in the initial phase of crisis. Thus, our work may provide a model of a general, system-level adaptation mechanism to environmental changes.

Categorizing biases in high-confidence high-throughput protein-protein interaction data sets

We characterized and evaluated the functional attributes of three yeast high-confidence protein-protein interaction data sets derived from affinity purification/mass spectrometry, protein-fragment complementation assay, and yeast two-hybrid experiments. The interacting proteins retrieved from these data sets formed distinct, partially overlapping sets with different protein-protein interaction characteristics. These differences were primarily a function of the deployed experimental technologies used to recover these interactions. This affected the total coverage of interactions and was especially evident in the recovery of interactions among different functional classes of proteins. We found that the interaction data obtained by the yeast two-hybrid method was the least biased toward any particular functional characterization. In contrast, interacting proteins in the affinity purification/mass spectrometry and protein-fragment complementation assay data sets were over- and under-represented among distinct and different functional categories. We delineated how these differences affected protein complex organization in the network of interactions, in particular for strongly interacting complexes (e.g. RNA and protein synthesis) versus weak and transient interacting complexes (e.g. protein transport). We quantified methodological differences in detecting protein interactions from larger protein complexes, in the correlation of protein abundance among interacting proteins, and in their connectivity of essential proteins. In the latter case, we showed that minimizing inherent methodology biases removed many of the ambiguous conclusions about protein essentiality and protein connectivity. We used these findings to rationalize how biological insights obtained by analyzing data sets originating from different sources sometimes do not agree or may even contradict each other. An important corollary of this work was that discrepancies in biological insights did not necessarily imply that one detection methodology was better or worse, but rather that, to a large extent, the insights reflected the methodological biases themselves. Consequently, interpreting the protein interaction data within their experimental or cellular context provided the best avenue for overcoming biases and inferring biological knowledge.

ppiTrim: constructing non-redundant and up-to-date interactomes

Robust advances in interactome analysis demand comprehensive, non-redundant and consistently annotated data sets. By non-redundant, we mean that the accounting of evidence for every interaction should be faithful: each independent experimental support is counted exactly once, no more, no less. While many interactions are shared among public repositories, none of them contains the complete known interactome for any model organism. In addition, the annotations of the same experimental result by different repositories often disagree. This brings up the issue of which annotation to keep while consolidating evidences that are the same. The iRefIndex database, including interactions from most popular repositories with a standardized protein nomenclature, represents a significant advance in all aspects, especially in comprehensiveness. However, iRefIndex aims to maintain all information/annotation from original sources and requires users to perform additional processing to fully achieve the aforementioned goals. Another issue has to do with protein complexes. Some databases represent experimentally observed complexes as interactions with more than two participants, while others expand them into binary interactions using spoke or matrix model. To avoid untested interaction information buildup, it is preferable to replace the expanded protein complexes, either from spoke or matrix models, with a flat list of complex members.

To address these issues and to achieve our goals, we have developed ppiTrim, a script that processes iRefIndex to produce non-redundant, consistently annotated data sets of physical interactions. Our script proceeds in three stages: mapping all interactants to gene identifiers and removing all undesired raw interactions, deflating potentially expanded complexes, and reconciling for each interaction the annotation labels among different source databases. As an illustration, we have processed the three largest organismal data sets: yeast, human and fruitfly. While ppiTrim can resolve most apparent conflicts between different labelings, we also discovered some unresolvable disagreements mostly resulting from different annotation policies among repositories.

Database URL:http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads/ppiTrim.html

A comparison of the functional modules identified from time course and static PPI network data

BACKGROUND

Cellular systems are highly dynamic and responsive to cues from the environment. Cellular function and response patterns to external stimuli are regulated by biological networks. A protein-protein interaction (PPI) network with static connectivity is dynamic in the sense that the nodes implement so-called functional activities that evolve in time. The shift from static to dynamic network analysis is essential for further understanding of molecular systems.

RESULTS

In this paper, Time Course Protein Interaction Networks (TC-PINs) are reconstructed by incorporating time series gene expression into PPI networks. Then, a clustering algorithm is used to create functional modules from three kinds of networks: the TC-PINs, a static PPI network and a pseudorandom network. For the functional modules from the TC-PINs, repetitive modules and modules contained within bigger modules are removed. Finally, matching and GO enrichment analyses are performed to compare the functional modules detected from those networks.

CONCLUSIONS

The comparative analyses show that the functional modules from the TC-PINs have much more significant biological meaning than those from static PPI networks. Moreover, it implies that many studies on static PPI networks can be done on the TC-PINs and accordingly, the experimental results are much more satisfactory. The 36 PPI networks corresponding to 36 time points, identified as part of this study, and other materials are available at http://bioinfo.csu.edu.cn/txw/TC-PINs.

Comparing signaling networks between normal and transformed hepatocytes using discrete logical models

Substantial effort in recent years has been devoted to constructing and analyzing large-scale gene and protein networks on the basis of "omic" data and literature mining. These interaction graphs provide valuable insight into the topologies of complex biological networks but are rarely context specific and cannot be used to predict the responses of cell signaling proteins to specific ligands or drugs. Conversely, traditional approaches to analyzing cell signaling are narrow in scope and cannot easily make use of network-level data. Here, we combine network analysis and functional experimentation by using a hybrid approach in which graphs are converted into simple mathematical models that can be trained against biochemical data. Specifically, we created Boolean logic models of immediate-early signaling in liver cells by training a literature-based prior knowledge network against biochemical data obtained from primary human hepatocytes and 4 hepatocellular carcinoma cell lines exposed to combinations of cytokines and small-molecule kinase inhibitors. Distinct families of models were recovered for each cell type, and these families clustered topologically into normal and diseased sets.

Brain Hodotopy: From Esoteric Concept to Practical Surgical Applications

BACKGROUND:

The traditional neurosurgical approach to cerebral lesions is based on the classic view of a rigid brain organization in fixed “eloquent” areas. However, this method is brought into discussion by the conceptual and methodological advances in neurosciences that provide a more dynamic representation of the anatomo-functional distribution of the human central nervous system (CNS).

OBJECTIVE AND METHODS:

We review the relevant literature concerning the main features of the modern CNS representation and their implications in neurosurgical practice.

RESULTS:

The CNS is an integrated, wide, plastic network made up of cortical functional epicenters, “topic organization,” connected by both short-local and large-scale white matter fibers, ie, “hodological organization.” According to this model, called hodotopic, brain function results from parallel streams of information dynamically modulated within an interactive, multimodal, and widely distributed circuit. The application of this framework, which can be studied by combining preoperative, intraoperative, and postoperative mapping techniques, enables the neurosurgeon exploration of the individual anatomo-functional architecture, including neurocognitive and emotional aspects. Thus, it is possible to adapt the surgical approach specifically to each patient and to each lesion according to the individual organization. Several experiences demonstrate the possibility of removing regions traditionally considered inoperable without inducing permanent deficits and the potential use of these areas as a safe passage to deeper territories.

CONCLUSION:

We advocate the more systematic integration of a hodotopical view of the CNS to improve the surgical indications and planning for brain lesions, with the goal of optimizing both the extent of resection and functional outcome.

From in silico target prediction to multi-target drug design: current databases, methods and applications

Given the tremendous growth of bioactivity databases, the use of computational tools to predict protein targets of small molecules has been gaining importance in recent years. Applications span a wide range, from the 'designed polypharmacology' of compounds to mode-of-action analysis. In this review, we firstly survey databases that can be used for ligand-based target prediction and which have grown tremendously in size in the past. We furthermore outline methods for target prediction that exist, both based on the knowledge of bioactivities from the ligand side and methods that can be applied in situations when a protein structure is known. Applications of successful in silico target identification attempts are discussed in detail, which were based partly or in whole on computational target predictions in the first instance. This includes the authors' own experience using target prediction tools, in this case considering phenotypic antibacterial screens and the analysis of high-throughput screening data. Finally, we will conclude with the prospective application of databases to not only predict, retrospectively, the protein targets of a small molecule, but also how to design ligands with desired polypharmacology in a prospective manner.

On the role of electrostatics in protein-protein interactions

The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

Prediction of protein-protein interactions: unifying evolution and structure at protein interfaces

The vast majority of the chores in the living cell involve protein-protein interactions. Providing details of protein interactions at the residue level and incorporating them into protein interaction networks are crucial toward the elucidation of a dynamic picture of cells. Despite the rapid increase in the number of structurally known protein complexes, we are still far away from a complete network. Given experimental limitations, computational modeling of protein interactions is a prerequisite to proceed on the way to complete structural networks. In this work, we focus on the question 'how do proteins interact?' rather than 'which proteins interact?' and we review structure-based protein-protein interaction prediction approaches. As a sample approach for modeling protein interactions, PRISM is detailed which combines structural similarity and evolutionary conservation in protein interfaces to infer structures of complexes in the protein interaction network. This will ultimately help us to understand the role of protein interfaces in predicting bound conformations.

IsoBase: a database of functionally related proteins across PPI networks

We describe IsoBase, a database identifying functionally related proteins, across five major eukaryotic model organisms: Saccharomyces cerevisiae, Drosophila melanogaster, Caenorhabditis elegans, Mus musculus and Homo Sapiens. Nearly all existing algorithms for orthology detection are based on sequence comparison. Although these have been successful in orthology prediction to some extent, we seek to go beyond these methods by the integration of sequence data and protein–protein interaction (PPI) networks to help in identifying true functionally related proteins. With that motivation, we introduce IsoBase, the first publicly available ortholog database that focuses on functionally related proteins. The groupings were computed using the IsoRankN algorithm that uses spectral methods to combine sequence and PPI data and produce clusters of functionally related proteins. These clusters compare favorably with those from existing approaches: proteins within an IsoBase cluster are more likely to share similar Gene Ontology (GO) annotation. A total of 48 120 proteins were clustered into 12 693 functionally related groups. The IsoBase database may be browsed for functionally related proteins across two or more species and may also be queried by accession numbers, species-specific identifiers, gene name or keyword. The database is freely available for download at http://isobase.csail.mit.edu/.

Dynamics of cullin-RING ubiquitin ligase network revealed by systematic quantitative proteomics

Dynamic reorganization of signaling systems frequently accompanies pathway perturbations, yet quantitative studies of network remodeling by pathway stimuli are lacking. Here, we report the development of a quantitative proteomics platform centered on multiplex absolute quantification (AQUA) technology to elucidate the architecture of the cullin-RING ubiquitin ligase (CRL) network and to evaluate current models of dynamic CRL remodeling. Current models suggest that CRL complexes are controlled by cycles of CRL deneddylation and CAND1 binding. Contrary to expectations, acute CRL inhibition with MLN4924, an inhibitor of the NEDD8-activating enzyme, does not result in a global reorganization of the CRL network. Examination of CRL complex stoichiometry reveals that, independent of cullin neddylation, a large fraction of cullins are assembled with adaptor modules, whereas only a small fraction are associated with CAND1. These studies suggest an alternative model of CRL dynamicity where the abundance of adaptor modules, rather than cycles of neddylation and CAND1 binding, drives CRL network organization.

Network biology: a direct approach to study biological function

In this paper we discuss the dualism of gene networks and their role in systems biology. We argue that gene networks (1) can serve as a conceptual framework, forming a fundamental level of a phenomenological description, and (2) are a means to represent and analyze data. The latter point does not only allow a systems analysis but is even amenable for a direct approach to study biological function. Here we focus on the clarity of our main arguments and conceptual meaning of gene networks, rather than the causal inference of gene networks from data. WIREs Syst Biol Med 2011 3 379-391 DOI: 10.1002/wsbm.134 For further resources related to this article, please visit the WIREs website.

Stability in flux: community structure in dynamic networks

The structure of many biological, social and technological systems can usefully be described in terms of complex networks. Although often portrayed as fixed in time, such networks are inherently dynamic, as the edges that join nodes are cut and rewired, and nodes themselves update their states. Understanding the structure of these networks requires us to understand the dynamic processes that create, maintain and modify them. Here, we build upon existing models of coevolving networks to characterize how dynamic behaviour at the level of individual nodes generates stable aggregate behaviours. We focus particularly on the dynamics of groups of nodes formed endogenously by nodes that share similar properties (represented as node state) and demonstrate that, under certain conditions, network modularity based on state compares well with network modularity based on topology. We show that if nodes rewire their edges based on fixed node states, the network modularity reaches a stable equilibrium which we quantify analytically. Furthermore, if node state is not fixed, but can be adopted from neighbouring nodes, the distribution of group sizes reaches a dynamic equilibrium, which remains stable even as the composition and identity of the groups change. These results show that dynamic networks can maintain the stable community structure that has been observed in many social and biological systems.

Cross-platform microarray data normalisation for regulatory network inference

BACKGROUND

Inferring Gene Regulatory Networks (GRNs) from time course microarray data suffers from the dimensionality problem created by the short length of available time series compared to the large number of genes in the network. To overcome this, data integration from diverse sources is mandatory. Microarray data from different sources and platforms are publicly available, but integration is not straightforward, due to platform and experimental differences.

METHODS

We analyse here different normalisation approaches for microarray data integration, in the context of reverse engineering of GRN quantitative models. We introduce two preprocessing approaches based on existing normalisation techniques and provide a comprehensive comparison of normalised datasets.

CONCLUSIONS

Results identify a method based on a combination of Loess normalisation and iterative K-means as best for time series normalisation for this problem.

Networks inferred from biochemical data reveal profound differences in toll-like receptor and inflammatory signaling between normal and transformed hepatocytes

Systematic study of cell signaling networks increasingly involves high throughput proteomics, transcriptional profiling, and automated literature mining with the aim of assembling large scale interaction networks. In contrast, functional analysis of cell signaling usually focuses on a much smaller sets of proteins and eschews computation but focuses directly on cellular responses to environment and perturbation. We sought to combine these two traditions by collecting cell response measures on a reasonably large scale and then attempting to infer differences in network topology between two cell types. Human hepatocytes and hepatocellular carcinoma cell lines were exposed to inducers of inflammation, innate immunity, and proliferation in the presence and absence of small molecule drugs, and multiplex biochemical measurement was then performed on intra- and extracellular signaling molecules. We uncovered major differences between primary and transformed hepatocytes with respect to the engagement of toll-like receptor and NF-κB-dependent secretion of chemokines and cytokines that prime and attract immune cells. Overall, our results serve as a proof of principle for an approach to network analysis that is systematic, comparative, and biochemically focused. More specifically, our data support the hypothesis that hepatocellular carcinoma cells down-regulate normal inflammatory and immune responses to avoid immune editing.

Network integration meets network dynamics

Molecular interaction networks provide a window on the workings of the cell. However, combining various types of networks into one coherent large-scale dynamic model remains a formidable challenge. A recent paper in BMC Systems Biology describes a promising step in this direction.