MD-Miner: a network-based approach for personalized drug repositioning

Background

Due to advances in next generation sequencing technologies and corresponding reductions in cost, it is now attainable to investigate genome-wide gene expression and variants at a patient-level, so as to better understand and anticipate heterogeneous responses to therapy. Consequently, it is feasible to inform personalized drug treatment decisions using personal genomics data. However, these efforts are limited due to a lack of reliable computational approaches for predicting effective drugs for individual patients. The reverse gene set enrichment analysis (i.e., connectivity mapping) approach and its variants have been widely and successfully used for drug prediction. However, the performance of these methods is limited by undefined mechanism of action (MoA) of drugs and reliance on cohorts of patients rather than personalized predictions for individual patients.

Results

In this study, we have developed and evaluated a computational approach, known as Mechanism and Drug Miner (MD-Miner), using a network-based computational approach to predict effective drugs and reveal potential drug mechanisms of action at the level of signaling pathways. Specifically, the patient-specific signaling network is constructed by integrating known disease associated genes with patient-derived gene expression profiles. In parallel, a drug mechanism of action network is constructed by integrating drug targets and z-score profiles of drug-induced gene expression (pre vs. post-drug treatment). Potentially effective candidate drugs are prioritized according to the number of common genes between the patient-specific dysfunctional signaling network and drug MoA network. We evaluated the MD-Miner method on the PC-3 prostate cancer cell line, and showed that it significantly improved the success rate of discovering effective drugs compared with the random selection, and could provide insight into potential mechanisms of action.

Conclusions

This work provides a signaling network-based drug repositioning approach. Compared with the reverse gene signature based drug repositioning approaches, the proposed method can provide clues of mechanism of action in terms of signaling transduction networks.

     

Analysis of bypass signaling in EGFR pathway and profiling of bypass genes for predicting response to anticancer EGFR tyrosine kinase inhibitors

Some drugs, such as anticancer EGFR tyrosine kinase inhibitors, elicit markedly different clinical response rates due to differences in drug bypass signaling as well as genetic variations of drug target and downstream drug-resistant genes. The profiles of these bypass signaling are expected to be useful for improved drug response prediction, which have not been systematically explored previously. In this work, we searched and analyzed 16 literature-reported EGFR tyrosine kinase inhibitor bypass signaling routes in the EGFR pathway, which include 5 compensatory routes of EGFR transactivation by another receptor, and 11 alternative routes activated by another receptor. These 16 routes are reportedly regulated by 11 bypass genes. Their expression profiles together with the mutational, amplification and expression profiles of EGFR and 4 downstream drug-resistant genes, were used as new sets of biomarkers for identifying 53 NSCLC cell-lines sensitive or resistant to EGFR tyrosine kinase inhibitors gefitinib, erlotinib and lapatinib. The collective profiles of all 16 genes distinguish sensitive and resistant cell-lines are better than those of individual genes and the combined EGFR and downstream drug resistant genes, and their derived cell-line response rates are consistent with the reported clinical response rates of the three drugs. The usefulness of cell-line data for drug response studies was further analyzed by comparing the expression profiles of EGFR and bypass genes in NSCLC cell-lines and patient samples, and by using a machine learning feature selection method for selecting drug response biomarkers. Our study suggested that the profiles of drug bypass signaling are highly useful for improved drug response prediction.     

Large-scale elucidation of drug response pathways in humans

Elucidating signaling pathways is a fundamental step in understanding cellular processes and developing new therapeutic strategies. Here we introduce a method for the large-scale elucidation of signaling pathways involved in cellular response to drugs. Combining drug targets, drug response expression profiles, and the human physical interaction network, we infer 99 human drug response pathways and study their properties. Based on the newly inferred pathways, we develop a pathway-based drug-drug similarity measure and compare it to two common, gold standard drug-drug similarity measures. Remarkably, our measure provides better correspondence to these gold standards than similarity measures that are based on associations between drugs and known pathways, or on drug-specific gene expression profiles. It further improves the prediction of drug side effects and indications, elucidating specific response pathways that may be associated with these drug properties. Supplementary Material for this article is available at www.liebertonline.com/cmb.     

