Unraveling druggable cancer-driving proteins and targeted drugs using artificial intelligence and multi-omics analyses

The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .

B- and T-Lymphocyte Attenuator in Systemic Lupus Erythematosus Disease Pathogenesis

B- and T-lymphocyte attenuator (BTLA; CD272) is an immunoglobulin superfamily member and part of a family of checkpoint inhibitory receptors that negatively regulate immune cell activation. The natural ligand for BTLA is herpes virus entry mediator (HVEM; TNFRSF14), and binding of HVEM to BTLA leads to attenuation of lymphocyte activation. In this study, we evaluated the role of BTLA and HVEM expression in the pathogenesis of systemic lupus erythematosus (SLE), a multisystem autoimmune disease. Peripheral blood mononuclear cells from healthy volunteers (N = 7) were evaluated by mass cytometry by time-of-flight to establish baseline expression of BTLA and HVEM on human lymphocytes compared with patients with SLE during a self-reported flare (N = 5). High levels of BTLA protein were observed on B cells, CD4+, and CD8+ T cells, and plasmacytoid dendritic cells in healthy participants. HVEM protein levels were lower in patients with SLE compared with healthy participants, while BTLA levels were similar between SLE and healthy groups. Correlations of BTLA-HVEM hub genes' expression with patient and disease characteristics were also analyzed using whole blood gene expression data from patients with SLE (N = 1,760) and compared with healthy participants (N = 60). HVEM, being one of the SLE-associated genes, showed an exceptionally strong negative association with disease activity. Several other genes in the BTLA-HVEM signaling network were strongly (negative or positive) correlated, while BTLA had a low association with disease activity. Collectively, these data provide a clinical rationale for targeting BTLA with an agonist in SLE patients with low HVEM expression.

Innovative target mining stratagems to navigate drug repurposing endeavours

The conventional theory linking a single gene with a particular disease and a specific drug contributes to the dwindling success rates of traditional drug discovery. This requires a substantial shift focussing on contemporary drug design or drug repurposing, which entails linking multiple genes to diverse physiological or pathological pathways and drugs. Lately, drug repurposing, the art of discovering new/unlabelled indications for existing drugs or candidates in clinical trials, is gaining attention owing to its success rates. The rate-limiting phase of this strategy lies in target identification, which is generally driven through disease-centric and/or drug-centric approaches. The disease-centric approach is based on exploration of crucial biomolecules such as genes or proteins underlying pathological cascades of the disease of interest. Investigating these pathological interplays aids in the identification of potential drug targets that can be leveraged for novel therapeutic interventions. The drug-centric approach involves various strategies such as exploring the mechanism of adverse drug reactions that can unearth potential targets, as these untoward reactions might be considered desirable therapeutic actions in other disease conditions. Currently, artificial intelligence is an emerging robust tool that can be used to translate the aforementioned intricate biological networks to render interpretable data for extracting precise molecular targets. Integration of multiple approaches, big data analytics, and clinical corroboration are essential for successful target mining. This chapter highlights the contemporary strategies steering target identification and diverse frameworks for drug repurposing. These strategies are illustrated through case studies curated from recent drug repurposing research inclined towards neurodegenerative diseases, cancer, infections, immunological, and cardiovascular disorders.

PWN: enhanced random walk on a warped network for disease target prioritization

BACKGROUND

Extracting meaningful information from unbiased high-throughput data has been a challenge in diverse areas. Specifically, in the early stages of drug discovery, a considerable amount of data was generated to understand disease biology when identifying disease targets. Several random walk-based approaches have been applied to solve this problem, but they still have limitations. Therefore, we suggest a new method that enhances the effectiveness of high-throughput data analysis with random walks.

RESULTS

We developed a new random walk-based algorithm named prioritization with a warped network (PWN), which employs a warped network to achieve enhanced performance. Network warping is based on both internal and external features: graph curvature and prior knowledge.

CONCLUSIONS

We showed that these compositive features synergistically increased the resulting performance when applied to random walk algorithms, which led to PWN consistently achieving the best performance among several other known methods. Furthermore, we performed subsequent experiments to analyze the characteristics of PWN.

Systemic Efficacy of Sirolimus via the ERBB Signaling Pathway in Breast Cancer

Rapamycin, also known as sirolimus, inhibits the mTOR pathway in complex diseases such as cancer, and its downstream targets are ribosomal S6 kinases (RPS6K). Sirolimus is involved in regulating cell growth and cell survival through roles such as the mediation of epidermal growth factor signaling. However, the systemic efficacy of sirolimus in pathway regulation is unclear. The purpose of this study is to determine systemic drug efficacy using computational methods and drug-induced datasets. We suggest a computational method using gene expression datasets induced by sirolimus and an inverse algorithm that simultaneously identifies parameters referring to gene–gene interactions. We downloaded two sirolimus-induced microarray gene expression datasets and used a computational method to obtain the most enriched pathway, then adopted an inverse algorithm to discover the gene–gene interactions of that pathway. In the results, RPS6KB1 was a target gene of sirolimus and was associated with genes in the pathway. The common gene interactions from two datasets were a hub gene, RPS6KB1, and 10 related genes (AKT3, CBLC, MAP2K7, NRG1/2, PAK3, PIK3CD/G, PRKCG, and SHC3) in the epidermal growth factor (ERBB) signaling pathway.

Integrative system biology and mathematical modeling of genetic networks identifies shared biomarkers for obesity and diabetes

Obesity and type 2 and diabetes mellitus (T2D) are two dual epidemics whose shared genetic pathological mechanisms are still far from being fully understood. Therefore, this study is aimed at discovering key genes, molecular mechanisms, and new drug targets for obesity and T2D by analyzing the genome wide gene expression data with different computational biology approaches. In this study, the RNA-sequencing data of isolated primary human adipocytes from individuals who are lean, obese, and T2D was analyzed by an integrated framework consisting of gene expression, protein interaction network (PIN), tissue specificity, and druggability approaches. Our findings show a total of 1932 unique differentially expressed genes (DEGs) across the diabetes versus obese group comparison (p≤0.05). The PIN analysis of these 1932 DEGs identified 190 high centrality network (HCN) genes, which were annotated against 3367 GO terms and functional pathways, like response to insulin signaling, phosphorylation, lipid metabolism, glucose metabolism, etc. (p≤0.05). By applying additional PIN and topological parameters to 190 HCN genes, we further mapped 25 high confidence genes, functionally connected with diabetes and obesity traits. Interestingly, ERBB2, FN1, FYN, HSPA1A, HBA1, and ITGB1 genes were found to be tractable by small chemicals, antibodies, and/or enzyme molecules. In conclusion, our study highlights the potential of computational biology methods in correlating expression data to topological parameters, functional relationships, and druggability characteristics of the candidate genes involved in complex metabolic disorders with a common etiological basis.

