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.

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.

Discovery of an Anion-Dependent Farnesyltransferase Inhibitor from a Phenotypic Screen.

By employing a phenotypic screen, a set of compounds, exemplified by 1, were identified which potentiate the ability of histone deacetylase inhibitor vorinostat to reverse HIV latency. Proteome enrichment followed by quantitative mass spectrometric analysis employing a modified analogue of 1 as affinity bait identified farnesyl transferase (FTase) as the primary interacting protein in cell lysates. This ligand-FTase binding interaction was confirmed via X-ray crystallography and temperature dependent fluorescence studies, despite 1 lacking structural and binding similarity to known FTase inhibitors. Although multiple lines of evidence established the binding interaction, these ligands exhibited minimal inhibitory activity in a cell-free biochemical FTase inhibition assay. Subsequent modification of the biochemical assay by increasing anion concentration demonstrated FTase inhibitory activity in this novel class. We propose 1 binds together with the anion in the active site to inhibit farnesyl transferase. Implications for phenotypic screening deconvolution and HIV reactivation are discussed.

Molecular Profiling Reveals a Common Metabolic Signature of Tissue Fibrosis.

Fibrosis, or the accumulation of extracellular matrix, is a common feature of many chronic diseases. To interrogate core molecular pathways underlying fibrosis, we cross-examine human primary cells from various tissues treated with TGF-ß, as well as kidney and liver fibrosis models. Transcriptome analyses reveal that genes involved in fatty acid oxidation are significantly perturbed. Furthermore, mitochondrial dysfunction and acylcarnitine accumulation are found in fibrotic tissues. Substantial downregulation of the PGC1α gene is evident in both in vitro and in vivo fibrosis models, suggesting a common node of metabolic signature for tissue fibrosis. In order to identify suppressors of fibrosis, we carry out a compound library phenotypic screen and identify AMPK and PPAR as highly enriched targets. We further show that pharmacological treatment of MK-8722 (AMPK activator) and MK-4074 (ACC inhibitor) reduce fibrosis in vivo. Altogether, our work demonstrate that metabolic defect is integral to TGF-ß signaling and fibrosis.

CHEMGENIE: integration of chemogenomics data for applications in chemical biology.

Increasing amounts of biological data are accumulating in the pharmaceutical industry and academic institutions. However, data does not equal actionable information, and guidelines for appropriate data capture, harmonization, integration, mining, and visualization need to be established to fully harness its potential. Here, we describe ongoing efforts at Merck & Co. to structure data in the area of chemogenomics. We are integrating complementary data from both internal and external data sources into one chemogenomics database (Chemical Genetic Interaction Enterprise; CHEMGENIE). Here, we demonstrate how this well-curated database facilitates compound set design, tool compound selection, target deconvolution in phenotypic screening, and predictive model building.

Linking High-Throughput Screens to Identify MoAs and Novel Inhibitors of Mycobacterium tuberculosis Dihydrofolate Reductase.

Though phenotypic and target-based high-throughput screening approaches have been employed to discover new antibiotics, the identification of promising therapeutic candidates remains challenging. Each approach provides different information, and understanding their results can provide hypotheses for a mechanism of action (MoA) and reveal actionable chemical matter. Here, we describe a framework for identifying efficacy targets of bioactive compounds. High throughput biophysical profiling against a broad range of targets coupled with machine learning was employed to identify chemical features with predicted efficacy targets for a given phenotypic screen. We validate the approach on data from a set of 55 000 compounds in 24 historical internal antibacterial phenotypic screens and 636 bacterial targets screened in high-throughput biophysical binding assays. Models were built to reveal the relationships between phenotype, target, and chemotype, which recapitulated mechanisms for known antibacterials. We also prospectively identified novel inhibitors of dihydrofolate reductase with nanomolar antibacterial efficacy against Mycobacterium tuberculosis. Molecular modeling provided structural insight into target-ligand interactions underlying selective killing activity toward mycobacteria over human cells.

Prospective Assessment of Virtual Screening Heuristics Derived Using a Novel Fusion Score.

