Predicting protein targets for drug-like compounds using transcriptomics

Polypharmacology prediction: the long road toward comprehensively anticipating small-molecule selectivity to de-risk drug discovery

INTRODUCTION

Small molecules often bind to multiple targets, a behavior termed polypharmacology. Anticipating polypharmacology is essential for drug discovery since unknown off-targets can modulate safety and efficacy - profoundly affecting drug discovery success. Unfortunately, experimental methods to assess selectivity present significant limitations and drugs still fail in the clinic due to unanticipated off-targets. Computational methods are a cost-effective, complementary approach to predict polypharmacology.

AREAS COVERED

This review aims to provide a comprehensive overview of the state of polypharmacology prediction and discuss its strengths and limitations, covering both classical cheminformatics methods and bioinformatic approaches. The authors review available data sources, paying close attention to their different coverage. The authors then discuss major algorithms grouped by the types of data that they exploit using selected examples.

EXPERT OPINION

Polypharmacology prediction has made impressive progress over the last decades and contributed to identify many off-targets. However, data incompleteness currently limits most approaches to comprehensively predict selectivity. Moreover, our limited agreement on model assessment challenges the identification of the best algorithms - which at present show modest performance in prospective real-world applications. Despite these limitations, the exponential increase of multidisciplinary Big Data and AI hold much potential to better polypharmacology prediction and de-risk drug discovery.

Systematic computational strategies for identifying protein targets and lead discovery

Computational algorithms and tools have retrenched the drug discovery and development timeline. The applicability of computational approaches has gained immense relevance owing to the dramatic surge in the structural information of biomacromolecules and their heteromolecular complexes. Computational methods are now extensively used in identifying new protein targets, druggability assessment, pharmacophore mapping, molecular docking, the virtual screening of lead molecules, bioactivity prediction, molecular dynamics of protein-ligand complexes, affinity prediction, and for designing better ligands. Herein, we provide an overview of salient components of recently reported computational drug-discovery workflows that includes algorithms, tools, and databases for protein target identification and optimized ligand selection.

Deep representation learning of chemical-induced transcriptional profile for phenotype-based drug discovery

Artificial intelligence transforms drug discovery, with phenotype-based approaches emerging as a promising alternative to target-based methods, overcoming limitations like lack of well-defined targets. While chemical-induced transcriptional profiles offer a comprehensive view of drug mechanisms, inherent noise often obscures the true signal, hindering their potential for meaningful insights. Here, we highlight the development of TranSiGen, a deep generative model employing self-supervised representation learning. TranSiGen analyzes basal cell gene expression and molecular structures to reconstruct chemical-induced transcriptional profiles with high accuracy. By capturing both cellular and compound information, TranSiGen-derived representations demonstrate efficacy in diverse downstream tasks like ligand-based virtual screening, drug response prediction, and phenotype-based drug repurposing. Notably, in vitro validation of TranSiGen's application in pancreatic cancer drug discovery highlights its potential for identifying effective compounds. We envisage that integrating TranSiGen into the drug discovery and mechanism research holds significant promise for advancing biomedicine.

Towards Precision Ecotoxicology: Leveraging Evolutionary Conservation of Pharmaceutical and Personal Care Product Targets to Understand Adverse Outcomes Across Species and Life Stages

Translation of environmental science to the practice aims to protect biodiversity and ecosystem services, and our future ability to do so relies on the development of a precision ecotoxicology approach wherein we leverage the genetics and informatics of species to better understand and manage the risks of global pollution. A little over a decade ago, a workshop focusing on the risks of pharmaceuticals and personal care products (PPCPs) in the environment identified a priority research question, "What can be learned about the evolutionary conservation of PPCP targets across species and life stages in the context of potential adverse outcomes and effects?" We review the activities in this area over the past decade, consider prospects of more recent developments, and identify future research needs to develop next-generation approaches for PPCPs and other global chemicals and waste challenges. Environ Toxicol Chem 2024;43:526-536. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Drug target prediction through deep learning functional representation of gene signatures

Many machine learning applications in bioinformatics currently rely on matching gene identities when analyzing input gene signatures and fail to take advantage of preexisting knowledge about gene functions. To further enable comparative analysis of OMICS datasets, including target deconvolution and mechanism of action studies, we develop an approach that represents gene signatures projected onto their biological functions, instead of their identities, similar to how the word2vec technique works in natural language processing. We develop the Functional Representation of Gene Signatures (FRoGS) approach by training a deep learning model and demonstrate that its application to the Broad Institute's L1000 datasets results in more effective compound-target predictions than models based on gene identities alone. By integrating additional pharmacological activity data sources, FRoGS significantly increases the number of high-quality compound-target predictions relative to existing approaches, many of which are supported by in silico and/or experimental evidence. These results underscore the general utility of FRoGS in machine learning-based bioinformatics applications. Prediction networks pre-equipped with the knowledge of gene functions may help uncover new relationships among gene signatures acquired by large-scale OMICs studies on compounds, cell types, disease models, and patient cohorts.