GuiltyTargets: Prioritization of Novel Therapeutic Targets With Network Representation Learning

The majority of clinical trials fail due to low efficacy of investigated drugs, often resulting from a poor choice of target protein. Existing computational approaches aim to support target selection either via genetic evidence or by putting potential targets into the context of a disease specific network reconstruction. The purpose of this work was to investigate whether network representation learning techniques could be used to allow for a machine learning based prioritization of putative targets. We propose a novel target prioritization approach, GuiltyTargets, which relies on attributed network representation learning of a genome-wide protein-protein interaction network annotated with disease-specific differential gene expression and uses positive-unlabeled (PU) machine learning for candidate ranking. We evaluated our approach on 12 datasets from six diseases of different type (cancer, metabolic, neurodegenerative) within a 10 times repeated 5-fold stratified cross-validation and achieved AUROC values between 0.92 - 0.97, significantly outperforming previous approaches that relied on manually engineered topological features. Moreover, we showed that GuiltyTargets allows for target repositioning across related disease areas. An application of GuiltyTargets to Alzheimer's disease resulted in a number of highly ranked candidates that are currently discussed as targets in the literature. Interestingly, one (COMT) is also the target of an approved drug (Tolcapone) for Parkinson's disease, highlighting the potential for target repositioning with our method. The GuiltyTargets Python package is available on PyPI and all code used for analysis can be found under the MIT License at https://github.com/GuiltyTargets. Attributed network representation learning techniques provide an interesting approach to effectively leverage the existing knowledge about the molecular mechanisms in different diseases. In this work, the combination with positive-unlabeled learning for target prioritization demonstrated a clear superiority compared to classical feature engineering approaches. Our work highlights the potential of attributed network representation learning for target prioritization. Given the overarching relevance of networks in computational biology we believe that attributed network representation learning techniques could have a broader impact in the future.

Predicting mechanism of action of cellular perturbations with pathway activity signatures

MOTIVATION

Misregulation of signaling pathway activity is etiologic for many human diseases, and modulating activity of signaling pathways is often the preferred therapeutic strategy. Understanding the mechanism of action (MOA) of bioactive chemicals in terms of targeted signaling pathways is the essential first step in evaluating their therapeutic potential. Changes in signaling pathway activity are often not reflected in changes in expression of pathway genes which makes MOA inferences from transcriptional signatures (TSeses) a difficult problem.

RESULTS

We developed a new computational method for implicating pathway targets of bioactive chemicals and other cellular perturbations by integrated analysis of pathway network topology, the Library of Integrated Network-based Cellular Signature TSes of genetic perturbations of pathway genes and the TS of the perturbation. Our methodology accurately predicts signaling pathways targeted by the perturbation when current pathway analysis approaches utilizing only the TS of the perturbation fail.

AVAILABILITY AND IMPLEMENTATION

Open source R package paslincs is available at https://github.com/uc-bd2k/paslincs.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TREAP: A New Topological Approach to Drug Target Inference

Over 50% of drugs fail in stage 3 clinical trials, many because of a poor understanding of the drug's mechanisms of action (MoA). A better comprehension of drug MoA will significantly improve research and development (R&D). Current proposed algorithms, such as ProTINA and DeMAND, can be overly complex. Additionally, they are unable to predict whether the drug-induced gene expression or the topology of the networks used to model gene regulation primarily impacts accurate drug target inference. In this work, we evaluate how network and gene expression data affect ProTINA's accuracy. We find that network topology predominantly determines the accuracy of ProTINA's predictions. We further show that the size of an interaction network and/or selecting cell-specific networks has a limited effect on accuracy. We then demonstrate that a specific network topology measure, betweenness, can be used to improve drug target prediction. Based on these results, we create a new algorithm, TREAP, that combines betweenness values and adjusted p-values for target inference. TREAP offers an alternative approach to drug target inference and is advantageous because it is not computationally demanding, provides easy-to-interpret results, and is often more accurate at predicting drug targets than current state-of-the-art approaches.

A New Computational Approach to Evaluating Systemic Gene–Gene Interactions in a Pathway Affected by Drug LY294002

In this study, we investigate how drugs systemically affect genes via pathways by integrating information from interactions between chemical compounds and molecular expression datasets, and from pathway information such as gene sets using mathematical models. First, we adopt drug-induced gene expression datasets; then, employ gene set enrichment analysis tools for selecting candidate enrichment pathways; and lastly, implement the inverse algorithm package for identifying gene–gene regulatory networks in a pathway. We tested LY294002-induced datasets of the MCF7 breast cancer cell lines, and found a CELL CYCLE pathway with 101 genes, ERBB signaling pathway consisting of 82 genes, and MTOR pathway consisting of 45 genes. We consider two interactions: quantity strength depending on number of interactions, and quality strength depending on weight of interaction as positive (+) and negative (−) interactions. Our methods revealed ANAPC1-CDK6 (−0.412) and ORC2L- CHEK1(0.951) for the CELL CYCLE pathway; INS-RPS6 (−3.125) and PRKAA2-PRKAA2 (+1.319) for the MTOR pathway; and CBLB-RPS6KB1 (−0.141), RPS6KB1-CBLC (+0.238) for the ERBB signaling pathway to be top quality interactions. Top quantity interactions discovered include 12; the CDC (−,+) gene family for the CELL CYCLE pathway, 20; PIK3 (−), 23; PIK3CG (+) for the MTOR pathway, 11; PAK (−), 10; PIK3 (+) for the ERBB signaling pathway.

Network hub-node prioritization of gene regulation with intra-network association

Background

