Reversal of cancer gene expression correlates with drug efficacy and reveals therapeutic targets

Combinatorial prediction of therapeutic perturbations using causally-inspired neural networks

Phenotype-driven approaches identify disease-counteracting compounds by analyzing the phenotypic signatures that distinguish diseased from healthy states. These approaches can guide the discovery of targeted perturbations, including small-molecule drugs and genetic interventions, that modulate disease phenotypes toward healthier states. Here, we introduce PDGrapher, a causally inspired graph neural network (GNN) designed to predict combinatorial perturbagens (sets of therapeutic targets) capable of reversing disease phenotypes. Unlike methods that learn how perturbations alter phenotypes, PDGrapher solves the inverse problem of directly predicting the perturbagens needed to achieve a desired response. PDGrapher is a GNN that embeds disease cell states into gene regulatory or protein-protein interaction networks, learns a latent representation of these states, and identifies the optimal combinatorial perturbations that most effectively shift the diseased state toward the desired treated state within that latent space. In experiments in nine cell lines with chemical perturbations, PDGrapher identified effective per-turbagens in up to 13.33% more test samples than competing methods and achieved a normalized discounted cumulative gain of up to 0.12 higher to classify therapeutic targets. It also demonstrated competitive performance on ten genetic perturbation datasets. A key advantage of PDGrapher is its direct prediction paradigm, in contrast to the indirect and computationally intensive models traditionally employed in phenotype-driven research. This approach accelerates training by up to 25 times compared to existing methods. PDGrapher provides a fast approach for identifying therapeutic perturbations and advancing phenotype-driven drug discovery.

Reversal gene expression assessment for drug repurposing, a case study of glioblastoma

BACKGROUND

Glioblastoma (GBM) is a rare brain cancer with an exceptionally high mortality rate, which illustrates the pressing demand for more effective therapeutic options. Despite considerable research efforts on GBM, its underlying biological mechanisms remain unclear. Furthermore, none of the United States Food and Drug Administration (FDA) approved drugs used for GBM deliver satisfactory survival improvement.

METHODS

This study presents a novel computational pipeline by utilizing gene expression data analysis for GBM for drug repurposing to address the challenges in rare disease drug development, particularly focusing on GBM. The GBM Gene Expression Profile (GGEP) was constructed with multi-omics data to identify drugs with reversal gene expression to GGEP from the Integrated Network-Based Cellular Signatures (iLINCS) database.

RESULTS

We prioritized the candidates via hierarchical clustering of their expression signatures and quantification of their reversal strength by calculating two self-defined indices based on the GGEP genes' log2 foldchange (LFC) that the drug candidates could induce. Among five prioritized candidates, in-vitro experiments validated Clofarabine and Ciclopirox as highly efficacious in selectively targeting GBM cancer cells.

CONCLUSIONS

The success of this study illustrated a promising avenue for accelerating drug development by uncovering underlying gene expression effect between drugs and diseases, which can be extended to other rare diseases and non-rare diseases.

HERB 2.0: an updated database integrating clinical and experimental evidence for traditional Chinese medicine

Clinical trials and meta-analyses are considered high-level medical evidence with solid credibility. However, such clinical evidence for traditional Chinese medicine (TCM) is scattered, requiring a unified entrance to navigate all available evaluations on TCM therapies under modern standards. Besides, novel experimental evidence has continuously accumulated for TCM since the publication of HERB 1.0. Therefore, we updated the HERB database to integrate four types of evidence for TCM: (i) we curated 8558 clinical trials and 8032 meta-analyses information for TCM and extracted clear clinical conclusions for 1941 clinical trials and 593 meta-analyses with companion supporting papers. (ii) we updated experimental evidence for TCM, increased the number of high-throughput experiments to 2231, and curated references to 6 644. We newly added high-throughput experiments for 376 diseases and evaluated all pairwise similarities among TCM herbs/ingredients/formulae, modern drugs and diseases. (iii) we provide an automatic analyzing interface for users to upload their gene expression profiles and map them to our curated datasets. (iv) we built knowledge graph representations of HERB entities and relationships to retrieve TCM knowledge better. In summary, HERB 2.0 represents rich data type, content, utilization, and visualization improvements to support TCM research and guide modern drug discovery. It is accessible through http://herb.ac.cn/v2 or http://47.92.70.12.

PerturBase: a comprehensive database for single-cell perturbation data analysis and visualization

