Reconciling multiple connectivity scores for drug repurposing

Translational regulation of SND1 governs endothelial homeostasis during stress

Translational control shapes the proteome and is particularly important in regulating gene expression under stress. A key source of endothelial stress is treatment with tyrosine kinase inhibitors (TKIs), which lowers cancer mortality but increases cardiovascular mortality. Using a human induced pluripotent stem cell-derived endothelial cell (hiPSC-EC) model of sunitinib-induced vascular dysfunction combined with ribosome profiling, we assessed the role of translational control in hiPSC-ECs in response to stress. We identified staphylococcal nuclease and tudor domain-containing protein 1 (SND1) as a sunitinib-dependent translationally repressed gene. SND1 translational repression was mediated by the mTORC1/4E-BP1 pathway. SND1 inhibition led to endothelial dysfunction, whereas SND1 OE protected against sunitinib-induced endothelial dysfunction. Mechanistically, SND1 transcriptionally regulated UBE2N, an E2-conjugating enzyme that mediates K63-linked ubiquitination. UBE2N along with the E3 ligases RNF8 and RNF168 regulated the DNA damage repair response pathway to mitigate the deleterious effects of sunitinib. In silico analysis of FDA-approved drugs led to the identification of an ACE inhibitor, ramipril, that protected against sunitinib-induced vascular dysfunction in vitro and in vivo, all while preserving the efficacy of cancer therapy. Our study established a central role for translational control of SND1 in sunitinib-induced endothelial dysfunction that could potentially be therapeutically targeted to reduce sunitinib-induced vascular toxicity.

A review on drug repurposing applicable to obesity

Obesity is a major public health concern and burden on individuals and healthcare systems. Due to the challenges and limitations of lifestyle adjustments, it is advisable to consider pharmacological treatment for people affected by obesity. However, the side effects and limited efficacy of available drugs make the obesity drug market far from sufficient. Drug repurposing involves identifying new applications for existing drugs and offers some advantages over traditional drug development approaches including lower costs and shorter development timelines. This review aims to provide an overview of drug repurposing for anti-obesity medications, including the rationale for repurposing, the challenges and approaches, and the potential drugs that are being investigated for repurposing. Through advanced computational techniques, researchers can unlock the potential of repurposed drugs to tackle the global obesity epidemic. Further research, clinical trials, and collaborative efforts are essential to fully explore and leverage the potential of drug repurposing in the fight against obesity.

Molecular Clustering of Metabolic Dysfunction-Associated Steatotic Liver Disease Based on Transcriptome Analysis

Background: Metabolic dysfunction-associated steatotic liver disease (MASLD) is a complex metabolic disorder with a diverse spectrum. This study aimed to classify patients with MASLD into molecular subtypes based on the underlying pathophysiology. Methods: We performed high-throughput RNA sequencing on 164 liver tissue samples from healthy controls and patients with MASLD. The clustering was based on individual genes or pathways that showed high variation across the samples. Second, the clustering was based on single-sample gene set enrichment analysis. Results: Optimal homogeneity was achieved by dividing the samples into four clusters (one healthy control and three MASLD clusters I-III) based on the top genes or pathways with differences across the samples. No significant differences were observed in waist circumference, blood pressure, glucose, alanine transaminase (ALT), or aspartate transferase (AST) levels between cluster I patients with MASLD and the healthy controls. Cluster I showed the clinical features of lean MASLD. Cluster III of MASLD demonstrated hypertension and a T2DM prevalence of 57.1% and 50.0%, respectively, with a significantly higher fibrosis burden (stage of fibrosis, liver stiffness, and FIB-4 value) than clusters I and II. Cluster III was associated with severe fibrosis and abnormal glucose homeostasis. In MASLD cluster I, the sphingolipid and GPCR pathways were upregulated, whereas the complement and phagocytosis pathways were downregulated. In MASLD cluster II, the cell cycle and NOTCH3 pathways increased, whereas the PI3K and insulin-related pathways decreased. In MASLD cluster III, increased activity occurred in the interleukin-2 and -4 and extracellular matrix pathways, coupled with decreased activity in the serotonin 2A and B pathways. Conclusions: MASLD can be divided into three distinct molecular phenotypes, wherein each is characterized by a specific molecular pathway.

