Machine Learning Models and Pathway Genome Data Base for Trypanosoma cruzi Drug Discovery

Current advances in drug discovery for Chagas disease

Chagas disease, also known as American trypanosomiasis, is one of the 17 neglected tropical diseases (NTDs) according to World Health Organization. It is estimated that 8-10 million people are infected worldwide, mainly in Latin America. Chagas disease is caused by the parasite Trypanosoma cruzi and is characterized by two phases: acute and chronic. The current therapy for Chagas disease is limited to drugs such as nifurtimox and benznidazole, which are effective in treating only the acute phase of the disease. In addition, several side effects ranging from hypersensitivity to bone marrow depression and peripheral polyneuropathy have been associated with these drugs. Therefore, the current challenge is to find new effective and safe drugs against this NTD. The aim of this review is to describe the advances in the medicinal chemistry of new anti-chagasic compounds reported in the literature in the last five years. We report promising prototypes for drug discovery identified through target-based and phenotype-based strategies and present some important targets for the development of new synthetic compounds.

Identification of Anti-Trypanosoma cruzi Lead Compounds with Putative Immunomodulatory Activity

In seeking substitutions for the current Chagas disease treatment, which has several relevant side effects, new therapeutic candidates have been extensively investigated. In this context, a balanced interaction between mediators of the host immune response seems to be a key element for therapeutic success, as a proinflammatory microenvironment modulated by interleukin-10 (IL-10) is shown to be relevant to potentiate anti-Trypanosoma cruzi drug activity. This study aimed to identify the potential immunomodulatory activities of the anti-T. cruzi K777, pyronaridine (PYR), and furazolidone (FUR) compounds in peripheral blood mononuclear cells (PBMC) from noninfected (NI) subjects and chronic Chagas disease (CD) patients. Our results showed low cytotoxicity to PBMC populations, with 50% cytotoxic concentrations (CC50) of 71.0 μM (K777), 9.0 μM (PYR), and greater than 20 μM (FUR). In addition, K777 showed no impact on the exposure index (EI) of phytohemagglutinin-stimulated leukocytes (PHA), while PYR and FUR treatments induced increased EI of monocytes and T lymphocytes at late stages of apoptosis in NI subjects. Moreover, K777 induced a more prominent proinflammatory response (tumor necrosis factor alpha-positive [TNF-α+] CD8+/CD4+, gamma interferon-positive [IFN-γ+] CD4+/CD8+ modulated by interleukin-10-positive [IL-10+] CD4+ T/CD8+ T) than did PYR (TNF-α+ CD8+, IL-10+ CD8+) and FUR (TNF-α+ CD8+, IL-10+ CD8+). Signature analysis of intracytoplasmic cytokines corroborated the proinflammatory/modulated (K777) and proinflammatory (PYR and FUR) profiles previously found. In conclusion, the lead compound K777 may induce beneficial changes in the immunological profile of patients presenting the chronic phase of Chagas disease and may contribute to a more effective therapy against the disease.

Computational Methods to Discover Compounds for the Treatment of Chagas Disease

Infectious diseases continue to be a major public health. Among these diseases, American trypanosomiasis or Chagas disease (CD) is a major cause of morbidity and death for millions of people in Latin America. The two drugs currently available for the treatment of CD have poor efficacy and major side effects. Thus, there is a pressing need to develop safe and effective drugs against this disease. Herein we review the diversity and coverage of chemical space of compounds tested as inhibitors of Trypanosoma cruzi, a parasite causing CD. We also review major molecular targets currently pursued to kill the parasite and recent computational approaches to identify inhibitors for such targets.

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.

Comparison of Deep Learning With Multiple Machine Learning Methods and Metrics Using Diverse Drug Discovery Data Sets

