SAMNet: a network-based approach to integrate multi-dimensional high throughput datasets

Integration of Transcriptomics, Proteomics and Loss-of-function Screening Reveals WEE1 as a Target for Combination with Dasatinib against Proneural Glioblastoma

Glioblastoma is characterized by a pronounced resistance to therapy with dismal prognosis. Transcriptomics classify glioblastoma into proneural (PN), mesenchymal (MES) and classical (CL) subtypes that show differential resistance to targeted therapies. The aim of this study was to provide a viable approach for identifying combination therapies in glioblastoma subtypes. Proteomics and phosphoproteomics were performed on dasatinib inhibited glioblastoma stem cells (GSCs) and complemented by an shRNA loss-of-function screen to identify genes whose knockdown sensitizes GSCs to dasatinib. Proteomics and screen data were computationally integrated with transcriptomic data using the SamNet 2.0 algorithm for network flow learning to reveal potential combination therapies in PN GSCs. In vitro viability assays and tumor spheroid models were used to verify the synergy of identified therapy. Further in vitro and TCGA RNA-Seq data analyses were utilized to provide a mechanistic explanation of these effects. Integration of data revealed the cell cycle protein WEE1 as a potential combination therapy target for PN GSCs. Validation experiments showed a robust synergistic effect through combination of dasatinib and the WEE1 inhibitor, MK-1775, in PN GSCs. Combined inhibition using dasatinib and MK-1775 propagated DNA damage in PN GCSs, with GCSs showing a differential subtype-driven pattern of expression of cell cycle genes in TCGA RNA-Seq data. The integration of proteomics, loss-of-function screens and transcriptomics confirmed WEE1 as a target for combination with dasatinib against PN GSCs. Utilizing this integrative approach could be of interest for studying resistance mechanisms and revealing combination therapy targets in further tumor entities.

Current and future directions in network biology

Summary

Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These stem from various factors, notably the growing complexity and volume of data together with the increased diversity of data types describing different tiers of biological organization. We discuss prevailing research directions in network biology, focusing on molecular/cellular networks but also on other biological network types such as biomedical knowledge graphs, patient similarity networks, brain networks, and social/contact networks relevant to disease spread. In more detail, we highlight areas of inference and comparison of biological networks, multimodal data integration and heterogeneous networks, higher-order network analysis, machine learning on networks, and network-based personalized medicine. Following the overview of recent breakthroughs across these five areas, we offer a perspective on future directions of network biology. Additionally, we discuss scientific communities, educational initiatives, and the importance of fostering diversity within the field. This article establishes a roadmap for an immediate and long-term vision for network biology.

Availability and implementation

Not applicable.

Selective CDK7 Inhibition Suppresses Cell Cycle Progression and MYC Signaling While Enhancing Apoptosis in Therapy-resistant Estrogen Receptor-positive Breast Cancer

PURPOSE

Resistance to endocrine therapy (ET) and CDK4/6 inhibitors (CDK4/6i) is a clinical challenge in estrogen receptor (ER)-positive (ER+) breast cancer. Cyclin-dependent kinase 7 (CDK7) is a candidate target in endocrine-resistant ER+ breast cancer models and selective CDK7 inhibitors (CDK7i) are in clinical development for the treatment of ER+ breast cancer. Nonetheless, the precise mechanisms responsible for the activity of CDK7i in ER+ breast cancer remain elusive. Herein, we sought to unravel these mechanisms.

EXPERIMENTAL DESIGN

We conducted multi-omic analyses in ER+ breast cancer models in vitro and in vivo, including models with different genetic backgrounds. We also performed genome-wide CRISPR/Cas9 knockout screens to identify potential therapeutic vulnerabilities in CDK4/6i-resistant models.

RESULTS

We found that the on-target antitumor effects of CDK7 inhibition in ER+ breast cancer are in part p53 dependent, and involve cell cycle inhibition and suppression of c-Myc. Moreover, CDK7 inhibition exhibited cytotoxic effects, distinctive from the cytostatic nature of ET and CDK4/6i. CDK7 inhibition resulted in suppression of ER phosphorylation at S118; however, long-term CDK7 inhibition resulted in increased ER signaling, supporting the combination of ET with a CDK7i. Finally, genome-wide CRISPR/Cas9 knockout screens identified CDK7 and MYC signaling as putative vulnerabilities in CDK4/6i resistance, and CDK7 inhibition effectively inhibited CDK4/6i-resistant models.