High-throughput screening (HTS) is a widespread method in early drug discovery for identifying promising chemical matter that modulates a target or phenotype of interest. Because HTS campaigns involve screening millions of compounds, it is often desirable to initiate screening with a subset of the full collection. Subsequently, virtual screening methods prioritize likely active compounds in the remaining collection in an iterative process. With this approach, orthogonal virtual screening methods are often applied, necessitating the prioritization of hits from different approaches. Here, we introduce a novel method of fusing these prioritizations and benchmark it prospectively on 17 screening campaigns using virtual screening methods in three descriptor spaces. We found that the fusion approach retrieves 15% to 65% more active chemical series than any single machine-learning method and that appropriately weighting contributions of similarity and machine-learning scoring techniques can increase enrichment by 1% to 19%. We also use fusion scoring to evaluate the tradeoff between screening more chemical matter initially in lieu of replicate samples to prevent false-positives and find that the former option leads to the retrieval of more active chemical series. These results represent guidelines that can increase the rate of identification of promising active compounds in future iterative screens.

Representing high throughput expression profiles via perturbation barcodes reveals compound targets.

High throughput mRNA expression profiling can be used to characterize the response of cell culture models to perturbations such as pharmacologic modulators and genetic perturbations. As profiling campaigns expand in scope, it is important to homogenize, summarize, and analyze the resulting data in a manner that captures significant biological signals in spite of various noise sources such as batch effects and stochastic variation. We used the L1000 platform for large-scale profiling of 978 representative genes across thousands of compound treatments. Here, a method is described that uses deep learning techniques to convert the expression changes of the landmark genes into a perturbation barcode that reveals important features of the underlying data, performing better than the raw data in revealing important biological insights. The barcode captures compound structure and target information, and predicts a compound's high throughput screening promiscuity, to a higher degree than the original data measurements, indicating that the approach uncovers underlying factors of the expression data that are otherwise entangled or masked by noise. Furthermore, we demonstrate that visualizations derived from the perturbation barcode can be used to more sensitively assign functions to unknown compounds through a guilt-by-association approach, which we use to predict and experimentally validate the activity of compounds on the MAPK pathway. The demonstrated application of deep metric learning to large-scale chemical genetics projects highlights the utility of this and related approaches to the extraction of insights and testable hypotheses from big, sometimes noisy data.

Iterative Focused Screening with Biological Fingerprints Identifies Selective Asc-1 Inhibitors Distinct from Traditional High Throughput Screening.

N-methyl-d-aspartate receptors (NMDARs) mediate glutamatergic signaling that is critical to cognitive processes in the central nervous system, and NMDAR hypofunction is thought to contribute to cognitive impairment observed in both schizophrenia and Alzheimer's disease. One approach to enhance the function of NMDAR is to increase the concentration of an NMDAR coagonist, such as glycine or d-serine, in the synaptic cleft. Inhibition of alanine-serine-cysteine transporter-1 (Asc-1), the primary transporter of d-serine, is attractive because the transporter is localized to neurons in brain regions critical to cognitive function, including the hippocampus and cortical layers III and IV, and is colocalized with d-serine and NMDARs. To identify novel Asc-1 inhibitors, two different screening approaches were performed with whole-cell amino acid uptake in heterologous cells stably expressing human Asc-1: (1) a high-throughput screen (HTS) of 3 M compounds measuring 35S l-cysteine uptake into cells attached to scintillation proximity assay beads in a 1536 well format and (2) an iterative focused screen (IFS) of a 45 000 compound diversity set using a 3H d-serine uptake assay with a liquid scintillation plate reader in a 384 well format. Critically important for both screening approaches was the implementation of counter screens to remove nonspecific inhibitors of radioactive amino acid uptake. Furthermore, a 15 000 compound expansion step incorporating both on- and off-target data into chemical and biological fingerprint-based models for selection of additional hits enabled the identification of novel Asc-1-selective chemical matter from the IFS that was not identified in the full-collection HTS.

Synthesis of Complex Druglike Molecules by the Use of Highly Functionalized Bench-Stable Organozinc Reagents.

The reactivity of a representative set of 17 organozinc pivalates with 18 polyfunctional druglike electrophiles (informers) in Negishi cross-coupling reactions was evaluated by high-throughput experimentation protocols. The high-fidelity scaleup of successful reactions in parallel enabled the isolation of sufficient material for biological testing, thus demonstrating the high value of these new solid zinc reagents in a drug-discovery setting and potentially for many other applications in chemistry. Principal component analysis (PCA) clearly defined the independent roles of the zincates and the informers toward druggable-space coverage.

Chemistry informer libraries: a chemoinformatics enabled approach to evaluate and advance synthetic methods.

