Using High-Throughput Screening Data To Discriminate Compounds with Single-Target Effects from Those with Side Effects

Preliminary insight on diarylpentanoids as potential antimalarials: In silico, in vitro pLDH and in vivo zebrafish toxicity assessment

Malaria remains a major public health problem worldwide, including in Southeast Asia. Chemotherapeutic agents such as chloroquine (CQ) are effective, but problems with drug resistance and toxicity have necessitated a continuous search for new effective antimalarial agents. Here we report on a virtual screening of ∼300 diarylpentanoids and derivatives, in search of potential Plasmodium falciparum lactate dehydrogenase (PfLDH) inhibitors with acceptable drug-like properties. Several molecules with binding affinities comparable to CQ were chosen for in vitro validation of antimalarial efficacy. Among them, MS33A, MS33C and MS34C are the most promising against CQ-sensitive (3D7) with EC50 values of 1.6, 2.5 and 3.1 μM, respectively. Meanwhile, MS87 (EC50 of 1.85 μM) shown the most active against the CQ-resistant Gombak A strain, and MS33A and MS33C the most effective P. knowlesi inhibitors (EC50 of 3.6 and 5.1 μM, respectively). The in vitro cytotoxicity of selected diarylpentanoids (MS33A, MS33C, MS34C and MS87) was tested on Vero mammalian cells to evaluate parasite selectivity (SI), showing moderate to low cytotoxicity (CC50 > 82 μM). In addition, MS87 exhibited a high SI and the lowest resistance index (RI), suggesting that MS87 may exert effective parasite inhibition with low resistance potential in the CQ-resistant P. falciparum strain. Furthermore, the in vivo toxicity of the molecules on early embryonic development, the cardiovascular system, heart rate, motor activity and apoptosis were assessed in a zebrafish animal model. The overall results indicate the preliminary potential of diarylpentanoids, which need further investigation for their development as new antimalarial agents.

Artificial intelligence and machine-learning approaches in structure and ligand-based discovery of drugs affecting central nervous system

CNS disorders are indications with a very high unmet medical needs, relatively smaller number of available drugs, and a subpar satisfaction level among patients and caregiver. Discovery of CNS drugs is extremely expensive affair with its own unique challenges leading to extremely high attrition rates and low efficiency. With explosion of data in information age, there is hardly any aspect of life that has not been touched by data driven technologies such as artificial intelligence (AI) and machine learning (ML). Drug discovery is no exception, emergence of big data via genomic, proteomic, biological, and chemical technologies has driven pharmaceutical giants to collaborate with AI oriented companies to revolutionise drug discovery, with the goal of increasing the efficiency of the process. In recent years many examples of innovative applications of AI and ML techniques in CNS drug discovery has been reported. Research on therapeutics for diseases such as schizophrenia, Alzheimer's and Parkinsonism has been provided with a new direction and thrust from these developments. AI and ML has been applied to both ligand-based and structure-based drug discovery and design of CNS therapeutics. In this review, we have summarised the general aspects of AI and ML from the perspective of drug discovery followed by a comprehensive coverage of the recent developments in the applications of AI/ML techniques in CNS drug discovery.

Big Data and Artificial Intelligence Modeling for Drug Discovery

Due to the massive data sets available for drug candidates, modern drug discovery has advanced to the big data era. Central to this shift is the development of artificial intelligence approaches to implementing innovative modeling based on the dynamic, heterogeneous, and large nature of drug data sets. As a result, recently developed artificial intelligence approaches such as deep learning and relevant modeling studies provide new solutions to efficacy and safety evaluations of drug candidates based on big data modeling and analysis. The resulting models provided deep insights into the continuum from chemical structure to in vitro, in vivo, and clinical outcomes. The relevant novel data mining, curation, and management techniques provided critical support to recent modeling studies. In summary, the new advancement of artificial intelligence in the big data era has paved the road to future rational drug development and optimization, which will have a significant impact on drug discovery procedures and, eventually, public health.

Comparing Multiple Machine Learning Algorithms and Metrics for Estrogen Receptor Binding Prediction

Many chemicals that disrupt endocrine function have been linked to a variety of adverse biological outcomes. However, screening for endocrine disruption using in vitro or in vivo approaches is costly and time-consuming. Computational methods, e.g., quantitative structure-activity relationship models, have become more reliable due to bigger training sets, increased computing power, and advanced machine learning algorithms, such as multilayered artificial neural networks. Machine learning models can be used to predict compounds for endocrine disrupting capabilities, such as binding to the estrogen receptor (ER), and allow for prioritization and further testing. In this work, an exhaustive comparison of multiple machine learning algorithms, chemical spaces, and evaluation metrics for ER binding was performed on public data sets curated using in-house cheminformatics software (Assay Central). Chemical features utilized in modeling consisted of binary fingerprints (ECFP6, FCFP6, ToxPrint, or MACCS keys) and continuous molecular descriptors from RDKit. Each feature set was subjected to classic machine learning algorithms (Bernoulli Naive Bayes, AdaBoost Decision Tree, Random Forest, Support Vector Machine) and Deep Neural Networks (DNN). Models were evaluated using a variety of metrics: recall, precision, F1-score, accuracy, area under the receiver operating characteristic curve, Cohen's Kappa, and Matthews correlation coefficient. For predicting compounds within the training set, DNN has an accuracy higher than that of other methods; however, in 5-fold cross validation and external test set predictions, DNN and most classic machine learning models perform similarly regardless of the data set or molecular descriptors used. We have also used the rank normalized scores as a performance-criteria for each machine learning method, and Random Forest performed best on the validation set when ranked by metric or by data sets. These results suggest classic machine learning algorithms may be sufficient to develop high quality predictive models of ER activity.

