A review of connectivity map and computational approaches in pharmacogenomics

Application of GWAS summary data and drug-induced gene expression profiles of neural progenitor cells in psychiatric drug prioritization analysis

Common psychiatric disorders constitute one of the most substantial healthcare burdens worldwide. However, drug development in psychiatry remains hampered partially due to the lack of approaches to estimating drugs that can simultaneously modulate the expression of a nontrivial fraction of disease susceptibility genes. We proposed a new drug prioritization strategy under the framework of our previously proposed phenotype-associated tissues estimation approach (DESE) by investigating the drugs' selective perturbation effect on disease susceptibility genes. Based on the genome-wide association study summary data and drug-induced gene expression profiles of neural progenitor cells, we applied this strategy to prioritize candidate drugs for schizophrenia, depression and bipolar I disorder and identified several known therapeutic drugs among the top-ranked drug candidates. Also, our results revealed that the disease susceptibility genes involved in the selective gene perturbation analysis were enriched with many biologically sensible function terms and interacted with known therapeutic drugs. Our results suggested that selective gene perturbation analysis could be a promising starting point to prioritize biologically sensible drug candidates under the "one drug, multiple targets" paradigm for the drug development of common psychiatric disorders.

DrugSK: A Stacked Ensemble Learning Framework for Predicting Drug Combinations of Multiple Diseases

Combination therapy is an important direction of continuous exploration in the field of medicine, with the core goals of improving treatment efficacy, reducing adverse reactions, and optimizing clinical outcomes. Machine learning technology holds great promise in improving the prediction of drug synergy combinations. However, most studies focus on single disease-oriented collaborative predictive models or involve excessive feature categories, making it challenging to predict the majority of new drugs. To address these challenges, the DrugSK comprehensive model was developed, which utilizes SMILES-BERT to extract structural information from 3492 drugs and trains on reactions from 48,756 drug combinations. DrugSK is an integrated learning model capable of predicting interactions among various drug categories. First, the primary learner is trained from the initial data set. Random forest, support vector machine, and XGboost model are selected as primary learners and logistic regression as secondary learners. A new data set is then "generated" to train level 2 learners, which can be thought of as a prediction for each model. Finally, the results are filtered using logistic regression. Furthermore, the combination of the new antibacterial drug Drafloxacin with other antibacterial agents was tested. The synergistic effect of Drafloxacin and Isavuconazonium in the fight against Candida albicans has been confirmed, providing enlightenment for the clinical treatment of skin infection. DrugSK's prediction is accurate in practical application and can also predict the probability of the outcome. In addition, the tendency of Drafloxacin and antifungal drugs to be synergistic was found. The development of DrugSK will provide a new blueprint for predicting drug combination synergies.

Computational drug discovery approaches identify mebendazole as a candidate treatment for autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is a rare genetic disorder characterised by numerous renal cysts, the progressive expansion of which can impact kidney function and lead eventually to renal failure. Tolvaptan is the only disease-modifying drug approved for the treatment of ADPKD, however its poor side effect and safety profile necessitates the need for the development of new therapeutics in this area. Using a combination of transcriptomic and machine learning computational drug discovery tools, we predicted that a number of existing drugs could have utility in the treatment of ADPKD, and subsequently validated several of these drug predictions in established models of disease. We determined that the anthelmintic mebendazole was a potent anti-cystic agent in human cellular and in vivo models of ADPKD, and is likely acting through the inhibition of microtubule polymerisation and protein kinase activity. These findings demonstrate the utility of combining computational approaches to identify and understand potential new treatments for traditionally underserved rare diseases.

Network-based drug repurposing for schizophrenia

Despite recent progress, the challenges in drug discovery for schizophrenia persist. However, computational drug repurposing has gained popularity as it leverages the wealth of expanding biomedical databases. Network analyses provide a comprehensive understanding of transcription factor (TF) regulatory effects through gene regulatory networks, which capture the interactions between TFs and target genes by integrating various lines of evidence. Using the PANDA algorithm, we examined the topological variances in TF-gene regulatory networks between individuals with schizophrenia and healthy controls. This algorithm incorporates binding motifs, protein interactions, and gene co-expression data. To identify these differences, we subtracted the edge weights of the healthy control network from those of the schizophrenia network. The resulting differential network was then analysed using the CLUEreg tool in the GRAND database. This tool employs differential network signatures to identify drugs that potentially target the gene signature associated with the disease. Our analysis utilised a large RNA-seq dataset comprising 532 post-mortem brain samples from the CommonMind project. We constructed co-expression gene regulatory networks for both schizophrenia cases and healthy control subjects, incorporating 15,831 genes and 413 overlapping TFs. Through drug repurposing, we identified 18 promising candidates for repurposing as potential treatments for schizophrenia. The analysis of TF-gene regulatory networks revealed that the TFs in schizophrenia predominantly regulate pathways associated with energy metabolism, immune response, cell adhesion, and thyroid hormone signalling. These pathways represent significant targets for therapeutic intervention. The identified drug repurposing candidates likely act through TF-targeted pathways. These promising candidates, particularly those with preclinical evidence such as rimonabant and kaempferol, warrant further investigation into their potential mechanisms of action and efficacy in alleviating the symptoms of schizophrenia.

Decoding cardiovascular risks: analyzing type 2 diabetes mellitus and ASCVD gene expression

Background

ASCVD is the primary cause of mortality in individuals with T2DM. A potential link between ASCVD and T2DM has been suggested, prompting further investigation.

Methods

We utilized linear and multivariate logistic regression, Wilcoxon test, and Spearman's correlation toanalyzethe interrelation between ASCVD and T2DM in NHANES data from 2001-2018.The Gene Expression Omnibus (GEO) database and Weighted Gene Co-expression Network Analysis (WGCNA) wereconducted to identify co-expression networks between ASCVD and T2DM. Hub genes were identified using LASSO regression analysis and further validated in two additional cohorts. Bioinformatics methods were employed for gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, along with the prediction of candidate small molecules.

Results

Our analysis of the NHANES dataset indicated a significant impact of blood glucose on lipid levels within diabetic cohort, suggesting that abnormal lipid metabolism is a critical factor in ASCVD development. Cross-phenotyping analysis revealed two pivotal genes, ABCC5 and WDR7, associated with both T2DM and ASCVD. Enrichment analyses demonstrated the intertwining of lipid metabolism in both conditions, encompassing adipocytokine signaling pathway, fatty acid degradation and metabolism, and the regulation of adipocyte lipolysis. Immune infiltration analysis underscored the involvement of immune processes in both diseases. Notably, RITA, ON-01910, doxercalciferol, and topiramate emerged as potential therapeutic agents for both T2DM and ASCVD, indicating their possible clinical significance.

Conclusion

Our findings pinpoint ABCC5 and WDR7 as new target genes between T2DM and ASCVD, with RITA, ON-01910, doxercalciferol, and topiramate highlighted as promising therapeutic agents.

Combinatorial prediction of therapeutic perturbations using causally-inspired neural networks

As an alternative to target-driven drug discovery, phenotype-driven approaches identify compounds that counteract the overall disease effects by analyzing phenotypic signatures. Our study introduces a novel approach to this field, aiming to expand the search space for new therapeutic agents. We introduce PDGrapher, a causally-inspired graph neural network model designed to predict arbitrary perturbagens - sets of therapeutic targets - capable of reversing disease effects. Unlike existing methods that learn responses to perturbations, PDGrapher solves the inverse problem, which is to infer the perturbagens necessary to achieve a specific response - i.e., directly predicting perturbagens by learning which perturbations elicit a desired response. Experiments across eight datasets of genetic and chemical perturbations show that PDGrapher successfully predicted effective perturbagens in up to 9% additional test samples and ranked therapeutic targets up to 35% higher than competing methods. A key innovation of PDGrapher is its direct prediction capability, which contrasts with the indirect, computationally intensive models traditionally used in phenotypedriven drug discovery that only predict changes in phenotypes due to perturbations. The direct approach enables PDGrapher to train up to 30 times faster, representing a significant leap in efficiency. Our results suggest that PDGrapher can advance phenotype-driven drug discovery, offering a fast and comprehensive approach to identifying therapeutically useful perturbations.

