Drugmonizome and Drugmonizome-ML: integration and abstraction of small molecule attributes for drug enrichment analysis and machine learning

DReAmocracy: A Method to Capitalise on Prior Drug Discovery Efforts to Highlight Candidate Drugs for Repurposing

In the area of drug research, several computational drug repurposing studies have highlighted candidate repurposed drugs, as well as clinical trial studies that have tested/are testing drugs in different phases. To the best of our knowledge, the aggregation of the proposed lists of drugs by previous studies has not been extensively exploited towards generating a dynamic reference matrix with enhanced resolution. To fill this knowledge gap, we performed weight-modulated majority voting of the modes of action, initial indications and targeted pathways of the drugs in a well-known repository, namely the Drug Repurposing Hub. Our method, DReAmocracy, exploits this pile of information and creates frequency tables and, finally, a disease suitability score for each drug from the selected library. As a testbed, we applied this method to a group of neurodegenerative diseases (Alzheimer's, Parkinson's, Huntington's disease and Multiple Sclerosis). A super-reference table with drug suitability scores has been created for all four neurodegenerative diseases and can be queried for any drug candidate against them. Top-scored drugs for Alzheimer's Disease include agomelatine, mirtazapine and vortioxetine; for Parkinson's Disease, they include apomorphine, pramipexole and lisuride; for Huntington's, they include chlorpromazine, fluphenazine and perphenazine; and for Multiple Sclerosis, they include zonisamide, disopyramide and priralfimide. Overall, DReAmocracy is a methodology that focuses on leveraging the existing drug-related experimental and/or computational knowledge rather than a predictive model for drug repurposing, offering a quantified aggregation of existing drug discovery results to (1) reveal trends in selected tracks of drug discovery research with increased resolution that includes modes of action, targeted pathways and initial indications for the investigated drugs and (2) score new candidate drugs for repurposing against a selected disease.

Expression of the checkpoint kinase BUB1 is a predictor of response to cancer therapies

The identification of clinically-relevant biomarkers is of upmost importance for the management of cancer, from diagnosis to treatment choices. We performed a pan-cancer analysis of the mitotic checkpoint budding uninhibited by benzimidazole 1 gene BUB1, in the attempt to ascertain its diagnostic and prognostic values, specifically in the context of drug response. BUB1 was found to be overexpressed in the majority of cancers, and particularly elevated in clinically aggressive molecular subtypes. Its expression was correlated with clinico-phenotypic features, notably tumour staging, size, invasion, hypoxia, and stemness. In terms of prognostic value, the expression of BUB1 bore differential clinical outcomes depending on the treatment administered in TCGA cancer cohorts, suggesting sensitivity or resistance, depending on the expression levels. We also integrated in vitro drug sensitivity data from public projects based on correlation between drug efficacy and BUB1 expression to produce a list of candidate compounds with differential responses according to BUB1 levels. Gene Ontology enrichment analyses revealed that BUB1 overexpression in cancer is associated with biological processes related to mitosis and chromosome segregation machinery, reflecting the mechanisms of action of drugs with a differential effect based on BUB1 expression.

DendroX: multi-level multi-cluster selection in dendrograms

BACKGROUND

Cluster heatmaps are widely used in biology and other fields to uncover clustering patterns in data matrices. Most cluster heatmap packages provide utility functions to divide the dendrograms at a certain level to obtain clusters, but it is often difficult to locate the appropriate cut in the dendrogram to obtain the clusters seen in the heatmap or computed by a statistical method. Multiple cuts are required if the clusters locate at different levels in the dendrogram.

RESULTS

We developed DendroX, a web app that provides interactive visualization of a dendrogram where users can divide the dendrogram at any level and in any number of clusters and pass the labels of the identified clusters for functional analysis. Helper functions are provided to extract linkage matrices from cluster heatmap objects in R or Python to serve as input to the app. A graphic user interface was also developed to help prepare input files for DendroX from data matrices stored in delimited text files. The app is scalable and has been tested on dendrograms with tens of thousands of leaf nodes. As a case study, we clustered the gene expression signatures of 297 bioactive chemical compounds in the LINCS L1000 dataset and visualized them in DendroX. Seventeen biologically meaningful clusters were identified based on the structure of the dendrogram and the expression patterns in the heatmap. We found that one of the clusters consisting of mostly naturally occurring compounds is not previously reported and has its members sharing broad anticancer, anti-inflammatory and antioxidant activities.

CONCLUSIONS

DendroX solves the problem of matching visually and computationally determined clusters in a cluster heatmap and helps users navigate among different parts of a dendrogram. The identification of a cluster of naturally occurring compounds with shared bioactivities implicates a convergence of biological effects through divergent mechanisms.

Predicting drug response from single-cell expression profiles of tumours

BACKGROUND

Intra-tumour heterogeneity (ITH) presents a significant obstacle in formulating effective treatment strategies in clinical practice. Single-cell RNA sequencing (scRNA-seq) has evolved as a powerful instrument for probing ITH at the transcriptional level, offering an unparalleled opportunity for therapeutic intervention.

RESULTS

Drug response prediction at the single-cell level is an emerging field of research that aims to improve the efficacy and precision of cancer treatments. Here, we introduce DREEP (Drug Response Estimation from single-cell Expression Profiles), a computational method that leverages publicly available pharmacogenomic screens from GDSC2, CTRP2, and PRISM and functional enrichment analysis to predict single-cell drug sensitivity from transcriptomic data. We validated DREEP extensively in vitro using several independent single-cell datasets with over 200 cancer cell lines and showed its accuracy and robustness. Additionally, we also applied DREEP to molecularly barcoded breast cancer cells and identified drugs that can selectively target specific cell populations.

CONCLUSIONS

DREEP provides an in silico framework to prioritize drugs from single-cell transcriptional profiles of tumours and thus helps in designing personalized treatment strategies and accelerating drug repurposing studies. DREEP is available at https://github.com/gambalab/DREEP .

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 .

A systematic review of artificial intelligence-based COVID-19 modeling on multimodal genetic information

This study systematically reviews the Artificial Intelligence (AI) methods developed to resolve the critical process of COVID-19 gene data analysis, including diagnosis, prognosis, biomarker discovery, drug responsiveness, and vaccine efficacy. This systematic review follows the guidelines of Preferred Reporting for Systematic Reviews and Meta-Analyses (PRISMA). We searched PubMed, Embase, Web of Science, and Scopus databases to identify the relevant articles from January 2020 to June 2022. It includes the published studies of AI-based COVID-19 gene modeling extracted through relevant keyword searches in academic databases. This study included 48 articles discussing AI-based genetic studies for several objectives. Ten articles confer about the COVID-19 gene modeling with computational tools, and five articles evaluated ML-based diagnosis with observed accuracy of 97% on SARS-CoV-2 classification. Gene-based prognosis study reviewed three articles and found host biomarkers detecting COVID-19 progression with 90% accuracy. Twelve manuscripts reviewed the prediction models with various genome analysis studies, nine articles examined the gene-based in silico drug discovery, and another nine investigated the AI-based vaccine development models. This study compiled the novel coronavirus gene biomarkers and targeted drugs identified through ML approaches from published clinical studies. This review provided sufficient evidence to delineate the potential of AI in analyzing complex gene information for COVID-19 modeling on multiple aspects like diagnosis, drug discovery, and disease dynamics. AI models entrenched a substantial positive impact by enhancing the efficiency of the healthcare system during the COVID-19 pandemic.

Repurposable Drugs That Interact with Steroid Responsive Gene Targets for Inner Ear Disease

Corticosteroids, oral or transtympanic, remain the mainstay for inner ear diseases characterized by hearing fluctuation or sudden changes in hearing, including sudden sensorineural hearing loss (SSNHL), Meniere's disease (MD), and autoimmune inner ear disease (AIED). Despite their use across these diseases, the rate of complete recovery remains low, and results across the literature demonstrates significant heterogeneity with respect to the effect of corticosteroids, suggesting a need to identify more efficacious treatment options. Previously, our group has cross-referenced steroid-responsive genes in the cochlea with published single-cell and single-nucleus transcriptome datasets to demonstrate that steroid-responsive differentially regulated genes are expressed in spiral ganglion neurons (SGN) and stria vascularis (SV) cell types. These differentially regulated genes represent potential druggable gene targets. We utilized multiple gene target databases (DrugBank, Pharos, and LINCS) to identify orally administered, FDA approved medications that potentially target these genes. We identified 42 candidate drugs that have been shown to interact with these genes, with an emphasis on safety profile, and tolerability. This study utilizes multiple databases to identify drugs that can target a number of druggable genes in otologic disorders that are commonly treated with steroids, providing a basis for establishing novel repurposing treatment trials.

SigCom LINCS: data and metadata search engine for a million gene expression signatures

Millions of transcriptome samples were generated by the Library of Integrated Network-based Cellular Signatures (LINCS) program. When these data are processed into searchable signatures along with signatures extracted from Genotype-Tissue Expression (GTEx) and Gene Expression Omnibus (GEO), connections between drugs, genes, pathways and diseases can be illuminated. SigCom LINCS is a webserver that serves over a million gene expression signatures processed, analyzed, and visualized from LINCS, GTEx, and GEO. SigCom LINCS is built with Signature Commons, a cloud-agnostic skeleton Data Commons with a focus on serving searchable signatures. SigCom LINCS provides a rapid signature similarity search for mimickers and reversers given sets of up and down genes, a gene set, a single gene, or any search term. Additionally, users of SigCom LINCS can perform a metadata search to find and analyze subsets of signatures and find information about genes and drugs. SigCom LINCS is findable, accessible, interoperable, and reusable (FAIR) with metadata linked to standard ontologies and vocabularies. In addition, all the data and signatures within SigCom LINCS are available via a well-documented API. In summary, SigCom LINCS, available at https://maayanlab.cloud/sigcom-lincs, is a rich webserver resource for accelerating drug and target discovery in systems pharmacology.

Gene and drug landing page aggregator

Motivation

Many biological and biomedical researchers commonly search for information about genes and drugs to gather knowledge from these resources. For the most part, such information is served as landing pages in disparate data repositories and web portals.

Results

The Gene and Drug Landing Page Aggregator (GDLPA) provides users with access to 50 gene-centric and 19 drug-centric repositories, enabling them to retrieve landing pages corresponding to their gene and drug queries. Bringing these resources together into one dashboard that directs users to the landing pages across many resources can help centralize gene- and drug-centric knowledge, as well as raise awareness of available resources that may be missed when using standard search engines. To demonstrate the utility of GDLPA, case studies for the gene klotho and the drug remdesivir were developed. The first case study highlights the potential role of klotho as a drug target for aging and kidney disease, while the second study gathers knowledge regarding approval, usage, and safety for remdesivir, the first approved coronavirus disease 2019 therapeutic. Finally, based on our experience, we provide guidelines for developing effective landing pages for genes and drugs.

Availability and implementation

GDLPA is open source and is available from: https://cfde-gene-pages.cloud/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

DrugShot: querying biomedical search terms to retrieve prioritized lists of small molecules

Background

PubMed contains millions of abstracts that co-mention terms that describe drugs with other biomedical terms such as genes or diseases. Unique opportunities exist for leveraging these co-mentions by integrating them with other drug-drug similarity resources such as the Library of Integrated Network-based Cellular Signatures (LINCS) L1000 signatures to develop novel hypotheses.

Results

DrugShot is a web-based server application and an Appyter that enables users to enter any biomedical search term into a simple input form to receive ranked lists of drugs and other small molecules based on their relevance to the search term. To produce ranked lists of small molecules, DrugShot cross-references returned PubMed identifiers (PMIDs) with DrugRIF or AutoRIF, which are curated resources of drug-PMID associations, to produce an associated small molecule list where each small molecule is ranked according to total co-mentions with the search term from shared PubMed IDs. Additionally, using two types of drug-drug similarity matrices, lists of small molecules are predicted to be associated with the search term. Such predictions are based on literature co-mentions and signature similarity from LINCS L1000 drug-induced gene expression profiles.

Conclusions

DrugShot prioritizes drugs and small molecules associated with biomedical search terms. In addition to listing known associations, DrugShot predicts additional drugs and small molecules related to any search term. Hence, DrugShot can be used to prioritize drugs and preclinical compounds for drug repurposing and suggest indications and adverse events for preclinical compounds. DrugShot is freely and openly available at: https://maayanlab.cloud/drugshot and https://appyters.maayanlab.cloud/#/DrugShot.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04590-5.

Getting Started with the IDG KMC Datasets and Tools

The Illuminating the Druggable Genome (IDG) consortium is a National Institutes of Health (NIH) Common Fund program designed to enhance our knowledge of under-studied proteins, more specifically, proteins unannotated within the three most commonly drug-targeted protein families: G-protein coupled receptors, ion channels, and protein kinases. Since 2014, the IDG Knowledge Management Center (IDG-KMC) has generated several open-access datasets and resources that jointly serve as a highly translational machine-learning-ready knowledgebase focused on human protein-coding genes and their products. The goal of the IDG-KMC is to develop comprehensive integrated knowledge for the druggable genome to illuminate the uncharacterized or poorly annotated portion of the druggable genome. The tools derived from the IDG-KMC provide either user-friendly visualizations or ways to impute the knowledge about potential targets using machine learning strategies. In the following protocols, we describe how to use each web-based tool to accelerate illumination in under-studied proteins. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Interacting with the Pharos user interface Basic Protocol 2: Accessing the data in Harmonizome Basic Protocol 3: The ARCHS4 resource Basic Protocol 4: Making predictions about gene function with PrismExp Basic Protocol 5: Using Geneshot to illuminate knowledge about under-studied targets Basic Protocol 6: Exploring under-studied targets with TIN-X Basic Protocol 7: Interacting with the DrugCentral user interface Basic Protocol 8: Estimating Anti-SARS-CoV-2 activities with DrugCentral REDIAL-2020 Basic Protocol 9: Drug Set Enrichment Analysis using Drugmonizome Basic Protocol 10: The Drugmonizome-ML Appyter Basic Protocol 11: The Harmonizome-ML Appyter Basic Protocol 12: GWAS target illumination with TIGA Basic Protocol 13: Prioritizing kinases for lists of proteins and phosphoproteins with KEA3 Basic Protocol 14: Converting PubMed searches to drug sets with the DrugShot Appyter.

Drugmonizome and Drugmonizome-ML: integration and abstraction of small molecule attributes for drug enrichment analysis and machine learning

Understanding the underlying molecular and structural similarities between seemingly heterogeneous sets of drugs can aid in identifying drug repurposing opportunities and assist in the discovery of novel properties of preclinical small molecules. A wealth of information about drug and small molecule structure, targets, indications and side effects; induced gene expression signatures; and other attributes are publicly available through web-based tools, databases and repositories. By processing, abstracting and aggregating information from these resources into drug set libraries, knowledge about novel properties of drugs and small molecules can be systematically imputed with machine learning. In addition, drug set libraries can be used as the underlying database for drug set enrichment analysis. Here, we present Drugmonizome, a database with a search engine for querying annotated sets of drugs and small molecules for performing drug set enrichment analysis. Utilizing the data within Drugmonizome, we also developed Drugmonizome-ML. Drugmonizome-ML enables users to construct customized machine learning pipelines using the drug set libraries from Drugmonizome. To demonstrate the utility of Drugmonizome, drug sets from 12 independent SARS-CoV-2 in vitro screens were subjected to consensus enrichment analysis. Despite the low overlap among these 12 independent in vitro screens, we identified common biological processes critical for blocking viral replication. To demonstrate Drugmonizome-ML, we constructed a machine learning pipeline to predict whether approved and preclinical drugs may induce peripheral neuropathy as a potential side effect. Overall, the Drugmonizome and Drugmonizome-ML resources provide rich and diverse knowledge about drugs and small molecules for direct systems pharmacology applications. Database URL: https://maayanlab.cloud/drugmonizome/.