Single-cell perturbation (scPerturbation) sequencing techniques, represented by single-cell genetic perturbation (e.g. Perturb-seq) and single-cell chemical perturbation (e.g. sci-Plex), result from the integration of single-cell toolkits with conventional bulk screening methods. These innovative sequencing techniques empower researchers to dissect perturbation effects in biological systems at an unprecedented resolution. Despite these advancements, a notable gap exists in the availability of a dedicated database for exploring scPerturbation data. To address this gap, we present PerturBase, the most comprehensive database designed for the analysis and visualization of scPerturbation data (http://www.perturbase.cn/). PerturBase curates 122 datasets from 46 publicly available studies, covering 115 single-modal and 7 multi-modal datasets that include 24 254 genetic and 230 chemical perturbations from approximately 5 million cells. The database, comprising the 'Dataset' and 'Perturbation' modules, provides insights into various results, encompassing quality control, denoising, differential gene expression analysis, functional analysis of perturbation effects and characterization of relationships between perturbations. All the datasets and results are presented on user-friendly, easy-to-browse web pages and can be visualized through intuitive and interactive plot and table formats. In summary, PerturBase stands as a pioneering, high-content database intended for searching, visualizing and analyzing scPerturbation datasets, contributing to a deeper understanding of perturbation effects.

Revolutionising Neurological Therapeutics: Investigating Drug Repurposing Strategies

Repurposing drugs (DR) has become a viable approach to hasten the search for cures for neurodegenerative diseases (NDs). This review examines different off-target and on-target drug discovery techniques and how they might be used to find possible treatments for non-diagnostic depressions. Off-target strategies look at the known or unknown side effects of currently approved drugs for repositioning, whereas on-target strategies connect disease pathways to targets that can be treated with drugs. The review highlights the potential of experimental and computational methodologies, such as machine learning, proteomic techniques, network and genomics-based approaches, and in silico screening, in uncovering new drug-disease correlations. It also looks at difficulties and failed attempts at drug repurposing for NDs, highlighting the necessity of exact and standardised procedures to increase success rates. This review's objectives are to address the purpose of drug repurposing in human disorders, particularly neurological diseases, and to provide an overview of repurposing candidates that are presently undergoing clinical trials for neurological conditions, along with any possible causes and early findings. We then include a list of drug repurposing strategies, restrictions, and difficulties for upcoming research.

Repurposing an epithelial sodium channel inhibitor as a therapy for murine and human skin inflammation

Inflammatory skin disease is characterized by a pathologic interplay between skin cells and immunocytes and can result in disfiguring cutaneous lesions and systemic inflammation. Immunosuppression is commonly used to target the inflammatory component; however, these drugs are often expensive and associated with side effects. To identify previously unidentified targets, we carried out a nonbiased informatics screen to identify drug compounds with an inverse transcriptional signature to keratinocyte inflammatory signals. Using psoriasis, a prototypic inflammatory skin disease, as a model, we used pharmacologic, transcriptomic, and proteomic characterization to find that benzamil, the benzyl derivative of the US Food and Drug Administration-approved diuretic amiloride, effectively reversed keratinocyte-driven inflammatory signaling. Through three independent mouse models of skin inflammation (Rac1G12V transgenic mice, topical imiquimod, and human skin xenografts from patients with psoriasis), we found that benzamil disrupted pathogenic interactions between the small GTPase Rac1 and its adaptor NCK1. This reduced STAT3 and NF-κB signaling and downstream cytokine production in keratinocytes. Genetic knockdown of sodium channels or pharmacological inhibition by benzamil prevented excess Rac1-NCK1 binding and limited proinflammatory signaling pathway activation in patient-derived keratinocytes without systemic immunosuppression. Both systemic and topical applications of benzamil were efficacious, suggesting that it may be a potential therapeutic avenue for treating skin inflammation.

Unlocking biological insights from differentially expressed genes: Concepts, methods, and future perspectives

BACKGROUND

Identifying differentially expressed genes (DEGs) is a core task of transcriptome analysis, as DEGs can reveal the molecular mechanisms underlying biological processes. However, interpreting the biological significance of large DEG lists is challenging. Currently, gene ontology, pathway enrichment and protein-protein interaction analysis are common strategies employed by biologists. Additionally, emerging analytical strategies/approaches (such as network module analysis, knowledge graph, drug repurposing, cell marker discovery, trajectory analysis, and cell communication analysis) have been proposed. Despite these advances, comprehensive guidelines for systematically and thoroughly mining the biological information within DEGs remain lacking.

AIM OF REVIEW

This review aims to provide an overview of essential concepts and methodologies for the biological interpretation of DEGs, enhancing the contextual understanding. It also addresses the current limitations and future perspectives of these approaches, highlighting their broad applications in deciphering the molecular mechanism of complex diseases and phenotypes. To assist users in extracting insights from extensive datasets, especially various DEG lists, we developed DEGMiner (https://www.ciblab.net/DEGMiner/), which integrates over 300 easily accessible databases and tools.

KEY SCIENTIFIC CONCEPTS OF REVIEW

This review offers strong support and guidance for exploring DEGs, and also will accelerate the discovery of hidden biological insights within genomes.

OrthologAL: A Shiny application for quality-aware humanization of non-human pre-clinical high-dimensional gene expression data

Single-cell and spatial transcriptomics provide unprecedented insight into the inner workings of disease. Pharmacotranscriptomic approaches are powerful tools that leverage gene expression data for drug repurposing and treatment discovery in many diseases. Multiple databases attempt to connect human cellular transcriptional responses to small molecules for use in transcriptome-based drug discovery efforts. However, pre-clinical research often requires in vivo experiments in non-human species, which makes capitalizing on such valuable resources difficult. To facilitate the application of pharmacotranscriptomic databases to pre-clinical research models and to facilitate human orthologous conversion of non-human transcriptomes, we introduce OrthologAL. OrthologAL leverages the BioMart database to access different gene sets from Ensembl, facilitating the interaction between these servers without needing user-generated code. Researchers can input their single-cell or other high-dimensional gene expression data from any species, and OrthologAL will output a human ortholog-converted dataset for download and use. To demonstrate the utility of this application, we characterized orthologous conversion in single-cell, single-nuclei, and spatial transcriptomic data derived from common pre-clinical models, including patient-derived orthotopic xenografts of medulloblastoma, and mouse and rat models of spinal cord injury. We show that OrthologAL can convert these data types efficiently to that of corresponding orthologs while preserving the dimensional architecture of the original non-human expression data. OrthologAL will be broadly useful for applying pre-clinical, high-dimensional transcriptomics data in functional small molecule predictions using existing human-annotated databases.

CDCM: a correlation-dependent connectivity map approach to rapidly screen drugs during outbreaks of infectious diseases

In the context of the global damage caused by coronavirus disease 2019 (COVID-19) and the emergence of the monkeypox virus (MPXV) outbreak as a public health emergency of international concern, research into methods that can rapidly test potential therapeutics during an outbreak of a new infectious disease is urgently needed. Computational drug discovery is an effective way to solve such problems. The existence of various large open databases has mitigated the time and resource consumption of traditional drug development and improved the speed of drug discovery. However, the diversity of cell lines used in various databases remains limited, and previous drug discovery methods are ineffective for cross-cell prediction. In this study, we propose a correlation-dependent connectivity map (CDCM) to achieve cross-cell predictions of drug similarity. The CDCM mainly identifies drug-drug or disease-drug relationships from the perspective of gene networks by exploring the correlation changes between genes and identifying similarities in the effects of drugs or diseases on gene expression. We validated the CDCM on multiple datasets and found that it performed well for drug identification across cell lines. A comparison with the Connectivity Map revealed that our method was more stable and performed better across different cell lines. In the application of the CDCM to COVID-19 and MPXV data, the predictions of potential therapeutic compounds for COVID-19 were consistent with several previous studies, and most of the predicted drugs were found to be experimentally effective against MPXV. This result confirms the practical value of the CDCM. With the ability to predict across cell lines, the CDCM outperforms the Connectivity Map, and it has wider application prospects and a reduced cost of use.

Identification of Key Immune and Cell Cycle Modules and Prognostic Genes for Glioma Patients through Transcriptome Analysis

BACKGROUND

Gliomas, the most prevalent type of primary brain tumor, stand out as one of the most aggressive and lethal types of human cancer.

METHODS & RESULTS

To uncover potential prognostic markers, we employed the weighted correlation network analysis (WGCNA) on the Chinese Glioma Genome Atlas (CGGA) 693 dataset to reveal four modules significantly associated with glioma clinical traits, primarily involved in immune function, cell cycle regulation, and ribosome biogenesis. Using the least absolute shrinkage and selection operator (LASSO) regression algorithm, we identified 11 key genes and developed a prognostic risk score model, which exhibits precise prognostic prediction in the CGGA 325 dataset. More importantly, we also validated the model in 12 glioma patients with overall survival (OS) ranging from 4 to 132 months using mRNA sequencing and immunohistochemical analysis. The analysis of immune infiltration revealed that patients with high-risk scores exhibit a heightened immune infiltration, particularly immune suppression cells, along with increased expression of immune checkpoints. Furthermore, we explored potentially effective drugs targeting 11 key genes for gliomas using the library of integrated network-based cellular signatures (LINCS) L1000 database, identifying that in vitro, both torin-1 and clofarabine exhibit promising anti-glioma activity and inhibitory effect on the cell cycle, a significant pathway enriched in the identified glioma modules.

CONCLUSIONS

In conclusion, our study provides valuable insights into molecular mechanisms and identifying potential therapeutic targets for gliomas.

Reversal Gene Expression Assessment for Drug Repurposing, a Case Study of Glioblastoma

Glioblastoma (GBM) is a rare brain cancer with an exceptionally high mortality rate, which illustrates the pressing demand for more effective therapeutic options. Despite considerable research efforts on GBM, its underlying biological mechanisms remain unclear. Furthermore, none of the United States Food and Drug Administration (FDA) approved drugs used for GBM deliver satisfactory survival improvement. This study presents a novel computational pipeline by utilizing gene expression data analysis for GBM for drug repurposing to address the challenges in rare disease drug development, particularly focusing on GBM. The GBM Gene Expression Profile (GGEP) was constructed with multi-omics data to identify drugs with reversal gene expression to GGEP from the Integrated Network-Based Cellular Signatures (iLINCS) database. We prioritized the candidates via hierarchical clustering of their expression signatures and quantification of their reversal strength by calculating two self-defined indices based on the GGEP genes' log2 foldchange (LFCs) that the drug candidates could induce. Among eight prioritized candidates, in-vitro experiments validated Clofarabine and Ciclopirox as highly efficacious in selectively targeting GBM cancer cells. The success of this study illustrated a promising avenue for accelerating drug development by uncovering underlying gene expression effect between drugs and diseases, which can be extended to other rare diseases and non-rare diseases.

Know your molecule: pharmacological characterization of drug candidates to enhance efficacy and reduce late-stage attrition

Despite advances in chemical, computational and biological sciences, the rate of attrition of drug candidates in clinical development is still high. A key point in the small-molecule discovery process that could provide opportunities to help address this challenge is the pharmacological characterization of hit and lead compounds, culminating in the selection of a drug candidate. Deeper characterization is increasingly important, because the 'quality' of drug efficacy, at least for G protein-coupled receptors (GPCRs), is now understood to be much more than activation of commonly evaluated pathways such as cAMP signalling, with many more 'efficacies' of ligands that could be harnessed therapeutically. Such characterization is being enabled by novel assays to characterize the complex behaviour of GPCRs, such as biased signalling and allosteric modulation, as well as advances in structural biology, such as cryo-electron microscopy. This article discusses key factors in the assessments of the pharmacology of hit and lead compounds in the context of GPCRs as a target class, highlighting opportunities to identify drug candidates with the potential to address limitations of current therapies and to improve the probability of them succeeding in clinical development.

Computer-aided drug discovery strategies for novel therapeutics for prostate cancer leveraging next-generating sequencing data

INTRODUCTION

Prostate cancer (PC) is the most common malignancy and accounts for a significant proportion of cancer deaths among men. Although initial therapy success can often be observed in patients diagnosed with localized PC, many patients eventually develop disease recurrence and metastasis. Without effective treatments, patients with aggressive PC display very poor survival. To curb the current high mortality rate, many investigations have been carried out to identify efficacious therapeutics. Compared to de novo drug designs, computational methods have been widely employed to offer actionable drug predictions in a fast and cost-efficient way. Particularly, powered by an increasing availability of next-generation sequencing molecular profiles from PC patients, computer-aided approaches can be tailored to screen for candidate drugs.

AREAS COVERED

Herein, the authors review the recent advances in computational methods for drug discovery utilizing molecular profiles from PC patients. Given the uniqueness in PC therapeutic needs, they discuss in detail the drug discovery goals of these studies, highlighting their translational values for clinically impactful drug nomination.

EXPERT OPINION

Evolving molecular profiling techniques may enable new perspectives for computer-aided approaches to offer drug candidates for different tumor microenvironments. With ongoing efforts to incorporate new compounds into large-scale high-throughput screens, the authors envision continued expansion of drug candidate pools.

Deep representation learning of chemical-induced transcriptional profile for phenotype-based drug discovery

Artificial intelligence transforms drug discovery, with phenotype-based approaches emerging as a promising alternative to target-based methods, overcoming limitations like lack of well-defined targets. While chemical-induced transcriptional profiles offer a comprehensive view of drug mechanisms, inherent noise often obscures the true signal, hindering their potential for meaningful insights. Here, we highlight the development of TranSiGen, a deep generative model employing self-supervised representation learning. TranSiGen analyzes basal cell gene expression and molecular structures to reconstruct chemical-induced transcriptional profiles with high accuracy. By capturing both cellular and compound information, TranSiGen-derived representations demonstrate efficacy in diverse downstream tasks like ligand-based virtual screening, drug response prediction, and phenotype-based drug repurposing. Notably, in vitro validation of TranSiGen's application in pancreatic cancer drug discovery highlights its potential for identifying effective compounds. We envisage that integrating TranSiGen into the drug discovery and mechanism research holds significant promise for advancing biomedicine.

Intravesical instillation-based mTOR-STAT3 dual targeting for bladder cancer treatment

BACKGROUND

Recent intravesical administration of adenoviral vectors, either as a single injection or in combination with immune checkpoint inhibitors, exemplified by cretostimogene grenadenorepvec and nadofaragene firadenovec, has demonstrated remarkable efficacy in clinical trials for non-muscle invasive bladder cancer. Despite their ability to induce an enhanced immune reaction within the lesion, the intracellular survival signaling of cancer cells has not been thoroughly addressed.

METHODS

An analysis of the prognostic data revealed a high probability of therapeutic efficacy with simultaneous inhibition of mTOR and STAT3. Considering the challenges of limited pharmaco-accessibility to the bladder due to its pathophysiological structure and the partially undruggable nature of target molecules, we designed a dual siRNA system targeting both mRNAs. Subsequently, this dual siRNA system was encoded into the adenovirus 5/3 (Ad 5/3) to enhance in vivo delivery efficiency.

RESULTS

Gene-targeting efficacy was assessed using cells isolated from xenografted tumors using a single-cell analysis system. Our strategy demonstrated a balanced downregulation of mTOR and STAT3 at the single-cell resolution, both in vitro and in vivo. This approach reduced tumor growth in bladder cancer xenograft and orthotopic animal experiments. In addition, increased infiltration of CD8+ T cells was observed in a humanized mouse model. We provided helpful and safe tissue distribution data for intravesical therapy of siRNAs coding adenoviruses.

CONCLUSIONS

The bi-specific siRNA strategy, encapsulated in an adenovirus, could be a promising tool to augment cancer treatment efficacy and overcome conventional therapy limitations associated with "undruggability." Hence, we propose that dual targeting of mTOR and STAT3 is an advantageous strategy for intravesical therapy using adenoviruses.

Drug repurposing for metabolic disorders: Scientific, technological and economic issues

Obesity, diabetes, and other metabolic disorders place a huge burden on both the physical health and financial well-being of the community. While the need for effective treatment of metabolic disorders remains urgent and the reality is that traditional drug development involves high costs and a very long time with many pre-clinical and clinical trials, the need for drug repurposing has emerged as a potential alternative. Scientific evidence has shown the anti-diabetic and anti-obesity effects of old drugs, which were initially utilized for the treatment of inflammation, depression, infections, and even cancers. The drug library used modern technological methods to conduct drug screening. Computational molecular docking, genome-wide association studies, or omics data mining are advantageous and unavoidable methods for drug repurposing. Drug repurposing offers a promising avenue for economic efficiency in healthcare, especially for less common metabolic diseases, despite the need for rigorous research and validation. In this chapter, we aim to explore the scientific, technological, and economic issues surrounding drug repurposing for metabolic disorders. We hope to shed light on the potential of this approach and the challenges that need to be addressed to make it a viable option in the treatment of metabolic disorders, especially in the future fight against metabolic disorders.

Signature reversion of three disease-associated gene signatures prioritizes cancer drug repurposing candidates

Drug repurposing is promising because approving a drug for a new indication requires fewer resources than approving a new drug. Signature reversion detects drug perturbations most inversely related to the disease-associated gene signature to identify drugs that may reverse that signature. We assessed the performance and biological relevance of three approaches for constructing disease-associated gene signatures (i.e., limma, DESeq2, and MultiPLIER) and prioritized the resulting drug repurposing candidates for four low-survival human cancers. Our results were enriched for candidates that had been used in clinical trials or performed well in the PRISM drug screen. Additionally, we found that pamidronate and nimodipine, drugs predicted to be efficacious against the brain tumor glioblastoma (GBM), inhibited the growth of a GBM cell line and cells isolated from a patient-derived xenograft (PDX). Our results demonstrate that by applying multiple disease-associated gene signature methods, we prioritized several drug repurposing candidates for low-survival cancers.

Antiretroviral Drug Repositioning for Glioblastoma

Outcomes for glioblastoma (GBM) remain poor despite standard-of-care treatments including surgical resection, radiation, and chemotherapy. Intratumoral heterogeneity contributes to treatment resistance and poor prognosis, thus demanding novel therapeutic approaches. Drug repositioning studies on antiretroviral therapy (ART) have shown promising potent antineoplastic effects in multiple cancers; however, its efficacy in GBM remains unclear. To better understand the pleiotropic anticancer effects of ART on GBM, we conducted a comprehensive drug repurposing analysis of ART in GBM to highlight its utility in translational neuro-oncology. To uncover the anticancer role of ART in GBM, we conducted a comprehensive bioinformatic and in vitro screen of antiretrovirals against glioblastoma. Using the DepMap repository and reversal of gene expression score, we conducted an unbiased screen of 16 antiretrovirals in 40 glioma cell lines to identify promising candidates for GBM drug repositioning. We utilized patient-derived neurospheres and glioma cell lines to assess neurosphere viability, proliferation, and stemness. Our in silico screen revealed that several ART drugs including reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) demonstrated marked anti-glioma activity with the capability of reversing the GBM disease signature. RTIs effectively decreased cell viability, GBM stem cell markers, and proliferation. Our study provides mechanistic and functional insight into the utility of ART repurposing for malignant gliomas, which supports the current literature. Given their safety profile, preclinical efficacy, and neuropenetrance, ARTs may be a promising adjuvant treatment for GBM.

Modulator of TMB-associated immune infiltration (MOTIF) predicts immunotherapy response and guides combination therapy

Patients with high tumor mutational burden (TMB) levels do not consistently respond to immune checkpoint inhibitors (ICIs), possibly because a high TMB level does not necessarily result in adequate infiltration of CD8+ T cells. Using bulk ribonucleic acid sequencing (RNA-seq) data from 9311 tumor samples across 30 cancer types, we developed a novel tool called the modulator of TMB-associated immune infiltration (MOTIF), which comprises genes that can determine the extent of CD8+ T cell infiltration prompted by a certain TMB level. We confirmed that MOTIF can accurately reflect the integrity and defects of the cancer-immunity cycle. By analyzing 84 human single-cell RNA-seq datasets from 32 types of solid tumors, we revealed that MOTIF can provide insights into the diverse roles of various cell types in the modulation of CD8+ T cell infiltration. Using pretreatment RNA-seq data from 13 ICI-treated cohorts, we validated the use of MOTIF in predicting CD8+ T cell infiltration and ICI efficacy. Among the components of MOTIF, we identified EMC3 as a negative regulator of CD8+ T cell infiltration, which was validated via in vivo studies. Additionally, MOTIF provided guidance for the potential combinations of programmed death 1 blockade with certain immunostimulatory drugs to facilitate CD8+ T cell infiltration and improve ICI efficacy.

Mapping drug biology to disease genetics to discover drug impacts on the human phenome

Motivation

Medications can have unexpected effects on disease, including not only harmful drug side effects, but also beneficial drug repurposing. These effects on disease may result from hidden influences of drugs on disease gene networks. Then, discovering how biological effects of drugs relate to disease biology can both provide insight into the mechanism of latent drug effects, and can help predict new effects.

Results

Here, we develop Draphnet, a model that integrates molecular data on 429 drugs and gene associations of nearly 200 common phenotypes to learn a network that explains drug effects on disease in terms of these molecular signals. We present evidence that our method can both predict drug effects, and can provide insight into the biology of unexpected drug effects on disease. Using Draphnet to map a drug's known molecular effects to downstream effects on the disease genome, we put forward disease genes impacted by drugs, and we suggest a new grouping of drugs based on shared effects on the disease genome. Our approach has multiple applications, including predicting drug uses and learning drug biology, with implications for personalized medicine.

Availability and implementation

Code to reproduce the analysis is available at https://github.com/RDMelamed/drug-phenome.

Connectivity mapping-based identification of pharmacological inhibitor targeting HDAC6 in aggressive pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) remains highly lethal due to limited therapeutic options and expensive/burdensome drug discovery processes. Utilizing genomic-data-driven Connectivity Mapping (CMAP) to identify a drug closer to real-world PC targeting may improve pancreatic cancer (PC) patient outcomes. Initially, we mapped CMAP data to gene expression from 106 PC patients, identifying nine negatively connected drugs. These drugs were further narrowed down using a similar analysis for PC cell lines, human tumoroids, and patient-derived xenografts datasets, where ISOX emerged as the most potent agent to target PC. We used human and mouse syngeneic PC cells, human and mouse tumoroids, and in vivo mice to assess the ability of ISOX alone and in combination with 5FU to inhibit tumor growth. Global transcriptomic and pathway analysis of the ISOX-LINCS signature identified HDAC 6/cMyc as the target axis for ISOX. Specifically, we discovered that genetic and pharmacological targeting of HDAC 6 affected non-histone protein cMyc acetylation, leading to cMyc instability, thereby disrupting PC growth and metastasis by affecting cancer stemness. Finally, KrasG12D harboring tumoroids and mice responded effectively against ISOX and 5FU treatment by enhancing survival and controlling metastasis incidence. Overall, our data validate ISOX as a new drug to treat advanced PC patients without toxicity to normal cells. Our study supports the clinical utility of ISOX along with 5FU in future PC clinical trials.