Lansoprazole inhibits the development of sessile serrated lesions by inducing G1 arrest via Skp2/p27 signaling pathway

BACKGROUND

Although the serrated-neoplasia pathway reportedly accounts for 15-30% of colorectal cancer (CRC), no studies on chemoprevention of sessile serrated lesions (SSLs) have been reported. We searched for effective compounds comprehensively from a large series of compounds by employing Connectivity Map (CMAP) analysis of SSL-specific gene expression profiles coupled with in vitro screening using SSL patient-derived organoids (PDOs), and validated their efficacy using a xenograft mouse model of SSL.

METHODS

We generated SSL-specific gene signatures based on DNA microarray data, and applied them to CMAP analysis with 1309 FDA-approved compounds to select candidate compounds. We evaluated their inhibitory effects on SSL-PDOs using a cell viability assay. SSL-PDOs were orthotopically transplanted into NOG mice for in vivo evaluation. The signal transduction pathway was evaluated by gene expression profile and protein expression analysis.

RESULTS

We identified 221 compounds by employing CMAP analysis of SSL-specific signatures, which should cancel the gene signatures, and narrowed them down to 17 compounds. Cell viability assay using SSL-PDOs identified lansoprazole as having the lowest IC50 value (47 µM) among 17 compounds. When SSL-PDO was orthotopically transplanted into murine intestinal tract, the tumor grew gradually. Administration of lansoprazole to mice inhibited the growth of SSL xenograft whereas the tumor in control mice treated with vehicle alone grew gradually over time. The Ki67 index in xenograft lesions from the lansoprazole group was significantly lower compared with the control group. Cell cycle analysis of SSL-PDOs treated with lansoprazole exhibited a significant increase in G1 phase cell population. Microarray and protein analysis revealed that lansoprazole downregulated Skp2 expression and upregulated p27 expression in SSL-PDOs.

CONCLUSIONS

Our data strongly suggest that lansoprazole is the most effective chemopreventive agent against SSL, and that lansoprazole induces G1 cell cycle arrest by downregulating Skp2 and upregulating p27 in SSL cells.

Interrogating colorectal cancer metastasis to liver: a search for clinically viable compounds and mechanistic insights in colorectal cancer Patient Derived Organoids

BACKGROUND

Approximately 20-50% of patients presenting with localized colorectal cancer progress to stage IV metastatic disease (mCRC) following initial treatment and this is a major prognostic determinant. Here, we have interrogated a heterogeneous set of primary colorectal cancer (CRC), liver CRC metastases and adjacent liver tissue to identify molecular determinants of the colon to liver spreading. Screening Food and Drug Administration (FDA) approved drugs for their ability to interfere with an identified colon to liver metastasis signature may help filling an unmet therapeutic need.

METHODS

RNA sequencing of primary colorectal cancer specimens vs adjacent liver tissue vs synchronous and asynchronous liver metastases. Pathways enrichment analyses. The Library of Integrated Network-based Cellular Signatures (LINCS)-based and Connectivity Map (CMAP)-mediated identification of FDA-approved compounds capable to interfere with a 22 gene signature from primary CRC and liver metastases. Testing the identified compounds on CRC-Patient Derived Organoid (PDO) cultures. Microscopy and Fluorescence Activated Cell Sorting (FACS) based analysis of the treated PDOs.

RESULTS

We have found that liver metastases acquire features of the adjacent liver tissue while partially losing those of the primary tumors they derived from. We have identified a 22-gene signature differentially expressed among primary tumors and metastases and validated in public databases. A pharmacogenomic screening for FDA-approved compounds capable of interfering with this signature has been performed. We have validated some of the identified representative compounds in CRC-Patient Derived Organoid cultures (PDOs) and found that pentoxyfilline and, to a minor extent, dexketoprofen and desloratadine, can variably interfere with number, size and viability of the CRC -PDOs in a patient-specific way. We explored the pentoxifylline mechanism of action and found that pentoxifylline treatment attenuated the 5-FU elicited increase of ALDHhigh cells by attenuating the IL-6 mediated STAT3 (tyr705) phosphorylation.

CONCLUSIONS

Pentoxifylline synergizes with 5-Fluorouracil (5-FU) in attenuating organoid formation. It does so by interfering with an IL-6-STAT3 axis leading to the emergence of chemoresistant ALDHhigh cell subpopulations in 5-FU treated PDOs. A larger cohort of CRC-PDOs will be required to validate and expand on the findings of this proof-of-concept study.

Identification of genetic mechanisms underlying lipid metabolism-mediated tumor immunity in head and neck squamous cell carcinoma

OBJECTIVE

To identify the genetic mechanisms underlying lipid metabolism-mediated tumor immunity in head and neck squamous carcinoma (HNSC).

MATERIALS AND METHODS

RNA sequencing data and clinical characteristics of HNSC patients were procured from The Cancer Genome Atlas (TCGA) database. Lipid metabolism-related genes were collected from KEGG and MSigDB databases. Immune cells and immune-related genes were obtained from the TISIDB database. The differentially expressed genes (DEGs) in HNSC were identified and weighted correlation network analysis (WGCNA) was performed to identify the significant gene modules. Lasso regression analysis was performed to identify hub genes. The differential gene expression pattern, diagnostic values, relationships with clinical features, prognostic values, relationships with tumor mutation burden (TMB), and signaling pathways involved, were each investigated.

RESULTS

One thousand six hundred sixty-eight DEGs were identified as dysregulated between HNSC tumor samples and healthy control head and neck samples. WGCNA analysis and Lasso regression analysis identified 8 hub genes, including 3 immune-related genes (PLA2G2D, TNFAIP8L2 and CYP27A1) and 5 lipid metabolism-related genes (FOXP3, IL21R, ITGAL, TRAF1 and WIPF1). Except CYP27A1, the other hub genes were upregulated in HNSC as compared with healthy control samples, and a low expression of these hub genes indicated a higher risk of death in HNSC. Except PLA2G2D, all other hub genes were significantly and negatively related with TMB in HNSC. The hub genes were implicated in several immune-related signaling pathways including T cell receptor signaling, Th17 cell differentiation, and natural killer (NK) cell mediated cytotoxicity.

CONCLUSION

Three immune genes (PLA2G2D, TNFAIP8L2, and CYP27A1) and immune-related pathways (T cell receptor signaling, Th17 cell differentiation, and natural killer (NK) cell mediated cytotoxicity) were predicted to play significant roles in the lipid metabolism-mediated tumor immunity in HNSC.

Identification of significant gene expression changes in multiple perturbation experiments using knockoffs

Large-scale multiple perturbation experiments have the potential to reveal a more detailed understanding of the molecular pathways that respond to genetic and environmental changes. A key question in these studies is which gene expression changes are important for the response to the perturbation. This problem is challenging because (i) the functional form of the nonlinear relationship between gene expression and the perturbation is unknown and (ii) identification of the most important genes is a high-dimensional variable selection problem. To deal with these challenges, we present here a method based on the model-X knockoffs framework and Deep Neural Networks to identify significant gene expression changes in multiple perturbation experiments. This approach makes no assumptions on the functional form of the dependence between the responses and the perturbations and it enjoys finite sample false discovery rate control for the selected set of important gene expression responses. We apply this approach to the Library of Integrated Network-Based Cellular Signature data sets which is a National Institutes of Health Common Fund program that catalogs how human cells globally respond to chemical, genetic and disease perturbations. We identified important genes whose expression is directly modulated in response to perturbation with anthracycline, vorinostat, trichostatin-a, geldanamycin and sirolimus. We compare the set of important genes that respond to these small molecules to identify co-responsive pathways. Identification of which genes respond to specific perturbation stressors can provide better understanding of the underlying mechanisms of disease and advance the identification of new drug targets.

Telomere-related genes as potential biomarkers to predict endometriosis and immune response: Development of a machine learning-based risk model

Introduction

Endometriosis (EM) is an aggressive, pleomorphic, and common gynecological disease. Its clinical presentation includes abnormal menstruation, dysmenorrhea, and infertility, which seriously affect the patient's quality of life. However, the pathogenesis underlying EM and associated regulatory genes are unknown.

Methods