To identify and prioritize the influential hub genes in a gene-set or biological pathway, most analyses rely on calculation of marginal effects or tests of statistical significance. These procedures may be inappropriate since hub nodes are common connection points and therefore may interact with other nodes more often than non-hub nodes do. Such dependence among gene nodes can be conjectured based on the topology of the pathway network or the correlation between them.

Results

Here we develop a pathway activity score incorporating the marginal (local) effects of gene nodes as well as intra-network affinity measures. This score summarizes the expression levels in a gene-set/pathway for each sample, with weights on local and network information, respectively. The score is next used to examine the impact of each node through a leave-one-out evaluation. To illustrate the procedure, two cancer studies, one involving RNA-Seq from breast cancer patients with high-grade ductal carcinoma in situ and one microarray expression data from ovarian cancer patients, are used to assess the performance of the procedure, and to compare with existing methods, both ones that do and do not take into consideration correlation and network information. The hub nodes identified by the proposed procedure in the two cancer studies are known influential genes; some have been included in standard treatments and some are currently considered in clinical trials for target therapy. The results from simulation studies show that when marginal effects are mild or weak, the proposed procedure can still identify causal nodes, whereas methods relying only on marginal effect size cannot.

Conclusions

The NetworkHub procedure proposed in this research can effectively utilize the network information in combination with local effects derived from marker values, and provide a useful and complementary list of recommendations for prioritizing causal hubs.

Applications of Network Analysis in Biomedicine

The abundance of high-throughput data and technical refinements in graph theories have allowed network analysis to become an effective approach for various medical fields. This chapter introduces co-expression, Bayesian, and regression-based network construction methods, which are the basis of network analysis. Various methods in network topology analysis are explained, along with their unique features and applications in biomedicine. Furthermore, we explain the role of network embedding in reducing the dimensionality of networks and outline several popular algorithms used by researchers today. Current literature has implemented different combinations of topology analysis and network embedding techniques, and we outline several studies in the fields of genetic-based disease prediction, drug-target identification, and multi-level omics integration.

Network-based method for drug target discovery at the isoform level

Identification of primary targets associated with phenotypes can facilitate exploration of the underlying molecular mechanisms of compounds and optimization of the structures of promising drugs. However, the literature reports limited effort to identify the target major isoform of a single known target gene. The majority of genes generate multiple transcripts that are translated into proteins that may carry out distinct and even opposing biological functions through alternative splicing. In addition, isoform expression is dynamic and varies depending on the developmental stage and cell type. To identify target major isoforms, we integrated a breast cancer type-specific isoform coexpression network with gene perturbation signatures in the MCF7 cell line in the Connectivity Map database using the ‘shortest path’ drug target prioritization method. We used a leukemia cancer network and differential expression data for drugs in the HL-60 cell line to test the robustness of the detection algorithm for target major isoforms. We further analyzed the properties of target major isoforms for each multi-isoform gene using pharmacogenomic datasets, proteomic data and the principal isoforms defined by the APPRIS and STRING datasets. Then, we tested our predictions for the most promising target major protein isoforms of DNMT1, MGEA5 and P4HB4 based on expression data and topological features in the coexpression network. Interestingly, these isoforms are not annotated as principal isoforms in APPRIS. Lastly, we tested the affinity of the target major isoform of MGEA5 for streptozocin through in silico docking. Our findings will pave the way for more effective and targeted therapies via studies of drug targets at the isoform level.

Network perturbation analysis of gene transcriptional profiles reveals protein targets and mechanism of action of drugs and influenza A viral infection

Genome-wide transcriptional profiling provides a global view of cellular state and how this state changes under different treatments (e.g. drugs) or conditions (e.g. healthy and diseased). Here, we present ProTINA (Protein Target Inference by Network Analysis), a network perturbation analysis method for inferring protein targets of compounds from gene transcriptional profiles. ProTINA uses a dynamic model of the cell-type specific protein-gene transcriptional regulation to infer network perturbations from steady state and time-series differential gene expression profiles. A candidate protein target is scored based on the gene network's dysregulation, including enhancement and attenuation of transcriptional regulatory activity of the protein on its downstream genes, caused by drug treatments. For benchmark datasets from three drug treatment studies, ProTINA was able to provide highly accurate protein target predictions and to reveal the mechanism of action of compounds with high sensitivity and specificity. Further, an application of ProTINA to gene expression profiles of influenza A viral infection led to new insights of the early events in the infection.

gene2drug: a computational tool for pathway-based rational drug repositioning

Motivation

Drug repositioning has been proposed as an effective shortcut to drug discovery. The availability of large collections of transcriptional responses to drugs enables computational approaches to drug repositioning directly based on measured molecular effects.

Results

We introduce a novel computational methodology for rational drug repositioning, which exploits the transcriptional responses following treatment with small molecule. Specifically, given a therapeutic target gene, a prioritization of potential effective drugs is obtained by assessing their impact on the transcription of genes in the pathway(s) including the target. We performed in silico validation and comparison with a state-of-art technique based on similar principles. We next performed experimental validation in two different real-case drug repositioning scenarios: (i) upregulation of the glutamate-pyruvate transaminase (GPT), which has been shown to induce reduction of oxalate levels in a mouse model of primary hyperoxaluria, and (ii) activation of the transcription factor TFEB, a master regulator of lysosomal biogenesis and autophagy, whose modulation may be beneficial in neurodegenerative disorders.

Availability and implementation

A web tool for Gene2drug is freely available at http://gene2drug.tigem.it. An R package is under development and can be obtained from https://github.com/franapoli/gep2pep.

Contact

dibernardo@tigem.it.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Comparative Study of Cluster Detection Algorithms in Protein-Protein Interaction for Drug Target Discovery and Drug Repurposing

