Enhancing the Small-Scale Screenable Biological Space beyond Known Chemogenomics Libraries with Gray Chemical Matter─Compounds with Novel Mechanisms from High-Throughput Screening Profiles

Phenotypic assays have become an established approach to drug discovery. Greater disease relevance is often achieved through cellular models with increased complexity and more detailed readouts, such as gene expression or advanced imaging. However, the intricate nature and cost of these assays impose limitations on their screening capacity, often restricting screens to well-characterized small compound sets such as chemogenomics libraries. Here, we outline a cheminformatics approach to identify a small set of compounds with likely novel mechanisms of action (MoAs), expanding the MoA search space for throughput limited phenotypic assays. Our approach is based on mining existing large-scale, phenotypic high-throughput screening (HTS) data. It enables the identification of chemotypes that exhibit selectivity across multiple cell-based assays, which are characterized by persistent and broad structure activity relationships (SAR). We validate the effectiveness of our approach in broad cellular profiling assays (Cell Painting, DRUG-seq, and Promotor Signature Profiling) and chemical proteomics experiments. These experiments revealed that the compounds behave similarly to known chemogenetic libraries, but with a notable bias toward novel protein targets. To foster collaboration and advance research in this area, we have curated a public set of such compounds based on the PubChem BioAssay dataset and made it available for use by the scientific community.

Compound Activity Prediction with Dose-Dependent Transcriptomic Profiles and Deep Learning

Predicting compound activity in assays is a long-standing challenge in drug discovery. Computational models based on compound-induced gene expression signatures from a single profiling assay have shown promise toward predicting compound activity in other, seemingly unrelated, assays. Applications of such models include predicting mechanisms-of-action (MoA) for phenotypic hits, identifying off-target activities, and identifying polypharmacologies. Here, we introduce transcriptomics-to-activity transformer (TAT) models that leverage gene expression profiles observed over compound treatment at multiple concentrations to predict the compound activity in other biochemical or cellular assays. We built TAT models based on gene expression data from a RASL-seq assay to predict the activity of 2692 compounds in 262 dose-response assays. We obtained useful models for 51% of the assays, as determined through a realistic held-out set. Prospectively, we experimentally validated the activity predictions of a TAT model in a malaria inhibition assay. With a 63% hit rate, TAT successfully identified several submicromolar malaria inhibitors. Our results thus demonstrate the potential of transcriptomic responses over compound concentration and the TAT modeling framework as a cost-efficient way to identify the bioactivities of promising compounds across many assays.

Pocket Crafter: a 3D generative modeling based workflow for the rapid generation of hit molecules in drug discovery

We present a user-friendly molecular generative pipeline called Pocket Crafter, specifically designed to facilitate hit finding activity in the drug discovery process. This workflow utilized a three-dimensional (3D) generative modeling method Pocket2Mol, for the de novo design of molecules in spatial perspective for the targeted protein structures, followed by filters for chemical-physical properties and drug-likeness, structure-activity relationship analysis, and clustering to generate top virtual hit scaffolds. In our WDR5 case study, we acquired a focused set of 2029 compounds after a targeted searching within Novartis archived library based on the virtual scaffolds. Subsequently, we experimentally profiled these compounds, resulting in a novel chemical scaffold series that demonstrated activity in biochemical and biophysical assays. Pocket Crafter successfully prototyped an effective end-to-end 3D generative chemistry-based workflow for the exploration of new chemical scaffolds, which represents a promising approach in early drug discovery for hit identification.

Activation of human STING by a molecular glue-like compound

Stimulator of interferon genes (STING) is a dimeric transmembrane adapter protein that plays a key role in the human innate immune response to infection and has been therapeutically exploited for its antitumor activity. The activation of STING requires its high-order oligomerization, which could be induced by binding of the endogenous ligand, cGAMP, to the cytosolic ligand-binding domain. Here we report the discovery through functional screens of a class of compounds, named NVS-STGs, that activate human STING. Our cryo-EM structures show that NVS-STG2 induces the high-order oligomerization of human STING by binding to a pocket between the transmembrane domains of the neighboring STING dimers, effectively acting as a molecular glue. Our functional assays showed that NVS-STG2 could elicit potent STING-mediated immune responses in cells and antitumor activities in animal models.

DRUG-seq Provides Unbiased Biological Activity Readouts for Neuroscience Drug Discovery

Unbiased transcriptomic RNA-seq data has provided deep insights into biological processes. However, its impact in drug discovery has been narrow given high costs and low throughput. Proof-of-concept studies with Digital RNA with pertUrbation of Genes (DRUG)-seq demonstrated the potential to address this gap. We extended the DRUG-seq platform by subjecting it to rigorous testing and by adding an open-source analysis pipeline. The results demonstrate high reproducibility and ability to resolve the mechanism(s) of action for a diverse set of compounds. Furthermore, we demonstrate how this data can be incorporated into a drug discovery project aiming to develop therapeutics for schizophrenia using human stem cell-derived neurons. We identified both an on-target activation signature, induced by a set of chemically distinct positive allosteric modulators of the N-methyl-d-aspartate (NMDA) receptor, and independent off-target effects. Overall, the protocol and open-source analysis pipeline are a step toward industrializing RNA-seq for high-complexity transcriptomics studies performed at a saturating scale.