Deciphering Drug Targets and Actions with Single-Cell and Spatial Resolution

Recent advances in chemical, molecular, and genetic approaches have provided us with an unprecedented capacity to identify drug-target interactions across the whole proteome and genome. Meanwhile, rapid developments of single-cell and spatial omics technologies are revolutionizing our understanding of the molecular architecture of biological systems. However, a significant gap remains in how we align our understanding of drug actions, traditionally based on molecular affinities, with the in vivo cellular and spatial tissue heterogeneity revealed by these newer techniques. Here, we review state-of-the-art methods for profiling drug-target interactions and emerging multiomics tools to delineate the tissue heterogeneity at single-cell resolution. Highlighting the recent technical advances enabling high-resolution, multiplexable in situ small-molecule drug imaging (clearing-assisted tissue click chemistry, or CATCH), we foresee the integration of single-cell and spatial omics platforms, data, and concepts into the future framework of defining and understanding in vivo drug-target interactions and mechanisms of actions.

CMAUP database update 2024: extended functional and association information of useful plants for biomedical research

Knowledge of the collective activities of individual plants together with the derived clinical effects and targeted disease associations is useful for plant-based biomedical research. To provide the information in complement to the established databases, we introduced a major update of CMAUP database, previously featured in NAR. This update includes (i) human transcriptomic changes overlapping with 1152 targets of 5765 individual plants, covering 74 diseases from 20 027 patient samples; (ii) clinical information for 185 individual plants in 691 clinical trials; (iii) drug development information for 4694 drug-producing plants with metabolites developed into approved or clinical trial drugs; (iv) plant and human disease associations (428 737 associations by target, 220 935 reversion of transcriptomic changes, 764 and 154121 associations by clinical trials of individual plants and plant ingredients); (v) the location of individual plants in the phylogenetic tree for navigating taxonomic neighbors, (vi) DNA barcodes of 3949 plants, (vii) predicted human oral bioavailability of plant ingredients by the established SwissADME and HobPre algorithm, (viii) 21-107% increase of CMAUP data over the previous version to cover 60 222 chemical ingredients, 7865 plants, 758 targets, 1399 diseases, 238 KEGG human pathways, 3013 gene ontologies and 1203 disease ontologies. CMAUP update version is freely accessible at https://bidd.group/CMAUP/index.html.

