Breaking the paradigm: Dr Insight empowers signature-free, enhanced drug repurposing

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.

DrugReSC: targeting disease-critical cell subpopulations with single-cell transcriptomic data for drug repurposing in cancer

The field of computational drug repurposing aims to uncover novel therapeutic applications for existing drugs through high-throughput data analysis. However, there is a scarcity of drug repurposing methods leveraging the cellular-level information provided by single-cell RNA sequencing data. To address this need, we propose DrugReSC, an innovative approach to drug repurposing utilizing single-cell RNA sequencing data, intending to target specific cell subpopulations critical to disease pathology. DrugReSC constructs a drug-by-cell matrix representing the transcriptional relationships between individual cells and drugs and utilizes permutation-based methods to assess drug contributions to cellular phenotypic changes. We demonstrate DrugReSC's superior performance compared to existing drug repurposing methods based on bulk or single-cell RNA sequencing data across multiple cancer case studies. In summary, DrugReSC offers a novel perspective on the utilization of single-cell sequencing data in drug repurposing methods, contributing to the advancement of precision medicine for cancer.

Novel insights into the progression and prognosis of the calpain family members in hepatocellular carcinoma: a comprehensive integrated analysis

Objectives: The goal of our bioinformatics study was to comprehensively analyze the association between the whole calpain family members and the progression and prognosis of hepatocellular carcinoma (HCC). Methods: The data were collected from The Cancer Genome Atlas (TCGA). The landscape of the gene expression, copy number variation (CNV), mutation, and DNA methylation of calpain members were analyzed. Clustering analysis was performed to stratify the calpain-related groups. The least absolute shrinkage and selection operator (LASSO)-based Cox model was used to select hub survival genes. Results: We found 14 out of 16 calpain members expressed differently between tumor and normal tissues of HCC. The clustering analyses revealed high- and low-risk calpain groups which had prognostic difference. We found the high-risk calpain group had higher B cell infiltration and higher expression of immune checkpoint genes HAVCR2, PDCD1, and TIGHT. The CMap analysis found that the histone deacetylase (HDAC) inhibitor trichostatin A and the PI3K-AKT-mTOR pathway inhibitors LY-294002 and wortmannin might have a therapeutic effect on the high-risk calpain group. The DEGs between calpain groups were identified. Subsequent univariate Cox analysis of each DEG and LASSO-based Cox model obtained a calpain-related prognostic signature. The risk score model of this signature showed good ability to predict the overall survival of HCC patients in TCGA datasets and external validation datasets from the Gene Expression Omnibus database and the International Cancer Genome Consortium database. Conclusion: We found that calpain family members were associated with the progression, prognosis, and drug response of HCC. Our results require further studies to confirm.

Informatics on Drug Repurposing for Breast Cancer

Moving a new drug from bench to bedside is a long and arduous process. The tactic of drug repurposing, which solves "new" diseases with "old" existing drugs, is more efficient and economical than conventional ab-initio way for drug development. Information technology has dramatically changed the paradigm of biomedical research in the new century, and drug repurposing studies have been significantly accelerated by implementing informatics techniques related to genomics, systems biology and biophysics during the past few years. A series of remarkable achievements in this field comes with the practical applications of in silico approaches including transcriptomic signature matching, gene-connection-based scanning, and simulated structure docking in repositioning drug therapies against breast cancer. In this review, we systematically curated these impressive accomplishments with summarization of the main findings on potentially repurposable drugs, and provide our insights into the current issues as well as future directions of the field. With the prospective improvement in reliability, the computer-assisted repurposing strategy will play a more critical role in drug research and development.

Drug mechanism enrichment analysis improves prioritization of therapeutics for repurposing

BACKGROUND

There is a pressing need for improved methods to identify effective therapeutics for diseases. Many computational approaches have been developed to repurpose existing drugs to meet this need. However, these tools often output long lists of candidate drugs that are difficult to interpret, and individual drug candidates may suffer from unknown off-target effects. We reasoned that an approach which aggregates information from multiple drugs that share a common mechanism of action (MOA) would increase on-target signal compared to evaluating drugs on an individual basis. In this study, we present drug mechanism enrichment analysis (DMEA), an adaptation of gene set enrichment analysis (GSEA), which groups drugs with shared MOAs to improve the prioritization of drug repurposing candidates.

RESULTS

First, we tested DMEA on simulated data and showed that it can sensitively and robustly identify an enriched drug MOA. Next, we used DMEA on three types of rank-ordered drug lists: (1) perturbagen signatures based on gene expression data, (2) drug sensitivity scores based on high-throughput cancer cell line screening, and (3) molecular classification scores of intrinsic and acquired drug resistance. In each case, DMEA detected the expected MOA as well as other relevant MOAs. Furthermore, the rankings of MOAs generated by DMEA were better than the original single-drug rankings in all tested data sets. Finally, in a drug discovery experiment, we identified potential senescence-inducing and senolytic drug MOAs for primary human mammary epithelial cells and then experimentally validated the senolytic effects of EGFR inhibitors.

CONCLUSIONS