An overview of recent advances and challenges in predicting compound-protein interaction (CPI)

Compound-protein interactions (CPIs) are critical in drug discovery for identifying therapeutic targets, drug side effects, and repurposing existing drugs. Machine learning (ML) algorithms have emerged as powerful tools for CPI prediction, offering notable advantages in cost-effectiveness and efficiency. This review provides an overview of recent advances in both structure-based and non-structure-based CPI prediction ML models, highlighting their performance and achievements. It also offers insights into CPI prediction-related datasets and evaluation benchmarks. Lastly, the article presents a comprehensive assessment of the current landscape of CPI prediction, elucidating the challenges faced and outlining emerging trends to advance the field.

On the correspondence between the transcriptomic response of a compound and its effects on its targets

Better understanding the transcriptomic response produced by a compound perturbing its targets can shed light on the underlying biological processes regulated by the compound. However, establishing the relationship between the induced transcriptomic response and the target of a compound is non-trivial, partly because targets are rarely differentially expressed. Therefore, connecting both modalities requires orthogonal information (e.g., pathway or functional information). Here, we present a comprehensive study aimed at exploring this relationship by leveraging thousands of transcriptomic experiments and target data for over 2000 compounds. Firstly, we confirm that compound-target information does not correlate as expected with the transcriptomic signatures induced by a compound. However, we reveal how the concordance between both modalities increases by connecting pathway and target information. Additionally, we investigate whether compounds that target the same proteins induce a similar transcriptomic response and conversely, whether compounds with similar transcriptomic responses share the same target proteins. While our findings suggest that this is generally not the case, we did observe that compounds with similar transcriptomic profiles are more likely to share at least one protein target and common therapeutic applications. Finally, we demonstrate how to exploit the relationship between both modalities for mechanism of action deconvolution by presenting a case scenario involving a few compound pairs with high similarity.

Benchmarking causal reasoning algorithms for gene expression-based compound mechanism of action analysis

BACKGROUND

Elucidating compound mechanism of action (MoA) is beneficial to drug discovery, but in practice often represents a significant challenge. Causal Reasoning approaches aim to address this situation by inferring dysregulated signalling proteins using transcriptomics data and biological networks; however, a comprehensive benchmarking of such approaches has not yet been reported. Here we benchmarked four causal reasoning algorithms (SigNet, CausalR, CausalR ScanR and CARNIVAL) with four networks (the smaller Omnipath network vs. 3 larger MetaBase™ networks), using LINCS L1000 and CMap microarray data, and assessed to what extent each factor dictated the successful recovery of direct targets and compound-associated signalling pathways in a benchmark dataset comprising 269 compounds. We additionally examined impact on performance in terms of the functions and roles of protein targets and their connectivity bias in the prior knowledge networks.

RESULTS

According to statistical analysis (negative binomial model), the combination of algorithm and network most significantly dictated the performance of causal reasoning algorithms, with the SigNet recovering the greatest number of direct targets. With respect to the recovery of signalling pathways, CARNIVAL with the Omnipath network was able to recover the most informative pathways containing compound targets, based on the Reactome pathway hierarchy. Additionally, CARNIVAL, SigNet and CausalR ScanR all outperformed baseline gene expression pathway enrichment results. We found no significant difference in performance between L1000 data or microarray data, even when limited to just 978 'landmark' genes. Notably, all causal reasoning algorithms also outperformed pathway recovery based on input DEGs, despite these often being used for pathway enrichment. Causal reasoning methods performance was somewhat correlated with connectivity and biological role of the targets.

CONCLUSIONS

Overall, we conclude that causal reasoning performs well at recovering signalling proteins related to compound MoA upstream from gene expression changes by leveraging prior knowledge networks, and that the choice of network and algorithm has a profound impact on the performance of causal reasoning algorithms. Based on the analyses presented here this is true for both microarray-based gene expression data as well as those based on the L1000 platform.

Mechanism of action deconvolution of the small-molecule pathological tau aggregation inhibitor Anle138b

BACKGROUND