Identification of hub genes to determine drug-disease correlation in breast carcinomas

The exact molecular mechanism underlying the heterogeneous drug response against breast carcinoma remains to be fully understood. It is urgently required to identify key genes that are intricately associated with varied clinical response of standard anti-cancer drugs, clinically used to treat breast cancer patients. In the present study, the utility of transcriptomic data of breast cancer patients in discerning the clinical drug response using machine learning-based approaches were evaluated. Here, a computational framework has been developed which can be used to identify key genes that can be linked with clinical drug response and progression of cancer, offering an immense opportunity to predict potential prognostic biomarkers and therapeutic targets. The framework concerned utilizes DeSeq2, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Cytoscape, and machine learning techniques to find these crucial genes. Total RNA extraction and qRT-PCR were performed to quantify relative expression of few hub genes selected from the networks. In our study, we have experimentally checked the expression of few key hub genes like APOA2, DLX5, APOC3, CAMK2B, and PAK6 that were predicted to play an immense role in breast cancer tumorigenesis and progression in response to anti-cancer drug Paclitaxel. However, further experimental validations will be required to get mechanistic insights of these genes in regulating the drug response and cancer progression which will likely to play pivotal role in cancer treatment and precision oncology.

Comprehensive Profiling and Therapeutic Insights into Differentially Expressed Genes in Hepatocellular Carcinoma

Background: Drug repurposing is a strategy that complements the conventional approach of developing new drugs. Hepatocellular carcinoma (HCC) is a highly prevalent type of liver cancer, necessitating an in-depth understanding of the underlying molecular alterations for improved treatment. Methods: We searched for a vast array of microarray experiments in addition to RNA-seq data. Through rigorous filtering processes, we have identified highly representative differentially expressed genes (DEGs) between tumor and non-tumor liver tissues and identified a distinct class of possible new candidate drugs. Results: Functional enrichment analysis revealed distinct biological processes associated with metal ions, including zinc, cadmium, and copper, potentially implicating chronic metal ion exposure in tumorigenesis. Conversely, up-regulated genes are associated with mitotic events and kinase activities, aligning with the relevance of kinases in HCC. To unravel the regulatory networks governing these DEGs, we employed topological analysis methods, identifying 25 hub genes and their regulatory transcription factors. In the pursuit of potential therapeutic options, we explored drug repurposing strategies based on computational approaches, analyzing their potential to reverse the expression patterns of key genes, including AURKA, CCNB1, CDK1, RRM2, and TOP2A. Potential therapeutic chemicals are alvocidib, AT-7519, kenpaullone, PHA-793887, JNJ-7706621, danusertibe, doxorubicin and analogues, mitoxantrone, podofilox, teniposide, and amonafide. Conclusion: This multi-omic study offers a comprehensive view of DEGs in HCC, shedding light on potential therapeutic targets and drug repurposing opportunities.