DMEA is a versatile bioinformatic tool that can improve the prioritization of candidates for drug repurposing. By grouping drugs with a shared MOA, DMEA increases on-target signal and reduces off-target effects compared to analysis of individual drugs. DMEA is publicly available as both a web application and an R package at https://belindabgarana.github.io/DMEA .

ASGARD is A Single-cell Guided Pipeline to Aid Repurposing of Drugs

Single-cell RNA sequencing technology has enabled in-depth analysis of intercellular heterogeneity in various diseases. However, its full potential for precision medicine has yet to be reached. Towards this, we propose A Single-cell Guided Pipeline to Aid Repurposing of Drugs (ASGARD) that defines a drug score to recommend drugs by considering all cell clusters to address the intercellular heterogeneity within each patient. ASGARD shows significantly better average accuracy on single-drug therapy compared to two bulk-cell-based drug repurposing methods. We also demonstrated that it performs considerably better than other cell cluster-level predicting methods. In addition, we validate ASGARD using the drug response prediction method TRANSACT with Triple-Negative-Breast-Cancer patient samples. We find that many top-ranked drugs are either approved by the Food and Drug Administration or in clinical trials treating corresponding diseases. In conclusion, ASGARD is a promising drug repurposing recommendation tool guided by single-cell RNA-seq for personalized medicine. ASGARD is free for educational use at https://github.com/lanagarmire/ASGARD .

A Comprehensive Analysis of the Effects of Key Mitophagy Genes on the Progression and Prognosis of Lung Adenocarcinoma

The aim of our study was to perform a comprehensive analysis of the gene expression, copy number variation (CNV) and mutation of key mitophagy genes in the progression and prognosis of lung adenocarcinoma (LUAD). We obtained the data from The Cancer Genome Atlas (TCGA). Clustering analysis was performed to stratify the mitophagy related groups. The least absolute shrinkage and selection operator (LASSO) based cox model was used to select hub survival genes. An independent validation cohort was retrieved from Gene Expression Omnibus database. We found 24 out of 27 mitophagy genes were aberrantly expressed between tumor and normal samples. CNV gains were associated with higher expression of mitophagy genes in 23 of 27 mitophagy genes. The clustering analysis identified high and low risk mitophagy groups with distinct survival differences. The high risk mitophagy groups had higher tumor mutation burden, stemness phenotype, total CNVs and lower CD4+ T cells infiltration. Drugs targeted to high risk mitophagy groups were identified including the PI3K/AKT/mTOR inhibitor, HDAC inhibitor and chemotherapy agents such as cisplatin and gemcitabine. In addition, the differentially expressed genes (DEGs) were identified between mitophagy groups. Further univariate Cox analysis of each DEG and subsequent LASSO-based Cox model revealed a mitophagy-related prognostic signature. The risk score model of this signature showed a strong ability to predict the overall survival of LUAD patients in training and validation datasets. In conclusion, the mitophagy genes played an important role in the progression and prognosis of LUAD, which might provide useful information for the treatment of LUAD.

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.

Analysis of transcriptomic responses to SARS-CoV-2 reveals plausible defective pathways responsible for increased susceptibility to infection and complications and helps to develop fast-track repositioning of drugs against COVID-19

BACKGROUND

To understand the transcriptomic response to SARS-CoV-2 infection, is of the utmost importance to design diagnostic tools predicting the severity of the infection.

METHODS

We have performed a deep sampling analysis of the viral transcriptomic data oriented towards drug repositioning. Using different samplers, the basic principle of this methodology the biological invariance, which means that the pathways altered by the disease, should be independent on the algorithm used to unravel them.

RESULTS

The transcriptomic analysis of the altered pathways, reveals a distinctive inflammatory response and potential side effects of infection. The virus replication causes, in some cases, acute respiratory distress syndrome in the lungs, and affects other organs such as heart, brain, and kidneys. Therefore, the repositioned drugs to fight COVID-19 should, not only target the interferon signalling pathway and the control of the inflammation, but also the altered genetic pathways related to the side effects of infection. We also show via Principal Component Analysis that the transcriptome signatures are different from influenza and RSV. The gene COL1A1, which controls collagen production, seems to play a key/vital role in the regulation of the immune system. Additionally, other small-scale signature genes appear to be involved in the development of other COVID-19 comorbidities.

CONCLUSIONS

Transcriptome-based drug repositioning offers possible fast-track antiviral therapy for COVID-19 patients. It calls for additional clinical studies using FDA approved drugs for patients with increased susceptibility to infection and with serious medical complications.

CPDR: An R Package of Recommending Personalized Drugs for Cancer Patients by Reversing the Individual’s Disease-Related Signature

Due to cancer heterogeneity, only some patients can benefit from drug therapy. The personalized drug usage is important for improving the treatment response rate of cancer patients. The value of the transcriptome of patients has been recently demonstrated in guiding personalized drug use, and the Connectivity Map (CMAP) is a reliable computational approach for drug recommendation. However, there is still no personalized drug recommendation tool based on transcriptomic profiles of patients and CMAP. To fill this gap, here, we proposed such a feasible workflow and a user-friendly R package—Cancer-Personalized Drug Recommendation (CPDR). CPDR has three features. 1) It identifies the individual disease signature by using the patient subgroup with transcriptomic profiles similar to those of the input patient. 2) Transcriptomic profile purification is supported for the subgroup with high infiltration of non-cancerous cells. 3) It supports in silico drug efficacy assessment using drug sensitivity data on cancer cell lines. We demonstrated the workflow of CPDR with the aid of a colorectal cancer dataset from GEO and performed the in silico validation of drug efficacy. We further assessed the performance of CPDR by a pancreatic cancer dataset with clinical response to gemcitabine. The results showed that CPDR can recommend promising therapeutic agents for the individual patient. The CPDR R package is available at https://github.com/AllenSpike/CPDR.