Machine learning methods have been applied to many data sets in pharmaceutical research for several decades. The relative ease and availability of fingerprint type molecular descriptors paired with Bayesian methods resulted in the widespread use of this approach for a diverse array of end points relevant to drug discovery. Deep learning is the latest machine learning algorithm attracting attention for many of pharmaceutical applications from docking to virtual screening. Deep learning is based on an artificial neural network with multiple hidden layers and has found considerable traction for many artificial intelligence applications. We have previously suggested the need for a comparison of different machine learning methods with deep learning across an array of varying data sets that is applicable to pharmaceutical research. End points relevant to pharmaceutical research include absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) properties, as well as activity against pathogens and drug discovery data sets. In this study, we have used data sets for solubility, probe-likeness, hERG, KCNQ1, bubonic plague, Chagas, tuberculosis, and malaria to compare different machine learning methods using FCFP6 fingerprints. These data sets represent whole cell screens, individual proteins, physicochemical properties as well as a data set with a complex end point. Our aim was to assess whether deep learning offered any improvement in testing when assessed using an array of metrics including AUC, F1 score, Cohen's kappa, Matthews correlation coefficient and others. Based on ranked normalized scores for the metrics or data sets Deep Neural Networks (DNN) ranked higher than SVM, which in turn was ranked higher than all the other machine learning methods. Visualizing these properties for training and test sets using radar type plots indicates when models are inferior or perhaps over trained. These results also suggest the need for assessing deep learning further using multiple metrics with much larger scale comparisons, prospective testing as well as assessment of different fingerprints and DNN architectures beyond those used.

Virtual screening and repositioning of inconclusive molecules of beta-lactamase Bioassays-A data mining approach

This study focuses on the best possible way forward in utilizing inconclusive molecules of PubChem bioassays AID 1332, AID 434987 and AID 434955, which are related to beta-lactamase inhibitors of Mycobacterium tuberculosis (Mtb). The inadequacy in the experimental methods that were observed during the invitro screening resulted in an inconclusive dataset. This could be due to certain moieties present within the molecules. In order to reconsider such molecules, insilico methods can be suggested in place of invitro methods For instance, datamining and medicinal chemistry methods: have been adopted to prioritise the inconclusive dataset into active or inactive molecules. These include the Random Forest algorithm for dataminning, Lilly MedChem rules for virtually screening out the promiscuity, and Self Organizing Maps (SOM) for clustering the active molecules and enlisting them for repositioning through the use of artificial neural networks. These repositioned molecules could then be prioritized for downstream drug discovery analysis.

Bioinformatics and Drug Discovery

Bioinformatic analysis can not only accelerate drug target identification and drug candidate screening and refinement, but also facilitate characterization of side effects and predict drug resistance. High-throughput data such as genomic, epigenetic, genome architecture, cistromic, transcriptomic, proteomic, and ribosome profiling data have all made significant contribution to mechanismbased drug discovery and drug repurposing. Accumulation of protein and RNA structures, as well as development of homology modeling and protein structure simulation, coupled with large structure databases of small molecules and metabolites, paved the way for more realistic protein-ligand docking experiments and more informative virtual screening. I present the conceptual framework that drives the collection of these high-throughput data, summarize the utility and potential of mining these data in drug discovery, outline a few inherent limitations in data and software mining these data, point out news ways to refine analysis of these diverse types of data, and highlight commonly used software and databases relevant to drug discovery.

Collaborative drug discovery for More Medicines for Tuberculosis (MM4TB)

Neglected disease drug discovery is generally poorly funded compared with major diseases and hence there is an increasing focus on collaboration and precompetitive efforts such as public-private partnerships (PPPs). The More Medicines for Tuberculosis (MM4TB) project is one such collaboration funded by the EU with the goal of discovering new drugs for tuberculosis. Collaborative Drug Discovery has provided a commercial web-based platform called CDD Vault which is a hosted collaborative solution for securely sharing diverse chemistry and biology data. Using CDD Vault alongside other commercial and free cheminformatics tools has enabled support of this and other large collaborative projects, aiding drug discovery efforts and fostering collaboration. We will describe CDD's efforts in assisting with the MM4TB project.

Machine Learning Model Analysis and Data Visualization with Small Molecules Tested in a Mouse Model of Mycobacterium tuberculosis Infection (2014-2015)

