Scoring functions for drug-effect similarity

A natural small molecule pinocembrin resists high-fat diet-induced obesity through GPR120-ERK1/2 pathway

Obesity is a widely concerned health problem. Mobilizing white adipose tissue and reducing fat synthesis are considered as effective strategies in the treatment of obesity. Here, using Connectivity Map (CMap) approach, we identified the pinocembrin (PB), a natural flavonoid primarily found in propolis, as a potential anti-obesity drug. Therefore, high-fat-diet (HFD) mice were randomly divided into two groups and fed a HFD or HFD with PB in this study. In vivo experiments showed that supplementation of PB reduced the body weight gain and ameliorated insulin resistance in HFD-induced mice. More importantly, PB did not cause side effect through detecting the levels of alanine transaminase (ALT), aspartate aminotransferase (AST), creatinine (CRE) and blood urea nitrogen (BUN) in serum of mice. Additionally, PB reduced expansion of white adipose tissue with upregulation of genes related lipolysis and downregulation of genes related lipogenesis. Furthermore, in vitro experiments revealed that PB treatment dose-dependently inhibited lipid droplet formation with upregulation of genes related lipolysis and downregulation of genes related lipogenesis. Molecular docking analysis combined with cellular thermal shift assay (CETSA) suggested that PB has a high affinity to the G protein-coupled receptor 120 (GPR120). Meanwhile, we confirmed that PB efficiently inhibited adipogenic differentiation of preadipocytes by directly binding to GPR120 and subsequently activating the downstream phosphorylation extracellular regulated kinase 1/2 (ERK1/2). Collectively, PB exerted anti-obesity effect through GPR120-ERK1/2 signaling pathway, providing a novel and promising natural drug for the treatment of obesity.

MOASL: Predicting drug mechanism of actions through similarity learning with transcriptomic signature

Understanding the mechanisms of actions (MOAs) of compounds is crucial in drug discovery. A common step in drug MOAs annotation is to query the dysregulated gene signatures induced by drugs in a reference library of pre-defined signatures. However, traditional similarity-based computational strategies face challenges when dealing with high-dimensional and noisy transcriptional signature data. To address this issue, we introduce MOASL (MOAs prediction via Similarity Learning), a novel approach that contrastive to learn similarity embeddings among signatures with shared MOAs automatically. We evaluated the accuracy of signature matching on various transcriptional activity score (TAS) datasets and individual cell lines by using MOASL. The results show MOASL achieved higher performance over several statistical and machine learning methods. Furthermore, we provided the rationale of our model by visualizing the signature annotation procedure. Using MOASL, the MOAs label of query signature could be conveniently defined by calculating the similarity between the query embedding and the reference embeddings. Finally, we applied MOASL to repurpose thousands of compounds as glucocorticoid receptor (GR) agonists, accurately identifying 8 out of the top 10 compounds. MOASL is conveniently accessible on GitHub at https://github.com/jianglikun/MOASL, empowering researchers and practitioners in the field of drug discovery to predict the MOAs of drug.

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.

Gene co-expression analyses of health(span) across multiple species

Health(span)-related gene clusters/modules were recently identified based on knowledge about the cross-species genetic basis of health, to interpret transcriptomic datasets describing health-related interventions. However, the cross-species comparison of health-related observations reveals a lot of heterogeneity, not least due to widely varying health(span) definitions and study designs, posing a challenge for the exploration of conserved healthspan modules and, specifically, their transfer across species. To improve the identification and exploration of conserved/transferable healthspan modules, here we apply an established workflow based on gene co-expression network analyses employing GEO/ArrayExpress data for human and animal models, and perform a comprehensive meta-study of the resulting modules related to health(span), yielding a small set of literature backed health(span) candidate genes. For each experiment, WGCNA (weighted gene correlation network analysis) was used to infer modules of genes which correlate in their expression with a 'health phenotype score' and to determine the most-connected (hub) genes (and their interactions) for each such module. After mapping these hub genes to their human orthologs, 12 health(span) genes were identified in at least two species (ACTN3, ANK1, MRPL18, MYL1, PAXIP1, PPP1CA, SCN3B, SDCBP, SKIV2L, TUBG1, TYROBP, WIPF1), for which enrichment analysis by g:profiler found an association with actin filament-based movement and associated organelles, as well as muscular structures. We conclude that a meta-study of hub genes from co-expression network analyses for the complex phenotype health(span), across multiple species, can yield molecular-mechanistic insights and can direct experimentalists to further investigate the contribution of individual genes and their interactions to health(span).

Integrated multi-omics approach revealed cellular senescence landscape

Cellular senescence is a complex multifactorial biological phenomenon that plays essential roles in aging, and aging-related diseases. During this process, the senescent cells undergo gene expression altering and chromatin structure remodeling. However, studies on the epigenetic landscape of senescence using integrated multi-omics approaches are limited. In this research, we performed ATAC-seq, RNA-seq and ChIP-seq on different senescent types to reveal the landscape of senescence and identify the prime regulatory elements. We also obtained 34 key genes and deduced that NAT1, PBX1 and RRM2, which interacted with each other, could be the potential markers of aging and aging-related diseases. In summary, our work provides the landscape to study accessibility dynamics and transcriptional regulations in cellular senescence. The application of this technique in different types of senescence allows us to identify the regulatory elements responsible for the substantial regulation of transcription, providing the insights into molecular mechanisms of senescence.

Identification of small molecules as novel anti-adipogenic compounds based on Connectivity Map