Telomere-related genes (TRGs) were uploaded from TelNet. RNA-sequencing (RNA-seq) data of EM patients were obtained from three datasets (GSE5108, GSE23339, and GSE25628) in the GEO database, and a random forest approach was used to identify telomere signature genes and build nomogram prediction models. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Gene Set Enrichment Analysis were used to identify the pathways involved in the action of the signature genes. Finally, the CAMP database was used to screen drugs for potential use in EM treatment.

Results

Fifteen total genes were screened as EM–telomere differentially expressed genes. Further screening by machine learning obtained six genes as characteristic predictive of EM. Immuno-infiltration analysis of the telomeric genes showed that expressions including macrophages and natural killer cells were significantly higher in cluster A. Further enrichment analysis showed that the differential genes were mainly enriched in biological pathways like cell cycle and extracellular matrix. Finally, the Connective Map database was used to screen 11 potential drugs for EM treatment.

Discussion

TRGs play a crucial role in EM development, and are associated with immune infiltration and act on multiple pathways, including the cell cycle. Telomere signature genes can be valuable predictive markers for EM.

Advancing Targeted Protein Degradation via Multiomics Profiling and Artificial Intelligence

Only around 20% of the human proteome is considered to be druggable with small-molecule antagonists. This leaves some of the most compelling therapeutic targets outside the reach of ligand discovery. The concept of targeted protein degradation (TPD) promises to overcome some of these limitations. In brief, TPD is dependent on small molecules that induce the proximity between a protein of interest (POI) and an E3 ubiquitin ligase, causing ubiquitination and degradation of the POI. In this perspective, we want to reflect on current challenges in the field, and discuss how advances in multiomics profiling, artificial intelligence, and machine learning (AI/ML) will be vital in overcoming them. The presented roadmap is discussed in the context of small-molecule degraders but is equally applicable for other emerging proximity-inducing modalities.

Specific gene module pair-based target identification and drug discovery

Identification of the biological targets of a compound is of paramount importance for the exploration of the mechanism of action of drugs and for the development of novel drugs. A concept of the Connectivity Map (CMap) was previously proposed to connect genes, drugs, and disease states based on the common gene-expression signatures. For a new query compound, the CMap-based method can infer its potential targets by searching similar drugs with known targets (reference drugs) and measuring the similarities into their specific transcriptional responses between the query compound and those reference drugs. However, the available methods are often inefficient due to the requirement of the reference drugs as a medium to link the query agent and targets. Here, we developed a general procedure to extract target-induced consensus gene modules from the transcriptional profiles induced by the treatment of perturbagens of a target. A specific transcriptional gene module pair (GMP) was automatically identified for each target and could be used as a direct target signature. Based on the GMPs, we built the target network and identified some target gene clusters with similar biological mechanisms. Moreover, a gene module pair-based target identification (GMPTI) approach was proposed to predict novel compound-target interactions. Using this method, we have discovered novel inhibitors for three PI3K pathway proteins PI3Kα/β/δ, including PU-H71, alvespimycin, reversine, astemizole, raloxifene HCl, and tamoxifen.

Knowledge Mapping of Drug Repositioning's Theme and Development

Background

In recent years, the emergence of new diseases and resistance to known diseases have led to increasing demand for new drugs. By means of bibliometric analysis, this paper studied the relevant articles on drug repositioning in recent years and analyzed the current research foci and trends.

Methodology

The Web of Science database was searched to collect all relevant literature on drug repositioning from 2001 to 2022. These data were imported into CiteSpace and bibliometric online analysis platforms for bibliometric analysis. The processed data and visualized images predict the development trends in the research field.

Results

The quality and quantity of articles published after 2011 have improved significantly, with 45 of them cited more than 100 times. Articles posted by journals from different countries have high citation values. Authors from other institutions have also collaborated to analyze drug rediscovery. Keywords found in the literature include molecular docking (N=223), virtual screening (N=170), drug discovery (N=126), machine learning (N=125), and drug-target interaction (N=68); these words represent the core content of drug repositioning.

Conclusion

The key focus of drug research and development is related to the discovery of new indications for drugs. Researchers are starting to retarget drugs after analyzing online databases and clinical trials. More and more drugs are being targeted at other diseases to treat more patients, based on saving money and time. It is worth noting that researchers need more financial and technical support to complete drug development.

Integrated genomic analysis to identify druggable targets for pancreatic cancer

According to the National Comprehensive Cancer Network and the American Society of Clinical Oncology, the standard treatment for pancreatic cancer (PC) is gemcitabine and fluorouracil. Other chemotherapeutic agents have been widely combined. However, drug resistance remains a huge challenge, leading to the ineffectiveness of cancer therapy. Therefore, we are trying to discover new treatments for PC by utilizing genomic information to identify PC-associated genes as well as drug target genes for drug repurposing. Genomic information from a public database, the cBio Cancer Genomics Portal, was employed to retrieve the somatic mutation genes of PC. Five functional annotations were applied to prioritize the PC risk genes: Kyoto Encyclopedia of Genes and Genomes; biological process; knockout mouse; Gene List Automatically Derived For You; and Gene Expression Omnibus Dataset. DrugBank database was utilized to extract PC drug targets. To narrow down the most promising drugs for PC, CMap Touchstone analysis was applied. Finally, ClinicalTrials.gov and a literature review were used to screen the potential drugs under clinical and preclinical investigation. Here, we extracted 895 PC-associated genes according to the cBioPortal database and prioritized them by using five functional annotations; 318 genes were assigned as biological PC risk genes. Further, 216 genes were druggable according to the DrugBank database. CMap Touchstone analysis indicated 13 candidate drugs for PC. Among those 13 drugs, 8 drugs are in the clinical trials, 2 drugs were supported by the preclinical studies, and 3 drugs are with no evidence status for PC. Importantly, we found that midostaurin (targeted PRKA) and fulvestrant (targeted ESR1) are promising candidate drugs for PC treatment based on the genomic-driven drug repurposing pipelines. In short, integrated analysis using a genomic information database demonstrated the viability for drug repurposing. We proposed two drugs (midostaurin and fulvestrant) as promising drugs for PC.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

A risk model developed based on necroptosis to assess progression for ischemic cardiomyopathy and identify possible therapeutic drugs

Object: Ischemic cardiomyopathy (ICM), with high morbidity and mortality, is the most common cause of heart failure. Cardiovascular remodeling secondary to chronic myocardial ischemia is the main cause of its progression. A recently identified type of programmed cell death called necroptosis is crucial in the development of various cardiovascular diseases. However, the function role of necroptosis in cardiac remodeling of ICM has not been elucidated. Our study aimed to screen for genes associated with necroptosis and construct a risk score to assess the progression and evaluate the prognosis of ICM patients, and further to search for potentially therapeutic drugs. Methods: The gene expression profiling was obtained from the GEO database. LASSO regression analysis was used to construct necroptosis-related gene signatures associated with ICM progression and prognosis. TF-gene and miRNA-gene networks were constructed to identify the regulatory targets of potential necroptosis-related signature genes. Pathway alterations in patients with high necroptosis-related score (NRS) were analyzed by GO, KEGG, GSEA analysis, and immune cell infiltration was estimated by ImmuCellAI analysis. CMap analysis was performed to screen potential small molecule compounds targeting patients with high NRS. Independent risk analyses were performed using nomograms. Results: Six necroptosis-related signature genes (STAT4, TNFSF10, CHMP5, CHMP18, JAK1, and CFLAR) were used to define the NRS, with areas under the ROC curves of 0.833, 0.765, and 0.75 for training test, test set, and validation set, respectively. Transcription factors FOXC1 and hsa-miR-124-3p miRNA may be regulators of signature genes. Patients with higher NRS have pathway enriched in fibrosis and metabolism and elevated nTreg cells. AZD-7762 may be an effective drug to improve the prognosis of patients with high NRS. A feature-based nomogram was constructed from which patients could derive clinical benefit. Conclusion: Our results reveal 6 necroptosis gene signatures that can evaluate the progression and prognosis of ICM with high clinical value, and identify potential targets that could help improve cardiovascular remodeling.

Navigating Transcriptomic Connectivity Mapping Workflows to Link Chemicals with Bioactivities