The renewed urgency to develop new treatments for Mycobacterium tuberculosis (Mtb) infection has resulted in large-scale phenotypic screening and thousands of new active compounds in vitro. The next challenge is to identify candidates to pursue in a mouse in vivo efficacy model as a step to predicting clinical efficacy. We previously analyzed over 70 years of this mouse in vivo efficacy data, which we used to generate and validate machine learning models. Curation of 60 additional small molecules with in vivo data published in 2014 and 2015 was undertaken to further test these models. This represents a much larger test set than for the previous models. Several computational approaches have now been applied to analyze these molecules and compare their molecular properties beyond those attempted previously. Our previous machine learning models have been updated, and a novel aspect has been added in the form of mouse liver microsomal half-life (MLM t_1/2) and in vitro-based Mtb models incorporating cytotoxicity data that were used to predict in vivo activity for comparison. Our best Mtbin vivo models possess fivefold ROC values > 0.7, sensitivity > 80%, and concordance > 60%, while the best specificity value is >40%. Use of an MLM t_1/2 Bayesian model affords comparable results for scoring the 60 compounds tested. Combining MLM stability and in vitroMtb models in a novel consensus workflow in the best cases has a positive predicted value (hit rate) > 77%. Our results indicate that Bayesian models constructed with literature in vivoMtb data generated by different laboratories in various mouse models can have predictive value and may be used alongside MLM t_1/2 and in vitro-based Mtb models to assist in selecting antitubercular compounds with desirable in vivo efficacy. We demonstrate for the first time that consensus models of any kind can be used to predict in vivo activity for Mtb. In addition, we describe a new clustering method for data visualization and apply this to the in vivo training and test data, ultimately making the method accessible in a mobile app.

Open drug discovery for the Zika virus

The Zika virus (ZIKV) outbreak in the Americas has caused global concern that we may be on the brink of a healthcare crisis. The lack of research on ZIKV in the over 60 years that we have known about it has left us with little in the way of starting points for drug discovery. Our response can build on previous efforts with virus outbreaks and lean heavily on work done on other flaviviruses such as dengue virus. We provide some suggestions of what might be possible and propose an open drug discovery effort that mobilizes global science efforts and provides leadership, which thus far has been lacking. We also provide a listing of potential resources and molecules that could be prioritized for testing as in vitro assays for ZIKV are developed. We propose also that in order to incentivize drug discovery, a neglected disease priority review voucher should be available to those who successfully develop an FDA approved treatment. Learning from the response to the ZIKV, the approaches to drug discovery used and the success and failures will be critical for future infectious disease outbreaks.

Open Source Bayesian Models. 3. Composite Models for Prediction of Binned Responses

Bayesian models constructed from structure-derived fingerprints have been a popular and useful method for drug discovery research when applied to bioactivity measurements that can be effectively classified as active or inactive. The results can be used to rank candidate structures according to their probability of activity, and this ranking benefits from the high degree of interpretability when structure-based fingerprints are used, making the results chemically intuitive. Besides selecting an activity threshold, building a Bayesian model is fast and requires few or no parameters or user intervention. The method also does not suffer from such acute overtraining problems as quantitative structure–activity relationships or quantitative structure–property relationships (QSAR/QSPR). This makes it an approach highly suitable for automated workflows that are independent of user expertise or prior knowledge of the training data. We now describe a new method for creating a composite group of Bayesian models to extend the method to work with multiple states, rather than just binary. Incoming activities are divided into bins, each covering a mutually exclusive range of activities. For each of these bins, a Bayesian model is created to model whether or not the compound belongs in the bin. Analyzing putative molecules using the composite model involves making a prediction for each bin and examining the relative likelihood for each assignment, for example, highest value wins. The method has been evaluated on a collection of hundreds of data sets extracted from ChEMBL v20 and validated data sets for ADME/Tox and bioactivity.

Machine learning models identify molecules active against the Ebola virus in vitro

