Integrating systems biology sources illuminates drug action

Drug repositioning using drug-disease vectors based on an integrated network

BACKGROUND

Diverse interactions occur between biomolecules, such as activation, inhibition, expression, or repression. However, previous network-based studies of drug repositioning have employed interaction on the binary protein-protein interaction (PPI) network without considering the characteristics of the interactions. Recently, some studies of drug repositioning using gene expression data found that associations between drug and disease genes are useful information for identifying novel drugs to treat diseases. However, the gene expression profiles for drugs and diseases are not always available. Although gene expression profiles of drugs and diseases are available, existing methods cannot use the drugs or diseases, when differentially expressed genes in the profiles are not included in their network.

RESULTS

We developed a novel method for identifying candidate indications of existing drugs considering types of interactions between biomolecules based on known drug-disease associations. To obtain associations between drug and disease genes, we constructed a directed network using protein interaction and gene regulation data obtained from various public databases providing diverse biological pathways. The network includes three types of edges depending on relationships between biomolecules. To quantify the association between a target gene and a disease gene, we explored the shortest paths from the target gene to the disease gene and calculated the types and weights of the shortest paths. For each drug-disease pair, we built a vector consisting of values for each disease gene influenced by the drug. Using the vectors and known drug-disease associations, we constructed classifiers to identify novel drugs for each disease.

CONCLUSION

We propose a method for exploring candidate drugs of diseases using associations between drugs and disease genes derived from a directed gene network instead of gene regulation data obtained from gene expression profiles. Compared to existing methods that require information on gene relationships and gene expression data, our method can be applied to a greater number of drugs and diseases. Furthermore, to validate our predictions, we compared the predictions with drug-disease pairs in clinical trials using the hypergeometric test, which showed significant results. Our method also showed better performance compared to existing methods for the area under the receiver operating characteristic curve (AUC).

Predictive Systems Toxicology

In this review we address to what extent computational techniques can augment our ability to predict toxicity. The first section provides a brief history of empirical observations on toxicity dating back to the dawn of Sumerian civilization. Interestingly, the concept of dose emerged very early on, leading up to the modern emphasis on kinetic properties, which in turn encodes the insight that toxicity is not solely a property of a compound but instead depends on the interaction with the host organism. The next logical step is the current conception of evaluating drugs from a personalized medicine point of view. We review recent work on integrating what could be referred to as classical pharmacokinetic analysis with emerging systems biology approaches incorporating multiple omics data. These systems approaches employ advanced statistical analytical data processing complemented with machine learning techniques and use both pharmacokinetic and omics data. We find that such integrated approaches not only provide improved predictions of toxicity but also enable mechanistic interpretations of the molecular mechanisms underpinning toxicity and drug resistance. We conclude the chapter by discussing some of the main challenges, such as how to balance the inherent tension between the predicitive capacity of models, which in practice amounts to constraining the number of features in the models versus allowing for rich mechanistic interpretability, i.e., equipping models with numerous molecular features. This challenge also requires patient-specific predictions on toxicity, which in turn requires proper stratification of patients as regards how they respond, with or without adverse toxic effects. In summary, the transformation of the ancient concept of dose is currently successfully operationalized using rich integrative data encoded in patient-specific models.

Cohort-specific imputation of gene expression improves prediction of warfarin dose for African Americans

Background

Genome-wide association studies are useful for discovering genotype–phenotype associations but are limited because they require large cohorts to identify a signal, which can be population-specific. Mapping genetic variation to genes improves power and allows the effects of both protein-coding variation as well as variation in expression to be combined into “gene level” effects.

Methods

Previous work has shown that warfarin dose can be predicted using information from genetic variation that affects protein-coding regions.

Here, we introduce a method that improves dose prediction by integrating tissue-specific gene expression. In particular, we use drug pathways and expression quantitative trait loci knowledge to impute gene expression—on the assumption that differential expression of key pathway genes may impact dose requirement. We focus on 116 genes from the pharmacokinetic and pharmacodynamic pathways of warfarin within training and validation sets comprising both European and African-descent individuals.

Results

We build gene-tissue signatures associated with warfarin dose in a cohort-specific manner and identify a signature of 11 gene-tissue pairs that significantly augments the International Warfarin Pharmacogenetics Consortium dosage-prediction algorithm in both populations.

Conclusions

Our results demonstrate that imputed expression can improve dose prediction and bridge population-specific compositions.

MATLAB code is available at https://github.com/assafgo/warfarin-cohort

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0495-0) contains supplementary material, which is available to authorized users.

Drug voyager: a computational platform for exploring unintended drug action