The interactions between drugs and their target proteins induce altered expression of genes involved in complex intracellular networks. The properties of these functional network modules are critical for the identification of drug targets, for drug repurposing, and for understanding the underlying mode of action of the drug. The topological modules generated by a computational approach are defined as functional clusters. However, the functions inferred for these topological modules extracted from a large-scale molecular interaction network, such as a protein–protein interaction (PPI) network, could differ depending on different cluster detection algorithms. Moreover, the dynamic gene expression profiles among tissues or cell types causes differential functional interaction patterns between the molecular components. Thus, the connections in the PPI network should be modified by the transcriptomic landscape of specific cell lines before producing topological clusters. Here, we systematically investigated the clusters of a cell-based PPI network by using four cluster detection algorithms. We subsequently compared the performance of these algorithms for target gene prediction, which integrates gene perturbation data with the cell-based PPI network using two drug target prioritization methods, shortest path and diffusion correlation. In addition, we validated the proportion of perturbed genes in clusters by finding candidate anti-breast cancer drugs and confirming our predictions using literature evidence and cases in the ClinicalTrials.gov. Our results indicate that the Walktrap (CW) clustering algorithm achieved the best performance overall in our comparative study.

Predicting protein targets for drug-like compounds using transcriptomics

An expanded chemical space is essential for improved identification of small molecules for emerging therapeutic targets. However, the identification of targets for novel compounds is biased towards the synthesis of known scaffolds that bind familiar protein families, limiting the exploration of chemical space. To change this paradigm, we validated a new pipeline that identifies small molecule-protein interactions and works even for compounds lacking similarity to known drugs. Based on differential mRNA profiles in multiple cell types exposed to drugs and in which gene knockdowns (KD) were conducted, we showed that drugs induce gene regulatory networks that correlate with those produced after silencing protein-coding genes. Next, we applied supervised machine learning to exploit drug-KD signature correlations and enriched our predictions using an orthogonal structure-based screen. As a proof-of-principle for this regimen, top-10/top-100 target prediction accuracies of 26% and 41%, respectively, were achieved on a validation of set 152 FDA-approved drugs and 3104 potential targets. We then predicted targets for 1680 compounds and validated chemical interactors with four targets that have proven difficult to chemically modulate, including non-covalent inhibitors of HRAS and KRAS. Importantly, drug-target interactions manifest as gene expression correlations between drug treatment and both target gene KD and KD of genes that act up- or down-stream of the target, even for relatively weak binders. These correlations provide new insights on the cellular response of disrupting protein interactions and highlight the complex genetic phenotypes of drug treatment. With further refinement, our pipeline may accelerate the identification and development of novel chemical classes by screening compound-target interactions.

Bioactive compounds often disrupt cellular gene expression in ways that are difficult to predict. While the correlation between a cellular response after treatment with a small molecule and the knockdown of its target protein should be simple to establish, in practice this goal has been difficult to achieve. The main challenges are that data are noisy, drugs are not intended to be active in all cell types, and signals from a bona fide target(s) may be obscured by correlations with knockdowns of other proteins in the same pathway(s). Here, we find that a random forest classification model can detect meaningful correlational patterns when gene expression profiles after compound treatment and gene knockdowns in four or more cell lines are compared. When this approach is enriched by a structure-based screen, novel drug-target interactions can be predicted. We then validated new ligand-protein interactions for four difficult targets. Although the initial compounds are not especially potent in vitro, they are capable of disrupting their target pathway in the cell to an extent that generates a significant and characteristic gene expression profile. Collectively, our studies provide insight on cell-level transcriptomic responses to pharmaceutical intervention and the use of these patterns for target identification. In addition, the method provides a novel drug discovery pipeline to test chemistries without a priori knowledge of their target(s).

JDINAC: joint density-based non-parametric differential interaction network analysis and classification using high-dimensional sparse omics data

Motivation

A complex disease is usually driven by a number of genes interwoven into networks, rather than a single gene product. Network comparison or differential network analysis has become an important means of revealing the underlying mechanism of pathogenesis and identifying clinical biomarkers for disease classification. Most studies, however, are limited to network correlations that mainly capture the linear relationship among genes, or rely on the assumption of a parametric probability distribution of gene measurements. They are restrictive in real application.

Results

We propose a new Joint density based non-parametric Differential Interaction Network Analysis and Classification (JDINAC) method to identify differential interaction patterns of network activation between two groups. At the same time, JDINAC uses the network biomarkers to build a classification model. The novelty of JDINAC lies in its potential to capture non-linear relations between molecular interactions using high-dimensional sparse data as well as to adjust confounding factors, without the need of the assumption of a parametric probability distribution of gene measurements. Simulation studies demonstrate that JDINAC provides more accurate differential network estimation and lower classification error than that achieved by other state-of-the-art methods. We apply JDINAC to a Breast Invasive Carcinoma dataset, which includes 114 patients who have both tumor and matched normal samples. The hub genes and differential interaction patterns identified were consistent with existing experimental studies. Furthermore, JDINAC discriminated the tumor and normal sample with high accuracy by virtue of the identified biomarkers. JDINAC provides a general framework for feature selection and classification using high-dimensional sparse omics data.

Availability and implementation

R scripts available at https://github.com/jijiadong/JDINAC.

Contact

lxie@iscb.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

Paeonol induces the apoptosis of the SGC‑7901 gastric cancer cell line by downregulating ERBB2 and inhibiting the NF‑κB signaling pathway

The purpose of the present study was to analyze the association between paeonol and the known genes related to gastric cancer (GC) using bioinformatics methods, and to investigate the role of paeonol in the potential impact on the nuclear factor‑κB (NF‑κB) signaling pathway, in order to provide a theoretical basis for further elucidating the effect of paeonol on cancer cells. Cell viability, morphology and apoptosis were detected using an MTT assay, an inverted microscope, and flow cytometry, respectively. The correlation between drugs and genes was analyzed using the Search Tool for Interactions of Chemicals (STITCH) gene‑drug interaction network. The expression levels of related mRNA and proteins were determined using reverse transcription‑quantitative polymerase chain reaction analysis and enzyme‑linked immunosorbent assay. The changes in protein expression were examined using western blot analysis. The correlation network between target genes directly affected by paeonol and known GC genes was determined by analyzing the association between the compounds and genes recorded in the STITCH database. The GC‑related epidermal growth factor receptor 2 (ERBB2) gene was at the core position of the paeonol interaction network and may be an important potential target gene for the effect of paeonol on cancer cells. The effect of paeonol on the viability of the SGC‑7901 GC cell line was detected using an MTT assay, which showed that the inhibitory effect occurred in a time‑ and dose‑dependent manner. The observations of cell morphology demonstrated that the cells were floating, abnormal in shape, had unclear boundaries and were sparse in arrangement following paeonol treatment. Flow cytometry indicated that paeonol significantly accelerated the apoptotic rate of the SGC‑7901 GC cells. The examination of clinical samples suggested that ERBB2 was expressed at a high level in GC samples, and was significantly downregulated following the addition of paeonol. The western blot analysis revealed that downregulating ERBB2 affected the activation of the NF‑κB signaling pathway, thereby upregulating the pro‑apoptotic factor B‑cell lymphoma‑associated X protein. Taken together, paeonol significantly downregulated ERBB2 and inhibited the activation of the NF‑κB signaling pathway, thereby inhibiting the proliferation of SGC‑7901 cells and inducing apoptosis.