CONCLUSIONS

Taken together, these findings support the clinical investigation of selective CDK7 inhibition combined with ET to overcome treatment resistance in ER+ breast cancer. In addition, our study highlights the potential of increased c-Myc activity and intact p53 as predictors of sensitivity to CDK7i-based treatments.

The use of machine learning to discover regulatory networks controlling biological systems

Biological systems are composed of a vast web of multiscale molecular interactors and interactions. High-throughput technologies, both bulk and single cell, now allow for investigation of the properties and quantities of these interactors. Computational algorithms and machine learning methods then provide the tools to derive meaningful insights from the resulting data sets. One such approach is graphical network modeling, which provides a computational framework to explicitly model the molecular interactions within and between the cells comprising biological systems. These graphical networks aim to describe a putative chain of cause and effect between interacting molecules. This feature allows for determination of key molecules in a biological process, accelerated generation of mechanistic hypotheses, and simulation of experimental outcomes. We review the computational concepts and applications of graphical network models across molecular scales for both intracellular and intercellular regulatory biology, examples of successful applications, and the future directions needed to overcome current limitations.

Network modeling of kinase inhibitor polypharmacology reveals pathways targeted in chemical screens

Small molecule screens are widely used to prioritize pharmaceutical development. However, determining the pathways targeted by these molecules is challenging, since the compounds are often promiscuous. We present a network strategy that takes into account the polypharmacology of small molecules in order to generate hypotheses for their broader mode of action. We report a screen for kinase inhibitors that increase the efficacy of gemcitabine, the first-line chemotherapy for pancreatic cancer. Eight kinase inhibitors emerge that are known to affect 201 kinases, of which only three kinases have been previously identified as modifiers of gemcitabine toxicity. In this work, we use the SAMNet algorithm to identify pathways linking these kinases and genetic modifiers of gemcitabine toxicity with transcriptional and epigenetic changes induced by gemcitabine that we measure using DNaseI-seq and RNA-seq. SAMNet uses a constrained optimization algorithm to connect genes from these complementary datasets through a small set of protein-protein and protein-DNA interactions. The resulting network recapitulates known pathways including DNA repair, cell proliferation and the epithelial-to-mesenchymal transition. We use the network to predict genes with important roles in the gemcitabine response, including six that have already been shown to modify gemcitabine efficacy in pancreatic cancer and ten novel candidates. Our work reveals the important role of polypharmacology in the activity of these chemosensitizing agents.

Integrated systems biology analysis of KSHV latent infection reveals viral induction and reliance on peroxisome mediated lipid metabolism

Kaposi’s Sarcoma associated Herpesvirus (KSHV), an oncogenic, human gamma-herpesvirus, is the etiological agent of Kaposi’s Sarcoma the most common tumor of AIDS patients world-wide. KSHV is predominantly latent in the main KS tumor cell, the spindle cell, a cell of endothelial origin. KSHV modulates numerous host cell-signaling pathways to activate endothelial cells including major metabolic pathways involved in lipid metabolism. To identify the underlying cellular mechanisms of KSHV alteration of host signaling and endothelial cell activation, we identified changes in the host proteome, phosphoproteome and transcriptome landscape following KSHV infection of endothelial cells. A Steiner forest algorithm was used to integrate the global data sets and, together with transcriptome based predicted transcription factor activity, cellular networks altered by latent KSHV were predicted. Several interesting pathways were identified, including peroxisome biogenesis. To validate the predictions, we showed that KSHV latent infection increases the number of peroxisomes per cell. Additionally, proteins involved in peroxisomal lipid metabolism of very long chain fatty acids, including ABCD3 and ACOX1, are required for the survival of latently infected cells. In summary, novel cellular pathways altered during herpesvirus latency that could not be predicted by a single systems biology platform, were identified by integrated proteomics and transcriptomics data analysis and when correlated with our metabolomics data revealed that peroxisome lipid metabolism is essential for KSHV latent infection of endothelial cells.

Kaposi’s Sarcoma herpesvirus (KSHV) is the etiologic agent of Kaposi’s Sarcoma, the most common tumor of AIDS patients. KSHV modulates host cell signaling and metabolism to maintain a life-long latent infection. To unravel the underlying cellular mechanisms modulated by KSHV, we used multiple global systems biology platforms to identify and integrate changes in both cellular protein expression and transcription following KSHV infection of endothelial cells, the relevant cell type for KS tumors. The analysis identified several interesting pathways including peroxisome biogenesis. Peroxisomes are small cytoplasmic organelles involved in redox reactions and lipid metabolism. KSHV latent infection increases the number of peroxisomes per cell and proteins involved in peroxisomal lipid metabolism are required for the survival of latently infected cells. In summary, through integration of multiple global systems biology analyses we were able to identify novel pathways that could not be predicted by one platform alone and found that lipid metabolism in a small cytoplasmic organelle is necessary for the survival of latent infection with a herpesvirus.