BACKGROUND

The dominant paradigm in understanding drug action focuses on the intended therapeutic effects and frequent adverse reactions. However, this approach may limit opportunities to grasp unintended drug actions, which can open up channels to repurpose existing drugs and identify rare adverse drug reactions. Advances in systems biology can be exploited to comprehensively understand pharmacodynamic actions, although proper frameworks to represent drug actions are still lacking.

RESULTS

We suggest a novel platform to construct a drug-specific pathway in which a molecular-level mechanism of action is formulated based on pharmacologic, pharmacogenomic, transcriptomic, and phenotypic data related to drug response ( http://databio.gachon.ac.kr/tools/ ). In this platform, an adoption of three conceptual levels imitating drug perturbation allows these pathways to be realistically rendered in comparison to those of other models. Furthermore, we propose a new method that exploits functional features of the drug-specific pathways to predict new indications as well as adverse reactions. For therapeutic uses, our predictions significantly overlapped with clinical trials and an up-to-date drug-disease association database. Also, our method outperforms existing methods with regard to classification of active compounds for cancers. For adverse reactions, our predictions were significantly enriched in an independent database derived from the Food and Drug Administration (FDA) Adverse Event Reporting System and meaningfully cover an Adverse Reaction Database provided by Health Canada. Lastly, we discuss several predictions for both therapeutic indications and side-effects through the published literature.

CONCLUSIONS

Our study addresses how we can computationally represent drug-signaling pathways to understand unintended drug actions and to facilitate drug discovery and screening.

DIFFERENTIAL PATHWAY DEPENDENCY DISCOVERY ASSOCIATED WITH DRUG RESPONSE ACROSS CANCER CELL LINES

The effort to personalize treatment plans for cancer patients involves the identification of drug treatments that can effectively target the disease while minimizing the likelihood of adverse reactions. In this study, the gene-expression profile of 810 cancer cell lines and their response data to 368 small molecules from the Cancer Therapeutics Research Portal (CTRP) are analyzed to identify pathways with significant rewiring between genes, or differential gene dependency, between sensitive and non-sensitive cell lines. Identified pathways and their corresponding differential dependency networks are further analyzed to discover essentiality and specificity mediators of cell line response to drugs/compounds. For analysis we use the previously published method EDDY (Evaluation of Differential DependencY). EDDY first constructs likelihood distributions of gene-dependency networks, aided by known genegene interaction, for two given conditions, for example, sensitive cell lines vs. non-sensitive cell lines. These sets of networks yield a divergence value between two distributions of network likelihoods that can be assessed for significance using permutation tests. Resulting differential dependency networks are then further analyzed to identify genes, termed mediators, which may play important roles in biological signaling in certain cell lines that are sensitive or non-sensitive to the drugs. Establishing statistical correspondence between compounds and mediators can improve understanding of known gene dependencies associated with drug response while also discovering new dependencies. Millions of compute hours resulted in thousands of these statistical discoveries. EDDY identified 8,811 statistically significant pathways leading to 26,822 compound-pathway-mediator triplets. By incorporating STITCH and STRING databases, we could construct evidence networks for 14,415 compound-pathway-mediator triplets for support. The results of this analysis are presented in a searchable website to aid researchers in studying potential molecular mechanisms underlying cells' drug response as well as in designing experiments for the purpose of personalized treatment regimens.

A Multi-scale Computational Platform to Mechanistically Assess the Effect of Genetic Variation on Drug Responses in Human Erythrocyte Metabolism

Progress in systems medicine brings promise to addressing patient heterogeneity and individualized therapies. Recently, genome-scale models of metabolism have been shown to provide insight into the mechanistic link between drug therapies and systems-level off-target effects while being expanded to explicitly include the three-dimensional structure of proteins. The integration of these molecular-level details, such as the physical, structural, and dynamical properties of proteins, notably expands the computational description of biochemical network-level properties and the possibility of understanding and predicting whole cell phenotypes. In this study, we present a multi-scale modeling framework that describes biological processes which range in scale from atomistic details to an entire metabolic network. Using this approach, we can understand how genetic variation, which impacts the structure and reactivity of a protein, influences both native and drug-induced metabolic states. As a proof-of-concept, we study three enzymes (catechol-O-methyltransferase, glucose-6-phosphate dehydrogenase, and glyceraldehyde-3-phosphate dehydrogenase) and their respective genetic variants which have clinically relevant associations. Using all-atom molecular dynamic simulations enables the sampling of long timescale conformational dynamics of the proteins (and their mutant variants) in complex with their respective native metabolites or drug molecules. We find that changes in a protein's structure due to a mutation influences protein binding affinity to metabolites and/or drug molecules, and inflicts large-scale changes in metabolism.

Vascular endothelial growth factor A as predictive marker for mTOR inhibition in relapsing high-grade serous ovarian cancer

Background

Development of resistance against first line drug therapy including cisplatin and paclitaxel in high-grade serous ovarian cancer (HGSOC) presents a major challenge. Identifying drug candidates breaking resistance, ideally combined with predictive biomarkers allowing precision use are needed for prolonging progression free survival of ovarian cancer patients.

Modeling of molecular processes driving drug resistance in tumor tissue further combined with mechanism of action of drugs provides a strategy for identification of candidate drugs and associated predictive biomarkers.

Results

Consolidation of transcriptomics profiles and biomedical literature mining results provides 1242 proteins linked with ovarian cancer drug resistance. Integrating this set on a protein interaction network followed by graph segmentation results in a molecular process model representation of drug resistant HGSOC embedding 409 proteins in 24 molecular processes. Utilizing independent transcriptomics profiles with follow-up data on progression free survival allows deriving molecular biomarker-based classifiers for predicting recurrence under first line therapy. Biomarkers of specific relevance are identified in a molecular process encapsulating TGF-beta, mTOR, Jak-STAT and Neurotrophin signaling. Mechanism of action molecular model representations of cisplatin and paclitaxel embed the very same signaling components, and specifically proteins afflicted with the activation status of the mTOR pathway become evident, including VEGFA. Analyzing mechanism of action interference of the mTOR inhibitor sirolimus shows specific impact on the drug resistance signature imposed by cisplatin and paclitaxel, further holding evidence for a synthetic lethal interaction to paclitaxel mechanism of action involving cyclin D1.

Conclusions

Stratifying drug resistant high grade serous ovarian cancer via VEGFA, and specifically treating with mTOR inhibitors in case of activation of the pathway may allow adding precision for overcoming resistance to first line therapy.

The acute impact of polyphenols from Hibiscus sabdariffa in metabolic homeostasis: an approach combining metabolomics and gene-expression analyses

We explored the acute multifunctional effects of polyphenols from Hibiscus sabdariffa in humans to assess possible consequences on the host's health. The expected dynamic response was studied using a combination of transcriptomics and metabolomics to integrate specific functional pathways through network-based methods and to generate hypotheses established by acute metabolic effects and/or modifications in the expression of relevant genes. Data were obtained from healthy male volunteers after 3 hours of ingestion of an aqueous Hibiscus sabdariffa extract. The data were compared with data obtained prior to the ingestion, and the overall findings suggest that these particular polyphenols had a simultaneous role in mitochondrial function, energy homeostasis and protection of the cardiovascular system. These findings suggest beneficial actions in inflammation, endothelial dysfunction, and oxidation, which are interrelated mechanisms. Among other effects, the activation of the heme oxygenase-biliverdin reductase axis, the systemic inhibition of the renin-angiotensin system, the inhibition of the angiotensin-converting enzyme, and several actions mirroring those of the peroxisome proliferator-activated receptor agonists further support this notion. We also found concordant findings in the serum of the participants, which include a decrease in cortisol levels and a significant increase in the active vasodilator metabolite of bradykinin (des-Arg(9)-bradykinin). Therefore, our data support the view that polyphenols from Hibiscus sabdariffa play a regulatory role in metabolic health and in the maintenance of blood pressure, thus implying a multi-faceted impact in metabolic and cardiovascular diseases.

Systems biology approaches for identifying adverse drug reactions and elucidating their underlying biological mechanisms

Small molecules are indispensable to modern medical therapy. However, their use may lead to unintended, negative medical outcomes commonly referred to as adverse drug reactions (ADRs). These effects vary widely in mechanism, severity, and populations affected, making ADR prediction and identification important public health concerns. Current methods rely on clinical trials and postmarket surveillance programs to find novel ADRs; however, clinical trials are limited by small sample size, whereas postmarket surveillance methods may be biased and inherently leave patients at risk until sufficient clinical evidence has been gathered. Systems pharmacology, an emerging interdisciplinary field combining network and chemical biology, provides important tools to uncover and understand ADRs and may mitigate the drawbacks of traditional methods. In particular, network analysis allows researchers to integrate heterogeneous data sources and quantify the interactions between biological and chemical entities. Recent work in this area has combined chemical, biological, and large-scale observational health data to predict ADRs in both individual patients and global populations. In this review, we explore the rapid expansion of systems pharmacology in the study of ADRs. We enumerate the existing methods and strategies and illustrate progress in the field with a model framework that incorporates crucial data elements, such as diet and comorbidities, known to modulate ADR risk. Using this framework, we highlight avenues of research that may currently be underexplored, representing opportunities for future work.

CUSP9* treatment protocol for recurrent glioblastoma: aprepitant, artesunate, auranofin, captopril, celecoxib, disulfiram, itraconazole, ritonavir, sertraline augmenting continuous low dose temozolomide

CUSP9 treatment protocol for recurrent glioblastoma was published one year ago. We now present a slight modification, designated CUSP9*. CUSP9* drugs--aprepitant, artesunate, auranofin, captopril, celecoxib, disulfiram, itraconazole, sertraline, ritonavir, are all widely approved by regulatory authorities, marketed for non-cancer indications. Each drug inhibits one or more important growth-enhancing pathways used by glioblastoma. By blocking survival paths, the aim is to render temozolomide, the current standard cytotoxic drug used in primary glioblastoma treatment, more effective. Although esthetically unpleasing to use so many drugs at once, the closely similar drugs of the original CUSP9 used together have been well-tolerated when given on a compassionate-use basis in the cases that have come to our attention so far. We expect similarly good tolerability for CUSP9*. The combined action of this suite of drugs blocks signaling at, or the activity of, AKT phosphorylation, aldehyde dehydrogenase, angiotensin converting enzyme, carbonic anhydrase -2,- 9, -12, cyclooxygenase-1 and -2, cathepsin B, Hedgehog, interleukin-6, 5-lipoxygenase, matrix metalloproteinase -2 and -9, mammalian target of rapamycin, neurokinin-1, p-gp efflux pump, thioredoxin reductase, tissue factor, 20 kDa translationally controlled tumor protein, and vascular endothelial growth factor. We believe that given the current prognosis after a glioblastoma has recurred, a trial of CUSP9* is warranted.

The potential of translational bioinformatics approaches for pharmacology research

The field of bioinformatics has allowed the interpretation of massive amounts of biological data, ushering in the era of 'omics' to biomedical research. Its potential impact on pharmacology research is enormous and it has shown some emerging successes. A full realization of this potential, however, requires standardized data annotation for large health record databases and molecular data resources. Improved standardization will further stimulate the development of system pharmacology models, using translational bioinformatics methods. This new translational bioinformatics paradigm is highly complementary to current pharmacological research fields, such as personalized medicine, pharmacoepidemiology and drug discovery. In this review, I illustrate the application of transformational bioinformatics to research in numerous pharmacology subdisciplines.

Ranking adverse drug reactions with crowdsourcing

BACKGROUND

There is no publicly available resource that provides the relative severity of adverse drug reactions (ADRs). Such a resource would be useful for several applications, including assessment of the risks and benefits of drugs and improvement of patient-centered care. It could also be used to triage predictions of drug adverse events.

OBJECTIVE

The intent of the study was to rank ADRs according to severity.

METHODS

We used Internet-based crowdsourcing to rank ADRs according to severity. We assigned 126,512 pairwise comparisons of ADRs to 2589 Amazon Mechanical Turk workers and used these comparisons to rank order 2929 ADRs.

RESULTS

There is good correlation (rho=.53) between the mortality rates associated with ADRs and their rank. Our ranking highlights severe drug-ADR predictions, such as cardiovascular ADRs for raloxifene and celecoxib. It also triages genes associated with severe ADRs such as epidermal growth-factor receptor (EGFR), associated with glioblastoma multiforme, and SCN1A, associated with epilepsy.

CONCLUSIONS

ADR ranking lays a first stepping stone in personalized drug risk assessment. Ranking of ADRs using crowdsourcing may have useful clinical and financial implications, and should be further investigated in the context of health care decision making.

Systems pharmacology augments drug safety surveillance

Small molecule drugs are the foundation of modern medical practice yet their use is limited by the onset of unexpected and severe adverse events (AEs). Regulatory agencies rely on post-marketing surveillance to monitor safety once drugs are approved for clinical use. Despite advances in pharmacovigilance methods that address issues of confounding bias, clinical data of AEs are inherently noisy. Systems pharmacology– the integration of systems biology and chemical genomics – can illuminate drug mechanisms of action. We hypothesize that these data can improve drug safety surveillance by highlighting drugs with a mechanistic connection to the target phenotype (enriching true positives) and filtering those that do not (depleting false positives). We present an algorithm, the modular assembly of drug safety subnetworks (MADSS), to combine systems pharmacology and pharmacovigilance data and significantly improve drug safety monitoring for four clinically relevant adverse drug reactions.