Screening new compounds for potential bioactivities against cellular targets is vital for drug discovery and chemical safety. Transcriptomics offers an efficient approach for assessing global gene expression changes, but interpreting chemical mechanisms from these data is often challenging. Connectivity mapping is a potential data-driven avenue for linking chemicals to mechanisms based on the observation that many biological processes are associated with unique gene expression signatures (gene signatures). However, mining the effects of a chemical on gene signatures for biological mechanisms is challenging because transcriptomic data contain thousands of noisy genes. New connectivity mapping approaches seeking to distinguish signal from noise continue to be developed, spurred by the promise of discovering chemical mechanisms, new drugs, and disease targets from burgeoning transcriptomic data. Here, we analyze these approaches in terms of diverse transcriptomic technologies, public databases, gene signatures, pattern-matching algorithms, and statistical evaluation criteria. To navigate the complexity of connectivity mapping, we propose a harmonized scheme to coherently organize and compare published workflows. We first standardize concepts underlying transcriptomic profiles and gene signatures based on various transcriptomic technologies such as microarrays, RNA-Seq, and L1000 and discuss the widely used data sources such as Gene Expression Omnibus, ArrayExpress, and MSigDB. Next, we generalize connectivity mapping as a pattern-matching task for finding similarity between a query (e.g., transcriptomic profile for new chemical) and a reference (e.g., gene signature of known target). Published pattern-matching approaches fall into two main categories: vector-based use metrics like correlation, Jaccard index, etc., and aggregation-based use parametric and nonparametric statistics (e.g., gene set enrichment analysis). The statistical methods for evaluating the performance of different approaches are described, along with comparisons reported in the literature on benchmark transcriptomic data sets. Lastly, we review connectivity mapping applications in toxicology and offer guidance on evaluating chemical-induced toxicity with concentration-response transcriptomic data. In addition to serving as a high-level guide and tutorial for understanding and implementing connectivity mapping workflows, we hope this review will stimulate new algorithms for evaluating chemical safety and drug discovery using transcriptomic data.

Estimating intraclonal heterogeneity and subpopulation changes from bulk expression profiles in CMap

Premnas is a computational framework that provides a new perspective to interpret perturbational data in LINC L1000 CMap by learning an ad hoc subpopulation representation from scRNA-seq and performing the digital cytometry to estimate the abundance of undetermined subpopulations.

The connectivity among signatures upon perturbations curated in the CMap library provides a valuable resource for understanding therapeutic pathways and biological processes associated with the drugs and diseases. However, because of the nature of bulk-level expression profiling by the L1000 assay, intraclonal heterogeneity and subpopulation compositional change that could contribute to the responses to perturbations are largely neglected, hampering the interpretability and reproducibility of the connections. In this work, we proposed a computational framework, Premnas, to estimate the abundance of undetermined subpopulations from L1000 profiles in CMap directly according to an ad hoc subpopulation representation learned from a well-normalized batch of single-cell RNA-seq datasets by the archetypal analysis. By recovering the information of subpopulation changes upon perturbation, the potentials of drug-resistant/susceptible subpopulations with CMap L1000 were further explored and examined. The proposed framework enables a new perspective to understand the connectivity among cellular signatures and expands the scope of the CMAP and other similar perturbation datasets limited by the bulk profiling technology.

Signatures of Co-Deregulated Genes and Their Transcriptional Regulators in Lung Cancer

Despite the significant progress made towards comprehending the deregulated signatures in lung cancer, these vary from study to study. We reanalyzed 25 studies from the Gene Expression Omnibus (GEO) to detect and annotate co-deregulated signatures in lung cancer and in single-gene or single-drug perturbation experiments. We aimed to decipher the networks that these co-deregulated genes (co-DEGs) form along with their upstream regulators. Differential expression and upstream regulators were computed using Characteristic Direction and Systems Biology tools, including GEO2Enrichr and X2K. Co-deregulated gene expression profiles were further validated across different molecular and immune subtypes in lung adenocarcinoma (TCGA-LUAD) and lung adenocarcinoma (TCGA-LUSC) datasets, as well as using immunohistochemistry data from the Human Protein Atlas, before being subjected to subsequent GO and KEGG enrichment analysis. The functional alterations of the co-upregulated genes in lung cancer were mostly related to immune response regulating the cell surface signaling pathway, in contrast to the co-downregulated genes, which were related to S-nitrosylation. Networks of hub proteins across the co-DEGs consisted of overlapping TFs (SOX2, MYC, KAT2A) and kinases (MAPK14, CSNK2A1 and CDKs). Furthermore, using Connectivity Map we highlighted putative repurposing drugs, including valproic acid, betonicine and astemizole. Similarly, we analyzed the co-DEG signatures in single-gene and single-drug perturbation experiments in lung cancer cell lines. In summary, we identified critical co-DEGs in lung cancer providing an innovative framework for their potential use in developing personalized therapeutic strategies.

Optimization of cancer immunotherapy through pyroptosis: A pyroptosis-related signature predicts survival benefit and potential synergy for immunotherapy in glioma

Background

Pyroptosis is a critical type of programmed cell death that is strongly associated with the regulation of tumor and immune cell functions. However, the role of pyroptosis in tumor progression and remodeling of the tumor microenvironment in gliomas has not been extensively studied. Thus, in this study, we aimed to establish a comprehensive pyroptosis-related signature and uncover its potential clinical application in gliomas.

Methods

The TCGA glioma cohort was obtained and divided into training and internal validation cohorts, while the CGGA glioma cohort was used as an external validation cohort. Unsupervised consensus clustering was performed to identify pyroptosis-related expression patterns. A Cox regression analysis was performed to establish a pyroptosis-related risk signature. Real-time quantitative PCR was performed to analyze the expression of signature genes in glioma tissues. Immune infiltration was analyzed and validated by immunohistochemical staining. The expression patterns of signature genes in different cell types were analyzed using single-cell RNA sequencing data. Finally, therapeutic responses to chemotherapy, immunotherapy, and potential small-molecule inhibitors were investigated.

Results

Patients with glioma were stratified into clusters 1 and 2 based on the expression patterns of pyroptosis-related genes. Cluster 2 showed a longer overall (P<0.001) and progression-free survival time (P<0.001) than Cluster 1. CD8+ T cell enrichment was observed in Cluster 1. A pyroptosis-related risk signature (PRRS) was then established. The high PRRS group showed a significantly poorer prognosis than the low PRRS group in the training cohort (P<0.001), with validation in the internal and external validation cohorts. Immunohistochemical staining demonstrated that CD8+ T cells were enriched in high PRRS glioma tissues. PRRS genes also showed cell-specific expression in tumor and immune cells. Moreover, the high PRRS risk group showed higher temozolomide sensitivity and increased response to anti-PD1 treatment in a glioblastoma immunotherapy cohort. Finally, Bcl-2 inhibitors were screened as candidates for adjunct immunotherapy of gliomas.

Conclusion

The pyroptosis-related signature established in this study can be used to reliably predict clinical outcomes and immunotherapy responses in glioma patients. The correlation between the pyroptosis signature and the tumor immune microenvironment may be used to further guide the sensitization of glioma patients to immunotherapy.

Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) are popular approaches to test for association between imputed gene expression levels and traits of interest. Here, we propose an integrative method PUMICE (Prediction Using Models Informed by Chromatin conformations and Epigenomics) to integrate 3D genomic and epigenomic data with expression quantitative trait loci (eQTL) to more accurately predict gene expressions. PUMICE helps define and prioritize regions that harbor cis-regulatory variants, which outperforms competing methods. We further describe an extension to our method PUMICE +, which jointly combines TWAS results from single- and multi-tissue models. Across 79 traits, PUMICE + identifies 22% more independent novel genes and increases median chi-square statistics values at known loci by 35% compared to the second-best method, as well as achieves the narrowest credible interval size. Lastly, we perform computational drug repurposing and confirm that PUMICE + outperforms other TWAS methods.

Connectivity Mapping Using a Novel sv2a Loss-of-Function Zebrafish Epilepsy Model as a Powerful Strategy for Anti-epileptic Drug Discovery