Major new advances in synthetic chemistry methods are typically reported using simple, non-standardized reaction substrates, and reaction failures are rarely documented. This makes the evaluation and choice of a synthetic method difficult. We report a standardized complex molecule diagnostic approach using collections of relevant drug-like molecules which we call chemistry informer libraries. With this approach, all chemistry results, successes and failures, can be documented to compare and evolve synthetic methods. To aid in the visualization of chemistry results in drug-like physicochemical space we have used an informatics methodology termed principal component analysis. We have validated this method using palladium- and copper-catalyzed reactions, including Suzuki-Miyaura, cyanation and Buchwald-Hartwig amination.

Fragment library design: using cheminformatics and expert chemists to fill gaps in existing fragment libraries.

Fragment based screening (FBS) has emerged as a mainstream lead discovery strategy in academia, biotechnology start-ups, and large pharma. As a prerequisite of FBS, a structurally diverse library of fragments is desirable in order to identify chemical matter that will interact with the range of diverse target classes that are prosecuted in contemporary screening campaigns. In addition, it is also desirable to offer synthetically amenable starting points to increase the probability of a successful fragment evolution through medicinal chemistry. Herein we describe a method to identify biologically relevant chemical substructures that are missing from an existing fragment library (chemical gaps), and organize these chemical gaps hierarchically so that medicinal chemists can efficiently navigate the prioritized chemical space and subsequently select purchasable fragments for inclusion in an enhanced fragment library.

Large scale meta-analysis of fragment-based screening campaigns: privileged fragments and complementary technologies.

A first step in fragment-based drug discovery (FBDD) often entails a fragment-based screen (FBS) to identify fragment "hits." However, the integration of conflicting results from orthogonal screens remains a challenge. Here we present a meta-analysis of 35 fragment-based campaigns at Novartis, which employed a generic 1400-fragment library against diverse target families using various biophysical and biochemical techniques. By statistically interrogating the multidimensional FBS data, we sought to investigate three questions: (1) What makes a fragment amenable for FBS? (2) How do hits from different fragment screening technologies and target classes compare with each other? (3) What is the best way to pair FBS assay technologies? In doing so, we identified substructures that were privileged for specific target classes, as well as fragments that were privileged for authentic activity against many targets. We also revealed some of the discrepancies between technologies. Finally, we uncovered a simple rule of thumb in screening strategy: when choosing two technologies for a campaign, pairing a biochemical and biophysical screen tends to yield the greatest coverage of authentic hits.

Stitched α-helical peptides via bis ring-closing metathesis.

Conformationally stabilized α-helical peptides are capable of inhibiting disease-relevant intracellular or extracellular protein-protein interactions in vivo. We have previously reported that the employment of ring-closing metathesis to introduce a single all-hydrocarbon staple along one face of an α-helical peptide greatly increases α-helical content, binding affinity to a target protein, cell penetration through active transport, and resistance to proteolytic degradation. In an effort to improve upon this technology for stabilizing a peptide in a bioactive α-helical conformation, we report the discovery of an efficient and selective bis ring-closing metathesis reaction leading to peptides bearing multiple contiguous staples connected by a central spiro ring junction. Circular dichroism spectroscopy, NMR, and computational analyses have been used to investigate the conformation of these "stitched" peptides, which are shown to exhibit remarkable thermal stabilities. Likewise, trypsin proteolysis assays confirm the achievement of a structural rigidity unmatched by peptides bearing a single staple. Furthermore, fluorescence-activated cell sorting (FACS) and confocal microscopy assays demonstrate that stitched peptides display superior cell penetrating ability compared to their stapled counterparts, suggesting that this technology may be useful not only in the context of enhancing the drug-like properties of α-helical peptides but also in producing potent agents for the intracellular delivery of proteins and oligonucleotides.

Efficient search of chemical space: navigating from fragments to structurally diverse chemotypes.

We introduce a novel strategy to sample bioactive chemical space, which follows-up on hits from fragment campaigns without the need for a crystal structure. Our results strongly suggest that screening a few hundred or thousand fragments can substantially improve the selection of small-molecule screening subsets. By combining fragment-based screening with virtual fragment linking and HTS fingerprints, we have developed an effective strategy not only to expand from low-affinity hits to potent compounds but also to hop in chemical space to substantially novel chemotypes. In benchmark calculations, our approach accessed subsets of compounds that were substantially enriched in chemically diverse hit compounds for various activity classes. Overall, half of the hits in the screening collection were found by screening only 10% of the library. Furthermore, a prospective application led to the discovery of two structurally novel histone deacetylase 4 inhibitors.

Inside the mind of a medicinal chemist: the role of human bias in compound prioritization during drug discovery.