Molecular mechanisms underlying COPD-muscle dysfunction unveiled through a systems medicine approach

MOTIVATION

Skeletal muscle dysfunction is a systemic effect in one-third of patients with chronic obstructive pulmonary disease (COPD), characterized by high reactive-oxygen-species (ROS) production and abnormal endurance training-induced adaptive changes. However, the role of ROS in COPD remains unclear, not least because of the lack of appropriate tools to study multifactorial diseases.

RESULTS

We describe a discrete model-driven method combining mechanistic and probabilistic approaches to decipher the role of ROS on the activity state of skeletal muscle regulatory network, assessed before and after an 8-week endurance training program in COPD patients and healthy subjects. In COPD, our computational analysis indicates abnormal training-induced regulatory responses leading to defective tissue remodeling and abnormal energy metabolism. Moreover, we identified tnf, insr, inha and myc as key regulators of abnormal training-induced adaptations in COPD. The tnf-insr pair was identified as a promising target for therapeutic interventions. Our work sheds new light on skeletal muscle dysfunction in COPD, opening new avenues for cost-effective therapies. It overcomes limitations of previous computational approaches showing high potential for the study of other multi-factorial diseases such as diabetes or cancer.

CONTACT

jroca@clinic.ub.es or martacascante@ub.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Pathway-based network modeling finds hidden genes in shRNA screen for regulators of acute lymphoblastic leukemia

Data integration stands to improve interpretation of RNAi screens which, as a result of off-target effects, typically yield numerous gene hits of which only a few validate. These off-target effects can result from seed matches to unintended gene targets (reagent-based) or cellular pathways, which can compensate for gene perturbations (biology-based). We focus on the biology-based effects and use network modeling tools to discover pathways de novo around RNAi hits. By looking at hits in a functional context, we can uncover novel biology not identified from any individual 'omics measurement. We leverage multiple 'omic measurements using the Simultaneous Analysis of Multiple Networks (SAMNet) computational framework to model a genome scale shRNA screen investigating Acute Lymphoblastic Leukemia (ALL) progression in vivo. Our network model is enriched for cellular processes associated with hematopoietic differentiation and homeostasis even though none of the individual 'omic sets showed this enrichment. The model identifies genes associated with the TGF-beta pathway and predicts a role in ALL progression for many genes without this functional annotation. We further experimentally validate the hidden genes - Wwp1, a ubiquitin ligase, and Hgs, a multi-vesicular body associated protein - for their role in ALL progression. Our ALL pathway model includes genes with roles in multiple types of leukemia and roles in hematological development. We identify a tumor suppressor role for Wwp1 in ALL progression. This work demonstrates that network integration approaches can compensate for off-target effects, and that these methods can uncover novel biology retroactively on existing screening data. We anticipate that this framework will be valuable to multiple functional genomic technologies - siRNA, shRNA, and CRISPR - generally, and will improve the utility of functional genomic studies.

Network Modeling Identifies Patient-specific Pathways in Glioblastoma

Glioblastoma is the most aggressive type of malignant human brain tumor. Molecular profiling experiments have revealed that these tumors are extremely heterogeneous. This heterogeneity is one of the principal challenges for developing targeted therapies. We hypothesize that despite the diverse molecular profiles, it might still be possible to identify common signaling changes that could be targeted in some or all tumors. Using a network modeling approach, we reconstruct the altered signaling pathways from tumor-specific phosphoproteomic data and known protein-protein interactions. We then develop a network-based strategy for identifying tumor specific proteins and pathways that were predicted by the models but not directly observed in the experiments. Among these hidden targets, we show that the ERK activator kinase1 (MEK1) displays increased phosphorylation in all tumors. By contrast, protein numb homolog (NUMB) is present only in the subset of the tumors that are the most invasive. Additionally, increased S100A4 is associated with only one of the tumors. Overall, our results demonstrate that despite the heterogeneity of the proteomic data, network models can identify common or tumor specific pathway-level changes. These results represent an important proof of principle that can improve the target selection process for tumor specific treatments.