Simplified antimalarial therapeutic monitoring: using the day-7 drug level?

The blood concentration profiles of most antimalarial drugs vary considerably between patients. The interpretation of antimalarial drug trials evaluating efficacy and effectiveness would be improved considerably if the exposure of the infecting parasite population to the antimalarial drug treatment could be measured. Artemisinin combination treatments are now recommended as first-line drugs for the treatment of falciparum malaria. Measurement of the blood, serum or plasma concentration of the slowly eliminated partner antimalarial drug on day 7 of follow-up is simpler and might be a better determinant of therapeutic response than the area under the concentration-time curve. Measurement of the day-7 drug level should be considered as a routine part of antimalarial drug trials.     

Therapeutic resistance resulting from mutations in Raf/MEK/ERK and PI3K/PTEN/Akt/mTOR signaling pathways

Chemotherapy remains a commonly used therapeutic approach for many cancers. Indeed chemotherapy is relatively effective for treatment of certain cancers and it may be the only therapy (besides radiotherapy) that is appropriate for certain cancers. However, a common problem with chemotherapy is the development of drug resistance. Many studies on the mechanisms of drug resistance concentrated on the expression of membrane transporters and how they could be aberrantly regulated in drug resistant cells. Attempts were made to isolate specific inhibitors which could be used to treat drug resistant patients. Unfortunately most of these drug transporter inhibitors have not proven effective for therapy. Recently the possibilities of more specific, targeted therapies have sparked the interest of clinical and basic researchers as approaches to kill cancer cells. However, there are also problems associated with these targeted therapies. Two key signaling pathways involved in the regulation of cell growth are the Ras/Raf/MEK/ERK and PI3K/PTEN/Akt/mTOR pathways. Dysregulated signaling through these pathways is often the result of genetic alterations in critical components in these pathways as well as mutations in upstream growth factor receptors. Furthermore, these pathways may be activated by chemotherapeutic drugs and ionizing radiation. This review documents how their abnormal expression can contribute to drug resistance as well as resistance to targeted therapy. This review will discuss in detail PTEN regulation as this is a critical tumor suppressor gene frequently dysregulated in human cancer which contributes to therapy resistance. Controlling the expression of these pathways could improve cancer therapy and ameliorate human health.     

Role of TGF-β signaling regulatory microRNAs in the pathogenesis of colorectal cancer

Transforming growth factor β (TGF-β) modulates tumor progression by regulating cell proliferation, apoptosis, metastasis, angiogenesis, and drug resistance. Biological and pharmacological agonists/antagonists, the interplay between intracellular signaling pathways, and microRNAs (miRNAs) control the initiation and activation of the TGF-β signaling pathway. It has been proposed that the expression profiles of tumor suppressor and oncogenic TGF-β miRNAs may be used for the classification, diagnosis, and prognosis of human malignancies. Deregulated miRNAs and aberrant activation of TGF-β signaling are frequently found in human colorectal cancers (CRCs), but a little is known about their mechanisms of action in the development and progression of colorectal carcinoma. This review summarizes the current knowledge of the role of TGF-β signaling regulatory miRNAs in the pathogenesis of CRC for a better understanding and hence better management of this disease.     

Integrated analysis of numerous heterogeneous gene expression profiles for detecting robust disease-specific biomarkers and proposing drug targets

Genome-wide expression profiling has revolutionized biomedical research; vast amounts of expression data from numerous studies of many diseases are now available. Making the best use of this resource in order to better understand disease processes and treatment remains an open challenge. In particular, disease biomarkers detected in case–control studies suffer from low reliability and are only weakly reproducible. Here, we present a systematic integrative analysis methodology to overcome these shortcomings. We assembled and manually curated more than 14 000 expression profiles spanning 48 diseases and 18 expression platforms. We show that when studying a particular disease, judicious utilization of profiles from other diseases and information on disease hierarchy improves classification quality, avoids overoptimistic evaluation of that quality, and enhances disease-specific biomarker discovery. This approach yielded specific biomarkers for 24 of the analyzed diseases. We demonstrate how to combine these biomarkers with large-scale interaction, mutation and drug target data, forming a highly valuable disease summary that suggests novel directions in disease understanding and drug repurposing. Our analysis also estimates the number of samples required to reach a desired level of biomarker stability. This methodology can greatly improve the exploitation of the mountain of expression profiles for better disease analysis.     