Cell size homeostasis is maintained by CDK4-dependent activation of p38 MAPK

While molecules that promote the growth of animal cells have been identified, it remains unclear how such signals are orchestrated to determine a characteristic target size for different cell types. It is increasingly clear that cell size is determined by size checkpoints-mechanisms that restrict the cell cycle progression of cells that are smaller than their target size. Previously, we described a p38 MAPK-dependent cell size checkpoint mechanism whereby p38 is selectively activated and prevents cell cycle progression in cells that are smaller than a given target size. In this study, we show that the specific target size required for inactivation of p38 and transition through the cell cycle is determined by CDK4 activity. Our data suggest a model whereby p38 and CDK4 cooperate analogously to the function of a thermostat: while p38 senses irregularities in size, CDK4 corresponds to the thermostat dial that sets the target size.

CYP27A1-dependent anti-melanoma activity of limonoid natural products targets mitochondrial metabolism

Three limonoid natural products with selective anti-proliferative activity against BRAF(V600E) and NRAS(Q61K)-mutation-dependent melanoma cell lines were identified. Differential transcriptome analysis revealed dependency of compound activity on expression of the mitochondrial cytochrome P450 oxidase CYP27A1, a transcriptional target of melanogenesis-associated transcription factor (MITF). We determined that CYP27A1 activity is necessary for the generation of a reactive metabolite that proceeds to inhibit cellular proliferation. A genome-wide small interfering RNA screen in combination with chemical proteomics experiments revealed gene-drug functional epistasis, suggesting that these compounds target mitochondrial biogenesis and inhibit tumor bioenergetics through a covalent mechanism. Our work suggests a strategy for melanoma-specific targeting by exploiting the expression of MITF target gene CYP27A1 and inhibiting mitochondrial oxidative phosphorylation in BRAF mutant melanomas.

Contemporary Techniques for Target Deconvolution and Mode of Action Elucidation

The elucidation of the cellular efficacy target and mechanism of action of a screening hit remain key steps in phenotypic drug discovery. A large number of experimental and in silico approaches have been introduced to address these questions and are being discussed in this chapter with a focus on recent developments. In addition to practical considerations such as throughput and technological requirements, these approaches differ conceptually in the specific compound characteristic that they are focusing on, including physical and functional interactions, cellular response patterns as well as structural features. As a result, different approaches often provide complementary information and we describe a multipronged strategy that is frequently key to successful identification of the efficacy target but also other epistatic nodes and off-targets that together shape the overall cellular effect of a bioactive compound.

Systematic Chemogenetic Library Assembly