Network-Based Interpretation of Diverse High-Throughput Datasets through the Omics Integrator Software Package

High-throughput, ‘omic’ methods provide sensitive measures of biological responses to perturbations. However, inherent biases in high-throughput assays make it difficult to interpret experiments in which more than one type of data is collected. In this work, we introduce Omics Integrator, a software package that takes a variety of ‘omic’ data as input and identifies putative underlying molecular pathways. The approach applies advanced network optimization algorithms to a network of thousands of molecular interactions to find high-confidence, interpretable subnetworks that best explain the data. These subnetworks connect changes observed in gene expression, protein abundance or other global assays to proteins that may not have been measured in the screens due to inherent bias or noise in measurement. This approach reveals unannotated molecular pathways that would not be detectable by searching pathway databases. Omics Integrator also provides an elegant framework to incorporate not only positive data, but also negative evidence. Incorporating negative evidence allows Omics Integrator to avoid unexpressed genes and avoid being biased toward highly-studied hub proteins, except when they are strongly implicated by the data. The software is comprised of two individual tools, Garnet and Forest, that can be run together or independently to allow a user to perform advanced integration of multiple types of high-throughput data as well as create condition-specific subnetworks of protein interactions that best connect the observed changes in various datasets. It is available at http://fraenkel.mit.edu/omicsintegrator and on GitHub at https://github.com/fraenkel-lab/OmicsIntegrator.

Elucidating MicroRNA Regulatory Networks Using Transcriptional, Post-transcriptional, and Histone Modification Measurements

MicroRNAs (miRNAs) regulate diverse biological processes by repressing mRNAs, but their modest effects on direct targets, together with their participation in larger regulatory networks, make it challenging to delineate miRNA-mediated effects. Here, we describe an approach to characterizing miRNA-regulatory networks by systematically profiling transcriptional, post-transcriptional and epigenetic activity in a pair of isogenic murine fibroblast cell lines with and without Dicer expression. By RNA sequencing (RNA-seq) and CLIP (crosslinking followed by immunoprecipitation) sequencing (CLIP-seq), we found that most of the changes induced by global miRNA loss occur at the level of transcription. We then introduced a network modeling approach that integrated these data with epigenetic data to identify specific miRNA-regulated transcription factors that explain the impact of miRNA perturbation on gene expression. In total, we demonstrate that combining multiple genome-wide datasets spanning diverse regulatory modes enables accurate delineation of the downstream miRNA-regulated transcriptional network and establishes a model for studying similar networks in other systems.

Large-scale models of signal propagation in human cells derived from discovery phosphoproteomic data

Mass spectrometry is widely used to probe the proteome and its modifications in an untargeted manner, with unrivalled coverage. Applied to phosphoproteomics, it has tremendous potential to interrogate phospho-signalling and its therapeutic implications. However, this task is complicated by issues of undersampling of the phosphoproteome and challenges stemming from its high-content but low-sample-throughput nature. Hence, methods using such data to reconstruct signalling networks have been limited to restricted data sets and insights (for example, groups of kinases likely to be active in a sample). We propose a new method to handle high-content discovery phosphoproteomics data on perturbation by putting it in the context of kinase/phosphatase-substrate knowledge, from which we derive and train logic models. We show, on a data set obtained through perturbations of cancer cells with small-molecule inhibitors, that this method can study the targets and effects of kinase inhibitors, and reconcile insights obtained from multiple data sets, a common issue with these data.

Phosphoproteomics can offer significant insight into cell signalling and how signalling is modified in response to perturbations. Here the authors develop a new tool for the analysis of high-content phosphoproteomics in the context of kinase/phosphatase-substrate knowledge, which is used to train logic models.

Reconstruction of the temporal signaling network in Salmonella-infected human cells

