Learning from the data: mining of large high-throughput screening databases

Identifying de-NEDDylation inhibitors: Virtual high-throughput screens targeting SENP8

Protein modification can have far-reaching effects. NEDDylation, a protein modification process with the protein NEDD8, stabilizes and modifies how the targeted protein interacts with other proteins. Its role in system regulation makes it a prime therapeutic target, and virtual high-throughput screening has already identified new NEDD8 inhibitors. SENP8 matures the NEDD8 proenzyme into the active form and regulates NEDDylation by removing NEDD8 from over-NEDDylated proteins. In this work, SENP8 inhibitor candidates were identified in two rounds of virtual high-throughput screening. Of the ten candidates identified in the first round of screening, four were active in validation experiments to yield an experimental hit rate of 40%. Of the five candidates identified in the second round of screening, one was active in validation experiments to yield an experimental hit rate of 20%. Results indicate virtual high-throughput screening improved hit rates over traditional high-throughput screening. The SENP8 inhibitor candidates can be used to interrogate the NEDDylation regulation mechanism.

Virtual high-throughput screens identifying hPK-M2 inhibitors: Exploration of model extrapolation

Glycolysis with PK-M2 occurs typically in anaerobic conditions and atypically in aerobic conditions, which is known as the Warburg effect. The Warburg effect is found in many oncogenic situations and is believed to provide energy and biomass for oncogenesis to persist. The work presented targets human PK-M2 (hPK-M2) in a virtual high-throughput screen to identify new inhibitors and leads for further study. In the initial screen, one of the 12 candidates selected for experimental validation showed biological activity (hit-rate = 8.13%). In the second screen with retrained models, six of 11 candidates selected for experimental validation showed biological activity (hit-rate: 54.5%). Additionally, four different scaffolds were identified for further analysis when examining the tested candidates and compounds in the training data. Finally, extrapolation was necessary to identify a sufficient number of candidates to test in the second screen. Examination of the results suggested stepwise extrapolation to maximize efficiency.

Identifying new clotting factor XIa inhibitors in virtual high-throughput screens using PCA-GA-SVM models and signature

Blood Clotting Factor XI is an important actor in the clotting mechanism: it activates downstream zymogen involved in the clotting process. It can be targeted for activation or inhibition depending on treatment goals to enhance or inhibit clotting. In terms of antithrombosis treatment, Factor XI has emerged as a promising target to focus on. In this work, an iterative virtual high-throughput screening pipeline was proposed that can supplement current efforts to find inhibitors. The first iteration identified 11 compounds to test with 3 active for a hit-rate of 27.3%. The second iteration of the pipeline identified another 11 compounds to test with 7 active for a hit-rate of 63.6%. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1553-1565, 2018.

Pharmaceutical Machine Learning: Virtual High-Throughput Screens Identifying Promising and Economical Small Molecule Inhibitors of Complement Factor C1s

When excessively activated, C1 is insufficiently regulated, which results in tissue damage. Such tissue damage causes the complement system to become further activated to remove the resulting tissue damage, and a vicious cycle of activation/tissue damage occurs. Current Food and Drug Administration approved treatments include supplemental recombinant C1 inhibitor, but these are extremely costly and a more economical solution is desired. In our work, we have utilized an existing data set of 136 compounds that have been previously tested for activity against C1. Using these compounds and the activity data, we have created models using principal component analysis, genetic algorithm, and support vector machine approaches to characterize activity. The models were then utilized to virtually screen the 72 million compound PubChem repository. This first round of virtual high-throughput screening identified many economical and promising inhibitor candidates, a subset of which was tested to validate their biological activity. These results were used to retrain the models and rescreen PubChem in a second round vHTS. Hit rates for the first round vHTS were 57%, while hit rates for the second round vHTS were 50%. Additional structure⁻property analysis was performed on the active and inactive compounds to identify interesting scaffolds for further investigation.

Data Mining and Computational Modeling of High-Throughput Screening Datasets