A review of connectivity map and computational approaches in pharmacogenomics

Large-scale perturbation databases, such as Connectivity Map (CMap) or Library of Integrated Network-based Cellular Signatures (LINCS), provide enormous opportunities for computational pharmacogenomics and drug design. A reason for this is that in contrast to classical pharmacology focusing at one target at a time, the transcriptomics profiles provided by CMap and LINCS open the door for systems biology approaches on the pathway and network level. In this article, we provide a review of recent developments in computational pharmacogenomics with respect to CMap and LINCS and related applications.

Multi-target drug repositioning by bipartite block-wise sparse multi-task learning

Background

Finding potential drug targets is a crucial step in drug discovery and development. Recently, resources such as the Library of Integrated Network-Based Cellular Signatures (LINCS) L1000 database provide gene expression profiles induced by various chemical and genetic perturbations and thereby make it possible to analyze the relationship between compounds and gene targets at a genome-wide scale. Current approaches for comparing the expression profiles are based on pairwise connectivity mapping analysis. However, this method makes the simple assumption that the effect of a drug treatment is similar to knocking down its single target gene. Since many compounds can bind multiple targets, the pairwise mapping ignores the combined effects of multiple targets, and therefore fails to detect many potential targets of the compounds.

Results

We propose an algorithm to find sets of gene knock-downs that induce gene expression changes similar to a drug treatment. Assuming that the effects of gene knock-downs are additive, we propose a novel bipartite block-wise sparse multi-task learning model with super-graph structure (BBSS-MTL) for multi-target drug repositioning that overcomes the restrictive assumptions of connectivity mapping analysis.

Conclusions

The proposed method BBSS-MTL is more accurate for predicting potential drug targets than the simple pairwise connectivity mapping analysis on five datasets generated from different cancer cell lines.

Availability

The code can be obtained at http://gr.xjtu.edu.cn/web/liminli/codes.

COL3A1 and SNAP91: novel glioblastoma markers with diagnostic and prognostic value

Although patients with glioblastoma (GBM) have grave prognosis, significant variability in patient outcome is observed. This study aims to identify novel targets for GBM diagnosis and therapy. Microarray data (GSE4290, GSE7696, and GSE4412) obtained from the Gene Expression Omnibus was used to identify the differentially expressed genes (DEGs) by significant analysis of microarray (SAM). Intersection of the identified DEGs for each profile revealed 46 DEGs in GBM. A subset of common DEGs were validated by real-time reverse transcription quantitative PCR (qPCR). The prognostic value of some of the markers was also studied. We determined that RRM2 and COL3A1 were increased and directly correlated with glioma grade, while SH3GL2 and SNAP91 were decreased in GBM and inversely correlated with glioma grade. Kaplan-Meir analysis of GSE7696 revealed that COL3A1 and SNAP91 correlated with survival, suggesting that COL3A1 and SNAP91 may be suitable biomarkers for diagnostic or therapeutic strategies for GBM.

Big Data in Drug Discovery

Interpretation of Big Data in the drug discovery community should enhance project timelines and reduce clinical attrition through improved early decision making. The issues we encounter start with the sheer volume of data and how we first ingest it before building an infrastructure to house it to make use of the data in an efficient and productive way. There are many problems associated with the data itself including general reproducibility, but often, it is the context surrounding an experiment that is critical to success. Help, in the form of artificial intelligence (AI), is required to understand and translate the context. On the back of natural language processing pipelines, AI is also used to prospectively generate new hypotheses by linking data together. We explain Big Data from the context of biology, chemistry and clinical trials, showcasing some of the impressive public domain sources and initiatives now available for interrogation.

Unveiling network-based functional features through integration of gene expression into protein networks

Decoding health and disease phenotypes is one of the fundamental objectives in biomedicine. Whereas high-throughput omics approaches are available, it is evident that any single omics approach might not be adequate to capture the complexity of phenotypes. Therefore, integrated multi-omics approaches have been used to unravel genotype-phenotype relationships such as global regulatory mechanisms and complex metabolic networks in different eukaryotic organisms. Some of the progress and challenges associated with integrated omics studies have been reviewed previously in comprehensive studies. In this work, we highlight and review the progress, challenges and advantages associated with emerging approaches, integrating gene expression and protein-protein interaction networks to unravel network-based functional features. This includes identifying disease related genes, gene prioritization, clustering protein interactions, developing the modules, extract active subnetworks and static protein complexes or dynamic/temporal protein complexes. We also discuss how these approaches contribute to our understanding of the biology of complex traits and diseases. This article is part of a Special Issue entitled: Cardiac adaptations to obesity, diabetes and insulin resistance, edited by Professors Jan F.C. Glatz, Jason R.B. Dyck and Christine Des Rosiers.

Approaches to target tractability assessment - a practical perspective

The assessment of the suitability of novel targets to intervention by different modalities, e.g. small molecules or antibodies, is increasingly seen as important in helping to select the most progressable targets at the outset of a drug discovery project. This perspective considers differing aspects of tractability and how it can be assessed using in silico and experimental approaches. We also share some of our experiences in using these approaches.

In silico prediction of novel therapeutic targets using gene-disease association data

BACKGROUND

Target identification and validation is a pressing challenge in the pharmaceutical industry, with many of the programmes that fail for efficacy reasons showing poor association between the drug target and the disease. Computational prediction of successful targets could have a considerable impact on attrition rates in the drug discovery pipeline by significantly reducing the initial search space. Here, we explore whether gene-disease association data from the Open Targets platform is sufficient to predict therapeutic targets that are actively being pursued by pharmaceutical companies or are already on the market.