Chemogenetic libraries, collections of well-defined chemical probes, provide tremendous value to biomedical research but require substantial effort to ensure diversity as well as quality of the contents. We have assembled a chemogenetic library by data mining and crowdsourcing institutional expertise. We are sharing our approach, lessons learned, and disclosing our current collection of 4,185 compounds with their primary annotated gene targets (https://github.com/Novartis/MoaBox). This physical collection is regularly updated and used broadly both within Novartis and in collaboration with external partners.

Growth-restricting effects of siRNA transfections: a largely deterministic combination of off-target binding and hybridization-independent competition

Perturbation of gene expression by means of synthetic small interfering RNAs (siRNAs) is a powerful way to uncover gene function. However, siRNA technology suffers from sequence-specific off-target effects and from limitations in knock-down efficiency. In this study, we assess a further problem: unintended effects of siRNA transfections on cellular fitness/proliferation. We show that the nucleotide compositions of siRNAs at specific positions have reproducible growth-restricting effects on mammalian cells in culture. This is likely distinct from hybridization-dependent off-target effects, since each nucleotide residue is seen to be acting independently and additively. The effect is robust and reproducible across different siRNA libraries and also across various cell lines, including human and mouse cells. Analyzing the growth inhibition patterns in correlation to the nucleotide sequence of the siRNAs allowed us to build a predictor that can estimate growth-restricting effects for any arbitrary siRNA sequence. Competition experiments with co-transfected siRNAs further suggest that the growth-restricting effects might be linked to an oversaturation of the cellular miRNA machinery, thus disrupting endogenous miRNA functions at large. We caution that competition between siRNA molecules could complicate the interpretation of double-knockdown or epistasis experiments, and potential interactions with endogenous miRNAs can be a factor when assaying cell growth or viability phenotypes.

Discovery of a ZIP7 inhibitor from a Notch pathway screen

The identification of activating mutations in NOTCH1 in 50% of T cell acute lymphoblastic leukemia has generated interest in elucidating how these mutations contribute to oncogenic transformation and in targeting the pathway. A phenotypic screen identified compounds that interfere with trafficking of Notch and induce apoptosis via an endoplasmic reticulum (ER) stress mechanism. Target identification approaches revealed a role for SLC39A7 (ZIP7), a zinc transport family member, in governing Notch trafficking and signaling. Generation and sequencing of a compound-resistant cell line identified a V430E mutation in ZIP7 that confers transferable resistance to the compound NVS-ZP7-4. NVS-ZP7-4 altered zinc in the ER, and an analog of the compound photoaffinity labeled ZIP7 in cells, suggesting a direct interaction between the compound and ZIP7. NVS-ZP7-4 is the first reported chemical tool to probe the impact of modulating ER zinc levels and investigate ZIP7 as a novel druggable node in the Notch pathway.

Correction to Indolyl-Pyridinyl-Propenone-Induced Methuosis through the Inhibition of PIKFYVE

[This corrects the article DOI: 10.1021/acsomega.8b00202.].

Indolyl-Pyridinyl-Propenone-Induced Methuosis through the Inhibition of PIKFYVE

Methuosis is a form of nonapoptotic cell death characterized by the accumulation of macropinosome-derived vacuoles. Herein, we identify PIKFYVE, a class III phosphoinositide (PI) kinase, as the protein target responsible for the methuosis-inducing activity of indolyl-pyridinyl-propenones (3-(5-methoxy-2-methyl-1H-indol-3-yl)-1-(4-pyridinyl)-2-propen-1-one). We further characterize the effects of chemical substitutions at the 2- and 5-indolyl positions on cytoplasmic vacuolization and PIKFYVE binding and inhibitory activity. Our study provides a better understanding of the mechanism of methuosis-inducing indolyl-pyridinyl-propenones.

Size uniformity of animal cells is actively maintained by a p38 MAPK-dependent regulation of G1-length

Animal cells within a tissue typically display a striking regularity in their size. To date, the molecular mechanisms that control this uniformity are still unknown. We have previously shown that size uniformity in animal cells is promoted, in part, by size-dependent regulation of G1 length. To identify the molecular mechanisms underlying this process, we performed a large-scale small molecule screen and found that the p38 MAPK pathway is involved in coordinating cell size and cell cycle progression. Small cells display higher p38 activity and spend more time in G1 than larger cells. Inhibition of p38 MAPK leads to loss of the compensatory G1 length extension in small cells, resulting in faster proliferation, smaller cell size and increased size heterogeneity. We propose a model wherein the p38 pathway responds to changes in cell size and regulates G1 exit accordingly, to increase cell size uniformity.

Jenkins-CI, an Open-Source Continuous Integration System, as a Scientific Data and Image-Processing Platform

High-throughput screening generates large volumes of heterogeneous data that require a diverse set of computational tools for management, processing, and analysis. Building integrated, scalable, and robust computational workflows for such applications is challenging but highly valuable. Scientific data integration and pipelining facilitate standardized data processing, collaboration, and reuse of best practices. We describe how Jenkins-CI, an “off-the-shelf,” open-source, continuous integration system, is used to build pipelines for processing images and associated data from high-content screening (HCS). Jenkins-CI provides numerous plugins for standard compute tasks, and its design allows the quick integration of external scientific applications. Using Jenkins-CI, we integrated CellProfiler, an open-source image-processing platform, with various HCS utilities and a high-performance Linux cluster. The platform is web-accessible, facilitates access and sharing of high-performance compute resources, and automates previously cumbersome data and image-processing tasks. Imaging pipelines developed using the desktop CellProfiler client can be managed and shared through a centralized Jenkins-CI repository. Pipelines and managed data are annotated to facilitate collaboration and reuse. Limitations with Jenkins-CI (primarily around the user interface) were addressed through the selection of helper plugins from the Jenkins-CI community.

Data-Driven Derivation of an "Informer Compound Set" for Improved Selection of Active Compounds in High-Throughput Screening

Despite the usefulness of high-throughput screening (HTS) in drug discovery, for some systems, low assay throughput or high screening cost can prohibit the screening of large numbers of compounds. In such cases, iterative cycles of screening involving active learning (AL) are employed, creating the need for smaller "informer sets" that can be routinely screened to build predictive models for selecting compounds from the screening collection for follow-up screens. Here, we present a data-driven derivation of an informer compound set with improved predictivity of active compounds in HTS, and we validate its benefit over randomly selected training sets on 46 PubChem assays comprising at least 300,000 compounds and covering a wide range of assay biology. The informer compound set showed improvement in BEDROC(α = 100), PRAUC, and ROCAUC values averaged over all assays of 0.024, 0.014, and 0.016, respectively, compared to randomly selected training sets, all with paired t-test p-values <10(-15). A per-assay assessment showed that the BEDROC(α = 100), which is of particular relevance for early retrieval of actives, improved for 38 out of 46 assays, increasing the success rate of smaller follow-up screens. Overall, we showed that an informer set derived from historical HTS activity data can be employed for routine small-scale exploratory screening in an assay-agnostic fashion. This approach led to a consistent improvement in hit rates in follow-up screens without compromising scaffold retrieval. The informer set is adjustable in size depending on the number of compounds one intends to screen, as performance gains are realized for sets with more than 3,000 compounds, and this set is therefore applicable to a variety of situations. Finally, our results indicate that random sampling may not adequately cover descriptor space, drawing attention to the importance of the composition of the training set for predicting actives.

Evidence-Based and Quantitative Prioritization of Tool Compounds in Phenotypic Drug Discovery

The use of potent and selective chemical tools with well-defined targets can help elucidate biological processes driving phenotypes in phenotypic screens. However, identification of selective compounds en masse to create targeted screening sets is non-trivial. A systematic approach is needed to prioritize probes, which prevents the repeated use of published but unselective compounds. Here we performed a meta-analysis of integrated large-scale, heterogeneous bioactivity data to create an evidence-based, quantitative metric to systematically rank tool compounds for targets. Our tool score (TS) was then tested on hundreds of compounds by assessing their activity profiles in a panel of 41 cell-based pathway assays. We demonstrate that high-TS tools show more reliably selective phenotypic profiles than lower-TS compounds. Additionally we highlight frequently tested compounds that are non-selective tools and distinguish target family polypharmacology from cross-family promiscuity. TS can therefore be used to prioritize compounds from heterogeneous databases for phenotypic screening.

Identifying compound efficacy targets in phenotypic drug discovery

The identification of the efficacy target(s) for hits from phenotypic compound screens remains a key step to progress compounds into drug development. In addition to efficacy targets, the characterization of epistatic proteins influencing compound activity often facilitates the elucidation of the underlying mechanism of action; and, further, early determination of off-targets that cause potentially unwanted secondary phenotypes helps in assessing potential liabilities. This short review discusses the most important technologies currently available for characterizing the direct and indirect target space of bioactive compounds following phenotypic screening. We present a comprehensive strategy employing complementary approaches to balance individual technology strengths and weaknesses.

Using information from historical high-throughput screens to predict active compounds

Modern high-throughput screening (HTS) is a well-established approach for hit finding in drug discovery that is routinely employed in the pharmaceutical industry to screen more than a million compounds within a few weeks. However, as the industry shifts to more disease-relevant but more complex phenotypic screens, the focus has moved to piloting smaller but smarter chemically/biologically diverse subsets followed by an expansion around hit compounds. One standard method for doing this is to train a machine-learning (ML) model with the chemical fingerprints of the tested subset of molecules and then select the next compounds based on the predictions of this model. An alternative approach would be to take advantage of the wealth of bioactivity information contained in older (full-deck) screens using so-called HTS fingerprints, where each element of the fingerprint corresponds to the outcome of a particular assay, as input to machine-learning algorithms. We constructed HTS fingerprints using two collections of data: 93 in-house assays and 95 publicly available assays from PubChem. For each source, an additional set of 51 and 46 assays, respectively, was collected for testing. Three different ML methods, random forest (RF), logistic regression (LR), and naïve Bayes (NB), were investigated for both the HTS fingerprint and a chemical fingerprint, Morgan2. RF was found to be best suited for learning from HTS fingerprints yielding area under the receiver operating characteristic curve (AUC) values >0.8 for 78% of the internal assays and enrichment factors at 5% (EF(5%)) >10 for 55% of the assays. The RF(HTS-fp) generally outperformed the LR trained with Morgan2, which was the best ML method for the chemical fingerprint, for the majority of assays. In addition, HTS fingerprints were found to retrieve more diverse chemotypes. Combining the two models through heterogeneous classifier fusion led to a similar or better performance than the best individual model for all assays. Further validation using a pair of in-house assays and data from a confirmatory screen--including a prospective set of around 2000 compounds selected based on our approach--confirmed the good performance. Thus, the combination of machine-learning with HTS fingerprints and chemical fingerprints utilizes information from both domains and presents a very promising approach for hit expansion, leading to more hits. The source code used with the public data is provided.

Causal Network Models for Predicting Compound Targets and Driving Pathways in Cancer

Gene-expression data are often used to infer pathways regulating transcriptional responses. For example, differentially expressed genes (DEGs) induced by compound treatment can help characterize hits from phenotypic screens, either by correlation with known drug signatures or by pathway enrichment. Pathway enrichment is, however, typically computed with DEGs rather than "upstream" nodes that are potentially causal of "downstream" changes. Here, we present graph-based models to predict causal targets from compound-microarray data. We test several approaches to traversing network topology, and show that a consensus minimum-rank score (SigNet) beat individual methods and could highly rank compound targets among all network nodes. In addition, larger, less canonical networks outperformed linear canonical interactions. Importantly, pathway enrichment using causal nodes rather than DEGs recovers relevant pathways more often. To further validate our approach, we used integrated data sets from the Cancer Genome Atlas to identify driving pathways in triple-negative breast cancer. Critical pathways were uncovered, including the epidermal growth factor receptor 2-phosphatidylinositide 3-kinase-AKT-MAPK growth pathway andATR-p53-BRCA DNA damage pathway, in addition to unexpected pathways, such as TGF-WNT cytoskeleton remodeling, IL12-induced interferon gamma production, and TNFR-IAP (inhibitor of apoptosis) apoptosis; the latter was validated by pooled small hairpin RNA profiling in cancer cells. Overall, our approach can bridge transcriptional profiles to compound targets and driving pathways in cancer.

Translation of off-target effects: prediction of ADRs by integrated experimental and computational approach

Adverse drug reactions (ADRs) are associated with most drugs, often discovered late in drug development and sometimes only during extended course of clinical use.

The multidimensional perturbation value: a single metric to measure similarity and activity of treatments in high-throughput multidimensional screens

Screens using high-throughput, information-rich technologies such as microarrays, high-content screening (HCS), and next-generation sequencing (NGS) have become increasingly widespread. Compared with single-readout assays, these methods produce a more comprehensive picture of the effects of screened treatments. However, interpreting such multidimensional readouts is challenging. Univariate statistics such as t-tests and Z-factors cannot easily be applied to multidimensional profiles, leaving no obvious way to answer common screening questions such as "Is treatment X active in this assay?" and "Is treatment X different from (or equivalent to) treatment Y?" We have developed a simple, straightforward metric, the multidimensional perturbation value (mp-value), which can be used to answer these questions. Here, we demonstrate application of the mp-value to three data sets: a multiplexed gene expression screen of compounds and genomic reagents, a microarray-based gene expression screen of compounds, and an HCS compound screen. In all data sets, active treatments were successfully identified using the mp-value, and simulations and follow-up analyses supported the mp-value's statistical and biological validity. We believe the mp-value represents a promising way to simplify the analysis of multidimensional data while taking full advantage of its richness.

Rethinking molecular similarity: comparing compounds on the basis of biological activity

Since the advent of high-throughput screening (HTS), there has been an urgent need for methods that facilitate the interrogation of large-scale chemical biology data to build a mode of action (MoA) hypothesis. This can be done either prior to the HTS by subset design of compounds with known MoA or post HTS by data annotation and mining. To enable this process, we developed a tool that compares compounds solely on the basis of their bioactivity: the chemical biological descriptor "high-throughput screening fingerprint" (HTS-FP). In the current embodiment, data are aggregated from 195 biochemical and cell-based assays developed at Novartis and can be used to identify bioactivity relationships among the in-house collection comprising ~1.5 million compounds. We demonstrate the value of the HTS-FP for virtual screening and in particular scaffold hopping. HTS-FP outperforms state of the art methods in several aspects, retrieving bioactive compounds with remarkable chemical dissimilarity to a probe structure. We also apply HTS-FP for the design of screening subsets in HTS. Using retrospective data, we show that a biodiverse selection of plates performs significantly better than a chemically diverse selection of plates, both in terms of number of hits and diversity of chemotypes retrieved. This is also true in the case of hit expansion predictions using HTS-FP similarity. Sets of compounds clustered with HTS-FP are biologically meaningful, in the sense that these clusters enrich for genes and gene ontology (GO) terms, showing that compounds that are bioactively similar also tend to target proteins that operate together in the cell. HTS-FP are valuable not only because of their predictive power but mainly because they relate compounds solely on the basis of bioactivity, harnessing the accumulated knowledge of a high-throughput screening facility toward the understanding of how compounds interact with the proteome.

Large-Scale QSAR in Target Prediction and Phenotypic HTS Assessment

The advent of in silico compound target prediction offers a potential paradigm shift in how large compound collections are understood and used strategically in high-throughput screens (HTS). Specifically, phenotypic HTS hits may be annotated both with known targets and predicted targets using large-scale QSAR models, enabling a more sophisticated hit assessment. Efforts in massive bioactivity data integration and standardization is empowering such compound-target annotations. These approaches differ fundamentally from the traditional role of QSAR in lead optimization and binding affinity predictions to global, probabilistic target predictions for thousands of human proteins.

Determination of minimal transcriptional signatures of compounds for target prediction

The identification of molecular target and mechanism of action of compounds is a key hurdle in drug discovery. Multiplexed techniques for bead-based expression profiling allow the measurement of transcriptional signatures of compound-treated cells in high-throughput mode. Such profiles can be used to gain insight into compounds' mode of action and the protein targets they are modulating. Through the proxy of target prediction from such gene signatures we explored important aspects of the use of transcriptional profiles to capture biological variability of perturbed cellular assays. We found that signatures derived from expression data and signatures derived from biological interaction networks performed equally well, and we showed that gene signatures can be optimised using a genetic algorithm. Gene signatures of approximately 128 genes seemed to be most generic, capturing a maximum of the perturbation inflicted on cells through compound treatment. Moreover, we found evidence for oxidative phosphorylation to be one of the most general ways to capture compound perturbation.

From in silico target prediction to multi-target drug design: current databases, methods and applications

Given the tremendous growth of bioactivity databases, the use of computational tools to predict protein targets of small molecules has been gaining importance in recent years. Applications span a wide range, from the 'designed polypharmacology' of compounds to mode-of-action analysis. In this review, we firstly survey databases that can be used for ligand-based target prediction and which have grown tremendously in size in the past. We furthermore outline methods for target prediction that exist, both based on the knowledge of bioactivities from the ligand side and methods that can be applied in situations when a protein structure is known. Applications of successful in silico target identification attempts are discussed in detail, which were based partly or in whole on computational target predictions in the first instance. This includes the authors' own experience using target prediction tools, in this case considering phenotypic antibacterial screens and the analysis of high-throughput screening data. Finally, we will conclude with the prospective application of databases to not only predict, retrospectively, the protein targets of a small molecule, but also how to design ligands with desired polypharmacology in a prospective manner.

A lead discovery strategy driven by a comprehensive analysis of proteases in the peptide substrate space

We present here a comprehensive analysis of proteases in the peptide substrate space and demonstrate its applicability for lead discovery. Aligned octapeptide substrates of 498 proteases taken from the MEROPS peptidase database were used for the in silico analysis. A multiple-category naïve Bayes model, trained on the two-dimensional chemical features of the substrates, was able to classify the substrates of 365 (73%) proteases and elucidate statistically significant chemical features for each of their specific substrate positions. The positional awareness of the method allows us to identify the most similar substrate positions between proteases. Our analysis reveals that proteases from different families, based on the traditional classification (aspartic, cysteine, serine, and metallo), could have substrates that differ at the cleavage site (P1-P1') but are similar away from it. Caspase-3 (cysteine protease) and granzyme B (serine protease) are previously known examples of cross-family neighbors identified by this method. To assess whether peptide substrate similarity between unrelated proteases could reliably translate into the discovery of low molecular weight synthetic inhibitors, a lead discovery strategy was tested on two other cross-family neighbors--namely cathepsin L2 and matrix metallo proteinase 9, and calpain 1 and pepsin A. For both these pairs, a naïve Bayes classifier model trained on inhibitors of one protease could successfully enrich those of its neighbor from a different family and vice versa, indicating that this approach could be prospectively applied to lead discovery for a novel protease target with no known synthetic inhibitors.

Plate-Based Diversity Selection Based on Empirical HTS Data to Enhance the Number of Hits and Their Chemical Diversity

Typically, screening collections of pharmaceutical companies contain more than a million compounds today. However, for certain high-throughput screening (HTS) campaigns, constraints posed by the assay throughput and/or the reagent costs make it impractical to screen the entire deck. Therefore, it is desirable to effectively screen subsets of the collection based on a hypothesis or a diversity selection. How to select compound subsets is a subject of ongoing debate. The authors present an approach based on extended connectivity fingerprints to carry out diversity selection on a per plate basis (instead of a per compound basis). HTS data from 35 Novartis screens spanning 5 target classes were investigated to assess the performance of this approach. The analysis shows that selecting a fingerprint-diverse subset of 250K compounds, representing 20% of the screening deck, would have achieved significantly higher hit rates for 86% of the screens. This measure also outperforms the Murcko scaffold-based plate selection described previously, where only 49% of the screens showed similar improvements. Strikingly, the 2-fold improvement in average hit rates observed for 3 of 5 target classes in the data set indicates a target bias of the plate (and thus compound) selection method. Even though the diverse subset selection lacks any target hypothesis, its application shows significantly better results for some targets—namely, G-protein-coupled receptors, proteases, and protein-protein interactions—but not for kinase and pathway screens. The synthetic origin of the compounds in the diverse subset appears to influence the screening hit rates. Natural products were the most diverse compound class, with significantly higher hit rates compared to the compounds from the traditional synthetic and combinatorial libraries. These results offer empirical guidelines for plate-based diversity selection to enhance hit rates, based on target class and the library type being screened. ( Journal of Biomolecular Screening 2009:690-699)

How similar are similarity searching methods? A principal component analysis of molecular descriptor space

Different molecular descriptors capture different aspects of molecular structures, but this effect has not yet been quantified systematically on a large scale. In this work, we calculate the similarity of 37 descriptors by repeatedly selecting query compounds and ranking the rest of the database. Euclidean distances between the rank-ordering of different descriptors are calculated to determine descriptor (as opposed to compound) similarity, followed by PCA for visualization. Four broad descriptor classes are identified, which are circular fingerprints; circular fingerprints considering counts; path-based and keyed fingerprints; and pharmacophoric descriptors. Descriptor behavior is much more defined by those four classes than the particular parametrization. Using counts instead of the presence/absence of fingerprints significantly changes descriptor behavior, which is crucial for performance of topological autocorrelation vectors, but not circular fingerprints. Four-point pharmacophores (piDAPH4) surprisingly lead to much higher retrieval rates than three-point pharmacophores (28.21% vs 19.15%) but still similar rank-ordering of compounds (retrieval of similar actives). Looking into individual rankings, circular fingerprints seem more appropriate than path-based fingerprints if complex ring systems or branching patterns are present; count-based fingerprints could be more suitable in databases with a large number of repeated subunits (amide bonds, sugar rings, terpenes). Information-based selection of diverse fingerprints for consensus scoring (ECFP4/TGD fingerprints) led only to marginal improvement over single fingerprint results. While it seems to be nontrivial to exploit orthogonal descriptor behavior to improve retrieval rates in consensus virtual screening, those descriptors still each retrieve different actives which corroborates the strategy of employing diverse descriptors individually in prospective virtual screening settings.

Chemogenomic analysis of safety profiling data

Understanding the safety of newly developed compounds is a key task in each early drug discovery project. In early stages, pharmaceutical companies address this task by using so-called preclinical safety profiling, in which compounds are screened in inexpensive large-scale assays to understand possible liabilities. This process generates a large amount of binding data on various compounds against a panel of targets - usually thousands or tens of thousands of compounds profiled against approximately 100 different targets. This data matrix is highly valuable and elicits further analysis. After briefly introducing the nature of safety profiling data, we describe several computational methods used internally at Novartis to analyze it. We showcase protocols that can be used to understand compound promiscuity on a chemical structure level and protocols to evaluate the promiscuity of targets used in safety profiling. We also describe a method to quickly determine the chemical similarity of compounds active against different targets. Next, it is shown what protocols can be used to evaluate global chemical similarity of targets. The above approaches can be used either to optimize the composition of a panel of targets or to better understand certain toxicities. Finally, we will explain a simple method to elucidate hidden patterns in safety profiling data.

Knowledge-based virtual screening: application to the MDM4/p53 protein-protein interaction

Chemogenomics knowledge-based drug discovery approaches aim to extract the knowledge gained from one target and to apply it for the discovery of ligands and hopefully drugs of a new target which is related to the parent target by homology or conserved molecular recognition. Herein, we demonstrate the potential of knowledge-based virtual screening by applying it to the MDM4-p53 protein-protein interaction where the MDM2-p53 protein-protein interaction constitutes the parent reference system; both systems are potentially relevant to cancer therapy. We show that a combination of virtual screening methods, including homology based similarity searching, QSAR (Quantitative Structure-Activity Relationship) methods, HTD (High Throughput Docking), and UNITY pharmacophore searching provide a successful approach to the discovery of inhibitors. The virtual screening hit list is of the magnitude of 50,000 compounds picked from the corporate compound library of approximately 1.2 million compounds. Emphasis is placed on the facts that such campaigns are only feasible because of the now existing HTCP (High throughput Cherry-Picking) automation systems in combination with robust MTS (Medium Throughput Screening) fluorescence-based assays. Given that the MDM2-p53 system constitutes the reference system, it is not surprising that significantly more and stronger hits are found for this interaction compared to the MDM4-p53 system. Novel, selective and dual hits are discovered for both systems. A hit rate analysis will be provided compared to the full HTS (High-throughput Screening).

Analysis of pharmacology data and the prediction of adverse drug reactions and off-target effects from chemical structure

Preclinical Safety Pharmacology (PSP) attempts to anticipate adverse drug reactions (ADRs) during early phases of drug discovery by testing compounds in simple, in vitro binding assays (that is, preclinical profiling). The selection of PSP targets is based largely on circumstantial evidence of their contribution to known clinical ADRs, inferred from findings in clinical trials, animal experiments, and molecular studies going back more than forty years. In this work we explore PSP chemical space and its relevance for the prediction of adverse drug reactions. Firstly, in silico (computational) Bayesian models for 70 PSP-related targets were built, which are able to detect 93% of the ligands binding at IC(50) < or = 10 microM at an overall correct classification rate of about 94%. Secondly, employing the World Drug Index (WDI), a model for adverse drug reactions was built directly based on normalized side-effect annotations in the WDI, which does not require any underlying functional knowledge. This is, to our knowledge, the first attempt to predict adverse drug reactions across hundreds of categories from chemical structure alone. On average 90% of the adverse drug reactions observed with known, clinically used compounds were detected, an overall correct classification rate of 92%. Drugs withdrawn from the market (Rapacuronium, Suprofen) were tested in the model and their predicted ADRs align well with known ADRs. The analysis was repeated for acetylsalicylic acid and Benperidol which are still on the market. Importantly, features of the models are interpretable and back-projectable to chemical structure, raising the possibility of rationally engineering out adverse effects. By combining PSP and ADR models new hypotheses linking targets and adverse effects can be proposed and examples for the opioid mu and the muscarinic M2 receptors, as well as for cyclooxygenase-1 are presented. It is hoped that the generation of predictive models for adverse drug reactions is able to help support early SAR to accelerate drug discovery and decrease late stage attrition in drug discovery projects. In addition, models such as the ones presented here can be used for compound profiling in all development stages.

Which aspects of HTS are empirically correlated with downstream success?

High-throughput screening (HTS) is a well-established hit-finding approach used in the pharmaceutical industry. In this article, recent experience at Novartis with respect to factors influencing the success of HTS campaigns is discussed. An inherent measure of HTS quality could be defined by the assay Z and Z' factors, the number of hits and their biological potencies; however, such measures of quality do not always correlate with the advancement of hits to the later stages of drug discovery. Also, for many target classes, such as kinases, it is easy to identify hits, but, as a result of selectivity, intellectual property and other issues, the projects do not result in lead declarations. In this article, HTS success is defined as the fraction of HTS campaigns that advance into the later stages of drug discovery, and the major influencing factors are examined. Interestingly, screening compounds in individual wells or in mixtures did not have a major impact on the HTS success and, equally interesting, there was no difference in the progression rates of biochemical and cell-based assays. Particular target types, assay technologies, structure-activity relationships and powder availability had a much greater impact on success as defined above. In addition, significant mutual dependencies can be observed - while one assay format works well with one target type, this situation might be completely reversed for a combination of the same readout technology with a different target type. The results and opinions presented here should be regarded as groundwork, and a plethora of factors that influence the fate of a project, such as biophysical measurements, chemical attractiveness of the hits, strategic reasons and safety pharmacology, are not covered here. Nonetheless, it is hoped that this information will be used industry-wide to improve success rates in terms of hits progressing into exploratory chemistry and beyond. The support that can be obtained from new in silico approaches to phase transitions are also described, along with the gaps they are designed to fill.

Chemogenomic data analysis: prediction of small-molecule targets and the advent of biological fingerprint

Chemogenomics comprises a systematic relationship between targets and ligands that are used as target modulators in living systems such as cells or organisms. In recent years, data on small molecule-bioactivity relationships have become increasingly available, and consequently so have the number of approaches used to translate bioactivity data into knowledge. This review will focus on two aspects of chemogenomics. Firstly, in cases such as cell-based screens, the question of which target(s) a compound is modulating in order to cause the observed phenotype is crucial. In silico target prediction tools can suggest likely biological targets of small molecules via data mining in target-annotated chemical databases. This review presents some of the current tools available for this task and shows some sample applications relevant to a pharmaceutical industry setting. These applications are the prediction of false-positives in cell-based reporter gene assays, the prediction of targets by linking bioassay data with protein domain annotations, and the direct prediction of adverse reactions. Secondly, in recent years a shift from structure-derived chemical descriptors to biological descriptors has occurred. Here, the effect of a compound on a number of biological endpoints is used to make predictions about other properties, such as putative targets, associated adverse reactions, and pathways modulated by the compound. This review further summarizes these "performance" descriptors and their applications, focusing on gene expression profiles and high-content screening data. The advent of such biological fingerprints suggests that the field of drug discovery is currently at a crossroads, where single target bioassay results are supplanted by multidimensional biological fingerprints that reflect a new awareness of biological networks and polypharmacology.

Blocking binding of Bacillus thuringiensis Cry1Aa to Bombyx mori cadherin receptor results in only a minor reduction of toxicity

Background

Bacillus thuringiensis Cry1Aa insecticidal protein is the most active known B. thuringiensis toxin against the forest insect pest Lymantria dispar (gypsy moth), unfortunately it is also highly toxic against the non-target insect Bombyx mori (silk worm).

Results

Surface exposed hydrophobic residues over domains II and III were targeted for site-directed mutagenesis. Substitution of a phenylalanine residue (F328) by alanine reduced binding to the Bombyx mori cadherin by 23-fold, reduced biological activity against B. mori by 4-fold, while retaining activity against Lymantria dispar.

Conclusion

The results identify a novel receptor-binding epitope and demonstrate that virtual elimination of binding to cadherin BR-175 does not completely remove toxicity in the case of B. mori.

Understanding false positives in reporter gene assays: in silico chemogenomics approaches to prioritize cell-based HTS data

High throughput screening (HTS) data is often noisy, containing both false positives and negatives. Thus, careful triaging and prioritization of the primary hit list can save time and money by identifying potential false positives before incurring the expense of followup. Of particular concern are cell-based reporter gene assays (RGAs) where the number of hits may be prohibitively high to be scrutinized manually for weeding out erroneous data. Based on statistical models built from chemical structures of 650 000 compounds tested in RGAs, we created "frequent hitter" models that make it possible to prioritize potential false positives. Furthermore, we followed up the frequent hitter evaluation with chemical structure based in silico target predictions to hypothesize a mechanism for the observed "off target" response. It was observed that the predicted cellular targets for the frequent hitters were known to be associated with undesirable effects such as cytotoxicity. More specifically, the most frequently predicted targets relate to apoptosis and cell differentiation, including kinases, topoisomerases, and protein phosphatases. The mechanism-based frequent hitter hypothesis was tested using 160 additional druglike compounds predicted by the model to be nonspecific actives in RGAs. This validation was successful (showing a 50% hit rate compared to a normal hit rate as low as 2%), and it demonstrates the power of computational models toward understanding complex relations between chemical structure and biological function.

Modeling Promiscuity Based on in vitro Safety Pharmacology Profiling Data

This study describes a method for mining and modeling binding data obtained from a large panel of targets (in vitro safety pharmacology) to distinguish differences between promiscuous and selective compounds. Two naïve Bayes models for promiscuity and selectivity were generated and validated on a test set as well as publicly available drug databases. The model shows a higher score (lower promiscuity) for marketed drugs than for compounds in early development or compounds that failed during clinical development. Such models can be used in triaging high-throughput screening data or for lead optimization.