Global Mapping of Traditional Chinese Medicine into Bioactivity Space and Pathways Annotation Improves Mechanistic Understanding and Discovers Relationships between Therapeutic Action (Sub)classes

Traditional Chinese medicine (TCM) still needs more scientific rationale to be proven for it to be accepted further in the West. We are now in the position to propose computational hypotheses for the mode-of-actions (MOAs) of 45 TCM therapeutic action (sub)classes from in silico target prediction algorithms, whose target was later annotated with Kyoto Encyclopedia of Genes and Genomes pathway, and to discover the relationship between them by generating a hierarchical clustering. The results of 10,749 TCM compounds showed 183 enriched targets and 99 enriched pathways from Estimation Score ≤ 0 and ≥ 5% of compounds/targets in a (sub)class. The MOA of a (sub)class was established from supporting literature. Overall, the most frequent top three enriched targets/pathways were immune-related targets such as tyrosine-protein phosphatase nonreceptor type 2 (PTPN2) and digestive system such as mineral absorption. We found two major protein families, G-protein coupled receptor (GPCR), and protein kinase family contributed to the diversity of the bioactivity space, while digestive system was consistently annotated pathway motif, which agreed with the important treatment principle of TCM, “the foundation of acquired constitution” that includes spleen and stomach. In short, the TCM (sub)classes, in many cases share similar targets/pathways despite having different indications.

Big data in chemical toxicity research: the use of high-throughput screening assays to identify potential toxicants

High-throughput screening (HTS) assays that measure the in vitro toxicity of environmental compounds have been widely applied as an alternative to in vivo animal tests of chemical toxicity. Current HTS studies provide the community with rich toxicology information that has the potential to be integrated into toxicity research. The available in vitro toxicity data is updated daily in structured formats (e.g., deposited into PubChem and other data-sharing web portals) or in an unstructured way (papers, laboratory reports, toxicity Web site updates, etc.). The information derived from the current toxicity data is so large and complex that it becomes difficult to process using available database management tools or traditional data processing applications. For this reason, it is necessary to develop a big data approach when conducting modern chemical toxicity research. In vitro data for a compound, obtained from meaningful bioassays, can be viewed as a response profile that gives detailed information about the compound's ability to affect relevant biological proteins/receptors. This information is critical for the evaluation of complex bioactivities (e.g., animal toxicities) and grows rapidly as big data in toxicology communities. This review focuses mainly on the existing structured in vitro data (e.g., PubChem data sets) as response profiles for compounds of environmental interest (e.g., potential human/animal toxicants). Potential modeling and mining tools to use the current big data pool in chemical toxicity research are also described.

A high-content subtractive screen for selecting small molecules affecting internalization of GPCRs

G-protein-coupled receptors (GPCRs) are pivotal in cellular responses to the environment and are common drug targets. Identification of selective small molecules acting on single GPCRs is complicated by the shared machinery coupling signal transduction to physiology. Here, we demonstrate a high-content screen using a panel of GPCR assays to identify receptor selective molecules acting within the kinase/phosphatase inhibitor family. A collection of 88 kinase and phosphatase inhibitors was screened against seven agonist-induced GPCR internalization cell models as well as transferrin uptake in human embryonic kidney cells. Molecules acting on a single receptor were identified through excluding pan-specific compounds affecting housekeeping endocytosis or disrupting internalization of multiple receptors. We identified compounds acting on a sole GPCR from activities in a broad range of chemical structures that could not be easily sorted by conventional means. Selective analysis can therefore rapidly select compounds selectively affecting GPCR activity with specificity to one receptor class through high-content screening.

Quantifying structure and performance diversity for sets of small molecules comprising small-molecule screening collections

Using a diverse collection of small molecules we recently found that compound sets from different sources (commercial; academic; natural) have different protein-binding behaviors, and these behaviors correlate with trends in stereochemical complexity for these compound sets. These results lend insight into structural features that synthetic chemists might target when synthesizing screening collections for biological discovery. We report extensive characterization of structural properties and diversity of biological performance for these compounds and expand comparative analyses to include physicochemical properties and three-dimensional shapes of predicted conformers. The results highlight additional similarities and differences between the sets, but also the dependence of such comparisons on the choice of molecular descriptors. Using a protein-binding dataset, we introduce an information-theoretic measure to assess diversity of performance with a constraint on specificity. Rather than relying on finding individual active compounds, this measure allows rational judgment of compound subsets as groups. We also apply this measure to publicly available data from ChemBank for the same compound sets across a diverse group of functional assays. We find that performance diversity of compound sets is relatively stable across a range of property values as judged by this measure, both in protein-binding studies and functional assays. Because building screening collections with improved performance depends on efficient use of synthetic organic chemistry resources, these studies illustrate an important quantitative framework to help prioritize choices made in building such collections.