METHODS

To test our hypothesis, we train four different classifiers (a random forest, a support vector machine, a neural network and a gradient boosting machine) on partially labelled data and evaluate their performance using nested cross-validation and testing on an independent set. We then select the best performing model and use it to make predictions on more than 15,000 genes. Finally, we validate our predictions by mining the scientific literature for proposed therapeutic targets.

RESULTS

We observe that the data types with the best predictive power are animal models showing a disease-relevant phenotype, differential expression in diseased tissue and genetic association with the disease under investigation. On a test set, the neural network classifier achieves over 71% accuracy with an AUC of 0.76 when predicting therapeutic targets in a semi-supervised learning setting. We use this model to gain insights into current and failed programmes and to predict 1431 novel targets, of which a highly significant proportion has been independently proposed in the literature.

CONCLUSIONS

Our in silico approach shows that data linking genes and diseases is sufficient to predict novel therapeutic targets effectively and confirms that this type of evidence is essential for formulating or strengthening hypotheses in the target discovery process. Ultimately, more rapid and automated target prioritisation holds the potential to reduce both the costs and the development times associated with bringing new medicines to patients.

Global vision of druggability issues: applications and perspectives

During the preliminary stage of a drug discovery project, the lack of druggability information and poor target selection are the main causes of frequent failures. Elaborating on accurate computational druggability prediction methods is a requirement for prioritizing target selection, designing new drugs and avoiding side effects. In this review, we describe a survey of recently reported druggability prediction methods mainly based on networks, statistical pocket druggability predictions and virtual screening. An application for a frequent mutation of p53 tumor suppressor is presented, illustrating the complementarity of druggability prediction approaches, the remaining challenges and potential new drug development perspectives.

PPIC, EMP3 and CHI3L1 Are Novel Prognostic Markers for High Grade Glioma

Current treatment methods for patients diagnosed with gliomas have shown limited success. This is partly due to the lack of prognostic genes available to accurately predict disease outcomes. The aim of this study was to investigate novel prognostic genes based on the molecular profile of tumor samples and their correlation with clinical parameters. In the current study, microarray data (GSE4412 and GSE7696) downloaded from Gene Expression Omnibus were used to identify differentially expressed prognostic genes (DEPGs) by significant analysis of microarray (SAM) between long-term survivors (>2 years) and short-term survivors (≤2 years). DEPGs generated from these two datasets were intersected to obtain a list of common DEPGs. The expression of a subset of common DEPGs was then independently validated by real-time reverse transcription quantitative PCR (qPCR). Survival value of the common DEPGs was validated using known survival data from the GSE4412 and TCGA dataset. After intersecting DEPGs generated from the above two datasets, three genes were identified which may potentially be used to determine glioma patient prognosis. Independent validation with glioma patients tissue (n = 70) and normal brain tissue (n = 19) found PPIC, EMP3 and CHI3L1 were up-regulated in glioma tissue. Survival value validation showed that the three genes correlated with patient survival by Kaplan-Meir analysis, including grades, age and therapy.

A powerful weighted statistic for detecting group differences of directed biological networks

Complex disease is largely determined by a number of biomolecules interwoven into networks, rather than a single biomolecule. Different physiological conditions such as cases and controls may manifest as different networks. Statistical comparison between biological networks can provide not only new insight into the disease mechanism but statistical guidance for drug development. However, the methods developed in previous studies are inadequate to capture the changes in both the nodes and edges, and often ignore the network structure. In this study, we present a powerful weighted statistical test for group differences of directed biological networks, which is independent of the network attributes and can capture the changes in both the nodes and edges, as well as simultaneously accounting for the network structure through putting more weights on the difference of nodes locating on relatively more important position. Simulation studies illustrate that this method had better performance than previous ones under various sample sizes and network structures. One application to GWAS of leprosy successfully identifies the specific gene interaction network contributing to leprosy. Another real data analysis significantly identifies a new biological network, which is related to acute myeloid leukemia. One potential network responsible for lung cancer has also been significantly detected. The source R code is available on our website.

Transcriptional response networks for elucidating mechanisms of action of multitargeted agents

Drug discovery is moving away from the single target-based approach towards harnessing the potential of polypharmacological agents that modulate the activity of multiple nodes in the complex networks of deregulations underlying disease phenotypes. Computational network pharmacology methods that use systems-level drug-response phenotypes, such as those originating from genome-wide transcriptomic profiles, have proved particularly effective for elucidating the mechanisms of action of multitargeted compounds. Here, we show, via the case study of the natural product pinosylvin, how the combination of two complementary network-based methods can provide novel, unexpected mechanistic insights. This case study also illustrates that elucidating the mechanism of action of multitargeted natural products through transcriptional response-based approaches is a challenging endeavor, often requiring multiple computational-experimental iterations.

A powerful score-based statistical test for group difference in weighted biological networks

BACKGROUND

Complex disease is largely determined by a number of biomolecules interwoven into networks, rather than a single biomolecule. A key but inadequately addressed issue is how to test possible differences of the networks between two groups. Group-level comparison of network properties may shed light on underlying disease mechanisms and benefit the design of drug targets for complex diseases. We therefore proposed a powerful score-based statistic to detect group difference in weighted networks, which simultaneously capture the vertex changes and edge changes.

RESULTS

Simulation studies indicated that the proposed network difference measure (NetDifM) was stable and outperformed other methods existed, under various sample sizes and network topology structure. One application to real data about GWAS of leprosy successfully identified the specific gene interaction network contributing to leprosy. For additional gene expression data of ovarian cancer, two candidate subnetworks, PI3K-AKT and Notch signaling pathways, were considered and identified respectively.

CONCLUSIONS

The proposed method, accounting for the vertex changes and edge changes simultaneously, is valid and powerful to capture the group difference of biological networks.

Prediction of Druggable Proteins Using Machine Learning and Systems Biology: A Mini-Review