Integrating population-level and cell-based signatures for drug repositioning

Drug repositioning presents a streamlined and cost-efficient way to expand the range of therapeutic possibilities. Furthermore, drugs with genetic evidence are more likely to progress successfully through clinical trials towards FDA approval. Exploiting these developments, single gene-based drug repositioning methods have been implemented, but approaches leveraging the entire spectrum of molecular signatures are critically underexplored. Most multi-gene-based approaches rely on differential gene expression (DGE) analysis, which is prone to identify the molecular consequence of disease and renders causal inference challenging. We propose a framework TReD (Transcriptome-informed Reversal Distance) that integrates population-level disease signatures robust to reverse causality and cell-based drug-induced transcriptome response profiles. TReD embeds the disease signature and drug profile in a high-dimensional normed space, quantifying the reversal potential of candidate drugs in a disease-related cell screen assay. The robustness is ensured by evaluation in additional cell screens. For an application, we implement the framework to identify potential drugs against COVID-19. Taking transcriptome-wide association study (TWAS) results from four relevant tissues and three DGE results as disease features, we identify 37 drugs showing potential reversal roles in at least four of the seven disease signatures. Notably, over 70% (27/37) of the drugs have been linked to COVID-19 from other studies, and among them, eight drugs are supported by ongoing/completed clinical trials. For example, TReD identifies the well-studied JAK1/JAK2 inhibitor baricitinib, the first FDA-approved immunomodulatory treatment for COVID-19. Novel potential candidates, including enzastaurin, a selective inhibitor of PKC-beta which can be activated by SARS-CoV-2, are also identified. In summary, we propose a comprehensive genetics-anchored framework integrating population-level signatures and cell-based screens that can accelerate the search for new therapeutic strategies.

Applying a Gene Reversal Rate Computational Methodology to Identify Drugs for a Rare Cancer: Inflammatory Breast Cancer

The aim of this study was to utilize a computational methodology based on Gene Reversal Rate (GRR) scoring to repurpose existing drugs for a rare and understudied cancer: inflammatory breast cancer (IBC). This method uses IBC-related gene expression signatures (GES) and drug-induced gene expression profiles from the LINCS database to calculate a GRR score for each candidate drug, and is based on the idea that a compound that can counteract gene expression changes of a disease may have potential therapeutic applications for that disease. Genes related to IBC with associated differential expression data (265 up-regulated and 122 down-regulated) were collated from PubMed-indexed publications. Drug-induced gene expression profiles were downloaded from the LINCS database and candidate drugs to treat IBC were predicted using their GRR scores. Thirty-two (32) drug perturbations that could potentially reverse the pre-compiled list of 297 IBC genes were obtained using the LINCS Canvas Browser (LCB) analysis. Binary combinations of the 32 perturbations were assessed computationally to identify combined perturbations with the highest GRR scores, and resulted in 131 combinations with GRR greater than 80%, that reverse up to 264 of the 297 genes in the IBC-GES. The top 35 combinations involve 20 unique individual drug perturbations, and 19 potential drug candidates. A comprehensive literature search confirmed 17 of the 19 known drugs as having either anti-cancer or anti-inflammatory activities. AZD-7545, BMS-754807, and nimesulide target known IBC relevant genes: PDK, Met, and COX, respectively. AG-14361, butalbital, and clobenpropit are known to be functionally relevant in DNA damage, cell cycle, and apoptosis, respectively. These findings support the use of the GRR approach to identify drug candidates and potential combination therapies that could be used to treat rare diseases such as IBC.

COIMMR: a computational framework to reveal the contribution of herbal ingredients against human cancer via immune microenvironment and metabolic reprogramming

Immune evasion and metabolism reprogramming have been regarded as two vital hallmarks of the mechanism of carcinogenesis. Thus, targeting the immune microenvironment and the reprogrammed metabolic processes will aid in developing novel anti-cancer drugs. In recent decades, herbal medicine has been widely utilized to treat cancer through the modulation of the immune microenvironment and reprogrammed metabolic processes. However, labor-based herbal ingredient screening is time consuming, laborious and costly. Luckily, some computational approaches have been proposed to screen candidates for drug discovery rapidly. Yet, it has been challenging to develop methods to screen drug candidates exclusively targeting specific pathways, especially for herbal ingredients which exert anti-cancer effects by multiple targets, multiple pathways and synergistic ways. Meanwhile, currently employed approaches cannot quantify the contribution of the specific pathway to the overall curative effect of herbal ingredients. Hence, to address this problem, this study proposes a new computational framework to infer the contribution of the immune microenvironment and metabolic reprogramming (COIMMR) in herbal ingredients against human cancer and specifically screen herbal ingredients targeting the immune microenvironment and metabolic reprogramming. Finally, COIMMR was applied to identify isoliquiritigenin that specifically regulates the T cells in stomach adenocarcinoma and cephaelin hydrochloride that specifically targets metabolic reprogramming in low-grade glioma. The in silico results were further verified using in vitro experiments. Taken together, our approach opens new possibilities for repositioning drugs targeting immune and metabolic dysfunction in human cancer and provides new insights for drug development in other diseases. COIMMR is available at https://github.com/LYN2323/COIMMR.

Germline mechanisms of immunotherapy toxicities in the era of genome-wide association studies

Cancer immunotherapy has revolutionized the treatment of advanced cancers and is quickly becoming an option for early-stage disease. By reactivating the host immune system, immunotherapy harnesses patients' innate defenses to eradicate the tumor. By putatively similar mechanisms, immunotherapy can also substantially increase the risk of toxicities or immune-related adverse events (irAEs). Severe irAEs can lead to hospitalization, treatment discontinuation, lifelong immune complications, or even death. Many irAEs present with similar symptoms to heritable autoimmune diseases, suggesting that germline genetics may contribute to their onset. Recently, genome-wide association studies (GWAS) of irAEs have identified common germline associations and putative mechanisms, lending support to this hypothesis. A wide range of well-established GWAS methods can potentially be harnessed to understand the etiology of irAEs specifically and immunotherapy outcomes broadly. This review summarizes current findings regarding germline effects on immunotherapy outcomes and discusses opportunities and challenges for leveraging germline genetics to understand, predict, and treat irAEs.

Endometriosis in the era of precision medicine and impact on sexual and reproductive health across the lifespan and in diverse populations

Endometriosis is a common estrogen-dependent disorder wherein uterine lining tissue (endometrium) is found mainly in the pelvis where it causes inflammation, chronic pelvic pain, pain with intercourse and menses, and infertility. Recent evidence also supports a systemic inflammatory component that underlies associated co-morbidities, e.g., migraines and cardiovascular and autoimmune diseases. Genetics and environment contribute significantly to disease risk, and with the explosion of omics technologies, underlying mechanisms of symptoms are increasingly being elucidated, although novel and effective therapeutics for pain and infertility have lagged behind these advances. Moreover, there are stark disparities in diagnosis, access to care, and treatment among persons of color and transgender/nonbinary identity, socioeconomically disadvantaged populations, and adolescents, and a disturbing low awareness among health care providers, policymakers, and the lay public about endometriosis, which, if left undiagnosed and under-treated can lead to significant fibrosis, infertility, depression, and markedly diminished quality of life. This review summarizes endometriosis epidemiology, compelling evidence for its pathogenesis, mechanisms underlying its pathophysiology in the age of precision medicine, recent biomarker discovery, novel therapeutic approaches, and issues around reproductive justice for marginalized populations with this disorder spanning the past 100 years. As we enter the next revolution in health care and biomedical research, with rich molecular and clinical datasets, single-cell omics, and population-level data, endometriosis is well positioned to benefit from data-driven research leveraging computational and artificial intelligence approaches integrating data and predicting disease risk, diagnosis, response to medical and surgical therapies, and prognosis for recurrence.

Transcriptomic clock predicts vascular changes of prodromal diabetic retinopathy