Recent Advances in Cancer Drug Discovery Through the Use of Phenotypic Reporter Systems, Connectivity Mapping, and Pooled CRISPR Screening

Multi-omic approaches offer an unprecedented overview of the development, plasticity, and resistance of cancer. However, the translation from anti-cancer compounds identified in vitro to clinically active drugs have a notoriously low success rate. Here, we review how technical advances in cell culture, robotics, computational biology, and development of reporter systems have transformed drug discovery, enabling screening approaches tailored to clinically relevant functional readouts (e.g., bypassing drug resistance). Illustrating with selected examples of “success stories,” we describe the process of phenotype-based high-throughput drug screening to target malignant cells or the immune system. Second, we describe computational approaches that link transcriptomic profiling of cancers with existing pharmaceutical compounds to accelerate drug repurposing. Finally, we review how CRISPR-based screening can be applied for the discovery of mechanisms of drug resistance and sensitization. Overall, we explore how the complementary strengths of each of these approaches allow them to transform the paradigm of pre-clinical drug development.

Identification of HOX signatures contributing to oral cancer phenotype

The role of evolutionarily conserved homeobox-containing HOX genes as transcriptional regulators in the developmental specification of organisms is well known. The contribution of HOX genes involvement in oral cancer phenotype has yet to be fully ascertained. TCGA-HNSC HTSeq-counts and clinical data were retrieved from the GDC portal for oral cavity neoplasms. GEO datasets (GSE72627, GSE30784, GSE37991) were accessed and analyzed using GEO2R. Differential HOX gene expression was profiled using the DESeq2 R package with a log2 fold change cut-off (- 1 and + 1) and Benjamini-Hochberg p-adjusted value at ≤ 0.01. Gene set over-representation analysis and semantic analysis associated with the disease ontology was performed using the ClusterProfiler R package, and pathway over-representation analysis was performed using IMPaLa. HOX protein interaction network was constructed using the Pathfind R package. HOX phenotype associations were performed using Mammalian Phenotype Ontology, Human Phenotype Ontology, PhenGenI associations, Jensen tissues, and OMIM entries. Drug connectivity mapping was carried out with Dr. Insight R package. HOXA2 was upregulated in oral dysplasia but silenced during tumor progression. Loss of HOXB2 expression was consistent in the potentially malignant oral lesions as well as in the primary tumor. HOXA7, HOXA10, HOXB7, HOXC6, HOXC10, HOXD10, and HOXD11 were consistently upregulated from premalignancy to malignancy and were notably associated with risk factors. Overrepresentation analysis suggested HOXA10 was involved in the transcriptional misregulation contributing to the oral cancer phenotype. HOX genes subnetwork analysis showed crucial interactions with cell cycle regulators, growth responsive elements, and proto-oncogenes. Phenotype associations specific to the oral region involving HOX genes provide intrinsic cues to tumor development. The 5' HOX genes were aberrantly upregulated during oral carcinogenesis reflecting their posterior prevalence.

The Diagnostic and Prognostic Values of HOXA Gene Family in Kidney Clear Cell Renal Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) is one of the most common cancers with high mortality worldwide. As members of the homeobox (HOX) family, homeobox-A (HOXA) genes have been reported to play an increasingly important role in tumorigenesis and the progression of multiple cancers. However, limited studies have investigated the potential diagnostic and prognostic roles of HOXA genes in KIRC. In this research, we explored the expression pattern of the HOXA gene family in KIRC progression by differential analysis of expression profiles from The Cancer Genome Atlas (TCGA). By using univariate Cox analysis and lasso regression analysis, we comprehensively evaluated the prognostic value of HOXA genes and eventually identified a prognostic risk model consisting of five HOXA genes (HOXA2, HOXA3, HOXA7, HOXA11, and HOXA13). The risk model was further validated as a novel independent prognostic factor for KIRC patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses, and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, to explore the potential mechanism of tumorigenesis and clinical application of KIRC, we also developed the HOXA-based competing endogenous RNA (ceRNA) regulatory network and machine learning classification model. Valproic acid and tretinoin were predicted to be the most promising small molecules to adjuvant treatment of KIRC by mining the CMAP and DGIdb drug database. Subsequently, pathway and functional enrichment analyses provided us with new ways to search for a possible mechanism of action of drugs. Taken together, our study demonstrated the nonnegligible role of HOXA genes in KIRC and constructed an effective prognostic and diagnostic model, which offers novel insights into KIRC prognosis.

Considerations and challenges for sex-aware drug repurposing