The emergence of -omics technologies has allowed the collection of vast amounts of data on biological systems. Although, the pace of such collection has been exponential, the impact of these data remains small on many critical biomedical applications such as drug development. Limited resources, high costs, and low hit-to-lead ratio have led researchers to search for more cost effective methodologies. A possible alternative is to incorporate computational methods of potential drug target prediction early during drug discovery workflow. Computational methods based on systems approaches have the advantage of taking into account the global properties of a molecule not limited to its sequence, structure or function. Machine learning techniques are powerful tools that can extract relevant information from massive and noisy data sets. In recent years the scientific community has explored the combined power of these fields to propose increasingly accurate and low cost methods to propose interesting drug targets. In this mini-review, we describe promising approaches based on the simultaneous use of systems biology and machine learning to access gene and protein druggability. Moreover, we discuss the state-of-the-art of this emerging and interdisciplinary field, discussing data sources, algorithms and the performance of the different methodologies. Finally, we indicate interesting avenues of research and some remaining open challenges.

Drug target prioritization by perturbed gene expression and network information

Drugs bind to their target proteins, which interact with downstream effectors and ultimately perturb the transcriptome of a cancer cell. These perturbations reveal information about their source, i.e., drugs’ targets. Here, we investigate whether these perturbations and protein interaction networks can uncover drug targets and key pathways. We performed the first systematic analysis of over 500 drugs from the Connectivity Map. First, we show that the gene expression of drug targets is usually not significantly affected by the drug perturbation. Hence, expression changes after drug treatment on their own are not sufficient to identify drug targets. However, ranking of candidate drug targets by network topological measures prioritizes the targets. We introduce a novel measure, local radiality, which combines perturbed genes and functional interaction network information. The new measure outperforms other methods in target prioritization and proposes cancer-specific pathways from drugs to affected genes for the first time. Local radiality identifies more diverse targets with fewer neighbors and possibly less side effects.

Network pharmacology applications to map the unexplored target space and therapeutic potential of natural products

It is widely accepted that drug discovery often requires a systems-level polypharmacology approach to tackle problems such as lack of efficacy and emerging resistance of single-targeted compounds. Network pharmacology approaches are increasingly being developed and applied to find new therapeutic opportunities and to re-purpose approved drugs. However, these recent advances have been relatively slow to be translated into the field of natural products. Here, we argue that a network pharmacology approach would enable an effective mapping of the yet unexplored target space of natural products, hence providing a systematic means to extend the druggable space of proteins implicated in various complex diseases. We give an overview of the key network pharmacology concepts and recent experimental-computational approaches that have been successfully applied to natural product research, including unbiased elucidation of mechanisms of action as well as systematic prediction of effective therapeutic combinations. We focus specifically on anticancer applications that use in vivo and in vitro functional phenotypic measurements, such as genome-wide transcriptomic response profiles, which enable a global modelling of the multi-target activity at the level of the biological pathways and interaction networks. We also provide representative examples of other disease applications, databases and tools as well as existing and emerging resources, which may prove useful for future natural product research. Finally, we offer our personal view of the current limitations, prospective developments and open questions in this exciting field.

Identification of drug candidates and repurposing opportunities through compound-target interaction networks

INTRODUCTION

System-wide identification of both on- and off-targets of chemical probes provides improved understanding of their therapeutic potential and possible adverse effects, thereby accelerating and de-risking drug discovery process. Given the high costs of experimental profiling of the complete target space of drug-like compounds, computational models offer systematic means for guiding these mapping efforts. These models suggest the most potent interactions for further experimental or pre-clinical evaluation both in cell line models and in patient-derived material.

AREAS COVERED

The authors focus here on network-based machine learning models and their use in the prediction of novel compound-target interactions both in target-based and phenotype-based drug discovery applications. While currently being used mainly in complementing the experimentally mapped compound-target networks for drug repurposing applications, such as extending the target space of already approved drugs, these network pharmacology approaches may also suggest completely unexpected and novel investigational probes for drug development.

EXPERT OPINION

Although the studies reviewed here have already demonstrated that network-centric modeling approaches have the potential to identify candidate compounds and selective targets in disease networks, many challenges still remain. In particular, these challenges include how to incorporate the cellular context and genetic background into the disease networks to enable more stratified and selective target predictions, as well as how to make the prediction models more realistic for the practical drug discovery and therapeutic applications.

Identifying the causative proteins of similar side effect pairs to explore the common molecular basis of these side effects

Drug side effects, or adverse drug reactions (ADRs), have become a major public health concern and often cause drug development failure and withdrawal. Some ADRs always occur concomitantly. Therefore, identifying these ADRs and their common molecular basis can better promote their prevention and treatment. In this paper we predicted the potential proteins for ADR pairs with similar mechanisms based on three layers of information: (i) the drug co-occurrence between a pair of ADRs; (ii) the correlation between a protein and an ADR pair based on the co-occurrence of drugs and (iii) the interaction between these proteins within the protein-protein interaction (PPI) network. The methods of randomization and functional annotation are used to investigate and analyze the relation between causative proteins and similar ADR pairs. The prediction accuracy of the relation between similar ADR pairs and related proteins reached 80%, and it increases with the number of drugs shared by the ADR pairs. From the ADR network made of single ADRs from predicted similar ADR pairs, we found that some ADRs are involved in multiple ADR pairs. The functional analysis of these ADR-related proteins suggests that a similar molecular basis is shared by multiple ADR pairs containing the same ADR, and these ADR pairs are almost caused by the same drug sets. The results of this study are reliable and provide a theoretical basis for the better prevention and treatment of ADRs that always occur concomitantly.

Galahad: a web server for drug effect analysis from gene expression