We are now seeing the benefit of investments made over the last decade in high-throughput screening (HTS) that is resulting in large structure activity datasets entering public and open databases such as ChEMBL and PubChem. The growth of academic HTS screening centers and the increasing move to academia for early stage drug discovery suggests a great need for the informatics tools and methods to mine such data and learn from it. Collaborative Drug Discovery, Inc. (CDD) has developed a number of tools for storing, mining, securely and selectively sharing, as well as learning from such HTS data. We present a new web based data mining and visualization module directly within the CDD Vault platform for high-throughput drug discovery data that makes use of a novel technology stack following modern reactive design principles. We also describe CDD Models within the CDD Vault platform that enables researchers to share models, share predictions from models, and create models from distributed, heterogeneous data. Our system is built on top of the Collaborative Drug Discovery Vault Activity and Registration data repository ecosystem which allows users to manipulate and visualize thousands of molecules in real time. This can be performed in any browser on any platform. In this chapter we present examples of its use with public datasets in CDD Vault. Such approaches can complement other cheminformatics tools, whether open source or commercial, in providing approaches for data mining and modeling of HTS data.

Large-scale bioactivity analysis of the small-molecule assayed proteome

This study presents an analysis of the small molecule bioactivity profiles across large quantities of diverse protein families represented in PubChem BioAssay. We compared the bioactivity profiles of FDA approved drugs to non-FDA approved compounds, and report several distinct patterns characteristic of the approved drugs. We found that a large fraction of the previously reported higher target promiscuity among FDA approved compounds, compared to non-FDA approved bioactives, was frequently due to cross-reactivity within rather than across protein families. We identified 804 potentially novel protein target candidates for FDA approved drugs, as well as 901 potentially novel target candidates with active non-FDA approved compounds, but no FDA approved drugs with activity against these targets. We also identified 486348 potentially novel compounds active against the same targets as FDA approved drugs, as well as 153402 potentially novel compounds active against targets without active FDA approved drugs. By quantifying the agreement among replicated screens, we estimated that more than half of these novel outcomes are reproducible. Using biclustering, we identified many dense clusters of FDA approved drugs with enriched activity against a common set of protein targets. We also report the distribution of compound promiscuity using a Bayesian statistical model, and report the sensitivity and specificity of two common methods for identifying promiscuous compounds. Aggregator assays exhibited greater accuracy in identifying highly promiscuous compounds, while PAINS substructures were able to identify a much larger set of "middle range" promiscuous compounds. Additionally, we report a large number of promiscuous compounds not identified as aggregators or PAINS. In summary, the results of this study represent a rich reference for selecting novel drug and target protein candidates, as well as for eliminating candidate compounds with unselective activities.

bioassayR: Cross-Target Analysis of Small Molecule Bioactivity

Despite a large and rapidly growing body of small molecule bioactivity screens available in the public domain, systematic leverage of the data to assess target druggability and compound selectivity has been confounded by a lack of suitable cross-target analysis software. We have developed bioassayR, a computational tool that enables simultaneous analysis of thousands of bioassay experiments performed over a diverse set of compounds and biological targets. Unique features include support for large-scale cross-target analyses of both public and custom bioassays, generation of high throughput screening fingerprints (HTSFPs), and an optional preloaded database that provides access to a substantial portion of publicly available bioactivity data. bioassayR is implemented as an open-source R/Bioconductor package available from https://bioconductor.org/packages/bioassayR/.

An Automatic Quality Control Pipeline for High-Throughput Screening Hit Identification

The correction or removal of signal errors in high-throughput screening (HTS) data is critical to the identification of high-quality lead candidates. Although a number of strategies have been previously developed to correct systematic errors and to remove screening artifacts, they are not universally effective and still require fair amount of human intervention. We introduce a fully automated quality control (QC) pipeline that can correct generic interplate systematic errors and remove intraplate random artifacts. The new pipeline was first applied to ~100 large-scale historical HTS assays; in silico analysis showed auto-QC led to a noticeably stronger structure-activity relationship. The method was further tested in several independent HTS runs, where QC results were sampled for experimental validation. Significantly increased hit confirmation rates were obtained after the QC steps, confirming that the proposed method was effective in enriching true-positive hits. An implementation of the algorithm is available to the screening community.

Profiling the kinome for drug discovery

The human kinome is made up of 518 distinctive serine/threonine and tyrosine kinases, which are key components of virtually every mammalian signal transduction pathway. Consequently, kinases provide a compelling target family for the development of small molecule inhibitors, which could be used as tools to delineate the mechanism of action for biological processes and potentially be used as therapeutics to treat human diseases such as cancer. A myriad of recent technological advances have accelerated our understanding of kinome function, its relationship to tumorigenic development, and have contributed to the progression of small molecule kinase inhibitors into the clinic. Essential to the continued growth of the field are informatics tools that can assist in interpreting disparate and voluminous data sets and correctly guide decision making processes. These advances are expected to have a dramatic impact on kinase drug development and clinical diagnoses and treatment in the near future.:

Chemical information matters: an e-Research perspective on information and data sharing in the chemical sciences

Recently, a number of organisations have called for open access to scientific information and especially to the data obtained from publicly funded research, among which the Royal Society report and the European Commission press release are particularly notable. It has long been accepted that building research on the foundations laid by other scientists is both effective and efficient. Regrettably, some disciplines, chemistry being one, have been slow to recognise the value of sharing and have thus been reluctant to curate their data and information in preparation for exchanging it. The very significant increases in both the volume and the complexity of the datasets produced has encouraged the expansion of e-Research, and stimulated the development of methodologies for managing, organising, and analysing "big data". We review the evolution of cheminformatics, the amalgam of chemistry, computer science, and information technology, and assess the wider e-Science and e-Research perspective. Chemical information does matter, as do matters of communicating data and collaborating with data. For chemistry, unique identifiers, structure representations, and property descriptors are essential to the activities of sharing and exchange. Open science entails the sharing of more than mere facts: for example, the publication of negative outcomes can facilitate better understanding of which synthetic routes to choose, an aspiration of the Dial-a-Molecule Grand Challenge. The protagonists of open notebook science go even further and exchange their thoughts and plans. We consider the concepts of preservation, curation, provenance, discovery, and access in the context of the research lifecycle, and then focus on the role of metadata, particularly the ontologies on which the emerging chemical Semantic Web will depend. Among our conclusions, we present our choice of the "grand challenges" for the preservation and sharing of chemical information.

High-throughput and in silico screenings in drug discovery

BACKGROUND

In the current situation of weak drug pipelines, impending patent expiration of several blockbuster drugs, industry consolidation and changing business models that target special diseases like cancer, diabetes, Alzheimer's and obesity, the pharmaceutical industry is under intense pressure to generate a strong drug pipeline distinguished by better productivity, diversity and cost effectiveness. The goal is discovering high-quality leads in the initial stages of the development cycle, to minimize the costs associated with failures at later ones.

OBJECTIVE

Thus, there is a great amount of interest in further developing and optimizing high-throughput screening and in silico screening, the two methods responsible for generating most of the lead compounds. Although high-throughput screening is the predominant starting point for discovery programs, in silico methods have gradually made inroads by their more rational approach, to expedite the drug discovery and development process.

CONCLUSION

Modern drug discovery strategies include both methods in tandem or in an iterative way. This review primarily provides a succinct overview and comparison of experimental and in silico screening techniques, selected case studies where both methods were used in concert to investigate their performance and complementary nature and a statement on the developments in experimental and in silico approaches in the near future.

Metal impurities cause false positives in high-throughput screening campaigns

Organic impurities in compound libraries are known to often cause false-positive signals in screening campaigns for new leads, but organic impurities do not fully account for all false-positive results. We discovered inorganic impurities in our screening library that can also cause positive signals for a variety of targets and/or readout systems, including biochemical and biosensor assays. We investigated in depth the example of zinc for a specific project and in retrospect in various HTS screens at Roche and propose a straightforward counter screen using the chelator TPEN to rule out inhibition caused by zinc.

Imaging of Plasmodium liver stages to drive next-generation antimalarial drug discovery

Most malaria drug development focuses on parasite stages detected in red blood cells, even though, to achieve eradication, next-generation drugs active against both erythrocytic and exo-erythrocytic forms would be preferable. We applied a multifactorial approach to a set of >4000 commercially available compounds with previously demonstrated blood-stage activity (median inhibitory concentration < 1 micromolar) and identified chemical scaffolds with potent activity against both forms. From this screen, we identified an imidazolopiperazine scaffold series that was highly enriched among compounds active against Plasmodium liver stages. The orally bioavailable lead imidazolopiperazine confers complete causal prophylactic protection (15 milligrams/kilogram) in rodent models of malaria and shows potent in vivo blood-stage therapeutic activity. The open-source chemical tools resulting from our effort provide starting points for future drug discovery programs, as well as opportunities for researchers to investigate the biology of exo-erythrocytic forms.

The Interweaving of Cheminformatics and HTS

The aim of this chapter is to describe the stages of early drug discovery that can be assisted by techniques commonly used in the field of cheminformatics. In fact, cheminformatics tools can be applied all the way from the design of compound libraries and the analysis of HTS results, to the discovery of functional relationships between compounds and their targets.

Label-free screening assays: a strategy for finding better drug candidates