Diabetic retinopathy is a common complication of long-term diabetes and that could lead to vision loss. Unfortunately, early diabetic retinopathy remains poorly understood. There is no effective way to prevent or treat early diabetic retinopathy until patients develop later stages of diabetic retinopathy. Elevated acellular capillary density is considered a reliable quantitative trait present in the early development of retinopathy. Hence, in this study, we interrogated whole retinal vascular transcriptomic changes via a Nile rat model to better understand the early pathogenesis of diabetic retinopathy. We uncovered the complexity of associations between acellular capillary density and the joint factors of blood glucose, diet, and sex, which was modeled through a Bayesian network. Using segmented regressions, we have identified different gene expression patterns and enriched Gene Ontology (GO) terms associated with acellular capillary density increasing. We developed a random forest regression model based on expression patterns of 14 genes to predict the acellular capillary density. Since acellular capillary density is a reliable quantitative trait in early diabetic retinopathy, and thus our model can be used as a transcriptomic clock to measure the severity of the progression of early retinopathy. We also identified NVP-TAE684, geldanamycin, and NVP-AUY922 as the top three potential drugs which can potentially attenuate the early DR. Although we need more in vivo studies in the future to support our re-purposed drugs, we have provided a data-driven approach to drug discovery.

Network neighborhood operates as a drug repositioning method for cancer treatment

Computational drug repositioning approaches are important, as they cost less compared to the traditional drug development processes. This study proposes a novel network-based drug repositioning approach, which computes similarities between disease-causing genes and drug-affected genes in a network topology to suggest candidate drugs with highest similarity scores. This new method aims to identify better treatment options by integrating systems biology approaches. It uses a protein-protein interaction network that is the main topology to compute a similarity score between candidate drugs and disease-causing genes. The disease-causing genes were mapped on this network structure. Transcriptome profiles of drug candidates were taken from the LINCS project and mapped individually on the network structure. The similarity of these two networks was calculated by different network neighborhood metrics, including Adamic-Adar, PageRank and neighborhood scoring. The proposed approach identifies the best candidates by choosing the drugs with significant similarity scores. The method was experimented on melanoma, colorectal, and prostate cancers. Several candidate drugs were predicted by applying AUC values of 0.6 or higher. Some of the predictions were approved by clinical phase trials or other in-vivo studies found in literature. The proposed drug repositioning approach would suggest better treatment options with integration of functional information between genes and transcriptome level effects of drug perturbations and diseases.

PATH-SURVEYOR: pathway level survival enquiry for immuno-oncology and drug repurposing

Pathway-level survival analysis offers the opportunity to examine molecular pathways and immune signatures that influence patient outcomes. However, available survival analysis algorithms are limited in pathway-level function and lack a streamlined analytical process. Here we present a comprehensive pathway-level survival analysis suite, PATH-SURVEYOR, which includes a Shiny user interface with extensive features for systematic exploration of pathways and covariates in a Cox proportional-hazard model. Moreover, our framework offers an integrative strategy for performing Hazard Ratio ranked Gene Set Enrichment Analysis and pathway clustering. As an example, we applied our tool in a combined cohort of melanoma patients treated with checkpoint inhibition (ICI) and identified several immune populations and biomarkers predictive of ICI efficacy. We also analyzed gene expression data of pediatric acute myeloid leukemia (AML) and performed an inverse association of drug targets with the patient's clinical endpoint. Our analysis derived several drug targets in high-risk KMT2A-fusion-positive patients, which were then validated in AML cell lines in the Genomics of Drug Sensitivity database. Altogether, the tool offers a comprehensive suite for pathway-level survival analysis and a user interface for exploring drug targets, molecular features, and immune populations at different resolutions.

Transcriptomics-driven drug repositioning for the treatment of diabetic foot ulcer

Diabetic foot ulcers (DFUs) are a common complication of diabetes and can lead to severe disability and even amputation. Despite advances in treatment, there is currently no cure for DFUs and available drugs for treatment are limited. This study aimed to identify new candidate drugs and repurpose existing drugs to treat DFUs based on transcriptomics analysis. A total of 31 differentially expressed genes (DEGs) were identified and used to prioritize the biological risk genes for DFUs. Further investigation using the database DGIdb revealed 12 druggable target genes among 50 biological DFU risk genes, corresponding to 31 drugs. Interestingly, we highlighted that two drugs (urokinase and lidocaine) are under clinical investigation for DFU and 29 drugs are potential candidates to be repurposed for DFU therapy. The top 5 potential biomarkers for DFU from our findings are IL6ST, CXCL9, IL1R1, CXCR2, and IL10. This study highlights IL1R1 as a highly promising biomarker for DFU due to its high systemic score in functional annotations, that can be targeted with an existing drug, Anakinra. Our study proposed that the integration of transcriptomic and bioinformatic-based approaches has the potential to drive drug repurposing for DFUs. Further research will further examine the mechanisms by which targeting IL1R1 can be used to treat DFU.

Computational drug repositioning for the identification of new agents to sensitize drug-resistant breast tumors across treatments and receptor subtypes

Introduction

Drug resistance is a major obstacle in cancer treatment and can involve a variety of different factors. Identifying effective therapies for drug resistant tumors is integral for improving patient outcomes.

Methods

In this study, we applied a computational drug repositioning approach to identify potential agents to sensitize primary drug resistant breast cancers. We extracted drug resistance profiles from the I-SPY 2 TRIAL, a neoadjuvant trial for early stage breast cancer, by comparing gene expression profiles of responder and non-responder patients stratified into treatments within HR/HER2 receptor subtypes, yielding 17 treatment-subtype pairs. We then used a rank-based pattern-matching strategy to identify compounds in the Connectivity Map, a database of cell line derived drug perturbation profiles, that can reverse these signatures in a breast cancer cell line. We hypothesize that reversing these drug resistance signatures will sensitize tumors to treatment and prolong survival.

Results

We found that few individual genes are shared among the drug resistance profiles of different agents. At the pathway level, however, we found enrichment of immune pathways in the responders in 8 treatments within the HR+HER2+, HR+HER2-, and HR-HER2- receptor subtypes. We also found enrichment of estrogen response pathways in the non-responders in 10 treatments primarily within the hormone receptor positive subtypes. Although most of our drug predictions are unique to treatment arms and receptor subtypes, our drug repositioning pipeline identified the estrogen receptor antagonist fulvestrant as a compound that can potentially reverse resistance across 13/17 of the treatments and receptor subtypes including HR+ and triple negative. While fulvestrant showed limited efficacy when tested in a panel of 5 paclitaxel resistant breast cancer cell lines, it did increase drug response in combination with paclitaxel in HCC-1937, a triple negative breast cancer cell line.

Conclusion

We applied a computational drug repurposing approach to identify potential agents to sensitize drug resistant breast cancers in the I-SPY 2 TRIAL. We identified fulvestrant as a potential drug hit and showed that it increased response in a paclitaxel-resistant triple negative breast cancer cell line, HCC-1937, when treated in combination with paclitaxel.

Towards a Synthesis of the Non-Genetic and Genetic Views of Cancer in Understanding Pancreatic Ductal Adenocarcinoma Initiation and Prevention

While much of the research in oncogenesis and cancer therapy has focused on mutations in key cancer driver genes, more recent work suggests a complementary non-genetic paradigm. This paradigm focuses on how transcriptional and phenotypic heterogeneity, even in clonally derived cells, can create sub-populations associated with oncogenesis, metastasis, and therapy resistance. We discuss this complementary paradigm in the context of pancreatic ductal adenocarcinoma. A better understanding of cellular transcriptional heterogeneity and its association with oncogenesis can lead to more effective therapies that prevent tumor initiation and slow progression.

The Role of WNT Pathway Mutations in Cancer Development and an Overview of Therapeutic Options

It is well established that mutations in the canonical WNT-signalling pathway play a major role in various cancers. Critical to developing new therapeutic strategies is understanding which cancers are driven by WNT pathway activation and at what level these mutations occur within the pathway. Some cancers harbour mutations in genes whose protein products operate at the receptor level of the WNT pathway. For instance, tumours with RNF43 or RSPO mutations, still require exogenous WNT ligands to drive WNT signalling (ligand-dependent mutations). Conversely, mutations within the cytoplasmic segment of the Wnt pathway, such as in APC and CTNNB1, lead to constitutive WNT pathway activation even in the absence of WNT ligands (ligand-independent). Here, we review the predominant driving mutations found in cancer that lead to WNT pathway activation, as well as explore some of the therapeutic interventions currently available against tumours harbouring either ligand-dependent or ligand-independent mutations. Finally, we discuss a potentially new therapeutic avenue by targeting the translational apparatus downstream from WNT signalling.