A key histopathological hallmark of Alzheimer's disease (AD) is the presence of neurofibrillary tangles of aggregated microtubule-associated protein tau in neurons. Anle138b is a small molecule which has previously shown efficacy in mice in reducing tau aggregates and rescuing AD disease phenotypes.

METHODS

In this work, we employed bioinformatics analysis-including pathway enrichment and causal reasoning-of an in vitro tauopathy model. The model consisted of cultured rat cortical neurons either unseeded or seeded with tau aggregates derived from human AD patients, both of which were treated with Anle138b to generate hypotheses for its mode of action. In parallel, we used a collection of human target prediction models to predict direct targets of Anle138b based on its chemical structure.

RESULTS

Combining the different approaches, we found evidence supporting the hypothesis that the action of Anle138b involves several processes which are key to AD progression, including cholesterol homeostasis and neuroinflammation. On the pathway level, we found significantly enriched pathways related to these two processes including those entitled "Superpathway of cholesterol biosynthesis" and "Granulocyte adhesion and diapedesis". With causal reasoning, we inferred differential activity of SREBF1/2 (involved in cholesterol regulation) and mediators of the inflammatory response such as NFKB1 and RELA. Notably, our findings were also observed in Anle138b-treated unseeded neurons, meaning that the inferred processes are independent of tau pathology and thus represent the direct action of the compound in the cellular system. Through structure-based ligand-target prediction, we predicted the intracellular cholesterol carrier NPC1 as well as NF-κB subunits as potential targets of Anle138b, with structurally similar compounds in the model training set known to target the same proteins.

CONCLUSIONS

This study has generated feasible hypotheses for the potential mechanism of action of Anle138b, which will enable the development of future molecular interventions aiming to reduce tau pathology in AD patients.

Computational drug repurposing by exploiting large-scale gene expression data: Strategy, methods and applications

De novo drug development is an extremely complex, time-consuming and costly task. Urgent needs for therapies of various diseases have greatly accelerated searches for more effective drug development methods. Luckily, drug repurposing provides a new and effective perspective on disease treatment. Rapidly increased large-scale transcriptome data paints a detailed prospect of gene expression during disease onset and thus has received wide attention in the field of computational drug repurposing. However, how to efficiently mine transcriptome data and identify new indications for old drugs remains a critical challenge. This review discussed the irreplaceable role of transcriptome data in computational drug repurposing and summarized some representative databases, tools and strategies. More importantly, it proposed a practical guideline through establishing the correspondence between three gene expression data types and five strategies, which would facilitate researchers to adopt appropriate strategies to deeply mine large-scale transcriptome data and discover more effective therapies.

Evolution of precision oncology-guided treatment paradigms

Cancer treatment is gradually evolving from the classical use of nonspecific cytotoxic drugs targeting generic mechanisms of cell growth and proliferation. Instead, new "patient-specific treatment paradigms" that are based on an individual patient's tumor-specific molecular features are emerging, and these include "druggable" genomic alterations such as oncogenic driver mutations, downstream activities of cancer-signaling pathways, and the expression of specific genes involved in tumorigenesis and cancer progression. This evolving landscape of making evidence-based treatment decisions forms the foundation of precision oncology, which aims to deliver "the right drug, to the right patient and at the right time". The long-term vision for this approach is to maximize the treatment efficacy while minimizing exposure to ineffective therapy and reducing co-morbidity-related side effects. Successful clinical translation and implementation of this vision have the potential to revolutionize treatment paradigms from predominantly reactive, to more evidence-based, proactive and predictive care. In this article, we review the past and current approaches in precision oncology, and describe their remarkable power and limitations. We also speculate on the evolution of newly emerging methodologies of the future that can be used to address some of the key challenges associated with the existing paradigms. This article is categorized under: Cancer > Genetics/Genomics/Epigenetics Cancer > Molecular and Cellular Physiology Cancer > Computational Models.

How to Design Peptides

Novel design of proteins to target receptors for treatment or tissue augmentation has come to the fore owing to advancements in computing power, modeling frameworks, and translational successes. Shorter proteins, or peptides, can offer combinatorial synergies with dendrimer, polymer, or other peptide carriers for enhanced local signaling, which larger proteins may sterically hinder. Here, we present a generalized method for designing a novel peptide. We first show how to create a script protocol that can be used to iteratively optimize and screen novel peptide sequences for binding a target protein. We present a step-by-step introduction to utilizing file repositories, data bases, and the Rosetta software suite. RosettaScripts, an .xml interface that allows for sequential functions to be performed, is used to order the functions for repeatable performance. These strategies may lead to more groups venturing into computational design, which may result in synergies from artificial intelligence/machine learning (AI/ML) to phage display and screening. Importantly, the beginner is expected to be able to design their first peptide ligand and begin their journey in peptide drug discovery. Generally, these peptides potentially could be used to interact with any enzyme or receptor, for example, in the study of chemokines and their interactions with glycosoaminoglycans and their receptors.

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.

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

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