Sex differences are essential factors in disease etiology and manifestation in many diseases such as cardiovascular disease, cancer, and neurodegeneration [33]. The biological influence of sex differences (including genomic, epigenetic, hormonal, immunological, and metabolic differences between males and females) and the lack of biomedical studies considering sex differences in their study design has led to several policies. For example, the National Institute of Health's (NIH) sex as a biological variable (SABV) and Sex and Gender Equity in Research (SAGER) policies to motivate researchers to consider sex differences [204]. However, drug repurposing, a promising alternative to traditional drug discovery by identifying novel uses for FDA-approved drugs, lacks sex-aware methods that can improve the identification of drugs that have sex-specific responses [7, 11, 14, 33]. Sex-aware drug repurposing methods either select drug candidates that are more efficacious in one sex or deprioritize drug candidates based on if they are predicted to cause a sex-bias adverse event (SBAE), unintended therapeutic effects that are more likely to occur in one sex. Computational drug repurposing methods are encouraging approaches to develop for sex-aware drug repurposing because they can prioritize sex-specific drug candidates or SBAEs at lower cost and time than traditional drug discovery. Sex-aware methods currently exist for clinical, genomic, and transcriptomic information [1, 7, 155]. They have not expanded to other data types, such as DNA variation, which has been beneficial in other drug repurposing methods that do not consider sex [114]. Additionally, some sex-aware methods suffer from poorer performance because a disproportionate number of male and female samples are available to train computational methods [7]. However, there is development potential for several different categories (i.e., data mining, ligand binding predictions, molecular associations, and networks). Low-dimensional representations of molecular association and network approaches are also especially promising candidates for future sex-aware drug repurposing methodologies because they reduce the multiple hypothesis testing burden and capture sex-specific variation better than the other methods [151, 159]. Here we review how sex influences drug response, the current state of drug repurposing including with respect to sex-bias drug response, and how model organism study design choices influence drug repurposing validation.

A survey of optimal strategy for signature-based drug repositioning and an application to liver cancer

Pharmacologic perturbation projects, such as Connectivity Map (CMap) and Library of Integrated Network-based Cellular Signatures (LINCS), have produced many perturbed expression data, providing enormous opportunities for computational therapeutic discovery. However, there is no consensus on which methodologies and parameters are the most optimal to conduct such analysis. Aiming to fill this gap, new benchmarking standards were developed to quantitatively evaluate drug retrieval performance. Investigations of potential factors influencing drug retrieval were conducted based on these standards. As a result, we determined an optimal approach for LINCS data-based therapeutic discovery. With this approach, homoharringtonine (HHT) was identified to be a candidate agent with potential therapeutic and preventive effects on liver cancer. The antitumor and antifibrotic activity of HHT was validated experimentally using subcutaneous xenograft tumor model and carbon tetrachloride (CCL₄)-induced liver fibrosis model, demonstrating the reliability of the prediction results. In summary, our findings will not only impact the future applications of LINCS data but also offer new opportunities for therapeutic intervention of liver cancer.

A gene signature consisting of ubiquitin ligases and deubiquitinating enzymes of SKP2 is associated with clinical outcome in breast cancer

The ubiquitination of SKP2, an oncoprotein, is controlled by its E3 ligases, including APC/C^FZR1 and deubiquitinases such as USP10. We identified a two-gene signature for the ubiquitination of SKP2, consisting of the copy number of FZR1 compared to the copy number of USP10. The signature reflects the level of SKP2 activity, stratifying BC patients into two groups with significantly different protein levels of SKP2 ubiquitination substrate p27 (t-test p < 0.01) and recapitulating the expression patterns of SKP2 between tumor and normal tissue (Spearman’s ρ = 0.39.) The signature is also highly associated with clinical outcome in luminal BC but not other subtypes, characterizing patients into two groups with significantly different overall survival times (log-rank p = 0.006). In addition, it is dramatically associated with tumor grade (Chi-squared p = 6.7 × 10⁻³), stage (Chi-squared p = 1.6 × 10⁻¹¹), and the number of positive lymph nodes (negative binomial regression coefficient p = 2.0 × 10⁻³). Our study provides a rationale for targeting the SKP2 ubiquitination pathway in luminal BC and for further investigation of the use of ubiquitinase/deubiquitinase genes as prognosis and treatment biomarkers.

Combination of Radiosensitivity Gene Signature and PD-L1 Status Predicts Clinical Outcome of Patients With Locally Advanced Head and Neck Squamous Cell Carcinoma: A Study Based on The Cancer Genome Atlas Dataset

Aim: The aim of our study was to investigate the potential predictive value of the combination of radiosensitivity gene signature and PD-L1 expression for the prognosis of locally advanced head and neck squamous cell carcinoma (HNSCC).

Methods: The cohort was selected from The Cancer Genome Atlas (TCGA) and classified into the radiosensitive (RS) group and radioresistant (RR) group by a radiosensitivity-related gene signature. The cohort was also grouped as PD-L1-high or PD-L1-low based on PD-L1 mRNA expression. The least absolute shrinkage and selection operator (lasso)-based Cox model was used to select hub survival genes. An independent validation cohort was obtained from the Gene Expression Omnibus (GEO) database.