The search for small molecule inhibitors of Ebola virus (EBOV) has led to several high throughput screens over the past 3 years. These have identified a range of FDA-approved active pharmaceutical ingredients (APIs) with anti-EBOV activity in vitro and several of which are also active in a mouse infection model. There are millions of additional commercially-available molecules that could be screened for potential activities as anti-EBOV compounds. One way to prioritize compounds for testing is to generate computational models based on the high throughput screening data and then virtually screen compound libraries. In the current study, we have generated Bayesian machine learning models with viral pseudotype entry assay and the EBOV replication assay data. We have validated the models internally and externally. We have also used these models to computationally score the MicroSource library of drugs to select those likely to be potential inhibitors. Three of the highest scoring molecules that were not in the model training sets, quinacrine, pyronaridine and tilorone, were tested in vitro and had EC50 values of 350, 420 and 230 nM, respectively. Pyronaridine is a component of a combination therapy for malaria that was recently approved by the European Medicines Agency, which may make it more readily accessible for clinical testing. Like other known antimalarial drugs active against EBOV, it shares the 4-aminoquinoline scaffold. Tilorone, is an investigational antiviral agent that has shown a broad array of biological activities including cell growth inhibition in cancer cells, antifibrotic properties, α7 nicotinic receptor agonist activity, radioprotective activity and activation of hypoxia inducible factor-1. Quinacrine is an antimalarial but also has use as an anthelmintic. Our results suggest data sets with less than 1,000 molecules can produce validated machine learning models that can in turn be utilized to identify novel EBOV inhibitors in vitro.

Machine learning models identify molecules active against the Ebola virus in vitro

The search for small molecule inhibitors of Ebola virus (EBOV) has led to several high throughput screens over the past 3 years. These have identified a range of FDA-approved active pharmaceutical ingredients (APIs) with anti-EBOV activity in vitro and several of which are also active in a mouse infection model. There are millions of additional commercially-available molecules that could be screened for potential activities as anti-EBOV compounds. One way to prioritize compounds for testing is to generate computational models based on the high throughput screening data and then virtually screen compound libraries. In the current study, we have generated Bayesian machine learning models with viral pseudotype entry assay and the EBOV replication assay data. We have validated the models internally and externally. We have also used these models to computationally score the MicroSource library of drugs to select those likely to be potential inhibitors. Three of the highest scoring molecules that were not in the model training sets, quinacrine, pyronaridine and tilorone, were tested in vitro and had EC ₅₀ values of 350, 420 and 230 nM, respectively. Pyronaridine is a component of a combination therapy for malaria that was recently approved by the European Medicines Agency, which may make it more readily accessible for clinical testing. Like other known antimalarial drugs active against EBOV, it shares the 4-aminoquinoline scaffold. Tilorone, is an investigational antiviral agent that has shown a broad array of biological activities including cell growth inhibition in cancer cells, antifibrotic properties, α7 nicotinic receptor agonist activity, radioprotective activity and activation of hypoxia inducible factor-1. Quinacrine is an antimalarial but also has use as an anthelmintic. Our results suggest data sets with less than 1,000 molecules can produce validated machine learning models that can in turn be utilized to identify novel EBOV inhibitors in vitro.

Machine learning models identify molecules active against the Ebola virus in vitro

The search for small molecule inhibitors of Ebola virus (EBOV) has led to several high throughput screens over the past 3 years. These have identified a range of FDA-approved active pharmaceutical ingredients (APIs) with anti-EBOV activity in vitro and several of which are also active in a mouse infection model. There are millions of additional commercially-available molecules that could be screened for potential activities as anti-EBOV compounds. One way to prioritize compounds for testing is to generate computational models based on the high throughput screening data and then virtually screen compound libraries. In the current study, we have generated Bayesian machine learning models with viral pseudotype entry assay and the EBOV replication assay data. We have validated the models internally and externally. We have also used these models to computationally score the MicroSource library of drugs to select those likely to be potential inhibitors. Three of the highest scoring molecules that were not in the model training sets, quinacrine, pyronaridine and tilorone, were tested in vitro and had EC ₅₀ values of 350, 420 and 230 nM, respectively. Pyronaridine is a component of a combination therapy for malaria that was recently approved by the European Medicines Agency, which may make it more readily accessible for clinical testing. Like other known antimalarial drugs active against EBOV, it shares the 4-aminoquinoline scaffold. Tilorone, is an investigational antiviral agent that has shown a broad array of biological activities including cell growth inhibition in cancer cells, antifibrotic properties, α7 nicotinic receptor agonist activity, radioprotective activity and activation of hypoxia inducible factor-1. Quinacrine is an antimalarial but also has use as an anthelmintic. Our results suggest data sets with less than 1,000 molecules can produce validated machine learning models that can in turn be utilized to identify novel EBOV inhibitors in vitro.