RNA aptamers specific for transmembrane p24 trafficking protein 6 and Clusterin for the targeted delivery of imaging reagents and RNA therapeutics to human β cells

The ability to detect and target β cells in vivo can substantially refine how diabetes is studied and treated. However, the lack of specific probes still hampers a precise characterization of human β cell mass and the delivery of therapeutics in clinical settings. Here, we report the identification of two RNA aptamers that specifically and selectively recognize mouse and human β cells. The putative targets of the two aptamers are transmembrane p24 trafficking protein 6 (TMED6) and clusterin (CLUS). When given systemically in immune deficient mice, these aptamers recognize the human islet graft producing a fluorescent signal proportional to the number of human islets transplanted. These aptamers cross-react with endogenous mouse β cells and allow monitoring the rejection of mouse islet allografts. Finally, once conjugated to saRNA specific for X-linked inhibitor of apoptosis (XIAP), they can efficiently transfect non-dissociated human islets, prevent early graft loss, and improve the efficacy of human islet transplantation in immunodeficient in mice.

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.

Computational tools for analyzing single-cell data in pluripotent cell differentiation studies

Single-cell technologies are revolutionizing the ability of researchers to infer the causes and results of biological processes. Although several studies of pluripotent cell differentiation have recently utilized single-cell sequencing data, other aspects related to the optimization of differentiation protocols, their validation, robustness, and usage are still not taking full advantage of single-cell technologies. In this review, we focus on computational approaches for the analysis of single-cell omics and imaging data and discuss their use to address many of the major challenges involved in the development, validation, and use of cells obtained from pluripotent cell differentiation.

TRIOMPHE: Transcriptome-Based Inference and Generation of Molecules with Desired Phenotypes by Machine Learning

One of the most challenging tasks in the drug-discovery process is the efficient identification of small molecules with desired phenotypes. In this study, we propose a novel computational method for omics-based de novo drug design, which we call TRIOMPHE (transcriptome-based inference and generation of molecules with desired phenotypes). We investigated the correlation between chemically induced transcriptome profiles (reflecting cellular responses to compound treatment) and genetically perturbed transcriptome profiles (reflecting cellular responses to gene knock-down or gene overexpression of target proteins) in terms of ligand-target interactions. Subsequently, we developed novel machine learning methods to generate the chemical structures of new molecules with desired transcriptome profiles in the framework of a variational autoencoder. The use of desired transcriptome profiles enables the automatic design of molecules that are likely to have bioactivities for target proteins of interest. We showed that our methods can generate chemically valid molecules that are likely to have biological activities on 10 target proteins; moreover, they can outperform previous methods that had the same objective. Our omics-based structure generator is expected to be useful for the de novo design of drugs for a variety of target proteins.

Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: The Multi-Targeting Drug DREAM Challenge

A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ('polypharmacology'). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Novel FOXM1 inhibitor identified via gene network analysis induces autophagic FOXM1 degradation to overcome chemoresistance of human cancer cells

FOXM1 transcription factor is an oncogene and a master regulator of chemoresistance in multiple cancers. Pharmacological inhibition of FOXM1 is a promising approach but has proven to be challenging. We performed a network-centric transcriptomic analysis to identify a novel compound STL427944 that selectively suppresses FOXM1 by inducing the relocalization of nuclear FOXM1 protein to the cytoplasm and promoting its subsequent degradation by autophagosomes. Human cancer cells treated with STL427944 exhibit increased sensitivity to cytotoxic effects of conventional chemotherapeutic treatments (platinum-based agents, 5-fluorouracil, and taxanes). RNA-seq analysis of STL427944-induced gene expression changes revealed prominent suppression of gene signatures characteristic for FOXM1 and its downstream targets but no significant changes in other important regulatory pathways, thereby suggesting high selectivity of STL427944 toward the FOXM1 pathway. Collectively, the novel autophagy-dependent mode of FOXM1 suppression by STL427944 validates a unique pathway to overcome tumor chemoresistance and improve the efficacy of treatment with conventional cancer drugs.

Machine Learning Uses Chemo-Transcriptomic Profiles to Stratify Antimalarial Compounds With Similar Mode of Action