Results: We selected 288 locally advanced HNSCC patients from TCGA. The Kaplan–Meier method found that the RR and PD-L1-high group had a worse survival than others (p = 0.033). The differentially expressed gene (DEG) analysis identified 553 upregulated genes and 486 downregulated genes (p < 0.05, fold change >2) between the RR and PD-L1-high group and others. The univariate Cox analysis of each DEG and subsequent lasso-based Cox model revealed five hub survival genes (POU4F1, IL34, HLF, CBS, and RNF165). A further hub survival gene-based risk score model was constructed, which was validated by an external cohort. We observed that a higher risk score predicted a worse prognosis (p = 0.0013). The area under the receiver operating characteristic curve (AUC) plots showed that this risk score model had good prediction value (1-year AUC = 0.684, 2-year AUC = 0.702, and 3-year AUC = 0.688). Five different deconvolution methods all showed that the B cells were lower in the RR and PD-L1-high group (p < 0.05). Finally, connectivity mapping analysis showed that the histone deacetylase (HDAC) inhibitor trichostatin A might have the potential to reverse the phenotype of RR and PD-L1-high in locally advanced HNSCC (p < 0.05, false discovery rate <0.1).

Conclusion: The combination of 31-gene signature and the PD-L1 mRNA expression had a potential predictive value for the prognosis of locally advanced HNSCC who had RT. The B cells were lower in the RR and PD-L1-high group. The identified risk gene signature of locally advanced HNSCC and the potential therapeutic drug trichostatin A for the RR and PD-L1-high group are worth being further studied in a prospective homogenous cohort.

Drug repurposing for COVID-19: Approaches, challenges and promising candidates

Traditional drug development and discovery has not kept pace with threats from emerging and re-emerging diseases such as Ebola virus, MERS-CoV and more recently, SARS-CoV-2. Among other reasons, the exorbitant costs, high attrition rate and extensive periods of time from research to market approval are the primary contributing factors to the lag in recent traditional drug developmental activities. Due to these reasons, drug developers are starting to consider drug repurposing (or repositioning) as a viable alternative to the more traditional drug development process. Drug repurposing aims to find alternative uses of an approved or investigational drug outside of its original indication. The key advantages of this approach are that there is less developmental risk, and it is less time-consuming since the safety and pharmacological profile of the repurposed drug is already established. To that end, various approaches to drug repurposing are employed. Computational approaches make use of machine learning and algorithms to model disease and drug interaction, while experimental approaches involve a more traditional wet-lab experiments. This review would discuss in detail various ongoing drug repurposing strategies and approaches to combat the current COVID-19 pandemic, along with the advantages and the potential challenges.

Using predictive machine learning models for drug response simulation by calibrating patient-specific pathway signatures

The utility of pathway signatures lies in their capability to determine whether a specific pathway or biological process is dysregulated in a given patient. These signatures have been widely used in machine learning (ML) methods for a variety of applications including precision medicine, drug repurposing, and drug discovery. In this work, we leverage highly predictive ML models for drug response simulation in individual patients by calibrating the pathway activity scores of disease samples. Using these ML models and an intuitive scoring algorithm to modify the signatures of patients, we evaluate whether a given sample that was formerly classified as diseased, could be predicted as normal following drug treatment simulation. We then use this technique as a proxy for the identification of potential drug candidates. Furthermore, we demonstrate the ability of our methodology to successfully identify approved and clinically investigated drugs for four different cancers, outperforming six comparable state-of-the-art methods. We also show how this approach can deconvolute a drugs' mechanism of action and propose combination therapies. Taken together, our methodology could be promising to support clinical decision-making in personalized medicine by simulating a drugs' effect on a given patient.

ASGARD: A Single-cell Guided pipeline to Aid Repurposing of Drugs

Intercellular heterogeneity is a major obstacle to successful precision medicine. Single-cell RNA sequencing (scRNA-seq) technology has enabled in-depth analysis of intercellular heterogeneity in various diseases. However, its full potential for precision medicine has yet to be reached. Towards this, we propose a new drug recommendation system called: A Single-cell Guided Pipeline to Aid Repurposing of Drugs (ASGARD). ASGARD defines a novel drug score predicting drugs by considering all cell clusters to address the intercellular heterogeneity within each patient. We tested ASGARD on multiple diseases, including breast cancer, acute lymphoblastic leukemia, and coronavirus disease 2019 (COVID-19). On single-drug therapy, ASGARD shows significantly better average accuracy (AUC of 0.92) compared to two other bulk-cell-based drug repurposing methods (AUC of 0.80 and 0.76). It is also considerably better (AUC of 0.82) than other cell cluster level predicting methods (AUC of 0.67 and 0.55). In addition, ASGARD is also validated by the drug response prediction method TRANSACT with Triple-Negative-Breast-Cancer patient samples. Many top-ranked drugs are either approved by FDA or in clinical trials treating corresponding diseases. In silico cell-type specific drop-out experiments using triple-negative breast cancers show the importance of T cells in the tumor microenvironment in affecting drug predictions. In conclusion, ASGARD is a promising drug repurposing recommendation tool guided by single-cell RNA-seq for personalized medicine. ASGARD is free for educational use at https://github.com/lanagarmire/ASGARD.

Strategies for the treatment of breast cancer: from classical drugs to mathematical models