Synaptic vesicle glycoprotein 2A (SV2A) regulates action potential-dependent neurotransmitter release and is commonly known as the primary binding site of an approved anti-epileptic drug, levetiracetam. Although several rodent knockout models have demonstrated the importance of SV2A for functional neurotransmission, its precise physiological function and role in epilepsy pathophysiology remains to be elucidated. Here, we present a novel sv2a knockout model in zebrafish, a vertebrate with complementary advantages to rodents. We demonstrated that 6 days post fertilization homozygous sv2a^–/– mutant zebrafish larvae, but not sv2a^+/– and sv2a^+/+ larvae, displayed locomotor hyperactivity and spontaneous epileptiform discharges, however, no major brain malformations could be observed. A partial rescue of this epileptiform brain activity could be observed after treatment with two commonly used anti-epileptic drugs, valproic acid and, surprisingly, levetiracetam. This observation indicated that additional targets, besides Sv2a, maybe are involved in the protective effects of levetiracetam against epileptic seizures. Furthermore, a transcriptome analysis provided insights into the neuropathological processes underlying the observed epileptic phenotype. While gene expression profiling revealed only one differentially expressed gene (DEG) between wildtype and sv2a^+/– larvae, there were 4386 and 3535 DEGs between wildtype and sv2a^–/–, and sv2a^+/– and sv2a^–/– larvae, respectively. Pathway and gene ontology (GO) enrichment analysis between wildtype and sv2a^–/– larvae revealed several pathways and GO terms enriched amongst up- and down-regulated genes, including MAPK signaling, synaptic vesicle cycle, and extracellular matrix organization, all known to be involved in epileptogenesis and epilepsy. Importantly, we used the Connectivity map database to identify compounds with opposing gene signatures compared to the one observed in sv2a^–/– larvae, to finally rescue the epileptic phenotype. Two out of three selected compounds rescued electrographic discharges in sv2a^–/– larvae, while negative controls did not. Taken together, our results demonstrate that sv2a deficiency leads to increased seizure vulnerability and provide valuable insight into the functional importance of sv2a in the brain in general. Furthermore, we provided evidence that the concept of connectivity mapping represents an attractive and powerful approach in the discovery of novel compounds against epilepsy.

In silico Drug Screening Approach Using L1000-Based Connectivity Map and Its Application to COVID-19

Conventional drug screening methods search for a limited number of small molecules that directly interact with the target protein. This process can be slow, cumbersome and has driven the need for developing new drug screening approaches to counter rapidly emerging diseases such as COVID-19. We propose a pipeline for drug repurposing combining in silico drug candidate identification followed by in vitro characterization of these candidates. We first identified a gene target of interest, the entry receptor for the SARS-CoV-2 virus, angiotensin converting enzyme 2 (ACE2). Next, we employed a gene expression profile database, L1000-based Connectivity Map to query gene expression patterns in lung epithelial cells, which act as the primary site of SARS-CoV-2 infection. Using gene expression profiles from 5 different lung epithelial cell lines, we computationally identified 17 small molecules that were predicted to decrease ACE2 expression. We further performed a streamlined validation in the normal human epithelial cell line BEAS-2B to demonstrate that these compounds can indeed decrease ACE2 surface expression and to profile cell health and viability upon drug treatment. This proposed pipeline combining in silico drug compound identification and in vitro expression and viability characterization in relevant cell types can aid in the repurposing of FDA-approved drugs to combat rapidly emerging diseases.

Computational Methods for Drug Repurposing

The wealth of knowledge and multi-omics data available in drug research has allowed the rise of several computational methods in the drug discovery field, resulting in a novel and exciting strategy called drug repurposing. Drug repurposing consists in finding new applications for existing drugs. Numerous computational methods perform a high-level integration of different knowledge sources to facilitate the discovery of unknown mechanisms. In this chapter, we present a survey of data resources and computational tools available for drug repositioning.

Integrated Analysis of MATH-Based Subtypes Reveals a Novel Screening Strategy for Early-Stage Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is a frequently diagnosed cancer type, and many patients have already reached an advanced stage when diagnosed. Thus, it is crucial to develop a novel and efficient approach to diagnose and classify lung adenocarcinoma at an early stage. In our study, we combined in silico analysis and machine learning to develop a new five-gene–based diagnosis strategy, which was further verified in independent cohorts and in vitro experiments. Considering the heterogeneity in cancer, we used the MATH (mutant-allele tumor heterogeneity) algorithm to divide patients with early-stage LUAD into two groups (C1 and C2). Specifically, patients in C2 had lower intratumor heterogeneity and higher abundance of immune cells (including B cell, CD4 T cell, CD8 T cell, macrophage, dendritic cell, and neutrophil). In addition, patients in C2 had a higher likelihood of immunotherapy response and overall survival advantage than patients in C1. Combined drug sensitivity analysis (CTRP/PRISM/CMap/GDSC) revealed that BI-2536 might serve as a new therapeutic compound for patients in C1. In order to realize the application value of our study, we constructed the classifier (to classify early-stage LUAD patients into C1 or C2 groups) with multiple machine learning and bioinformatic analyses. The 21-gene–based classification model showed high accuracy and strong generalization ability, and it was verified in four independent validation cohorts. In summary, our research provided a new strategy for clinicians to make a quick preliminary assisting diagnosis of early-stage LUAD and make patient classification at the intratumor heterogeneity level. All data, codes, and study processes have been deposited to Github and are available online.

Applicability of Anticancer Drugs for the Triple-Negative Breast Cancer Based on Homologous Recombination Repair Deficiency

Triple-negative breast cancer (TNBC) is a highly aggressive disease with historically poor outcomes, primarily due to the lack of effective targeted therapies. Here, we established a drug sensitivity prediction model based on the homologous recombination deficiency (HRD) using 83 TNBC patients from TCGA. Through analyzing the effect of HRD status on response efficacy of anticancer drugs and elucidating its related mechanisms of action, we found rucaparib (PARP inhibitor) and doxorubicin (anthracycline) sensitive in HR-deficient patients, while paclitaxel sensitive in the HR-proficient. Further, we identified a HRD signature based on gene expression data and constructed a transcriptomic HRD score, for analyzing the functional association between anticancer drug perturbation and HRD. The results revealed that CHIR99021 (GSK3 inhibitor) and doxorubicin have similar expression perturbation patterns with HRD, and talazoparib (PARP inhibitor) could kill tumor cells by reversing the HRD activity. Genomic characteristics indicated that doxorubicin inhibited tumor cells growth by hindering the process of DNA damage repair, while the resistance of cisplatin was related to the activation of angiogenesis and epithelial-mesenchymal transition. The negative correlation of HRD signature score could interpret the association of doxorubicin pIC50 with worse chemotherapy response and shorter survival of TNBC patients. In summary, these findings explain the applicability of anticancer drugs in TNBC and underscore the importance of HRD in promoting personalized treatment development.

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.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Transcription Factor Activation Profiles (TFAP) identify compounds promoting differentiation of Acute Myeloid Leukemia cell lines

Repurposing of drugs for new therapeutic use has received considerable attention for its potential to limit time and cost of drug development. Here we present a new strategy to identify chemicals that are likely to promote a desired phenotype. We used data from the Connectivity Map (CMap) to produce a ranked list of drugs according to their potential to activate transcription factors that mediate myeloid differentiation of leukemic progenitor cells. To validate our strategy, we tested the in vitro differentiation potential of candidate compounds using the HL-60 human cell line as a myeloid differentiation model. Ten out of 22 compounds, which were ranked high in the inferred list, were confirmed to promote significant differentiation of HL-60. These compounds may be considered candidate for differentiation therapy. The method that we have developed is versatile and it can be adapted to different drug repurposing projects.

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.

Side effect prediction based on drug-induced gene expression profiles and random forest with iterative feature selection

One in every ten drug candidates fail in clinical trials mainly due to efficacy and safety related issues, despite in-depth preclinical testing. Even some of the approved drugs such as chemotherapeutics are notorious for their side effects that are burdensome on patients. In order to pave the way for new therapeutics with more tolerable side effects, the mechanisms underlying side effects need to be fully elucidated. In this work, we addressed the common side effects of chemotherapeutics, namely alopecia, diarrhea and edema. A strategy based on Random Forest algorithm unveiled an expression signature involving 40 genes that predicted these side effects with an accuracy of 89%. We further characterized the resulting signature and its association with the side effects using functional enrichment analysis and protein-protein interaction networks. This work contributes to the ongoing efforts in drug development for early identification of side effects to use the resources more effectively.