The rapid development of antimalarial resistance motivates the continued search for novel compounds with a mode of action (MoA) different to current antimalarials. Phenotypic screening has delivered thousands of promising hit compounds without prior knowledge of the compounds’ exact target or MoA. Whilst the latter is not initially required to progress a compound in a medicinal chemistry program, identifying the MoA early can accelerate hit prioritization, hit-to-lead optimization and preclinical combination studies in malaria research. The effects of drug treatment on a cell can be observed on systems level in changes in the transcriptome, proteome and metabolome. Machine learning (ML) algorithms are powerful tools able to deconvolute such complex chemically-induced transcriptional signatures to identify pathways on which a compound act and in this manner provide an indication of the MoA of a compound. In this study, we assessed different ML approaches for their ability to stratify antimalarial compounds based on varied chemically-induced transcriptional responses. We developed a rational gene selection approach that could identify predictive features for MoA to train and generate ML models. The best performing model could stratify compounds with similar MoA with a classification accuracy of 76.6 ± 6.4%. Moreover, only a limited set of 50 biomarkers was required to stratify compounds with similar MoA and define chemo-transcriptomic fingerprints for each compound. These fingerprints were unique for each compound and compounds with similar targets/MoA clustered together. The ML model was specific and sensitive enough to group new compounds into MoAs associated with their predicted target and was robust enough to be extended to also generate chemo-transcriptomic fingerprints for additional life cycle stages like immature gametocytes. This work therefore contributes a new strategy to rapidly, specifically and sensitively indicate the MoA of compounds based on chemo-transcriptomic fingerprints and holds promise to accelerate antimalarial drug discovery programs.

Scoring functions for drug-effect similarity

MOTIVATION

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

RESULTS

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

AVAILABILITY

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

CONTACT

steffen.moeller@uni-rostock.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Briefings in Bioinformatics online.

Drug-Target Interaction Prediction Based on Adversarial Bayesian Personalized Ranking

The prediction of drug-target interaction (DTI) is a key step in drug repositioning. In recent years, many studies have tried to use matrix factorization to predict DTI, but they only use known DTIs and ignore the features of drug and target expression profiles, resulting in limited prediction performance. In this study, we propose a new DTI prediction model named AdvB-DTI. Within this model, the features of drug and target expression profiles are associated with Adversarial Bayesian Personalized Ranking through matrix factorization. Firstly, according to the known drug-target relationships, a set of ternary partial order relationships is generated. Next, these partial order relationships are used to train the latent factor matrix of drugs and targets using the Adversarial Bayesian Personalized Ranking method, and the matrix factorization is improved by the features of drug and target expression profiles. Finally, the scores of drug-target pairs are achieved by the inner product of latent factors, and the DTI prediction is performed based on the score ranking. The proposed model effectively takes advantage of the idea of learning to rank to overcome the problem of data sparsity, and perturbation factors are introduced to make the model more robust. Experimental results show that our model could achieve a better DTI prediction performance.

Enhancing Chemogenomics with Predictive Pharmacology

One of the grand challenges in contemporary chemical biology is the generation of a probe for every member of the human proteome. Probe selection and optimization strategies typically rely on experimental bioactivity data to determine the potency and selectivity of candidate molecules. However, this approach is profoundly limited by the sparsity of the known data, the annotation bias often found in the literature, and the cost of physical screening. Recent advancements in predictive pharmacology, such as the application of multitask and transfer learning, as well as the use of biologically motivated, structure-agnostic features to characterize molecules, should serve to mitigate these issues. Computational modeling likely offers the only cost-effective approach to substantially increasing the bioactivity annotation density both on the local and global scale and thus, we argue, will need to make a substantial contribution if the ambitious goals of probing the human proteome are to be realized in the foreseeable future.

Learning to encode cellular responses to systematic perturbations with deep generative models

Cellular signaling systems play a vital role in maintaining homeostasis when a cell is exposed to different perturbations. Components of the systems are organized as hierarchical networks, and perturbing different components often leads to transcriptomic profiles that exhibit compositional statistical patterns. Mining such patterns to investigate how cellular signals are encoded is an important problem in systems biology, where artificial intelligence techniques can be of great assistance. Here, we investigated the capability of deep generative models (DGMs) to modeling signaling systems and learn representations of cellular states underlying transcriptomic responses to diverse perturbations. Specifically, we show that the variational autoencoder and the supervised vector-quantized variational autoencoder can accurately regenerate gene expression data in response to perturbagen treatments. The models can learn representations that reveal the relationships between different classes of perturbagens and enable mappings between drugs and their target genes. In summary, DGMs can adequately learn and depict how cellular signals are encoded. The resulting representations have broad applications, demonstrating the power of artificial intelligence in systems biology and precision medicine.