The last 10 years have seen advances in automation and high-throughput biochemistry in the drug-discovery arena. However, these advances have not led to improvements in drug-discovery success. Drug programs must find new ways to identify superior compounds. Advances in label-free assay technologies may provide advantages needed for improved drug discovery. In this article, we will discuss high-throughput MS, a technology that allows screening with native substrates and with targets inaccessible to standard assay formats. We will then discuss cell-based label-free biosensors, focusing on the increased information content available when using these platforms. We will conclude with speculation on the future and ways to obtain relevant biological information early in development to ensure the best compounds are promoted to medicinal chemistry campaigns.

Computational structure-activity relationship analysis of small-molecule agonists for human formyl peptide receptors

N-formyl peptide receptors (FPRs) are important in host defense. Because of the potential for FPRs as therapeutic targets, recent efforts have focused on identification of non-peptide agonists for two FPR subtypes, FPR1 and FPR2. Given that a number of specific small-molecule agonists have recently been identified, we hypothesized that computational structure-activity relationship (SAR) analysis of these molecules could provide new information regarding molecular features required for activity. We used a training set of 71 compounds, including 10 FPR1-specific agonists, 36 FPR2-specific agonists, and 25 non-active analogs. A sequence of (1) one-way analysis of variance selection, (2) cluster analysis, (3) linear discriminant analysis, and (4) classification tree analysis led to the derivation of SAR rules with high (95.8%) accuracy for correct classification of compounds. These SAR rules revealed key features distinguishing FPR1 versus FPR2 agonists. To verify predictive ability, we evaluated a test set of 17 additional FPR agonists, and found that the majority of these agonists (>94%) were classified correctly as agonists. This study represents the first successful application of classification tree methodology based on atom pairs to SAR analysis of FPR agonists. Importantly, these SAR rules represent a relatively simple classification approach for virtual screening of FPR1/FPR2 agonists.

New directions in library design and analysis

The high costs associated with high-throughput screening (HTS) coupled with the limited coverage and bias of current screening collections is such that diversity analysis continues to be an important criterion in lead generation. Whereas early approaches to diversity analysis were based on traditional descriptors such as two-dimensional fingerprints a recent emphasis has been on assessing scaffold coverage to ensure that a variety of different chemotypes are represented. Moreover, whether designing diverse or focused libraries, it is widely recognised that designs should aim to achieve a balance in a number of different properties and multiobjective optimisation provides an effective way of achieving such designs.

Genomics, systems biology and drug development for infectious diseases

Although a variety of drugs are available for many infectious diseases that predominantly affect the developing world reasons remain for continuing to search for new chemotherapeutics. First, the development of microbial resistance has made some of the most effective and inexpensive drug regimes unreliable and dangerous to use on severely ill patients. Second, many existing antimicrobial drugs show toxicity or are too expensive for countries where the per capita income is in the order of hundreds of dollars per year. In recognition of this, new publicly and privately financed drug discovery efforts have been established to identify and develop new therapies for diseases such as tuberculosis, malaria and AIDS. This in turn, has intensified the need for tools to facilitate drug identification for those microbes whose molecular biology is poorly understood, or which are difficult to grow in the laboratory. While much has been written about how functional genomics can be used to find novel protein targets for chemotherapeutics this review will concentrate on how genome-wide, systems biology approaches may be used following whole organism, cell-based screening to understand the mechanism of drug action or to identify biological targets of small molecules. Here we focus on protozoan parasites, however, many of the approaches can be applied to pathogenic bacteria or parasitic helminths, insects or disease-causing fungi.

In silico pharmacology for drug discovery: methods for virtual ligand screening and profiling

Pharmacology over the past 100 years has had a rich tradition of scientists with the ability to form qualitative or semi-quantitative relations between molecular structure and activity in cerebro. To test these hypotheses they have consistently used traditional pharmacology tools such as in vivo and in vitro models. Increasingly over the last decade however we have seen that computational (in silico) methods have been developed and applied to pharmacology hypothesis development and testing. These in silico methods include databases, quantitative structure-activity relationships, pharmacophores, homology models and other molecular modeling approaches, machine learning, data mining, network analysis tools and data analysis tools that use a computer. In silico methods are primarily used alongside the generation of in vitro data both to create the model and to test it. Such models have seen frequent use in the discovery and optimization of novel molecules with affinity to a target, the clarification of absorption, distribution, metabolism, excretion and toxicity properties as well as physicochemical characterization. The aim of this review is to illustrate some of the in silico methods for pharmacology that are used in drug discovery. Further applications of these methods to specific targets and their limitations will be discussed in the second accompanying part of this review.