Breast cancer is one of the most common cancers and generally affects women. It is a heterogeneous disease that presents different entities, different biological characteristics, and differentiated clinical behaviors. With this in mind, this literature review had as its main objective to analyze the path taken from the simple use of classical drugs to the application of mathematical models, which through the many ongoing studies, have been considered as one of the reliable strategies, explaining the reasons why chemotherapy is not always successful. Besides, the most commonly mentioned strategies are immunotherapy, which includes techniques and therapies such as the use of antibodies, cytokines, antitumor vaccines, oncolytic and genomic viruses, among others, and nanoparticles, including metallic, magnetic, polymeric, liposome, dendrimer, micelle, and others, as well as drug reuse, which is a process by which new therapeutic indications are found for existing and approved drugs. The most commonly used pharmacological categories are cardiac, antiparasitic, anthelmintic, antiviral, antibiotic, and others. For the efficient development of reused drugs, there must be a process of exchange of purposes, methods, and information already available, and for their better understanding, computational mathematical models are then used, of which the methods of blind search or screening, based on the target, knowledge, signature, pathway or network and the mechanism to which it is directed, stand out. To conclude it should be noted that these different strategies can be applied alone or in combination with each other always to improve breast cancer treatment.

Drug Repurposing in Oncology, an Attractive Opportunity for Novel Combinatorial Regimens

The costs of developing, validating and buying new drugs are dramatically increasing. On the other hand, sobering economies have difficulties in sustaining their healthcare systems, particularly in countries with an elderly population requiring increasing welfare. This conundrum requires immediate action, and a possible option is to study the large, already present arsenal of drugs approved and to use them for innovative therapies. This possibility is particularly interesting in oncology, where the complexity of the cancer genome dictates in most patients a multistep therapeutic approach. In this review, we discuss a) Computational approaches; b) preclinical models; c) currently ongoing or already published clinical trials in the drug repurposing field in oncology; and d) drug repurposing to overcome resistance to previous therapies.

Scoring functions for drug-effect similarity

MOTIVATION

The difficulty to find new drugs and bring them to the market has led to an increased interest to find new applications for known compounds. Biological samples from many disease contexts have been extensively profiled by transcriptomics, and, intuitively, this motivates to search for compounds with a reversing effect on the expression of characteristic disease genes. However, disease effects may be cell line-specific and also depend on other factors, such as genetics and environment. Transcription profile changes between healthy and diseased cells relate in complex ways to profile changes gathered from cell lines upon stimulation with a drug. Despite these differences, we expect that there will be some similarity in the gene regulatory networks at play in both situations. The challenge is to match transcriptomes for both diseases and drugs alike, even though the exact molecular pathology/pharmacogenomics may not be known.

RESULTS

We substitute the challenge to match a drug effect to a disease effect with the challenge to match a drug effect to the effect of the same drug at another concentration or in another cell line. This is welldefined, reproducible in vitro and in silico and extendable with external data. Based on the Connectivity Map (CMap) dataset, we combined 26 different similarity scores with six different heuristics to reduce the number of genes in the model. Such gene filters may also utilize external knowledge e.g. from biological networks. We found that no similarity score always outperforms all others for all drugs, but the Pearson correlation finds the same drug with the highest reliability. Results are improved by filtering for highly expressed genes and to a lesser degree for genes with large fold changes. Also a network-based reduction of contributing transcripts was beneficial, here implemented by the FocusHeuristics. We found no drop in prediction accuracy when reducing the whole transcriptome to the set of 1000 landmark genes of the CMap's successor project Library of Integrated Network-based Cellular Signatures. All source code to re-analyze and extend the CMap data, the source code of heuristics, filters and their evaluation are available to propel the development of new methods for drug repurposing.

AVAILABILITY

https://bitbucket.org/ibima/moldrugeffectsdb.

CONTACT

steffen.moeller@uni-rostock.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Briefings in Bioinformatics online.

Exploration of databases and methods supporting drug repurposing: a comprehensive survey

Drug development involves a deep understanding of the mechanisms of action and possible side effects of each drug, and sometimes results in the identification of new and unexpected uses for drugs, termed as drug repurposing. Both in case of serendipitous observations and systematic mechanistic explorations, confirmation of new indications for a drug requires hypothesis building around relevant drug-related data, such as molecular targets involved, and patient and cellular responses. These datasets are available in public repositories, but apart from sifting through the sheer amount of data imposing computational bottleneck, a major challenge is the difficulty in selecting which databases to use from an increasingly large number of available databases. The database selection is made harder by the lack of an overview of the types of data offered in each database. In order to alleviate these problems and to guide the end user through the drug repurposing efforts, we provide here a survey of 102 of the most promising and drug-relevant databases reported to date. We summarize the target coverage and types of data available in each database and provide several examples of how multi-database exploration can facilitate drug repurposing.

Drug Repositioning: New Approaches and Future Prospects for Life-Debilitating Diseases and the COVID-19 Pandemic Outbreak