181 Integrative analysis of gene expression and protein–protein interaction networks in Glioblastoma

Glioblastoma multiforme (GBM) is the most malignant of all the brain tumors with very low median survival time of one year, as per Central Brain Tumor Registry of the USA, 2001. Efforts are ongoing to understand this disease pathogenesis in complete details. Global gene expression changes in GBM pathogenesis have been studied by several groups using microarray technology (e.g. Carro et al., 2010). One of the many approaches to ‘understand the control mechanisms underlying the observed changes in the activity of a biological process’ (Cline et al., 2007) is integration of gene expression and protein–protein interactions (PPI) datasets. Among several examples, aberrant activation of Wnt/β-catenin signaling pathway as well as sonic hedgehog (SHH) signaling pathway is reported in GBMs (Klaus & Birchmeier, 2008). Further, these two pathways are also involved in proliferation and clonogenicity of glioma cancer stem cells (Li et al., 2009), which are thought to play a role in glioma initiation, proliferation, and invasion, and are one of the important points of intervention. Hedgehog–Gli1 signaling is also found to regulate the expression of stemness genes. In this paper, analyses of the relationship between the significant differential expression of these and other genes and the connectivity as well as topological features of a PPI network would be discussed. This way, genes potentially overlooked when relying solely on expression profiles may be identified which can be biologically relevant as possible drug target/s or disease biomarker/s.     

Gene expression profiling and its practice in drug development

The availability of sequenced genomes of human and many experimental animals necessitated the development of new technologies and powerful computational tools that are capable of exploiting these genomic data and ask intriguing questions about complex nature of biological processes. This gave impetus for developing whole genome approaches that can produce functional information of genes in the form of expression profiles and unscramble the relationships between variation in gene expression and the resulting physiological outcome. These profiles represent genetic fingerprints or catalogue of genes that characterize the cell or tissue being studied and provide a basis from which to begin an investigation of the underlying biology. Among the most powerful and versatile tools are high-density DNA microarrays to analyze the expression patterns of large numbers of genes across different tissues or within the same tissue under a variety of experimental conditions or even between species. The wide spread use of microarray technologies is generating large sets of data that is stimulating the development of better analytical tools so that functions can be predicted for novel genes. In this review, the authors discuss how these profiles are being used at various stages of the drug discovery process and help in the identification of new drug targets, predict the function of novel genes, and understand individual variability in response to drugs.     

Inferring pathway activity toward precise disease classification

The advent of microarray technology has made it possible to classify disease states based on gene expression profiles of patients. Typically, marker genes are selected by measuring the power of their expression profiles to discriminate among patients of different disease states. However, expression-based classification can be challenging in complex diseases due to factors such as cellular heterogeneity within a tissue sample and genetic heterogeneity across patients. A promising technique for coping with these challenges is to incorporate pathway information into the disease classification procedure in order to classify disease based on the activity of entire signaling pathways or protein complexes rather than on the expression levels of individual genes or proteins. We propose a new classification method based on pathway activities inferred for each patient. For each pathway, an activity level is summarized from the gene expression levels of its condition-responsive genes (CORGs), defined as the subset of genes in the pathway whose combined expression delivers optimal discriminative power for the disease phenotype. We show that classifiers using pathway activity achieve better performance than classifiers based on individual gene expression, for both simple and complex case-control studies including differentiation of perturbed from non-perturbed cells and subtyping of several different kinds of cancer. Moreover, the new method outperforms several previous approaches that use a static (i.e., non-conditional) definition of pathways. Within a pathway, the identified CORGs may facilitate the development of better diagnostic markers and the discovery of core alterations in human disease.