Prognostic biomarkers for predicting papillary thyroid carcinoma patients at high risk using nine genes of apoptotic pathway

Aberrant expressions of apoptotic genes have been associated with papillary thyroid carcinoma (PTC) in the past, however, their prognostic role and utility as biomarkers remains poorly understood. In this study, we analysed 505 PTC patients by employing Cox-PH regression techniques, prognostic index models and machine learning methods to elucidate the relationship between overall survival (OS) of PTC patients and 165 apoptosis related genes. It was observed that nine genes (ANXA1, TGFBR3, CLU, PSEN1, TNFRSF12A, GPX4, TIMP3, LEF1, BNIP3L) showed significant association with OS of PTC patients. Five out of nine genes were found to be positively correlated with OS of the patients, while the remaining four genes were negatively correlated. These genes were used for developing risk prediction models, which can be utilized to classify patients with a higher risk of death from the patients which have a good prognosis. Our voting-based model achieved highest performance (HR = 41.59, p = 3.36x10^-4, C = 0.84, logrank-p = 3.8x10^-8). The performance of voting-based model improved significantly when we used the age of patients with prognostic biomarker genes and achieved HR = 57.04 with p = 10⁻⁴ (C = 0.88, logrank-p = 1.44x10^-9). We also developed classification models that can classify high risk patients (survival ≤ 6 years) and low risk patients (survival > 6 years). Our best model achieved AUROC of 0.92. Further, the expression pattern of the prognostic genes was verified at mRNA level, which showed their differential expression between normal and PTC samples. Also, the immunostaining results from HPA validated these findings. Since these genes can also be used as potential therapeutic targets in PTC, we also identified potential drug molecules which could modulate their expression profile. The study briefly revealed the key prognostic biomarker genes in the apoptotic pathway whose altered expression is associated with PTC progression and aggressiveness. In addition to this, risk assessment models proposed here can help in efficient management of PTC patients.

Drug target inference by mining transcriptional data using a novel graph convolutional network framework

A fundamental challenge that arises in biomedicine is the need to characterize compounds in a relevant cellular context in order to reveal potential on-target or off-target effects. Recently, the fast accumulation of gene transcriptional profiling data provides us an unprecedented opportunity to explore the protein targets of chemical compounds from the perspective of cell transcriptomics and RNA biology. Here, we propose a novel Siamese spectral-based graph convolutional network (SSGCN) model for inferring the protein targets of chemical compounds from gene transcriptional profiles. Although the gene signature of a compound perturbation only provides indirect clues of the interacting targets, and the biological networks under different experiment conditions further complicate the situation, the SSGCN model was successfully trained to learn from known compound-target pairs by uncovering the hidden correlations between compound perturbation profiles and gene knockdown profiles. On a benchmark set and a large time-split validation dataset, the model achieved higher target inference accuracy as compared to previous methods such as Connectivity Map. Further experimental validations of prediction results highlight the practical usefulness of SSGCN in either inferring the interacting targets of compound, or reversely, in finding novel inhibitors of a given target of interest.

Investigating Pathogenetic Mechanisms of Alzheimer's Disease by Systems Biology Approaches for Drug Discovery

Alzheimer's disease (AD) is the most common cause of dementia, characterized by progressive cognitive decline and neurodegenerative disorder. Abnormal aggregations of intracellular neurofibrillary tangles (NFTs) and unusual accumulations of extracellular amyloid-β (Aβ) peptides are two important pathological features in AD brains. However, in spite of large-scale clinical studies and computational simulations, the molecular mechanisms of AD development and progression are still unclear. In this study, we divided all of the samples into two groups: early stage (Braak score I-III) and later stage (Braak score IV-VI). By big database mining, the candidate genetic and epigenetic networks (GEN) have been constructed. In order to find out the real GENs for two stages of AD, we performed systems identification and system order detection scheme to prune false positives with the help of corresponding microarray data. Applying the principal network projection (PNP) method, core GENs were extracted from real GENs based on the projection values. By the annotation of KEGG pathway, we could obtain core pathways from core GENs and investigate pathogenetic mechanisms for the early and later stage of AD, respectively. Consequently, according to pathogenetic mechanisms, several potential biomarkers are identified as drug targets for multiple-molecule drug design in the treatment of AD.

Advancing the use of genome-wide association studies for drug repurposing

Genome-wide association studies (GWAS) have revealed important biological insights into complex diseases, which are broadly expected to lead to the identification of new drug targets and opportunities for treatment. Drug development, however, remains hampered by the time taken and costs expended to achieve regulatory approval, leading many clinicians and researchers to consider alternative paths to more immediate clinical outcomes. In this Review, we explore approaches that leverage common variant genetics to identify opportunities for repurposing existing drugs, also known as drug repositioning. These approaches include the identification of compounds by linking individual loci to genes and pathways that can be pharmacologically modulated, transcriptome-wide association studies, gene-set association, causal inference by Mendelian randomization, and polygenic scoring.

Single-Cell RNA-Seq of Cisplatin-Treated Adult Stria Vascularis Identifies Cell Type-Specific Regulatory Networks and Novel Therapeutic Gene Targets

The endocochlear potential (EP) generated by the stria vascularis (SV) is necessary for hair cell mechanotransduction in the mammalian cochlea. We sought to create a model of EP dysfunction for the purposes of transcriptional analysis and treatment testing. By administering a single dose of cisplatin, a commonly prescribed cancer treatment drug with ototoxic side effects, to the adult mouse, we acutely disrupt EP generation. By combining these data with single cell RNA-sequencing findings, we identify transcriptional changes induced by cisplatin exposure, and by extension transcriptional changes accompanying EP reduction, in the major cell types of the SV. We use these data to identify gene regulatory networks unique to cisplatin treated SV, as well as the differentially expressed and druggable gene targets within those networks. Our results reconstruct transcriptional responses that occur in gene expression on the cellular level while identifying possible targets for interventions not only in cisplatin ototoxicity but also in EP dysfunction.

Evaluation of connectivity map shows limited reproducibility in drug repositioning

The Connectivity Map (CMap) is a popular resource designed for data-driven drug repositioning using a large transcriptomic compendium. However, evaluations of its performance are limited. We used two iterations of CMap (CMap 1 and 2) to assess their comparability and reliability. We queried CMap 2 with CMap 1-derived signatures, expecting CMap 2 would highly prioritize the queried compounds; the success rate was 17%. Analysis of previously published prioritizations yielded similar results. Low recall is caused by low differential expression (DE) reproducibility both between CMaps and within each CMap. DE strength was predictive of reproducibility, and is influenced by compound concentration and cell-line responsiveness. Reproducibility of CMap 2 sample expression levels was also lower than expected. We attempted to identify the "better" CMap by comparison with a third dataset, but they were mutually discordant. Our findings have implications for CMap usage and we suggest steps for investigators to limit false positives.

Reproducibility and robustness of high-throughput S1500+ transcriptomics on primary rat hepatocytes for chemical-induced hepatotoxicity assessment

•

TempO-Seq assays of rat hepatocytes in collagen sandwich are highly reproducible.

•

Gene expression analysis shows S1500+ is representative of the whole transcriptome.

•

Connectivity mapping shows consistency between TempO-Seq and Affymetrix data.

•

Gene set enrichment shows consistency between S1500+ and the whole transcriptome.

•

Gene set enrichment using hallmark gene sets informs hepatotoxicity.