A novel survival prediction signature outperforms PAM50 and artificial intelligence-based feature-selection methods

The robustness of a breast cancer gene signature, the super-proliferation set (SPS), is initially tested and investigated on breast cancer cell lines from the Cancer Cell Line Encyclopaedia (CCLE). Previously, SPS was derived via a meta-analysis of 47 independent breast cancer gene signatures, benchmarked on survival information from clinical data in the NKI dataset. Here, relying on the stability of cell line data and associative prior knowledge, we first demonstrate through Principal Component Analysis (PCA) that SPS prioritizes survival information over secondary subtype information, surpassing both PAM50 and Boruta, an artificial intelligence-based feature-selection algorithm, in this regard. We can also extract higher resolution 'progression' information using SPS, dividing survival outcomes into several clinically relevant stages ('good', 'intermediate', and 'bad) based on different quadrants of the PCA scatterplot. Furthermore, by transferring these 'progression' annotations onto independent clinical datasets, we demonstrate the generalisability of our method on actual patient data. Finally, via the characteristic genetic profiles of each quadrant/stage, we identified efficacious drugs using their gene reversal scores that can shift signatures across quadrants/stages, in a process known as gene signature reversal. This confirms the power of meta-analytical approaches for gene signature inference in breast cancer, as well as the clinical benefit in translating these inferences onto real-world patient data for more targeted therapies.

An Innovative Drug Repurposing Approach to Restrain Endometrial Cancer Metastatization

BACKGROUND

Endometrial cancer (EC) is the most common gynecologic tumor and the world's fourth most common cancer in women. Most patients respond to first-line treatments and have a low risk of recurrence, but refractory patients, and those with metastatic cancer at diagnosis, remain with no treatment options. Drug repurposing aims to discover new clinical indications for existing drugs with known safety profiles. It provides ready-to-use new therapeutic options for highly aggressive tumors for which standard protocols are ineffective, such as high-risk EC.

METHODS

Here, we aimed at defining new therapeutic opportunities for high-risk EC using an innovative and integrated computational drug repurposing approach.

RESULTS

We compared gene-expression profiles, from publicly available databases, of metastatic and non-metastatic EC patients being metastatization the most severe feature of EC aggressiveness. A comprehensive analysis of transcriptomic data through a two-arm approach was applied to obtain a robust prediction of drug candidates.

CONCLUSIONS

Some of the identified therapeutic agents are already successfully used in clinical practice to treat other types of tumors. This highlights the potential to repurpose them for EC and, therefore, the reliability of the proposed approach.

A comprehensive tool for tumor precision medicine with pharmaco-omics data analysis

Background: Cancer precision medicine is an effective strategy to fight cancers by bridging genomics and drug discovery to provide specific treatment for patients with different genetic characteristics. Although some public databases and modelling frameworks have been developed through studies on drug response, most of them only considered the ramifications of the drug on the cell line and the effects on the patient still require a huge amount of work to integrate data from various databases and calculations, especially concerning precision treatment. Furthermore, not only efficacy but also the adverse effects of drugs on patients should be taken into account during cancer treatment. However, the adverse effects as essential indicators of drug safety assessment are always neglected. Method: A holistic estimation explores various drugs' efficacy levels by calculating their potency both in reversing and enhancing cancer-associated gene expression change. And a method for bridging the gap between cell culture and living tissue estimates the effectiveness of a drug on individual patients through the mappings of various cell lines to each person according to their genetic mutation similarities. Result: We predicted the efficacy of FDA-recommended drugs, taking into account both efficacy and toxicity, and obtained consistent results. We also provided an intuitive and easy-to-use web server called DBPOM (http://www.dbpom.net/, a comprehensive database of pharmaco-omics for cancer precision medicine), which not only integrates the above methods but also provides calculation results on more than 10,000 small molecule compounds and drugs. As a one-stop web server, clinicians and drug researchers can also analyze the overall effect of a drug or a drug combination on cancer patients as well as the biological functions that they target. DBPOM is now public, free to use with no login requirement, and contains all the data and code. Conclusion: Both the positive and negative effects of drugs during precision treatment are essential for practical application of drugs. DBPOM based on the two effects will become a vital resource and analysis platform for drug development, drug mechanism studies and the discovery of new therapies.

Multiscale Analysis and Validation of Effective Drug Combinations Targeting Driver KRAS Mutations in Non-Small Cell Lung Cancer

Pharmacogenomics is a rapidly growing field with the goal of providing personalized care to every patient. Previously, we developed the Computational Analysis of Novel Drug Opportunities (CANDO) platform for multiscale therapeutic discovery to screen optimal compounds for any indication/disease by performing analytics on their interactions using large protein libraries. We implemented a comprehensive precision medicine drug discovery pipeline within the CANDO platform to determine which drugs are most likely to be effective against mutant phenotypes of non-small cell lung cancer (NSCLC) based on the supposition that drugs with similar interaction profiles (or signatures) will have similar behavior and therefore show synergistic effects. CANDO predicted that osimertinib, an EGFR inhibitor, is most likely to synergize with four KRAS inhibitors.Validation studies with cellular toxicity assays confirmed that osimertinib in combination with ARS-1620, a KRAS G12C inhibitor, and BAY-293, a pan-KRAS inhibitor, showed a synergistic effect on decreasing cellular proliferation by acting on mutant KRAS. Gene expression studies revealed that MAPK expression is strongly correlated with decreased cellular proliferation following treatment with KRAS inhibitor BAY-293, but not treatment with ARS-1620 or osimertinib. These results indicate that our precision medicine pipeline may be used to identify compounds capable of synergizing with inhibitors of KRAS G12C, and to assess their likelihood of becoming drugs by understanding their behavior at the proteomic/interactomic scales.

Signature-Based Computational Drug Repurposing for Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a late-onset fatal neurodegenerative disease characterized by progressive loss of the upper and lower motor neurons. There are currently limited approved drugs for the disorder, and for this reason the strategy of repositioning already approved therapeutics could exhibit a successful outcome. Herein, we used CMAP and L1000CDS2 databases which include gene expression profiles datasets (genomic signatures) to identify potent compounds and classes of compounds which reverse disease's signature which could in turn reverse its phenotype. ALS signature was obtained by comparing gene expression of muscle biopsy specimens between diseased and healthy individuals. Statistical analysis was conducted to explore differentially transcripts in patients' samples. Then, the list of upregulated and downregulated genes was used to query both databases in order to determine molecules which downregulate the genes which are upregulated by ALS and vice versa. These compounds, based on their chemical structure along with known treatments, were clustered to reveal drugs with novel and potentially more effective mode of action with most of them predicted to affect pathways heavily involved in ALS. This evidence suggests that these compounds are strong candidates for moving to the next phase of the drug repurposing pipeline which is in vitro and in vivo experimental evaluation.

Drug repositioning prediction for psoriasis using the adverse event reporting database

Introduction

Inverse signals produced from disproportional analyses using spontaneous drug adverse event reports can be used for drug repositioning purposes. The purpose of this study is to predict drug candidates using a computational method that integrates reported drug adverse event data, disease-specific gene expression profiles, and drug-induced gene expression profiles.

Methods

Drug and adverse events from 2015 through 2020 were downloaded from the United States Food and Drug Administration Adverse Event Reporting System (FAERS). The reporting odds ratio (ROR), information component (IC) and empirical Bayes geometric mean (EBGM) were used to calculate the inverse signals. Psoriasis was selected as the target disease. Disease specific gene expression profiles were obtained by the meta-analysis of the Gene Expression Omnibus (GEO). The reverse gene expression scores were calculated using the Library of Integrated Network-based Cellular Signatures (LINCS) and their correlations with the inverse signals were obtained.

Results

Reversal genes and the candidate compounds were identified. Additionally, these correlations were validated using the relationship between the reverse gene expression scores and the half-maximal inhibitory concentration (IC50) values from the Chemical European Molecular Biology Laboratory (ChEMBL).

Conclusion

Inverse signals produced from a disproportional analysis can be used for drug repositioning and to predict drug candidates against psoriasis.

Validation of transcriptome signature reversion for drug repurposing in oncology

Transcriptome signature reversion (TSR) has been extensively proposed and used to discover new indications for existing drugs (i.e. drug repositioning, drug repurposing) for various cancer types. TSR relies on the assumption that a drug that can revert gene expression changes induced by a disease back to original, i.e. healthy, levels is likely to be therapeutically active in treating the disease. Here, we aimed to validate the concept of TSR using the PRISM repurposing data set, which is-as of writing-the largest pharmacogenomic data set. The predictive utility of the TSR approach as it has currently been used appears to be much lower than previously reported and is completely nullified after the drug gene expression signatures are adjusted for the general anti-proliferative downstream effects of drug-induced decreased cell viability. Therefore, TSR mainly relies on generic anti-proliferative drug effects rather than on targeting cancer pathways specifically upregulated in tumor types.

Artificial intelligence to guide precision anticancer therapy with multitargeted kinase inhibitors

BACKGROUND

Vast amounts of rapidly accumulating biological data related to cancer and a remarkable progress in the field of artificial intelligence (AI) have paved the way for precision oncology. Our recent contribution to this area of research is CancerOmicsNet, an AI-based system to predict the therapeutic effects of multitargeted kinase inhibitors across various cancers. This approach was previously demonstrated to outperform other deep learning methods, graph kernel models, molecular docking, and drug binding pocket matching.

METHODS

CancerOmicsNet integrates multiple heterogeneous data by utilizing a deep graph learning model with sophisticated attention propagation mechanisms to extract highly predictive features from cancer-specific networks. The AI-based system was devised to provide more accurate and robust predictions than data-driven therapeutic discovery using gene signature reversion.

RESULTS

Selected CancerOmicsNet predictions obtained for "unseen" data are positively validated against the biomedical literature and by live-cell time course inhibition assays performed against breast, pancreatic, and prostate cancer cell lines. Encouragingly, six molecules exhibited dose-dependent antiproliferative activities, with pan-CDK inhibitor JNJ-7706621 and Src inhibitor PP1 being the most potent against the pancreatic cancer cell line Panc 04.03.

CONCLUSIONS

CancerOmicsNet is a promising AI-based platform to help guide the development of new approaches in precision oncology involving a variety of tumor types and therapeutics.

Drug repurposing using meta-analysis of gene expression in Alzheimer's disease

INTRODUCTION

Alzheimer's disease and other forms of dementia are disease that bring an increased global burden. However, the medicine developed to date remains limited. The purpose of this study is to predict drug repositioning candidates using a computational method that integrates gene expression profiles on Alzheimer's disease and compound-induced changes in gene expression levels.

METHODS

Gene expression data on Alzheimer's disease were obtained from the Gene Expression Omnibus (GEO) and we conducted a meta-analysis of their gene expression levels. The reverse scores of compound-induced gene expressions were computed based on the reversal relationship between disease and drug gene expression profiles.

RESULTS

Reversal genes and the candidate compounds were identified by the leave-one-out cross-validation procedure. Additionally, the half-maximal inhibitory concentration (IC50) values and the blood-brain barrier (BBB) permeability of candidate compounds were obtained from ChEMBL and PubChem, respectively.

CONCLUSION

New therapeutic target genes and drug candidates against Alzheimer's disease were identified by means of drug repositioning.

Bioinformatics roadmap for therapy selection in cancer genomics

Tumour heterogeneity is one of the main characteristics of cancer and can be categorised into inter- or intratumour heterogeneity. This heterogeneity has been revealed as one of the key causes of treatment failure and relapse. Precision oncology is an emerging field that seeks to design tailored treatments for each cancer patient according to epidemiological, clinical and omics data. This discipline relies on bioinformatics tools designed to compute scores to prioritise available drugs, with the aim of helping clinicians in treatment selection. In this review, we describe the current approaches for therapy selection depending on which type of tumour heterogeneity is being targeted and the available next-generation sequencing data. We cover intertumour heterogeneity studies and individual treatment selection using genomics variants, expression data or multi-omics strategies. We also describe intratumour dissection through clonal inference and single-cell transcriptomics, in each case providing bioinformatics tools for tailored treatment selection. Finally, we discuss how these therapy selection workflows could be integrated into the clinical practice.

Reversal of cancer gene expression identifies repurposed drugs for diffuse intrinsic pontine glioma

Diffuse intrinsic pontine glioma (DIPG) is an aggressive incurable brainstem tumor that targets young children. Complete resection is not possible, and chemotherapy and radiotherapy are currently only palliative. This study aimed to identify potential therapeutic agents using a computational pipeline to perform an in silico screen for novel drugs. We then tested the identified drugs against a panel of patient-derived DIPG cell lines. Using a systematic computational approach with publicly available databases of gene signature in DIPG patients and cancer cell lines treated with a library of clinically available drugs, we identified drug hits with the ability to reverse a DIPG gene signature to one that matches normal tissue background. The biological and molecular effects of drug treatment was analyzed by cell viability assay and RNA sequence. In vivo DIPG mouse model survival studies were also conducted. As a result, two of three identified drugs showed potency against the DIPG cell lines Triptolide and mycophenolate mofetil (MMF) demonstrated significant inhibition of cell viability in DIPG cell lines. Guanosine rescued reduced cell viability induced by MMF. In vivo, MMF treatment significantly inhibited tumor growth in subcutaneous xenograft mice models. In conclusion, we identified clinically available drugs with the ability to reverse DIPG gene signatures and anti-DIPG activity in vitro and in vivo. This novel approach can repurpose drugs and significantly decrease the cost and time normally required in drug discovery.

Deciphering COVID-19 host transcriptomic complexity and variations for therapeutic discovery against new variants

The molecular manifestations of host cells responding to SARS-CoV-2 and its evolving variants of infection are vastly different across the studied models and conditions, imposing challenges for host-based antiviral drug discovery. Based on the postulation that antiviral drugs tend to reverse the global host gene expression induced by viral infection, we retrospectively evaluated hundreds of signatures derived from 1,700 published host transcriptomic profiles of SARS/MERS/SARS-CoV-2 infection using an iterative data-driven approach. A few of these signatures could be reversed by known anti-SARS-CoV-2 inhibitors, suggesting the potential of extrapolating the biology for new variant research. We discovered IMD-0354 as a promising candidate to reverse the signatures globally with nanomolar IC50 against SARS-CoV-2 and its five variants. IMD-0354 stimulated type I interferon antiviral response, inhibited viral entry, and down-regulated hijacked proteins. This study demonstrates that the conserved coronavirus signatures and the transcriptomic reversal approach that leverages polypharmacological effects could guide new variant therapeutic discovery.

Therapeutic drug repositioning with special emphasis on neurodegenerative diseases: Threats and issues

Drug repositioning or repurposing is the process of discovering leading-edge indications for authorized or declined/abandoned molecules for use in different diseases. This approach revitalizes the traditional drug discovery method by revealing new therapeutic applications for existing drugs. There are numerous studies available that highlight the triumph of several drugs as repurposed therapeutics. For example, sildenafil to aspirin, thalidomide to adalimumab, and so on. Millions of people worldwide are affected by neurodegenerative diseases. According to a 2021 report, the Alzheimer's disease Association estimates that 6.2 million Americans are detected with Alzheimer's disease. By 2030, approximately 1.2 million people in the United States possibly acquire Parkinson's disease. Drugs that act on a single molecular target benefit people suffering from neurodegenerative diseases. Current pharmacological approaches, on the other hand, are constrained in their capacity to unquestionably alter the course of the disease and provide patients with inadequate and momentary benefits. Drug repositioning-based approaches appear to be very pertinent, expense- and time-reducing strategies for the enhancement of medicinal opportunities for such diseases in the current era. Kinase inhibitors, for example, which were developed for various oncology indications, demonstrated significant neuroprotective effects in neurodegenerative diseases. This review expounds on the classical and recent examples of drug repositioning at various stages of drug development, with a special focus on neurodegenerative disorders and the aspects of threats and issues viz. the regulatory, scientific, and economic aspects.