Salmonella enterica is a bacterial pathogen that usually infects its host through food sources. Translocation of the pathogen proteins into the host cells leads to changes in the signaling mechanism either by activating or inhibiting the host proteins. Given that the bacterial infection modifies the response network of the host, a more coherent view of the underlying biological processes and the signaling networks can be obtained by using a network modeling approach based on the reverse engineering principles. In this work, we have used a published temporal phosphoproteomic dataset of Salmonella-infected human cells and reconstructed the temporal signaling network of the human host by integrating the interactome and the phosphoproteomic dataset. We have combined two well-established network modeling frameworks, the Prize-collecting Steiner Forest (PCSF) approach and the Integer Linear Programming (ILP) based edge inference approach. The resulting network conserves the information on temporality, direction of interactions, while revealing hidden entities in the signaling, such as the SNARE binding, mTOR signaling, immune response, cytoskeleton organization, and apoptosis pathways. Targets of the Salmonella effectors in the host cells such as CDC42, RHOA, 14-3-3δ, Syntaxin family, Oxysterol-binding proteins were included in the reconstructed signaling network although they were not present in the initial phosphoproteomic data. We believe that integrated approaches, such as the one presented here, have a high potential for the identification of clinical targets in infectious diseases, especially in the Salmonella infections.

ONION: Functional Approach for Integration of Lipidomics and Transcriptomics Data

To date, the massive quantity of data generated by high-throughput techniques has not yet met bioinformatics treatment required to make full use of it. This is partially due to a mismatch in experimental and analytical study design but primarily due to a lack of adequate analytical approaches. When integrating multiple data types e.g. transcriptomics and metabolomics, multidimensional statistical methods are currently the techniques of choice. Typical statistical approaches, such as canonical correlation analysis (CCA), that are applied to find associations between metabolites and genes are failing due to small numbers of observations (e.g. conditions, diet etc.) in comparison to data size (number of genes, metabolites). Modifications designed to cope with this issue are not ideal due to the need to add simulated data resulting in a lack of p-value computation or by pruning of variables hence losing potentially valid information. Instead, our approach makes use of verified or putative molecular interactions or functional association to guide analysis. The workflow includes dividing of data sets to reach the expected data structure, statistical analysis within groups and interpretation of results. By applying pathway and network analysis, data obtained by various platforms are grouped with moderate stringency to avoid functional bias. As a consequence CCA and other multivariate models can be applied to calculate robust statistics and provide easy to interpret associations between metabolites and genes to leverage understanding of metabolic response. Effective integration of lipidomics and transcriptomics is demonstrated on publically available murine nutrigenomics data sets. We are able to demonstrate that our approach improves detection of genes related to lipid metabolism, in comparison to applying statistics alone. This is measured by increased percentage of explained variance (95% vs. 75–80%) and by identifying new metabolite-gene associations related to lipid metabolism.

ChiNet uncovers rewired transcription subnetworks in tolerant yeast for advanced biofuels conversion

Analysis of rewired upstream subnetworks impacting downstream differential gene expression aids the delineation of evolving molecular mechanisms. Cumulative statistics based on conventional differential correlation are limited for subnetwork rewiring analysis since rewiring is not necessarily equivalent to change in correlation coefficients. Here we present a computational method ChiNet to quantify subnetwork rewiring by statistical heterogeneity that enables detection of potential genotype changes causing altered transcription regulation in evolving organisms. Given a differentially expressed downstream gene set, ChiNet backtracks a rewired upstream subnetwork from a super-network including gene interactions known to occur under various molecular contexts. We benchmarked ChiNet for its high accuracy in distinguishing rewired artificial subnetworks, in silico yeast transcription-metabolic subnetworks, and rewired transcription subnetworks for Candida albicans versus Saccharomyces cerevisiae, against two differential-correlation based subnetwork rewiring approaches. Then, using transcriptome data from tolerant S. cerevisiae strain NRRL Y-50049 and a wild-type intolerant strain, ChiNet identified 44 metabolic pathways affected by rewired transcription subnetworks anchored to major adaptively activated transcription factor genes YAP1, RPN4, SFP1 and ROX1, in response to toxic chemical challenges involved in lignocellulose-to-biofuels conversion. These findings support the use of ChiNet in rewiring analysis of subnetworks where differential interaction patterns resulting from divergent nonlinear dynamics abound.

Multitask learning of signaling and regulatory networks with application to studying human response to flu

Reconstructing regulatory and signaling response networks is one of the major goals of systems biology. While several successful methods have been suggested for this task, some integrating large and diverse datasets, these methods have so far been applied to reconstruct a single response network at a time, even when studying and modeling related conditions. To improve network reconstruction we developed MT-SDREM, a multi-task learning method which jointly models networks for several related conditions. In MT-SDREM, parameters are jointly constrained across the networks while still allowing for condition-specific pathways and regulation. We formulate the multi-task learning problem and discuss methods for optimizing the joint target function. We applied MT-SDREM to reconstruct dynamic human response networks for three flu strains: H1N1, H5N1 and H3N2. Our multi-task learning method was able to identify known and novel factors and genes, improving upon prior methods that model each condition independently. The MT-SDREM networks were also better at identifying proteins whose removal affects viral load indicating that joint learning can still lead to accurate, condition-specific, networks. Supporting website with MT-SDREM implementation: http://sb.cs.cmu.edu/mtsdrem