Galahad (https://galahad.esat.kuleuven.be) is a web-based application for analysis of drug effects. It provides an intuitive interface to be used by anybody interested in leveraging microarray data to gain insights into the pharmacological effects of a drug, mainly identification of candidate targets, elucidation of mode of action and understanding of off-target effects. The core of Galahad is a network-based analysis method of gene expression. As an input, Galahad takes raw Affymetrix human microarray data from treatment versus control experiments and provides quality control and data exploration tools, as well as computation of differential expression. Alternatively, differential expression values can be uploaded directly. Using these differential expression values, drug target prioritization and both pathway and disease enrichment can be calculated and visualized. Drug target prioritization is based on the integration of the gene expression data with a functional protein association network. The web site is free and open to all and there is no login requirement.

Identification of structural features in chemicals associated with cancer drug response: a systematic data-driven analysis

Motivation: Analysis of relationships of drug structure to biological response is key to understanding off-target and unexpected drug effects, and for developing hypotheses on how to tailor drug therapies. New methods are required for integrated analyses of a large number of chemical features of drugs against the corresponding genome-wide responses of multiple cell models.

Results: In this article, we present the first comprehensive multi-set analysis on how the chemical structure of drugs impacts on genome-wide gene expression across several cancer cell lines [Connectivity Map (CMap) database]. The task is formulated as searching for drug response components across multiple cancers to reveal shared effects of drugs and the chemical features that may be responsible. The components can be computed with an extension of a recent approach called Group Factor Analysis. We identify 11 components that link the structural descriptors of drugs with specific gene expression responses observed in the three cell lines and identify structural groups that may be responsible for the responses. Our method quantitatively outperforms the limited earlier methods on CMap and identifies both the previously reported associations and several interesting novel findings, by taking into account multiple cell lines and advanced 3D structural descriptors. The novel observations include: previously unknown similarities in the effects induced by 15-delta prostaglandin J2 and HSP90 inhibitors, which are linked to the 3D descriptors of the drugs; and the induction by simvastatin of leukemia-specific response, resembling the effects of corticosteroids.

Availability and implementation: Source Code implementing the method is available at: http://research.ics.aalto.fi/mi/software/GFAsparse

Contact:suleiman.khan@aalto.fi or samuel.kaski@aalto.fi

Supplementary Information:Supplementary data are available at Bioinformatics online.

Gene co-expression network and function modules in three types of glioma

The aim of the present study was to identify the disease‑associated genes and their functions involved in the development of three types of glioma (astrocytoma, glioblastoma and oligodendroglioma) with DNA microarray technology, and to analyze their differences and correlations. First, the gene expression profile GSE4290 was downloaded from the Gene Expression Omnibus database, then the probe‑level data were pre‑processed and the differentially expressed genes (DEGs) were identified with limma package in R language. Gene functions of the selected DEGs were further analyzed with the Database for Annotation, Visualization and Integrated Discovery. After the co‑expression network of DEGs was constructed by Cytoscape, the functional modules were mined and enrichment analysis was performed, and then the similarities and differences between any two types of glioma were compared. A total of 1151 genes between normal and astrocytoma tissues, 684 genes between normal and malignant glioma tissues, and 551 genes between normal and oligodendroglioma tissues were filtered as DEGs, respectively. By constructing co‑expression networks of DEGs, a total of 77232, 455 and 987 interactions were involved in the differentially co‑expressed networks of astrocytoma, oligodendroglioma and glioblastoma, respectively. The functions of DEGs were consistent with the modules in astrocytoma, glioblastoma and oligodendroglioma, which were mainly enriched in neuron signal transmission, immune responses and synthesis of organic acids, respectively. Model functions of astrocytoma and glioblastoma were similar (mainly related with immune response), while the model functions of oligodendroglioma differed markedly from that of the other two types. The identification of the associations among these three types of glioma has potential clinical utility for improving the diagnosis of different types of glioma in the future. In addition, these results have marked significance in studying the underlying mechanisms, distinguishing between normal and cancer tissues, and examining novel therapeutic strategies for patients with glioma.

Pediatric perspective on pharmacogenomics

The advances in high-throughput genomic technologies have improved the understanding of disease pathophysiology and have allowed a better characterization of drug response and toxicity based on individual genetic make up. Pharmacogenomics is being recognized as a valid approach used to identify patients who are more likely to respond to medication, or those in whom there is a high probability of developing severe adverse drug reactions. An increasing number of pharmacogenomic studies are being published, most include only adults. A few studies have shown the impact of pharmacogenomics in pediatrics, highlighting a key difference between children and adults, which is the contribution of developmental changes to therapeutic responses across different age groups. This review focuses on pharmacogenomic research in pediatrics, providing examples from common pediatric conditions and emphasizing their developmental context.

Targeting molecular networks for drug research

The study of molecular networks has recently moved into the limelight of biomedical research. While it has certainly provided us with plenty of new insights into cellular mechanisms, the challenge now is how to modify or even restructure these networks. This is especially true for human diseases, which can be regarded as manifestations of distorted states of molecular networks. Of the possible interventions for altering networks, the use of drugs is presently the most feasible. In this mini-review, we present and discuss some exemplary approaches of how analysis of molecular interaction networks can contribute to pharmacology (e.g., by identifying new drug targets or prediction of drug side effects), as well as list pointers to relevant resources and software to guide future research. We also outline recent progress in the use of drugs for in vitro reprogramming of cells, which constitutes an example par excellence for altering molecular interaction networks with drugs.

Causal Network Models for Predicting Compound Targets and Driving Pathways in Cancer

Gene-expression data are often used to infer pathways regulating transcriptional responses. For example, differentially expressed genes (DEGs) induced by compound treatment can help characterize hits from phenotypic screens, either by correlation with known drug signatures or by pathway enrichment. Pathway enrichment is, however, typically computed with DEGs rather than "upstream" nodes that are potentially causal of "downstream" changes. Here, we present graph-based models to predict causal targets from compound-microarray data. We test several approaches to traversing network topology, and show that a consensus minimum-rank score (SigNet) beat individual methods and could highly rank compound targets among all network nodes. In addition, larger, less canonical networks outperformed linear canonical interactions. Importantly, pathway enrichment using causal nodes rather than DEGs recovers relevant pathways more often. To further validate our approach, we used integrated data sets from the Cancer Genome Atlas to identify driving pathways in triple-negative breast cancer. Critical pathways were uncovered, including the epidermal growth factor receptor 2-phosphatidylinositide 3-kinase-AKT-MAPK growth pathway andATR-p53-BRCA DNA damage pathway, in addition to unexpected pathways, such as TGF-WNT cytoskeleton remodeling, IL12-induced interferon gamma production, and TNFR-IAP (inhibitor of apoptosis) apoptosis; the latter was validated by pooled small hairpin RNA profiling in cancer cells. Overall, our approach can bridge transcriptional profiles to compound targets and driving pathways in cancer.