Rapid determination of nematode cell and organ susceptibility to toxic treatments

In research focused on the intestine of parasitic nematodes, we recently identified small molecule inhibitors toxic to intestinal cells of larval Ascaris suum (nematode intestinal toxins/toxicants; “NITs”). Some NITs had anthelmintic activity across the phylogenetic diversity of the Nematoda. The whole-worm motility inhibition assay quantified anthelmintic activity, but worm responses to NITs in relation to pathology or affected molecular pathways was not acquired. In this study we extended this research to more comprehensively determine in whole larval A. suum the cells, organ systems, molecular targets, and potential cellular pathways involved in mechanisms of toxicity leading to cell death. The experimental system utilized fluorescent nuclear probes (bisbenzimide, propidium iodide), NITs, an A. suum larval parasite culture system and transcriptional responses (RNA-seq) to NITs. The approach provides for rapid resolution of NIT-induced cell death among organ systems (e.g. intestine, excretory, esophagus, hypodermis and seam cells, and nervous), discriminates among NITs based on cell death profiles, and identifies cells and organ systems with the greatest NIT sensitivity (e.g. intestine and apparent neuronal cells adjacent to the nerve ring). Application was extended to identify cells and organs sensitive to several existing anthelmintics. This approach also resolved intestinal cell death and irreparable damage induced in adult A. suum by two NITs, establishing a new model to elucidate relevant pathologic mechanisms in adult worms. RNA-seq analysis resolved A. suum genes responsive to treatments with three NITs, identifying dihydroorotate dehydrogenase (uridine synthesis) and RAB GTPase(s) (vesicle transport) as potential targets/pathways leading to cell death. A set of genes induced by all three NITs tested suggest common stress or survival responses activated by NITs. Beyond the presented specific lines of research, elements of the overall experimental system presented in this study have broad application toward systematic development of new anthelmintics.

•

A unique rapid cell death assay was developed for parasitic nematodes.

•

Multiple drug-like molecules cause widespread cell death in many organs of A. suum.

•

Multiple cell and organ systems were validated as targets for anthelmintics.

•

Potential drug targets/pathways were implicated in activating cell death processes.

Universal Nature of Drug Treatment Responses in Drug-Tissue-Wide Model-Animal Experiments Using Tensor Decomposition-Based Unsupervised Feature Extraction

Gene expression profiles of tissues treated with drugs have recently been used to infer clinical outcomes. Although this method is often successful from the application point of view, gene expression altered by drugs is rarely analyzed in detail, because of the extremely large number of genes involved. Here, we applied tensor decomposition (TD)-based unsupervised feature extraction (FE) to the gene expression profiles of 24 mouse tissues treated with 15 drugs. TD-based unsupervised FE enabled identification of the common effects of 15 drugs including an interesting universal feature: these drugs affect genes in a gene-group-wide manner and were dependent on three tissue types (neuronal, muscular, and gastroenterological). For each tissue group, TD-based unsupervised FE enabled identification of a few tens to a few hundreds of genes affected by the drug treatment. These genes are distinctly expressed between drug treatments and controls as well as between tissues in individual tissue groups and other tissues. We also validated the assignment of genes to individual tissue groups using multiple enrichment analyses. We conclude that TD-based unsupervised FE is a promising method for integrated analysis of gene expression profiles from multiple tissues treated with multiple drugs in a completely unsupervised manner.

A systematic strategy for the investigation of vaccines and drugs targeting bacteria

Infectious and epidemic diseases induced by bacteria have historically caused great distress to people, and have even resulted in a large number of deaths worldwide. At present, many researchers are working on the discovery of viable drug and vaccine targets for bacteria through multiple methods, including the analyses of comparative subtractive genome, core genome, replication-related proteins, transcriptomics and riboswitches, which plays a significant part in the treatment of infectious and pandemic diseases. The 3D structures of the desired target proteins, drugs and epitopes can be predicted and modeled through target analysis. Meanwhile, molecular dynamics (MD) analysis of the constructed drug/epitope-protein complexes is an important standard for testing the suitability of these screened drugs and vaccines. Currently, target discovery, target analysis and MD analysis are integrated into a systematic set of drug and vaccine analysis strategy for bacteria. We hope that this comprehensive strategy will help in the design of high-performance vaccines and drugs.