To understand why some flu strains are more virulent than others, researchers attempt to profile and model the molecular human response to these strains and identify similarities and differences between the resulting models. So far, the modeling and analysis part has been done independently for each strain and the results contrasted in a post-processing step. Here we present a new method, termed MT-SDREM, that simultaneously models the response to all strains allowing us to identify both, the core response elements that are shared among the strains, and factors that are uniquely activated or repressed by individual strains. We applied this method to study the human response to three flu strains: H1N1, H3N2 and H5N1. As we show, the method was able to correctly identify several common and known factors regulating immune response to such strains and also identified unique factors for each of the strains. The models reconstructed by the simultaneous analysis method improved upon those generated by methods that model each strain response separately. Our joint models can be used to identify strain specific treatments as well as treatments that are likely to be effective against all three strains.

SAMNetWeb: identifying condition-specific networks linking signaling and transcription

MOTIVATION

High-throughput datasets such as genetic screens, mRNA expression assays and global phospho-proteomic experiments are often difficult to interpret due to inherent noise in each experimental system. Computational tools have improved interpretation of these datasets by enabling the identification of biological processes and pathways that are most likely to explain the measured results. These tools are primarily designed to analyse data from a single experiment (e.g. drug treatment versus control), creating a need for computational algorithms that can handle heterogeneous datasets across multiple experimental conditions at once.

SUMMARY

We introduce SAMNetWeb, a web-based tool that enables functional enrichment analysis and visualization of high-throughput datasets. SAMNetWeb can analyse two distinct data types (e.g. mRNA expression and global proteomics) simultaneously across multiple experimental systems to identify pathways activated in these experiments and then visualize the pathways in a single interaction network. Through the use of a multi-commodity flow based algorithm that requires each experiment 'share' underlying protein interactions, SAMNetWeb can identify distinct and common pathways across experiments.

AVAILABILITY AND IMPLEMENTATION

SAMNetWeb is freely available at http://fraenkel.mit.edu/samnetweb.

Sharing information to reconstruct patient-specific pathways in heterogeneous diseases

Advances in experimental techniques resulted in abundant genomic, transcriptomic, epigenomic, and proteomic data that have the potential to reveal critical drivers of human diseases. Complementary algorithmic developments enable researchers to map these data onto protein-protein interaction networks and infer which signaling pathways are perturbed by a disease. Despite this progress, integrating data across different biological samples or patients remains a substantial challenge because samples from the same disease can be extremely heterogeneous. Somatic mutations in cancer are an infamous example of this heterogeneity. Although the same signaling pathways may be disrupted in a cancer patient cohort, the distribution of mutations is long-tailed, and many driver mutations may only be detected in a small fraction of patients. We developed a computational approach to account for heterogeneous data when inferring signaling pathways by sharing information across the samples. Our technique builds upon the prize-collecting Steiner forest problem, a network optimization algorithm that extracts pathways from a protein-protein interaction network. We recover signaling pathways that are similar across all samples yet still reflect the unique characteristics of each biological sample. Leveraging data from related tumors improves our ability to recover the disrupted pathways and reveals patient-specific pathway perturbations in breast cancer.

Integrated network analyses for functional genomic studies in cancer

RNA-interference (RNAi) studies hold great promise for functional investigation of the significance of genetic variations and mutations, as well as potential synthetic lethalities, for understanding and treatment of cancer, yet technical and conceptual issues currently diminish the potential power of this approach. While numerous research groups are usefully employing this kind of functional genomic methodology to identify molecular mediators of disease severity, response, and resistance to treatment, findings are generally confounded by "off-target" effects. These effects arise from a variety of issues beyond non-specific reagent behavior, such as biological cross-talk and feedback processes so thus can occur even with specific perturbation. Interpreting RNAi results in a network framework instead of merely as individual "hits" or "targets" leverages contributions from all hit/target contributions to pathways via their relationships with other network nodes. This interpretation can ameliorate dependence upon individual reagent performance and increase confidence in biological validation. Here we provide background on RNAi studies in cancer applications, review key challenges with functional genomics, and motivate the use of network models grounded in pathway analyses.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.