Traditionally, drug discovery utilises a de novo design approach, which requires high cost and many years of drug development before it reaches the market. Novel drug development does not always account for orphan diseases, which have low demand and hence low-profit margins for drug developers. Recently, drug repositioning has gained recognition as an alternative approach that explores new avenues for pre-existing commercially approved or rejected drugs to treat diseases aside from the intended ones. Drug repositioning results in lower overall developmental expenses and risk assessments, as the efficacy and safety of the original drug have already been well accessed and approved by regulatory authorities. The greatest advantage of drug repositioning is that it breathes new life into the novel, rare, orphan, and resistant diseases, such as Cushing's syndrome, HIV infection, and pandemic outbreaks such as COVID-19. Repositioning existing drugs such as Hydroxychloroquine, Remdesivir, Ivermectin and Baricitinib shows good potential for COVID-19 treatment. This can crucially aid in resolving outbreaks in urgent times of need. This review discusses the past success in drug repositioning, the current technological advancement in the field, drug repositioning for personalised medicine and the ongoing research on newly emerging drugs under consideration for the COVID-19 treatment.

Prediction of repurposed drugs for treating lung injury in COVID-19

Background: Coronavirus disease (COVID-19) is an infectious disease discovered in 2019 and currently in outbreak across the world. Lung injury with severe respiratory failure is the leading cause of death in COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, there still lacks efficient treatment for COVID-19 induced lung injury and acute respiratory failure. Methods: Inhibition of angiotensin-converting enzyme 2 (ACE2) caused by the spike protein of SARS-CoV-2 is the most plausible mechanism of lung injury in COVID-19. We performed drug repositioning analysis to identify drug candidates that reverse gene expression pattern in L1000 lung cell line HCC515 treated with ACE2 inhibitor. We confirmed these drug candidates by similar bioinformatics analysis using lung tissues from patients deceased from COVID-19. We further investigated deregulated genes and pathways related to lung injury, as well as the gene-pathway-drug candidate relationships. Results: We propose two candidate drugs, COL-3 (a chemically modified tetracycline) and CGP-60474 (a cyclin-dependent kinase inhibitor), for treating lung injuries in COVID-19. Further bioinformatics analysis shows that 12 significantly enriched pathways (P-value <0.05) overlap between HCC515 cells treated with ACE2 inhibitor and human COVID-19 patient lung tissues. These include signaling pathways known to be associated with lung injury such as TNF signaling, MAPK signaling and chemokine signaling pathways. All 12 pathways are targeted in COL-3 treated HCC515 cells, in which genes such as RHOA, RAC2, FAS, CDC42 have reduced expression. CGP-60474 shares 11 of 12 pathways with COL-3 and common target genes such as RHOA. It also uniquely targets other genes related to lung injury, such as CALR and MMP14. Conclusions: This study shows that ACE2 inhibition is likely part of the mechanisms leading to lung injury in COVID-19, and that compounds such as COL-3 and CGP-60474 have potential as repurposed drugs for its treatment.

Prediction of repurposed drugs for treating lung injury in COVID-19

Background: Coronavirus disease (COVID-19) is an infectious disease discovered in 2019 and currently in outbreak across the world. Lung injury with severe respiratory failure is the leading cause of death in COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, there still lacks efficient treatment for COVID-19 induced lung injury and acute respiratory failure. Methods: Inhibition of angiotensin-converting enzyme 2 (ACE2) caused by the spike protein of SARS-CoV-2 is the most plausible mechanism of lung injury in COVID-19. We performed drug repositioning analysis to identify drug candidates that reverse gene expression pattern in L1000 lung cell line HCC515 treated with ACE2 inhibitor. We confirmed these drug candidates by similar bioinformatics analysis using lung tissues from patients deceased from COVID-19. We further investigated deregulated genes and pathways related to lung injury, as well as the gene-pathway-drug candidate relationships. Results: We propose two candidate drugs, COL-3 (a chemically modified tetracycline) and CGP-60474 (a cyclin-dependent kinase inhibitor), for treating lung injuries in COVID-19. Further bioinformatics analysis shows that 12 significantly enriched pathways (P-value <0.05) overlap between HCC515 cells treated with ACE2 inhibitor and human COVID-19 patient lung tissues. These include signaling pathways known to be associated with lung injury such as TNF signaling, MAPK signaling and chemokine signaling pathways. All 12 pathways are targeted in COL-3 treated HCC515 cells, in which genes such as RHOA, RAC2, FAS, CDC42 have reduced expression. CGP-60474 shares 11 of 12 pathways with COL-3 and common target genes such as RHOA. It also uniquely targets other genes related to lung injury, such as CALR and MMP14. Conclusions: This study shows that ACE2 inhibition is likely part of the mechanisms leading to lung injury in COVID-19, and that compounds such as COL-3 and CGP-60474 have potential as repurposed drugs for its treatment.

Robust Sampling of Defective Pathways in Alzheimer's Disease. Implications in Drug Repositioning