Combinatorial prediction of therapeutic perturbations using causally-inspired neural networks

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

Integrating pharmacogenomics and cheminformatics with diverse disease phenotypes for cell type-guided drug discovery

BACKGROUND

Large-scale pharmacogenomic resources, such as the Connectivity Map (CMap), have greatly assisted computational drug discovery. However, despite their widespread use, CMap-based methods have thus far been agnostic to the biological activity of drugs as well as to the genomic effects of drugs in multiple disease contexts. Here, we present a network-based statistical approach, Pathopticon, that uses CMap to build cell type-specific gene-drug perturbation networks and integrates these networks with cheminformatic data and diverse disease phenotypes to prioritize drugs in a cell type-dependent manner.

METHODS

We build cell type-specific gene-drug perturbation networks from CMap data using a statistical procedure we call Quantile-based Instance Z-score Consensus (QUIZ-C). Using these networks and a large-scale disease-gene network consisting of 569 disease signatures from the Enrichr database, we calculate Pathophenotypic Congruity Scores (PACOS) between input gene signatures and drug perturbation signatures and combine these scores with cheminformatic data from ChEMBL to prioritize drugs. We benchmark our approach by calculating area under the receiver operating characteristic curves (AUROC) for 73 gene sets from the Molecular Signatures Database (MSigDB) using target gene expression profiles from the Comparative Toxicogenomics Database (CTD). We validate the drugs predicted in our proofs-of-concept using real-time polymerase chain reaction (qPCR) experiments.

RESULTS

Cell type-specific gene-drug perturbation networks built using QUIZ-C are topologically distinct, reflecting the biological uniqueness of the cell lines in CMap, and are enriched in known drug targets. Pathopticon demonstrates a better prediction performance than solely cheminformatic measures as well as state-of-the-art network and deep learning-based methods. Top predictions made by Pathopticon have high chemical structural diversity, suggesting their potential for building compound libraries. In proof-of-concept applications on vascular diseases, we demonstrate that Pathopticon helps guide in vitro experiments by identifying pathways that are potentially regulated by the predicted therapeutic candidates.

CONCLUSIONS