Cell-based in vitro models coupled with high-throughput transcriptomics (HTTr) are increasingly utilized as alternative methods to animal-based toxicity testing. Here, using a panel of 14 chemicals with different risks of human drug-induced liver injury (DILI) and two dosing concentrations, we evaluated an HTTr platform comprised of collagen sandwich primary rat hepatocyte culture and the TempO-Seq surrogate S1500+ (ST) assay. First, the HTTr platform was found to exhibit high reproducibility between technical and biological replicates (r greater than 0.85). Connectivity mapping analysis further demonstrated a high level of inter-platform reproducibility between TempO-Seq data and Affymetrix GeneChip data from the Open TG-GATES project. Second, the TempO-Seq ST assay was shown to be a robust surrogate to the whole transcriptome (WT) assay in capturing chemical-induced changes in gene expression, as evident from correlation analysis, PCA and unsupervised hierarchical clustering. Gene set enrichment analysis (GSEA) using the Hallmark gene set collection also demonstrated consistency in enrichment scores between ST and WT assays. Lastly, unsupervised hierarchical clustering of hallmark enrichment scores broadly divided the samples into hepatotoxic, intermediate, and non-hepatotoxic groups. Xenobiotic metabolism, bile acid metabolism, apoptosis, p53 pathway, and coagulation were found to be the key hallmarks driving the clustering. Taken together, our results established the reproducibility and performance of collagen sandwich culture in combination with TempO-Seq S1500+ assay, and demonstrated the utility of GSEA using the hallmark gene set collection to identify potential hepatotoxicants for further validation.

The Complex Structure of the Pharmacological Drug-Disease Network

The complexity of drug–disease interactions is a process that has been explained in terms of the need for new drugs and the increasing cost of drug development, among other factors. Over the last years, diverse approaches have been explored to understand drug–disease relationships. Here, we construct a bipartite graph in terms of active ingredients and diseases based on thoroughly classified data from a recognized pharmacological website. We find that the connectivities between drugs (outgoing links) and diseases (incoming links) follow approximately a stretched-exponential function with different fitting parameters; for drugs, it is between exponential and power law functions, while for diseases, the behavior is purely exponential. The network projections, onto either drugs or diseases, reveal that the co-ocurrence of drugs (diseases) in common target diseases (drugs) lead to the appearance of connected components, which varies as the threshold number of common target diseases (drugs) is increased. The corresponding projections built from randomized versions of the original bipartite networks are considered to evaluate the differences. The heterogeneity of association at group level between active ingredients and diseases is evaluated in terms of the Shannon entropy and algorithmic complexity, revealing that higher levels of diversity are present for diseases compared to drugs. Finally, the robustness of the original bipartite network is evaluated in terms of most-connected nodes removal (direct attack) and random removal (random failures).

Predicting mechanism of action of novel compounds using compound structure and transcriptomic signature coembedding

MOTIVATION

Identifying mechanism of actions (MoA) of novel compounds is crucial in drug discovery. Careful understanding of MoA can avoid potential side effects of drug candidates. Efforts have been made to identify MoA using the transcriptomic signatures induced by compounds. However, these approaches fail to reveal MoAs in the absence of actual compound signatures.

RESULTS

We present MoAble, which predicts MoAs without requiring compound signatures. We train a deep learning-based coembedding model to map compound signatures and compound structure into the same embedding space. The model generates low-dimensional compound signature representation from the compound structures. To predict MoAs, pathway enrichment analysis is performed based on the connectivity between embedding vectors of compounds and those of genetic perturbation. Results show that MoAble is comparable to the methods that use actual compound signatures. We demonstrate that MoAble can be used to reveal MoAs of novel compounds without measuring compound signatures with the same prediction accuracy as that with measuring them.

AVAILABILITY AND IMPLEMENTATION

MoAble is available at https://github.com/dmis-lab/moable.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Prognostic Biomarker-Based Identification of Drugs for Managing the Treatment of Endometrial Cancer

INTRODUCTION

Uterine corpus endometrial carcinoma (UCEC) causes thousands of deaths per year. To improve the overall survival of patients with UCEC, there is a need to identify prognostic biomarkers and potential drugs.

OBJECTIVES

The aim of this study was twofold: the identification of prognostic gene signatures from expression profiles of pattern recognition receptor (PRR) genes and identification of the most effective existing drugs using the prognostic gene signature.

METHODS

This study was based on the expression profile of PRR genes of 541 patients with UCEC obtained from The Cancer Genome Atlas. Key prognostic signatures were identified using various approaches, including survival analysis, network, and clustering. Hub genes were identified by constructing a co-expression network. Representative genes were identified using k-means and k-medoids-based clustering. Univariate Cox proportional hazard (PH) analysis was used to identify survival-associated genes. 'cmap2' was used to identify potential drugs that can suppress/enhance the expression of prognostic genes.

RESULTS

Models were developed using hub genes and achieved a maximum hazard ratio (HR) of 1.37 (p = 0.294). Then, a clustering-based model was developed using seven genes (HR 9.14; p = 1.49 × 10^-12). Finally, a nine gene-based risk stratification model was developed (CLEC1B, CLEC3A, IRF7, CTSB, FCN1, RIPK2, NLRP10, NLRP9, and SARM1) and achieved HR 10.70; p = 1.1 × 10^-12. The performance of this model improved significantly in combination with the clinical stage and achieved HR 15.23; p = 2.21 × 10^-7. We also developed a model for predicting high-risk patients (survival ≤ 4.3 years) and achieved an area under the receiver operating characteristic curve (AUROC) of 0.86.

CONCLUSION

We identified potential immunotherapeutic agents based on prognostic gene signature: hexamethonium bromide and isoflupredone. Several novel candidate drugs were suggested, including human interferon-α-2b, paclitaxel, imiquimod, MESO-DAP1, and mifamurtide. These biomolecules and repurposed drugs may be utilised for prognosis and treatment for better survival.

Transcriptional drug repositioning and cheminformatics approach for differentiation therapy of leukaemia cells

Differentiation therapy is attracting increasing interest in cancer as it can be more specific than conventional chemotherapy approaches, and it has offered new treatment options for some cancer types, such as treating acute promyelocytic leukaemia (APL) by retinoic acid. However, there is a pressing need to identify additional molecules which act in this way, both in leukaemia and other cancer types. In this work, we hence developed a novel transcriptional drug repositioning approach, based on both bioinformatics and cheminformatics components, that enables selecting such compounds in a more informed manner. We have validated the approach for leukaemia cells, and retrospectively retinoic acid was successfully identified using our method. Prospectively, the anti-parasitic compound fenbendazole was tested in leukaemia cells, and we were able to show that it can induce the differentiation of leukaemia cells to granulocytes in low concentrations of 0.1 μM and within as short a time period as 3 days. This work hence provides a systematic and validated approach for identifying small molecules for differentiation therapy in cancer.

Transcriptomics-based drug repositioning pipeline identifies therapeutic candidates for COVID-19

The novel SARS-CoV-2 virus emerged in December 2019 and has few effective treatments. We applied a computational drug repositioning pipeline to SARS-CoV-2 differential gene expression signatures derived from publicly available data. We utilized three independent published studies to acquire or generate lists of differentially expressed genes between control and SARS-CoV-2-infected samples. Using a rank-based pattern matching strategy based on the Kolmogorov-Smirnov Statistic, the signatures were queried against drug profiles from Connectivity Map (CMap). We validated 16 of our top predicted hits in live SARS-CoV-2 antiviral assays in either Calu-3 or 293T-ACE2 cells. Validation experiments in human cell lines showed that 11 of the 16 compounds tested to date (including clofazimine, haloperidol and others) had measurable antiviral activity against SARS-CoV-2. These initial results are encouraging as we continue to work towards a further analysis of these predicted drugs as potential therapeutics for the treatment of COVID-19.

Connectivity Map Analysis Indicates PI3K/Akt/mTOR Inhibitors as Potential Anti-Hypoxia Drugs in Neuroblastoma

Neuroblastoma (NB) is one of the deadliest pediatric cancers, accounting for 15% of deaths in childhood. Hypoxia is a condition of low oxygen tension occurring in solid tumors and has an unfavorable prognostic factor for NB. In the present study, we aimed to identify novel promising drugs for NB treatment. Connectivity Map (CMap), an online resource for drug repurposing, was used to identify connections between hypoxia-modulated genes in NB tumors and compounds. Two sets of 34 and 21 genes up- and down-regulated between hypoxic and normoxic primary NB tumors, respectively, were analyzed with CMap. The analysis reported a significant negative connectivity score across nine cell lines for 19 compounds mainly belonging to the class of PI3K/Akt/mTOR inhibitors. The gene expression profiles of NB cells cultured under hypoxic conditions and treated with the mTORC complex inhibitor PP242, referred to as the Mohlin dataset, was used to validate the CMap findings. A heat map representation of hypoxia-modulated genes in the Mohlin dataset and the gene set enrichment analysis (GSEA) showed an opposite regulation of these genes in the set of NB cells treated with the mTORC inhibitor PP242. In conclusion, our analysis identified inhibitors of the PI3K/Akt/mTOR signaling pathway as novel candidate compounds to treat NB patients with hypoxic tumors and a poor prognosis.