Targeting the C-Terminal Domain Small Phosphatase 1

The human C-terminal domain small phosphatase 1 (CTDSP1/SCP1) is a protein phosphatase with a conserved catalytic site of DXDXT/V. CTDSP1’s major activity has been identified as dephosphorylation of the 5th Ser residue of the tandem heptad repeat of the RNA polymerase II C-terminal domain (RNAP II CTD). It is also implicated in various pivotal biological activities, such as acting as a driving factor in repressor element 1 (RE-1)-silencing transcription factor (REST) complex, which silences the neuronal genes in non-neuronal cells, G1/S phase transition, and osteoblast differentiation. Recent findings have denoted that negative regulation of CTDSP1 results in suppression of cancer invasion in neuroglioma cells. Several researchers have focused on the development of regulating materials of CTDSP1, due to the significant roles it has in various biological activities. In this review, we focused on this emerging target and explored the biological significance, challenges, and opportunities in targeting CTDSP1 from a drug designing perspective.

Exploring the Use of Compound-Induced Transcriptomic Data Generated From Cell Lines to Predict Compound Activity Toward Molecular Targets

Pharmaceutical or phytopharmaceutical molecules rely on the interaction with one or more specific molecular targets to induce their anticipated biological responses. Nonetheless, these compounds are also prone to interact with many other non-intended biological targets, also known as off-targets. Unfortunately, off-target identification is difficult and expensive. Consequently, QSAR models predicting the activity on a target have gained importance in drug discovery or in the de-risking of chemicals. However, a restricted number of targets are well characterized and hold enough data to build such in silico models. A good alternative to individual target evaluations is to use integrative evaluations such as transcriptomics obtained from compound-induced gene expression measurements derived from cell cultures. The advantage of these particular experiments is to capture the consequences of the interaction of compounds on many possible molecular targets and biological pathways, without having any constraints concerning the chemical space. In this work, we assessed the value of a large public dataset of compound-induced transcriptomic data, to predict compound activity on a selection of 69 molecular targets. We compared such descriptors with other QSAR descriptors, namely the Morgan fingerprints (similar to extended-connectivity fingerprints). Depending on the target, active compounds could show similar signatures in one or multiple cell lines, whether these active compounds shared similar or different chemical structures. Random forest models using gene expression signatures were able to perform similarly or better than counterpart models built with Morgan fingerprints for 25% of the target prediction tasks. These performances occurred mostly using signatures produced in cell lines showing similar signatures for active compounds toward the considered target. We show that compound-induced transcriptomic data could represent a great opportunity for target prediction, allowing to overcome the chemical space limitation of QSAR models.

Whole-Cell Phenotypic Screening of Medicines for Malaria Venture Pathogen Box Identifies Specific Inhibitors of Plasmodium falciparum Late-Stage Development and Egress

We report a systematic, cellular phenotype-based antimalarial screening of the Medicines for Malaria Venture Pathogen Box collection, which facilitated the identification of specific blockers of late-stage intraerythrocytic development of Plasmodium falciparum First, from standard growth inhibition assays, we identified 173 molecules with antimalarial activity (50% effective concentration [EC₅₀] ≤ 10 μM), which included 62 additional molecules over previously known antimalarial candidates from the Pathogen Box. We identified 90 molecules with EC₅₀ of ≤1 μM, which had significant effect on the ring-trophozoite transition, while 9 molecules inhibited the trophozoite-schizont transition and 21 molecules inhibited the schizont-ring transition (with ≥50% parasites failing to proceed to the next stage) at 1 μM. We therefore rescreened all 173 molecules and validated hits in microscopy to prioritize 12 hits as selective blockers of the schizont-ring transition. Seven of these molecules inhibited the calcium ionophore-induced egress of Toxoplasma gondii, a related apicomplexan parasite, suggesting that the inhibitors may be acting via a conserved mechanism which could be further exploited for target identification studies. We demonstrate that two molecules, MMV020670 and MMV026356, identified as schizont inhibitors in our screens, induce the fragmentation of DNA in merozoites, thereby impairing their ability to egress and invade. Further mechanistic studies would facilitate the therapeutic exploitation of these molecules as broadly active inhibitors targeting late-stage development and egress of apicomplexan parasites relevant to human health.

Tryptophan scanning mutagenesis as a way to mimic the compound-bound state and probe the selectivity of allosteric inhibitors in cells