Several physiological and pathological conditions such as aging, obesity, diabetes, anorexia nervosa are associated with increased adipogenesis in the bone marrow. A lack of effective drugs hinder the improved treatment for aberrant accumulation of bone marrow adipocytes. Given the higher costs, longer duration and sometimes lack of efficacy in drug discovery, computational and experimental strategies have been used to identify previously approved drugs for the treatment of diseases, also known as drug repurposing. Here, we describe the method of small molecule-prioritization by employing adipocyte-specific genes using the connectivity map (CMap). We then generated transcriptomic profiles using human mesenchymal stromal cells under adipogenic differentiation with the treatment of prioritized compounds, and identified emetine and kinetin-riboside to have a potent inhibitory effect on adipogenesis. Overall, we demonstrated a proof-of-concept method to identify repurposable drugs capable of inhibiting adipogenesis, using the Connectivity Map.

Integrated bioinformatic analysis reveals the underlying molecular mechanism of and potential drugs for pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a devastating cardiovascular disease without a clear mechanism or drugs for treatment. Therefore, it is crucial to reveal the underlying molecular mechanism and identify potential drugs for PAH. In this study, we first integrated three human lung tissue datasets (GSE113439, GSE53408, GSE117261) from GEO. A total of 151 differentially expressed genes (DEGs) were screened, followed by KEGG and GO enrichment analyses and PPI network construction. Five hub genes (CSF3R, NT5E, ANGPT2, FGF7, and CXCL9) were identified by Cytoscape (Cytohubba). GSEA and GSVA were performed for each hub gene to uncover the potential mechanism. Moreover, to repurpose known and therapeutic drugs, the CMap database was retrieved, and nine candidate compounds (lypressin, ruxolitinib, triclabendazole, L-BSO, tiaprofenic acid, AT-9283, QL-X-138, huperzine-a, and L-741742) with a high level of confidence were obtained. Then ruxolitinib was selected to perform molecular docking simulations with ANGPT2, FGF7, NT5E, CSF3R, JAK1, JAK2, JAK3, TYK2. A certain concentration of ruxolitinib could inhibit the proliferation and migration of rat pulmonary artery smooth muscle cells (rPASMCs) in vitro. Together, these analyses principally identified CSF3R, NT5E, ANGPT2, FGF7 and CXCL9 as candidate biomarkers of PAH, and ruxolitinib might exert promising therapeutic action for PAH.

Single-cell analyses of aging, inflammation and senescence

Single-cell gene expression (transcriptomics) data are becoming robust and abundant, and are increasingly used to track organisms along their life-course. This allows investigation into how aging affects cellular transcriptomes, and how changes in transcriptomes may underlie aging, including chronic inflammation (inflammaging), immunosenescence and cellular senescence. We compiled and tabulated aging-related single-cell datasets published to date, collected and discussed relevant findings, and inspected some of these datasets ourselves. We specifically note insights that cannot (or not easily) be based on bulk data. For example, in some datasets, the fraction of cells expressing p16 (CDKN2A), one of the most prominent markers of cellular senescence, was reported to increase, in addition to its upregulated mean expression over all cells. Moreover, we found evidence for inflammatory processes in most datasets, some of these driven by specific cells of the immune system. Further, single-cell data are specifically useful to investigate whether transcriptional heterogeneity (also called noise or variability) increases with age, and many (but not all) studies in our review report an increase in such heterogeneity. Finally, we demonstrate some stability of marker gene expression patterns across closely similar studies and suggest that single-cell experiments may hold the key to provide detailed insights whenever interventions (countering aging, inflammation, senescence, disease, etc.) are affecting cells depending on cell type.

Drug-drug similarity measure and its applications

Drug similarities play an important role in modern biology and medicine, as they help scientists gain deep insights into drugs' therapeutic mechanisms and conduct wet labs that may significantly improve the efficiency of drug research and development. Nowadays, a number of drug-related databases have been constructed, with which many methods have been developed for computing similarities between drugs for studying associations between drugs, human diseases, proteins (drug targets) and more. In this review, firstly, we briefly introduce the publicly available drug-related databases. Secondly, based on different drug features, interaction relationships and multimodal data, we summarize similarity calculation methods in details. Then, we discuss the applications of drug similarities in various biological and medical areas. Finally, we evaluate drug similarity calculation methods with common evaluation metrics to illustrate the important roles of drug similarity measures on different applications.

Suppression of the TGF-β pathway by a macrolide antibiotic decreases fibrotic responses by ocular fibroblasts in vitro

To elucidate and to inhibit post-surgical fibrotic processes after trabeculectomy in glaucoma therapy, we measured gene expression in a fibrotic cell culture model, based on transforming growth factor TGF-β induction in primary human tenon fibroblasts (hTFs), and used Connectivity Map (CMap) data for drug repositioning. We found that specific molecular mechanisms behind fibrosis are the upregulation of actins, the downregulation of CD34, and the upregulation of inflammatory cytokines such as IL6, IL11 and BMP6. The macrolide antibiotic Josamycin (JM) reverses these molecular mechanisms according to data from the CMap, and we thus tested JM as an inhibitor of fibrosis. JM was first tested for its toxic effects on hTFs, where it showed no influence on cell viability, but inhibited hTF proliferation in a concentration-dependent manner. We then demonstrated that JM suppresses the synthesis of extracellular matrix (ECM) components. In hTFs stimulated with TGF-β1, JM specifically inhibited α-smooth muslce actin expression, suggesting that it inhibits the transformation of fibroblasts into fibrotic myofibroblasts. In addition, a decrease of components of the ECM such as fibronectin, which is involved in in vivo scarring, was observed. We conclude that JM may be a promising candidate for the treatment of fibrosis after glaucoma filtration surgery or drainage device implantation in vivo.