We present the analysis of the defective genetic pathways of the Late-Onset Alzheimer’s Disease (LOAD) compared to the Mild Cognitive Impairment (MCI) and Healthy Controls (HC) using different sampling methodologies. These algorithms sample the uncertainty space that is intrinsic to any kind of highly underdetermined phenotype prediction problem, by looking for the minimum-scale signatures (header genes) corresponding to different random holdouts. The biological pathways can be identified performing posterior analysis of these signatures established via cross-validation holdouts and plugging the set of most frequently sampled genes into different ontological platforms. That way, the effect of helper genes, whose presence might be due to the high degree of under determinacy of these experiments and data noise, is reduced. Our results suggest that common pathways for Alzheimer’s disease and MCI are mainly related to viral mRNA translation, influenza viral RNA transcription and replication, gene expression, mitochondrial translation, and metabolism, with these results being highly consistent regardless of the comparative methods. The cross-validated predictive accuracies achieved for the LOAD and MCI discriminations were 84% and 81.5%, respectively. The difference between LOAD and MCI could not be clearly established (74% accuracy). The most discriminatory genes of the LOAD-MCI discrimination are associated with proteasome mediated degradation and G-protein signaling. Based on these findings we have also performed drug repositioning using Dr. Insight package, proposing the following different typologies of drugs: isoquinoline alkaloids, antitumor antibiotics, phosphoinositide 3-kinase PI3K, autophagy inhibitors, antagonists of the muscarinic acetylcholine receptor and histone deacetylase inhibitors. We believe that the potential clinical relevance of these findings should be further investigated and confirmed with other independent studies.

Prediction of repurposed drugs for treating lung injury in COVID-19

Coronavirus disease (COVID-19) is an infectious disease discovered in 2019 and currently in outbreak across the world. Lung injury with severe respiratory failure is the leading cause of death in COVID-19, brought by severe acute respiratory syndrome coronavirus 2 (SARS- CoV-2). However, there still lacks efficient treatment for COVID-19 induced lung injury and acute respiratory failure. Inhibition of Angiotensin-converting enzyme 2 (ACE2) caused by spike protein of SARS-CoV-2 is the most plausible mechanism of lung injury in COVID-19. We propose two candidate drugs, COL-3 (a chemically modified tetracycline) and CGP-60474 (a cyclin-dependent kinase inhibitor), for treating lung injuries in COVID-19, based on their abilities to reverse the gene expression patterns in HCC515 cells treated with ACE2 inhibitor and in human COVID-19 patient lung tissues. Further bioinformatics analysis shows that twelve significantly enriched pathways (P-value <0.05) overlap between HCC515 cells treated with ACE2 inhibitor and human COVID-19 patient lung tissues, including signaling pathways known to be associated with lung injury such as TNF signaling, MAPK signaling and Chemokine signaling pathways. All these twelve pathways are targeted in COL-3 treated HCC515 cells, in which genes such as RHOA, RAC2, FAS, CDC42 have reduced expression. CGP-60474 shares eleven of twelve pathways with COL-3 with common target genes such as RHOA. It also uniquely targets genes related to lung injury, such as CALR and MMP14. In summary, this study shows that ACE2 inhibition is likely part of the mechanisms leading to lung injury in COVID-19, and that compounds such as COL-3 and CGP-60474 have the potential as repurposed drugs for its treatment.

On the Role of Artificial Intelligence in Genomics to Enhance Precision Medicine

The complexity of orphan diseases, which are those that do not have an effective treatment, together with the high dimensionality of the genetic data used for their analysis and the high degree of uncertainty in the understanding of the mechanisms and genetic pathways which are involved in their development, motivate the use of advanced techniques of artificial intelligence and in-depth knowledge of molecular biology, which is crucial in order to find plausible solutions in drug design, including drug repositioning. Particularly, we show that the use of robust deep sampling methodologies of the altered genetics serves to obtain meaningful results and dramatically decreases the cost of research and development in drug design, influencing very positively the use of precision medicine and the outcomes in patients. The target-centric approach and the use of strong prior hypotheses that are not matched against reality (disease genetic data) are undoubtedly the cause of the high number of drug design failures and attrition rates. Sampling and prediction under uncertain conditions cannot be avoided in the development of precision medicine.

Analysis of defective pathways and drug repositioning in Multiple Sclerosis via machine learning approaches

BACKGROUND

Although some studies show that there could be a genetic predisposition to develop Multiple Sclerosis (MS), attempts to find genetic signatures related to MS diagnosis and development are extremely rare.

METHOD

We carried out a retrospective analysis of two different microarray datasets, using machine learning techniques to understand the defective pathways involved in this disease. We have modeled two data sets that are publicly accessible. The first was used to establish the list of most discriminatory genes; whereas, the second one was utilized for validation purposes.

RESULTS

The analysis provided a list of high discriminatory genes with predictive cross-validation accuracy higher than 95%, both in learning and in blind validation. The results were confirmed via the holdout sampler. The most discriminatory genes were related to the production of Hemoglobin. The biological processes involved were related to T-cell Receptor Signaling and co-stimulation, Interferon-Gamma Signaling and Antigen Processing and Presentation. Drug repositioning via CMAP methodologies highlighted the importance of Trichostatin A and other HDAC inhibitors.

CONCLUSIONS

The defective pathways suggest viral or bacterial infections as plausible mechanisms involved in MS development. The pathway analysis also confirmed coincidences with Epstein-Barr virus, Influenza A, Toxoplasmosis, Tuberculosis and Staphylococcus Aureus infections. Th17 Cell differentiation, and CD28 co-stimulation seemed to be crucial in the development of this disease. Furthermore, the additional knowledge provided by this analysis helps to identify new therapeutic targets.