Medicinal chemists' "intuition" is critical for success in modern drug discovery. Early in the discovery process, chemists select a subset of compounds for further research, often from many viable candidates. These decisions determine the success of a discovery campaign, and ultimately what kind of drugs are developed and marketed to the public. Surprisingly little is known about the cognitive aspects of chemists' decision-making when they prioritize compounds. We investigate 1) how and to what extent chemists simplify the problem of identifying promising compounds, 2) whether chemists agree with each other about the criteria used for such decisions, and 3) how accurately chemists report the criteria they use for these decisions. Chemists were surveyed and asked to select chemical fragments that they would be willing to develop into a lead compound from a set of ~4,000 available fragments. Based on each chemist's selections, computational classifiers were built to model each chemist's selection strategy. Results suggest that chemists greatly simplified the problem, typically using only 1-2 of many possible parameters when making their selections. Although chemists tended to use the same parameters to select compounds, differing value preferences for these parameters led to an overall lack of consensus in compound selections. Moreover, what little agreement there was among the chemists was largely in what fragments were undesirable. Furthermore, chemists were often unaware of the parameters (such as compound size) which were statistically significant in their selections, and overestimated the number of parameters they employed. A critical evaluation of the problem space faced by medicinal chemists and cognitive models of categorization were especially useful in understanding the low consensus between chemists.

Enforced presentation of an extrahelical guanine to the lesion recognition pocket of human 8-oxoguanine glycosylase, hOGG1.

A poorly understood aspect of DNA repair proteins is their ability to identify exceedingly rare sites of damage embedded in a large excess of nearly identical undamaged DNA, while catalyzing repair only at the damaged sites. Progress toward understanding this problem has been made by comparing the structures and biochemical behavior of these enzymes when they are presented with either a target lesion or a corresponding undamaged nucleobase. Trapping and analyzing such DNA-protein complexes is particularly difficult in the case of base extrusion DNA repair proteins because of the complexity of the repair reaction, which involves extrusion of the target base from DNA followed by its insertion into the active site where glycosidic bond cleavage is catalyzed. Here we report the structure of a human 8-oxoguanine (oxoG) DNA glycosylase, hOGG1, in which a normal guanine from DNA has been forcibly inserted into the enzyme active site. Although the interactions of the nucleobase with the active site are only subtly different for G versus oxoG, hOGG1 fails to catalyze excision of the normal nucleobase. This study demonstrates that even if hOGG1 mistakenly inserts a normal base into its active site, the enzyme can still reject it on the basis of catalytic incompatibility.

Mapping targetable sites on human telomerase RNA pseudoknot/template domain using 2'-OMe RNA-interacting polynucleotide (RIPtide) microarrays.

Most cellular RNAs engage in intrastrand base-pairing that gives rise to complex three-dimensional folds. This self-pairing presents an impediment toward binding of the RNA by nucleic acid-based ligands. An important step in the discovery of RNA-targeting ligands is therefore to identify those regions in a folded RNA that are accessible toward the nucleic acid-based ligand. Because the folding of RNA targets can involve interactions between nonadjacent regions and employ both Watson-Crick and non-Watson-Crick base-pairing, screening of candidate binder ensembles is typically necessary. Microarray-based screening approaches have shown great promise in this regard and have suggested that achieving complete sequence coverage would be a valuable attribute of a next generation system. Here, we report a custom microarray displaying a library of RNA-interacting polynucleotides comprising all possible 2'-OMe RNA sequences from 4- to 8-nucleotides in length. We demonstrate the utility of this array in identifying RNA-interacting polynucleotides that bind tightly and specifically to the highly conserved, functionally essential template/pseudoknot domain of human telomerase RNA and that inhibit telomerase function in vitro.

Structure of the stapled p53 peptide bound to Mdm2.

Mdm2 is a major negative regulator of the tumor suppressor p53 protein, a protein that plays a crucial role in maintaining genome integrity. Inactivation of p53 is the most prevalent defect in human cancers. Inhibitors of the Mdm2-p53 interaction that restore the functional p53 constitute potential nongenotoxic anticancer agents with a novel mode of action. We present here a 2.0 Å resolution structure of the Mdm2 protein with a bound stapled p53 peptide. Such peptides, which are conformationally and proteolytically stabilized with all-hydrocarbon staples, are an emerging class of biologics that are capable of disrupting protein-protein interactions and thus have broad therapeutic potential. The structure represents the first crystal structure of an i, i + 7 stapled peptide bound to its target and reveals that rather than acting solely as a passive conformational brace, a staple can intimately interact with the surface of a protein and augment the binding interface.