Reconciling multiple connectivity scores for drug repurposing

The basis of several recent methods for drug repurposing is the key principle that an efficacious drug will reverse the disease molecular ‘signature’ with minimal side effects. This principle was defined and popularized by the influential ‘connectivity map’ study in 2006 regarding reversal relationships between disease- and drug-induced gene expression profiles, quantified by a disease-drug ‘connectivity score.’ Over the past 15 years, several studies have proposed variations in calculating connectivity scores toward improving accuracy and robustness in light of massive growth in reference drug profiles. However, these variations have been formulated inconsistently using various notations and terminologies even though they are based on a common set of conceptual and statistical ideas. Therefore, we present a systematic reconciliation of multiple disease-drug similarity metrics (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$ES$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$css$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Sum$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Cosine$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XSum$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XCor$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XSpe$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XCos$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$EWCos$\end{document}) and connectivity scores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$CS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$RGES$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$NCS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$WCS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Tau$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$CSS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$EMUDRA$\end{document}) by defining them using consistent notation and terminology. In addition to providing clarity and deeper insights, this coherent definition of connectivity scores and their relationships provides a unified scheme that newer methods can adopt, enabling the computational drug-development community to compare and investigate different approaches easily. To facilitate the continuous and transparent integration of newer methods, this article will be available as a live document (https://jravilab.github.io/connectivity_scores) coupled with a GitHub repository (https://github.com/jravilab/connectivity_scores) that any researcher can build on and push changes to.

β-sitosterol reduces anxiety and synergizes with established anxiolytic drugs in mice

Anxiety and stress-related conditions represent a significant health burden in modern society. Unfortunately, most anxiolytic drugs are prone to side effects, limiting their long-term usage. Here, we employ a bioinformatics screen to identify drugs for repurposing as anxiolytics. Comparison of drug-induced gene-expression profiles with the hippocampal transcriptome of an importin α5 mutant mouse model with reduced anxiety identifies the hypocholesterolemic agent β-sitosterol as a promising candidate. β-sitosterol activity is validated by both intraperitoneal and oral application in mice, revealing it as the only clear anxiolytic from five closely related phytosterols. β-sitosterol injection reduces the effects of restraint stress, contextual fear memory, and c-Fos activation in the prefrontal cortex and dentate gyrus. Moreover, synergistic anxiolysis is observed when combining sub-efficacious doses of β-sitosterol with the SSRI fluoxetine. These preclinical findings support further development of β-sitosterol, either as a standalone anxiolytic or in combination with low-dose SSRIs.

Computational Drug Repositioning and Experimental Validation of Ivermectin in Treatment of Gastric Cancer

Objective: The aim of the present study was repositioning of ivermectin in treatment of gastric cancer (GC) by computational prediction based on gene expression profiles of human and mouse model of GC and validations with in silico, in vitro and in vivo approaches.

Methods: Computational drug repositioning was performed using connectivity map (cMap) and data/pathway mining with the Ingenuity Knowledge Base. Tissue samples of GC were collected from 16 patients and 57 mice for gene expression profiling. Additional seven independent datasets of gene expression of human GC from the TCGA database were used for validation. In silico testing was performed by constructing interaction networks of ivermectin and the downstream effects in targeted signaling pathways. In vitro testing was carried out in human GC cell lines (MKN74 and KATO-III). In vivo testing was performed in a transgenic mouse model of GC (INS-GAS mice).

Results: GC gene expression “signature” and data/pathway mining but not cMAP revealed nine molecular targets of ivermectin in both human and mouse GC associated with WNT/β-catenin signaling as well as cell proliferation pathways. In silico inhibition of the targets of ivermectin and concomitant activation of ivermectin led to the inhibition of WNT/β-catenin signaling pathway in “dose-depended” manner. In vitro, ivermectin inhibited cell proliferation in time- and concentration-depended manners, and cells were arrested in the G₁ phase at IC₅₀ and shifted to S phase arrest at >IC₅₀. In vivo, ivermectin reduced the tumor size which was associated with inactivation of WNT/β-catenin signaling and cell proliferation pathways and activation of cell death signaling pathways.

Conclusion: Ivermectin could be recognized as a repositioning candidate in treatment of gastric cancer.

Transcriptomics-based drug repositioning pipeline identifies therapeutic candidates for COVID-19

The novel SARS-CoV-2 virus emerged in December 2019 and has few effective treatments. We applied a computational drug repositioning pipeline to SARS-CoV-2 differential gene expression signatures derived from publicly available data. We utilized three independent published studies to acquire or generate lists of differentially expressed genes between control and SARS-CoV-2-infected samples. Using a rank-based pattern matching strategy based on the Kolmogorov-Smirnov Statistic, the signatures were queried against drug profiles from Connectivity Map (CMap). We validated sixteen of our top predicted hits in live SARS-CoV-2 antiviral assays in either Calu-3 or 293T-ACE2 cells. Validation experiments in human cell lines showed that 11 of the 16 compounds tested to date (including clofazimine, haloperidol and others) had measurable antiviral activity against SARS-CoV-2. These initial results are encouraging as we continue to work towards a further analysis of these predicted drugs as potential therapeutics for the treatment of COVID-19.

Repurposing novel therapeutic candidate drugs for coronavirus disease-19 based on protein-protein interaction network analysis

BACKGROUND

The coronavirus disease-19 (COVID-19) emerged in Wuhan, China and rapidly spread worldwide. Researchers are trying to find a way to treat this disease as soon as possible. The present study aimed to identify the genes involved in COVID-19 and find a new drug target therapy. Currently, there are no effective drugs targeting SARS-CoV-2, and meanwhile, drug discovery approaches are time-consuming and costly. To address this challenge, this study utilized a network-based drug repurposing strategy to rapidly identify potential drugs targeting SARS-CoV-2. To this end, seven potential drugs were proposed for COVID-19 treatment using protein-protein interaction (PPI) network analysis. First, 524 proteins in humans that have interaction with the SARS-CoV-2 virus were collected, and then the PPI network was reconstructed for these collected proteins. Next, the target miRNAs of the mentioned module genes were separately obtained from the miRWalk 2.0 database because of the important role of miRNAs in biological processes and were reported as an important clue for future analysis. Finally, the list of the drugs targeting module genes was obtained from the DGIDb database, and the drug-gene network was separately reconstructed for the obtained protein modules.

RESULTS

Based on the network analysis of the PPI network, seven clusters of proteins were specified as the complexes of proteins which are more associated with the SARS-CoV-2 virus. Moreover, seven therapeutic candidate drugs were identified to control gene regulation in COVID-19. PACLITAXEL, as the most potent therapeutic candidate drug and previously mentioned as a therapy for COVID-19, had four gene targets in two different modules. The other six candidate drugs, namely, BORTEZOMIB, CARBOPLATIN, CRIZOTINIB, CYTARABINE, DAUNORUBICIN, and VORINOSTAT, some of which were previously discovered to be efficient against COVID-19, had three gene targets in different modules. Eventually, CARBOPLATIN, CRIZOTINIB, and CYTARABINE drugs were found as novel potential drugs to be investigated as a therapy for COVID-19.

CONCLUSIONS

Our computational strategy for predicting repurposable candidate drugs against COVID-19 provides efficacious and rapid results for therapeutic purposes. However, further experimental analysis and testing such as clinical applicability, toxicity, and experimental validations are required to reach a more accurate and improved treatment. Our proposed complexes of proteins and associated miRNAs, along with discovered candidate drugs might be a starting point for further analysis by other researchers in this urgency of the COVID-19 pandemic.