Our network-based analytical framework integrating pharmacogenomics and cheminformatics (available at https://github.com/r-duh/Pathopticon ) provides a feasible blueprint for a cell type-specific drug discovery and repositioning platform with broad implications for the efficiency and success of drug development.

Transcriptomic analysis reveals Streptococcus agalactiae activation of oncogenic pathways in cervical adenocarcinoma

Cervical adenocarcinoma (AC), a subtype of uterine cervical cancer (CC), poses a challenge due to its resistance to therapy and poor prognosis compared with squamous cervical carcinoma. Streptococcus agalactiae [group B Streptococcus (GBS)], a Gram-positive coccus, has been associated with cervical intraepithelial neoplasia in CC. However, the underlying mechanism interaction between GBS and CC, particularly AC, remains elusive. Leveraging The Cancer Genome Atlas public data and time-series transcriptomic data, the present study investigated the interaction between GBS and AC, revealing activation of two pivotal pathways: 'MAPK signaling pathway' and 'mTORC1 signaling'. Western blotting, reverse transcription-quantitative PCR and cell viability assays were performed to validate the activation of these pathways and their role in promoting cancer cell proliferation. Subsequently, the present study evaluated the efficacy of two anticancer drugs targeting these pathways (binimetinib and ridaforolimus) in AC cell treatment. Binimetinib demonstrated a cytostatic effect, while ridaforolimus had a modest impact on HeLa cells after 48 h of treatment, as observed in both cell viability and cytotoxicity assays. The combination of binimetinib and ridaforolimus resulted in a significantly greater cytotoxic effect compared to binimetinib or ridaforolimus monotherapy, although the synergy score indicated an additive effect. In general, the MAPK and mTORC1 signaling pathways were identified as the main pathways associated with GBS and AC cells. The combination of binimetinib and ridaforolimus could be a potential AC treatment.

Signature Search Polestar: a comprehensive drug repurposing method evaluation assistant for customized oncogenic signature

SUMMARY

The burgeoning high-throughput technologies have led to a significant surge in the scale of pharmacotranscriptomic datasets, especially for oncology. Signature search methods (SSMs), utilizing oncogenic signatures formed by differentially expressed genes through sequencing, have been instrumental in anti-cancer drug screening and identifying mechanisms of action without relying on prior knowledge. However, various studies have found that different SSMs exhibit varying performance across pharmacotranscriptomic datasets. In addition, the size of the oncogenic signature can also significantly impact the result of drug repurposing. Therefore, finding the optimal SSMs and customized oncogenic signature for a specific disease remains a challenge. To address this, we introduce Signature Search Polestar (SSP), a webserver integrating the largest pharmacotranscriptomic datasets of anti-cancer drugs from LINCS L1000 with five state-of-the-art SSMs (XSum, CMap, GSEA, ZhangScore, XCos). SSP provides three main modules: Benchmark, Robustness, and Application. Benchmark uses two indices, Area Under the Curve and Enrichment Score, based on drug annotations to evaluate SSMs at different oncogenic signature sizes. Robustness, applicable when drug annotations are insufficient, uses a performance score based on drug self-retrieval for evaluation. Application provides three screening strategies, single method, SS_all, and SS_cross, allowing users to freely utilize optimal SSMs with tailored oncogenic signature for drug repurposing.

AVAILABILITY AND IMPLEMENTATION

SSP is free at https://web.biotcm.net/SSP/. The current version of SSP is archived in https://doi.org/10.6084/m9.figshare.26524741.v1, allowing users to directly use or customize their own SSP webserver.

Integrated transcriptomics- and structure-based drug repositioning identifies drugs with proteasome inhibitor properties

Computational pharmacogenomics can potentially identify new indications for already approved drugs and pinpoint compounds with similar mechanism-of-action. Here, we used an integrated drug repositioning approach based on transcriptomics data and structure-based virtual screening to identify compounds with gene signatures similar to three known proteasome inhibitors (PIs; bortezomib, MG-132, and MLN-2238). In vitro validation of candidate compounds was then performed to assess proteasomal proteolytic activity, accumulation of ubiquitinated proteins, cell viability, and drug-induced expression in A375 melanoma and MCF7 breast cancer cells. Using this approach, we identified six compounds with PI properties ((-)-kinetin-riboside, manumycin-A, puromycin dihydrochloride, resistomycin, tegaserod maleate, and thapsigargin). Although the docking scores pinpointed their ability to bind to the β5 subunit, our in vitro study revealed that these compounds inhibited the β1, β2, and β5 catalytic sites to some extent. As shown with bortezomib, only manumycin-A, puromycin dihydrochloride, and tegaserod maleate resulted in excessive accumulation of ubiquitinated proteins and elevated HMOX1 expression. Taken together, our integrated drug repositioning approach and subsequent in vitro validation studies identified six compounds demonstrating properties similar to proteasome inhibitors.

Optimizing in silico drug discovery: simulation of connected differential expression signatures and applications to benchmarking

BACKGROUND

We present a novel simulation method for generating connected differential expression signatures. Traditional methods have struggled with the lack of reliable benchmarking data and biases in drug-disease pair labeling, limiting the rigorous benchmarking of connectivity-based approaches.

OBJECTIVE

Our aim is to develop a simulation method based on a statistical framework that allows for adjustable levels of parametrization, especially the connectivity, to generate a pair of interconnected differential signatures. This could help to address the issue of benchmarking data availability for connectivity-based drug repurposing approaches.

METHODS

We first detailed the simulation process and how it reflected real biological variability and the interconnectedness of gene expression signatures. Then, we generated several datasets to enable the evaluation of different existing algorithms that compare differential expression signatures, providing insights into their performance and limitations.

RESULTS

Our findings demonstrate the ability of our simulation to produce realistic data, as evidenced by correlation analyses and the log2 fold-change distribution of deregulated genes. Benchmarking reveals that methods like extreme cosine similarity and Pearson correlation outperform others in identifying connected signatures.

CONCLUSION

Overall, our method provides a reliable tool for simulating differential expression signatures. The data simulated by our tool encompass a wide spectrum of possibilities to challenge and evaluate existing methods to estimate connectivity scores. This may represent a critical gap in connectivity-based drug repurposing research because reliable benchmarking data are essential for assessing and advancing in the development of new algorithms. The simulation tool is available as a R package (General Public License (GPL) license) at https://github.com/cgonzalez-gomez/cosimu.

PharmGWAS: a GWAS-based knowledgebase for drug repurposing

Leveraging genetics insights to promote drug repurposing has become a promising and active strategy in pharmacology. Indeed, among the 50 drugs approved by FDA in 2021, two-thirds have genetically supported evidence. In this regard, the increasing amount of widely available genome-wide association studies (GWAS) datasets have provided substantial opportunities for drug repurposing based on genetics discoveries. Here, we developed PharmGWAS, a comprehensive knowledgebase designed to identify candidate drugs through the integration of GWAS data. PharmGWAS focuses on novel connections between diseases and small-molecule compounds derived using a reverse relationship between the genetically-regulated expression signature and the drug-induced signature. Specifically, we collected and processed 1929 GWAS datasets across a diverse spectrum of diseases and 724 485 perturbation signatures pertaining to a substantial 33609 molecular compounds. To obtain reliable and robust predictions for the reverse connections, we implemented six distinct connectivity methods. In the current version, PharmGWAS deposits a total of 740 227 genetically-informed disease-drug pairs derived from drug-perturbation signatures, presenting a valuable and comprehensive catalog. Further equipped with its user-friendly web design, PharmGWAS is expected to greatly aid the discovery of novel drugs, the exploration of drug combination therapies and the identification of drug resistance or side effects. PharmGWAS is available at https://ngdc.cncb.ac.cn/pharmgwas.

Small molecule-mediated targeting of microRNAs for drug discovery: Experiments, computational techniques, and disease implications

Small molecules have been providing medical breakthroughs for human diseases for more than a century. Recently, identifying small molecule inhibitors that target microRNAs (miRNAs) has gained importance, despite the challenges posed by labour-intensive screening experiments and the significant efforts required for medicinal chemistry optimization. Numerous experimentally-verified cases have demonstrated the potential of miRNA-targeted small molecule inhibitors for disease treatment. This new approach is grounded in their posttranscriptional regulation of the expression of disease-associated genes. Reversing dysregulated gene expression using this mechanism may help control dysfunctional pathways. Furthermore, the ongoing improvement of algorithms has allowed for the integration of computational strategies built on top of laboratory-based data, facilitating a more precise and rational design and discovery of lead compounds. To complement the use of extensive pharmacogenomics data in prioritising potential drugs, our previous work introduced a computational approach based on only molecular sequences. Moreover, various computational tools for predicting molecular interactions in biological networks using similarity-based inference techniques have been accumulated in established studies. However, there are a limited number of comprehensive reviews covering both computational and experimental drug discovery processes. In this review, we outline a cohesive overview of both biological and computational applications in miRNA-targeted drug discovery, along with their disease implications and clinical significance. Finally, utilizing drug-target interaction (DTIs) data from DrugBank, we showcase the effectiveness of deep learning for obtaining the physicochemical characterization of DTIs.

DeepdlncUD: Predicting regulation types of small molecule inhibitors on modulating lncRNA expression by deep learning

Targeting lncRNAs by small molecules (SM-lncR) to alter their expression levels has emerged as an important therapeutic modality for disease treatment. To date, no computational tools have been dedicated to predicting small molecule-mediated upregulation or downregulation of lncRNA expression. Here, we introduce DeepdlncUD, which integrates predictions of nine deep learning algorithms together, to infer the regulation types of small molecules on modulating lncRNA expression. Through systematic optimization on a training set of 771 upregulation and 739 downregulation SM-lncR pairs, each encoding 1369 sequence, representational, and physiochemical features, this method outperforms a recently released program, DeepsmirUD, by achieving 0.674 in AUC (area under the receiver operating characteristic curve), 0.722 in AUCPR (area under the precision-recall curve), 0.681 in F1-score, and 0.516 in Jaccard Index on a test set of 222 SM-lncR pairs. By extracting 125 upregulation and 46 downregulation SM-lncR pairs that involve disease-associated lncRNAs, DeepdlncUD is shown to gain an accuracy of 0.700 in the pathological context. Using connectivity scores, around half of the small molecules are correctly estimated as drugs to treat lncRNA-regulated diseases. This tool can be run at a fast speed to assist the discovery of potential small molecule drugs of lncRNA targets on a large scale. DeepdlncUD is publicly available at https://github.com/2003100127/deepdlncud.

Cell-specific imputation of drug connectivity mapping with incomplete data

Drug repositioning allows expedited discovery of new applications for existing compounds, but re-screening vast compound libraries is often prohibitively expensive. "Connectivity mapping" is a process that links drugs to diseases by identifying compounds whose impact on expression in a collection of cells reverses the disease's impact on expression in disease-relevant tissues. The LINCS project has expanded the universe of compounds and cells for which data are available, but even with this effort, many clinically useful combinations are missing. To evaluate the possibility of repurposing drugs despite missing data, we compared collaborative filtering using either neighborhood-based or SVD imputation methods to two naive approaches via cross-validation. Methods were evaluated for their ability to predict drug connectivity despite missing data. Predictions improved when cell type was taken into account. Neighborhood collaborative filtering was the most successful method, with the best improvements in non-immortalized primary cells. We also explored which classes of compounds are most and least reliant on cell type for accurate imputation. We conclude that even for cells in which drug responses have not been fully characterized, it is possible to identify unassayed drugs that reverse in those cells the expression signatures observed in disease.

DeepsmirUD: Prediction of Regulatory Effects on microRNA Expression Mediated by Small Molecules Using Deep Learning

Aberrant miRNA expression has been associated with a large number of human diseases. Therefore, targeting miRNAs to regulate their expression levels has become an important therapy against diseases that stem from the dysfunction of pathways regulated by miRNAs. In recent years, small molecules have demonstrated enormous potential as drugs to regulate miRNA expression (i.e., SM-miR). A clear understanding of the mechanism of action of small molecules on the upregulation and downregulation of miRNA expression allows precise diagnosis and treatment of oncogenic pathways. However, outside of a slow and costly process of experimental determination, computational strategies to assist this on an ad hoc basis have yet to be formulated. In this work, we developed, to the best of our knowledge, the first cross-platform prediction tool, DeepsmirUD, to infer small-molecule-mediated regulatory effects on miRNA expression (i.e., upregulation or downregulation). This method is powered by 12 cutting-edge deep-learning frameworks and achieved AUC values of 0.843/0.984 and AUCPR values of 0.866/0.992 on two independent test datasets. With a complementarily constructed network inference approach based on similarity, we report a significantly improved accuracy of 0.813 in determining the regulatory effects of nearly 650 associated SM-miR relations, each formed with either novel small molecule or novel miRNA. By further integrating miRNA-cancer relationships, we established a database of potential pharmaceutical drugs from 1343 small molecules for 107 cancer diseases to understand the drug mechanisms of action and offer novel insight into drug repositioning. Furthermore, we have employed DeepsmirUD to predict the regulatory effects of a large number of high-confidence associated SM-miR relations. Taken together, our method shows promise to accelerate the development of potential miRNA targets and small molecule drugs.

Histone Deacetylase Inhibitors Restore Cancer Cell Sensitivity towards T Lymphocytes Mediated Cytotoxicity in Pancreatic Cancer

Despite medical advancements, the prognosis of pancreatic ductal adenocarcinoma (PDAC) has not improved significantly over the past 50 years. By utilising the large-scale genomic datasets available from the Australia Pancreatic Cancer Project (PACA-AU) and The Cancer Genomic Atlas Project (TCGA-PAAD), we studied the immunophenotype of PDAC in silico and identified that tumours with high cytotoxic T lymphocytes (CTL) killing activity were associated with favourable clinical outcomes. Using the STRING protein-protein interaction network analysis, the identified differentially expressed genes with low CTL killing activity were associated with TWIST/IL-6R, HDAC5, and EOMES signalling. Following Connectivity Map analysis, we identified 44 small molecules that could restore CTL sensitivity in the PDAC cells. Further high-throughput chemical library screening identified 133 inhibitors that effectively target both parental and CTL-resistant PDAC cells in vitro. Since CTL-resistant PDAC had a higher expression of histone proteins and its acetylated proteins compared to its parental cells, we further investigated the impact of histone deacetylase inhibitors (HDACi) on CTL-mediated cytotoxicity in PDAC cells in vitro, namely SW1990 and BxPC3. Further analyses revealed that givinostat and dacinostat were the two most potent HDAC inhibitors that restored CTL sensitivity in SW1990 and BxPC3 CTL-resistant cells. Through our in silico and in vitro studies, we demonstrate the novel role of HDAC inhibition in restoring CTL resistance and that combinations of HDACi with CTL may represent a promising therapeutic strategy, warranting its further detailed molecular mechanistic studies and animal studies before embarking on the clinical evaluation of these novel combined PDAC treatments.

Causal reasoning over knowledge graphs leveraging drug-perturbed and disease-specific transcriptomic signatures for drug discovery

Network-based approaches are becoming increasingly popular for drug discovery as they provide a systems-level overview of the mechanisms underlying disease pathophysiology. They have demonstrated significant early promise over other methods of biological data representation, such as in target discovery, side effect prediction and drug repurposing. In parallel, an explosion of -omics data for the deep characterization of biological systems routinely uncovers molecular signatures of disease for similar applications. Here, we present RPath, a novel algorithm that prioritizes drugs for a given disease by reasoning over causal paths in a knowledge graph (KG), guided by both drug-perturbed as well as disease-specific transcriptomic signatures. First, our approach identifies the causal paths that connect a drug to a particular disease. Next, it reasons over these paths to identify those that correlate with the transcriptional signatures observed in a drug-perturbation experiment, and anti-correlate to signatures observed in the disease of interest. The paths which match this signature profile are then proposed to represent the mechanism of action of the drug. We demonstrate how RPath consistently prioritizes clinically investigated drug-disease pairs on multiple datasets and KGs, achieving better performance over other similar methodologies. Furthermore, we present two case studies showing how one can deconvolute the predictions made by RPath as well as predict novel targets.

DREAM: a database of experimentally supported protein-coding RNAs and drug associations in human cancer

The Drug Response Gene Expression Associated Map, also referred as “DREAM” (http://bio-big-data.cn:8080/DREAM), is a manually curated database of experimentally supported protein-coding RNAs and drugs associations in human cancers. The current version of the DREAM documents 3048 entries about scientific literatures supported drug sensitivity or drug intervention related protein-coding RNAs from PubMed database and 195 high-throughput microarray data about drug sensitivity or drug intervention related protein-coding RNAs data from GEO database. Each entry in DREAM database contains detailed information on protein-coding RNA, drug, cancer, and other information including title, PubMed ID, journal, publish time. The DREAM database also provides some data visualization and online analysis services such as volcano plot, GO/KEGG enrichment function analysis, and novel drug discovery analysis. We hope the DREAM database should serve as a valuable resource for clinical practice and basic research, which could help researchers better understand the effects of protein-coding RNAs on drug response in human cancers.