A comparison of machine learning algorithms for chemical toxicity classification using a simulated multi-scale data model

BACKGROUND

Bioactivity profiling using high-throughput in vitro assays can reduce the cost and time required for toxicological screening of environmental chemicals and can also reduce the need for animal testing. Several public efforts are aimed at discovering patterns or classifiers in high-dimensional bioactivity space that predict tissue, organ or whole animal toxicological endpoints. Supervised machine learning is a powerful approach to discover combinatorial relationships in complex in vitro/in vivo datasets. We present a novel model to simulate complex chemical-toxicology data sets and use this model to evaluate the relative performance of different machine learning (ML) methods.

RESULTS

The classification performance of Artificial Neural Networks (ANN), K-Nearest Neighbors (KNN), Linear Discriminant Analysis (LDA), Naïve Bayes (NB), Recursive Partitioning and Regression Trees (RPART), and Support Vector Machines (SVM) in the presence and absence of filter-based feature selection was analyzed using K-way cross-validation testing and independent validation on simulated in vitro assay data sets with varying levels of model complexity, number of irrelevant features and measurement noise. While the prediction accuracy of all ML methods decreased as non-causal (irrelevant) features were added, some ML methods performed better than others. In the limit of using a large number of features, ANN and SVM were always in the top performing set of methods while RPART and KNN (k = 5) were always in the poorest performing set. The addition of measurement noise and irrelevant features decreased the classification accuracy of all ML methods, with LDA suffering the greatest performance degradation. LDA performance is especially sensitive to the use of feature selection. Filter-based feature selection generally improved performance, most strikingly for LDA.

CONCLUSION

We have developed a novel simulation model to evaluate machine learning methods for the analysis of data sets in which in vitro bioassay data is being used to predict in vivo chemical toxicology. From our analysis, we can recommend that several ML methods, most notably SVM and ANN, are good candidates for use in real world applications in this area.

eHiTS-to-VMD interface application. The search for tyrosine-tRNA ligase inhibitors

Owing to the recent development of virtual high-throughput screening (vHTS) and a vast number of compounds subjected to vHTS analyses, it has been essential to automate the processing of computational data required for the analysis and visualization of research results. Using the search for tyrosine-tRNA ligase inhibitors as an example, we present a computer application, an interface between eHiTS software for virtual high-throughput screening and VMD graphic software used to visualize calculation results.

A Cheminformatic Toolkit for Mining Biomedical Knowledge

PURPOSE

Cheminformatics can be broadly defined to encompass any activity related to the application of information technology to the study of properties, effects and uses of chemical agents. One of the most important current challenges in cheminformatics is to allow researchers to search databases of biomedical knowledge, using chemical structures as input.

MATERIALS AND METHODS

An important step towards this goal was the establishment of PubChem, an open, centralized database of small molecules accessible through the World Wide Web. While PubChem is primarily intended to serve as a repository for high throughput screening data from federally-funded screening centers and academic research laboratories, the major impact of PubChem could also reside in its ability to serve as a chemical gateway to biomedical databases such as PubMed.

CONCLUSION

This article will review cheminformatic tools that can be applied to facilitate annotation of PubChem through links to the scientific literature; to integrate PubChem with transcriptomic, proteomic, and metabolomic datasets; to incorporate results of numerical simulations of physiological systems into PubChem annotation; and ultimately, to translate data of chemical genomics screening efforts into information that will benefit biomedical researchers and physician scientists across all therapeutic areas.

The ToxCast program for prioritizing toxicity testing of environmental chemicals

The U.S. Environmental Protection Agency (EPA) is developing methods for utilizing computational chemistry, high-throughput screening (HTS), and various toxicogenomic technologies to predict potential for toxicity and prioritize limited testing resources toward chemicals that likely represent the greatest hazard to human health and the environment. This chemical prioritization research program, entitled "ToxCast," is being initiated with the purpose of developing the ability to forecast toxicity based on bioactivity profiling. The proof-of-concept phase of ToxCast will focus upon chemicals with an existing, rich toxicological database in order to provide an interpretive context for the ToxCast data. This set of several hundred reference chemicals will represent numerous structural classes and phenotypic outcomes, including tumorigens, developmental and reproductive toxicants, neurotoxicants, and immunotoxicants. The ToxCast program will evaluate chemical properties and bioactivity profiles across a broad spectrum of data domains: physical-chemical, predicted biological activities based on existing structure-activity models, biochemical properties based on HTS assays, cell-based phenotypic assays, and genomic and metabolomic analyses of cells. These data will be generated through a series of external contracts, along with collaborations across EPA, with the National Toxicology Program, and with the National Institutes of Health Chemical Genomics Center. The resulting multidimensional data set provides an informatics challenge requiring appropriate computational methods for integrating various chemical, biological, and toxicological data into profiles and models predicting toxicity.