Understanding the selectivity of a small molecule for its target(s) in cells is an important goal in chemical biology and drug discovery. One powerful way to address this question is with dominant negative (DN) mutants, in which an active site residue in the putative target is mutated. While powerful, this approach is less straightforward for allosteric sites. Here, we introduce tryptophan scanning mutagenesis as an expansion of this idea. As a test case, we focused on the challenging drug target, heat shock cognate protein 70 (Hsc70), and its allosteric inhibitor JG-98. Structure-based modelling predicted that mutating Y149W in human Hsc70 or Y145W in the bacterial ortholog DnaK would place an indole side chain into the allosteric pocket normally occupied by the compound. Indeed, we found that the tryptophan mutants acted as if they were engaged with JG-98. We then used DnaK Y145W to suggest that this protein may be an anti-bacterial target. Indeed, we found that DnaK inhibitors have minimum inhibitory concentration (MIC) values <0.125 μg mL-1 against several pathogens, including multidrug-resistant Staphylococcus aureus (MRSA) strains. We propose that tryptophan scanning mutagenesis may provide a distinct way to address the important problem of target engagement.

Retro-2 protects cells from ricin toxicity by inhibiting ASNA1-mediated ER targeting and insertion of tail-anchored proteins

The small molecule Retro-2 prevents ricin toxicity through a poorly-defined mechanism of action (MOA), which involves halting retrograde vesicle transport to the endoplasmic reticulum (ER). CRISPRi genetic interaction analysis revealed Retro-2 activity resembles disruption of the transmembrane domain recognition complex (TRC) pathway, which mediates post-translational ER-targeting and insertion of tail-anchored (TA) proteins, including SNAREs required for retrograde transport. Cell-based and in vitro assays show that Retro-2 blocks delivery of newly-synthesized TA-proteins to the ER-targeting factor ASNA1 (TRC40). An ASNA1 point mutant identified using CRISPR-mediated mutagenesis abolishes both the cytoprotective effect of Retro-2 against ricin and its inhibitory effect on ASNA1-mediated ER-targeting. Together, our work explains how Retro-2 prevents retrograde trafficking of toxins by inhibiting TA-protein targeting, describes a general CRISPR strategy for predicting the MOA of small molecules, and paves the way for drugging the TRC pathway to treat broad classes of viruses known to be inhibited by Retro-2.

Considerations for Strategic Use of High-Throughput Transcriptomics Chemical Screening Data in Regulatory Decisions

Recently, numerous organizations, including governmental regulatory agencies in the U.S. and abroad, have proposed using data from New Approach Methodologies (NAMs) for augmenting and increasing the pace of chemical assessments. NAMs are broadly defined as any technology, methodology, approach or combination thereof that can be used to provide information on chemical hazard and risk assessment that avoids the use of intact animals. High-throughput transcriptomics (HTTr) is a type of NAM that uses gene expression profiling as an endpoint for rapidly evaluating the effects of large numbers of chemicals on in vitro cell culture systems. As compared to targeted high-throughput screening (HTS) approaches that measure the effect of chemical X on target Y, HTTr is a non-targeted approach that allows researchers to more broadly characterize the integrated response of an intact biological system to chemicals that may affect a specific biological target or many biological targets under a defined set of treatment conditions (time, concentration, etc.). HTTr screening performed in concentration-response mode can provide potency estimates for the concentrations of chemicals that produce perturbations in cellular response pathways. Here, we discuss study design considerations for HTTr concentration-response screening and present a framework for the use of HTTr-based biological pathway-altering concentrations (BPACs) in a screening-level, risk-based chemical prioritization approach. The framework involves concentration-response modeling of HTTr data, mapping gene level responses to biological pathways, determination of BPACs, in vitro-to-in vivo extrapolation (IVIVE) and comparison to human exposure predictions.

A network-centric approach to drugging TNF-induced NF-κB signaling

Target-centric drug development strategies prioritize single-target potency in vitro and do not account for connectivity and multi-target effects within a signal transduction network. Here, we present a systems biology approach that combines transcriptomic and structural analyses with live-cell imaging to predict small molecule inhibitors of TNF-induced NF-κB signaling and elucidate the network response. We identify two first-in-class small molecules that inhibit the NF-κB signaling pathway by preventing the maturation of a rate-limiting multiprotein complex necessary for IKK activation. Our findings suggest that a network-centric drug discovery approach is a promising strategy to evaluate the impact of pharmacologic intervention in signaling.

Chemical perturbation of specific protein–protein interactions is notoriously difficult, yet necessary when complete inhibition of a signalling pathway is detrimental to the cell. Here, the authors use a systems approach and identify two first-in-class small molecules that specifically inhibit TNF-induced NF-κB activation.