Computational models to predict endocrine-disrupting chemical binding with androgen or oestrogen receptors

Screening of the Antagonistic Activity of Potential Bisphenol A Alternatives toward the Androgen Receptor Using Machine Learning and Molecular Dynamics Simulation

Over the past few decades, extensive research has indicated that exposure to bisphenol A (BPA) increases the health risks in humans. Toxicological studies have demonstrated that BPA can bind to the androgen receptor (AR), resulting in endocrine-disrupting effects. In recent investigations, many alternatives to BPA have been detected in various environmental media as major pollutants. However, related experimental evaluations of BPA alternatives have not been systematically implemented for the assessment of chemical safety and the effects of structural characteristics on the antagonistic activity of the AR. To promote the green development of BPA alternatives, high-throughput toxicological screening is fundamental for prioritizing chemical tests. Therefore, we proposed a hybrid deep learning architecture that combines molecular descriptors and molecular graphs to predict AR antagonistic activity. Compared to previous models, this hybrid architecture can extract substantial chemical information from various molecular representations to improve the model's generalization ability for BPA alternatives. Our predictions suggest that lignin-derivable bisguaiacols, as alternatives to BPA, are likely to be nonantagonist for AR compared to bisphenol analogues. Additionally, molecular dynamics (MD) simulations identified the dihydrotestosterone-bound pocket, rather than the surface, as the major binding site of bisphenol analogues. The conformational changes of key helix H12 from an agonistic to an antagonistic conformation can be evaluated qualitatively by accelerated MD simulations to explain the underlying mechanism. Overall, our computational study is helpful for toxicological screening of BPA alternatives and the design of environmentally friendly BPA alternatives.

Machine Learning Methods for Endocrine Disrupting Potential Identification Based on Single-Cell Data

Humans are continuously exposed to a variety of toxicants and chemicals which is exacerbated during and after environmental catastrophes such as floods, earthquakes, and hurricanes. The hazardous chemical mixtures generated during these events threaten the health and safety of humans and other living organisms. This necessitates the development of rapid decision-making tools to facilitate mitigating the adverse effects of exposure on the key modulators of the endocrine system, such as the estrogen receptor alpha (ERα), for example. The mechanistic stages of the estrogenic transcriptional activity can be measured with high content/high throughput microscopy-based biosensor assays at the single-cell level, which generates millions of object-based minable data points. By combining computational modeling and experimental analysis, we built a highly accurate data-driven classification framework to assess the endocrine disrupting potential of environmental compounds. The effects of these compounds on the ERα pathway are predicted as being receptor agonists or antagonists using the principal component analysis (PCA) projections of high throughput, high content image analysis descriptors. The framework also combines rigorous preprocessing steps and nonlinear machine learning algorithms, such as the Support Vector Machines and Random Forest classifiers, to develop highly accurate mathematical representations of the separation between ERα agonists and antagonists. The results show that Support Vector Machines classify the unseen chemicals correctly with more than 96% accuracy using the proposed framework, where the preprocessing and the PCA steps play a key role in suppressing experimental noise and unraveling hidden patterns in the dataset.

Water a major source of endocrine-disrupting chemicals: An overview on the occurrence, implications on human health and bioremediation strategies

Endocrine disrupting chemicals (EDCs) are toxic compounds that occur naturally or are the output of anthropogenic activities that negatively impact both humans and wildlife. A number of diseases are associated with these disruptors, including reproductive disorders, cardiovascular disorders, kidney disease, neurological disorders, autoimmune disorders, and cancer. Due to their integral role in pharmaceuticals and cosmetics, packaging companies, agro-industries, pesticides, and plasticizers, the scientific awareness on natural and artificial EDCs are increasing. As these xenobiotic compounds tend to bioaccumulate in body tissues and may also persist longer in the environment, the concentrations of these organic compounds may increase far from their original point of concentrations. Water remains as the major sources of how humans and animals are exposed to EDCs. However, these toxic compounds cannot be completely biodegraded nor bioremediated from the aqueous medium with conventional treatment strategies thereby requiring much more efficient strategies to combat EDC contamination. Recently, genetically engineered microorganism, genome editing, and the knowledge of protein and metabolic engineering has revolutionized the field of bioremediation thereby helping to breakdown EDCs effectively. This review shed lights on understanding the importance of aquatic mediums as a source of EDCs exposure. Furthermore, the review sheds light on the consequences of these EDCs on human health as well as highlights the importance of different remediation and bioremediation approaches. Particular attention is paid to the recent trends and perspectives in order to attain sustainable approaches to the bioremediation of EDCs. Additionally, rigorous restrictions to preclude the discharge of estrogenic chemicals into the environment should be followed in efforts to combat EDC pollution.

Ecotoxicological QSAR modelling of the acute toxicity of fused and non-fused polycyclic aromatic hydrocarbons (FNFPAHs) against two aquatic organisms: Consensus modelling and comparison with ECOSAR

Fused and non-fused polycyclic aromatic hydrocarbons (FNFPAHs) are a type of organic compounds widely occurring in the environment that pose a potential hazard to ecosystem and public health, and thus receive extensive attention from various regulatory agencies. Here, quantitative structure-activity relationship (QSAR) models were constructed to model the ecotoxicity of FNFPAHs against two aquatic species, Daphnia magna and Oncorhynchus mykiss. According to the stringent OECD guidelines, we used genetic algorithm (GA) plus multiple linear regression (MLR) approach to establish QSAR models of the two aquatic toxicity endpoints: D. magna (48 h LC₅₀) and O. mykiss (96 h LC₅₀). The models were established using simple 2D descriptors with explicit physicochemical significance and evaluated using various internal/external validation metrics. The results clearly show that both models are statistically robust (Q_LOO² = 0.7834 for D. magna and Q_LOO² = 0.8162 for O. mykiss), have good internal fitness (R² = 0.8159 for D. magna and R² = 0.8626 for O. mykiss and external predictive ability (D. magna: R_test² = 0.8259, Q_Fn² = 0.7640∼0.8140, CCC_test = 0.8972; O. mykiss:R_test² = 0.8077, Q_Fn² = 0.7615∼0.7722, CCC_test = 0.8910). To prove the predictive performance of the developed models, an additional comparison with the standard ECOSAR tool obviously shows that our models have lower RMSE values. Subsequently, we utilized the best models to predict the true external set compounds collected from the PPDB database to further fill the toxicity data gap. In addition, consensus models (CMs) that integrate all validated individual models (IMs) were more externally predictive than IMs, of which CM2 has the best prediction performance towards the two aquatic species. Overall, the models presented here could be used to evaluate unknown FNFPAHs inside the domain of applicability (AD), thus being very important for environmental risk assessment under current regulatory frameworks.

Predicting bioconcentration factor and estrogen receptor bioactivity of bisphenol a and its analogues in adult zebrafish by directed message passing neural networks

The bioconcentration factor (BCF) is a key parameter for bioavailability assessment of environmental pollutants in regulatory frameworks. The comparative toxicology and mechanism of action of congeners are also of concern. However, there are limitations to acquire them by conducting field and laboratory experiments while machinelearning is emerging as a promising predictive tool to fill the gap. In this study, the Direct Message Passing Neural Network (DMPNN) was applied to predict logBCFs of bisphenol A (BPA) and its four analogues (bisphenol AF (BPAF), bisphenol B (BPB), bisphenol F (BPF) and bisphenol S (BPS)). For the test set, the Pearson correlation coefficient (PCC) and mean square error (MSE) were 0.85 and 0.52 respectively, suggesting a good predictive performance. The predicted logBCFs values by the DMPNN ranging from 0.35 (BPS) to 2.14 (BPAF) coincided well with those by the classical EPI Suite (BCFBAF model). Besides, estrogen receptor α (ERα) bioactivity of these bisphenols was also predicted well by the DMPNN, with a probability of 97.0 % (BPB) to 99.7 % (BPAF), which was validated by the extent of vitellogenin (VTG) induction in male zebrafish as a biomarker except BPS. Thus, with little need for expert knowledge, DMPNN is confirmed to be a useful tool to accurately predict logBCF and screen for estrogenic activity from molecular structures. Moreover, a gender difference was noted in the changes of three endpoints (logBCF, ER binding affinity and VTG levels), the rank order of which was BPAF > BPB > BPA > BPF > BPS consistently, and abnormal amino acid metabolism is featured as an omics signature of abnormal hormone protein expression.

SApredictor: An Expert System for Screening Chemicals Against Structural Alerts

The rapid and accurate evaluation of chemical toxicity is of great significance for estimation of chemical safety. In the past decades, a great number of excellent computational models have been developed for chemical toxicity prediction. But most machine learning models tend to be “black box”, which bring about poor interpretability. In the present study, we focused on the identification and collection of structural alerts (SAs) responsible for a series of important toxicity endpoints. Then, we carried out effective storage of these structural alerts and developed a web-server named SApredictor (www.sapredictor.cn) for screening chemicals against structural alerts. People can quickly estimate the toxicity of chemicals with SApredictor, and the specific key substructures which cause the chemical toxicity will be intuitively displayed to provide valuable information for the structural optimization by medicinal chemists.

In Silico Prediction of Chemical Aquatic Toxicity by Multiple Machine Learning and Deep Learning Approaches

Fish is one of the model animals used to evaluate the adverse effects of a chemical exposed to the ecosystem. However, its low throughput and relevantly high expense make it impossible to test all new chemicals in manufacture. Hence using in silico models to prioritize compounds to be tested has been widely applied in environmental risk assessment and drug discovery. In this study, we constructed the local predictive models for four fish species, including bluegill sunfish, rainbow trout, fathead minnow, and sheepshead minnow, and the global models with all four fish data. A total of 1874 unique compounds with their labels, i.e. toxic (LC₅₀ < 10 ppm) or nontoxic were collected from ECOTOX and literature. Both conventional machine learning methods and the deep learning architecture, graph convolutional network (GCN), were used to build predictive models. The classification accuracy of the best local model for each fish species was higher than 0.83. For the global models, two strategies including consistency prediction and probability threshold were adopted to improve the predictive capability at the cost of limiting applicability domain. For 63% of compounds in domain, the accuracy was around 0.97. By comparison of the deep learning and machine learning methods, we found that the single-task GCN showed specific advantages in performance and multi-task GCN showed no advantages over the conventional machine learning methods. The data and models are available on GitHub (https://github.com/ChemPredict/ChemicalAquaticToxicity).

Graphene oxide enhanced the endocrine disrupting effects of bisphenol A in adult male zebrafish: Integrated deep learning and metabolomics studies

In our previous studies, it was found that graphene oxide (GO) reduced the endocrine disruption of bisphenol A (BPA) in zebrafish embryo and larvae, but through different mechanisms. In this study, adult male zebrafish were selected to further understand the interactions between GO and BPA considering that adult zebrafish have different uptake pathways and metabolism from embryo and larvae. BPA was predicted to bind with the estrogen receptor α (ERα) with a probability of 98.1% by training a directed-message passing deep neural network model, and was confirmed by molecular docking analysis. The results were in accordance with the significantly increased vitellogenin (VTG) and estradiol (E2) levels, while decreased testosterone (T) and follicle-stimulating hormone (FSH) levels in the adult male zebrafish after 7 d exposure to 500 μg/L BPA. Compared to BPA single exposure group, the presence of GO led to significantly lower T and FSH levels and fewer spermatozoa, indicating that GO enhanced the endocrine disruption effects of BPA in the adult zebrafish. Metabolomics analysis revealed that 5 μg/L BPA could elicit changes in the metabolome, and the responses were correlated with BPA concentrations. Metabolic pathway analysis revealed more disturbance was caused by the mixture of GO and BPA compared to BPA alone, including three additional pathways and stronger perturbations on carbohydrate, lipid, and amino acid metabolism, fortifying that GO exaggerated the toxic effects of BPA. This was opposite to the depression effect observed in zebrafish embryo and larvae, magnifying that the joint effects of exposure to nanomaterials and endocrine disrupting chemicals are relevant to the life stages of organisms.

Review of in silico studies dedicated to the nuclear receptor family: Therapeutic prospects and toxicological concerns

Being in the center of both therapeutic and toxicological concerns, NRs are widely studied for drug discovery application but also to unravel the potential toxicity of environmental compounds such as pesticides, cosmetics or additives. High throughput screening campaigns (HTS) are largely used to detect compounds able to interact with this protein family for both therapeutic and toxicological purposes. These methods lead to a large amount of data requiring the use of computational approaches for a robust and correct analysis and interpretation. The output data can be used to build predictive models to forecast the behavior of new chemicals based on their in vitro activities. This atrticle is a review of the studies published in the last decade and dedicated to NR ligands in silico prediction for both therapeutic and toxicological purposes. Over 100 articles concerning 14 NR subfamilies were carefully read and analyzed in order to retrieve the most commonly used computational methods to develop predictive models, to retrieve the databases deployed in the model building process and to pinpoint some of the limitations they faced.

Screening of potential oestrogen receptor α agonists in pesticides via in silico, in vitro and in vivo methods

In modern agricultural management, the use of pesticides is indispensable. Due to their massive use worldwide, pesticides represent a latent risk to both humans and the environment. In the present study, 1056 frequently used pesticides were screened for oestrogen receptor (ER) agonistic activity by using in silico methods. We found that 72 and 47 pesticides potentially have ER agonistic activity by the machine learning methods random forest (RF) and deep neural network (DNN), respectively. Among endocrine-disrupting chemicals (EDCs), 14 have been reported as EDCs or ER agonists by previous studies. We selected 3 reported and 7 previously unreported pesticides from 76 potential ER agonists to further assess ERα agonistic activity. All 10 selected pesticides exhibited ERα agonistic activity in human cells or zebrafish. In the dual-luciferase reporter gene assays, six pesticides exhibited ERα agonistic activity. Additionally, nine pesticides could induce mRNA expression of the pS2 and NRF1 genes in MCF-7 cells, and seven pesticides could induce mRNA expression of the vtg1 and vtg2 genes in zebrafish. Importantly, the remaining 48 out of 76 potential ER agonists, none of which have previously been reported to have endocrine-disrupting effects or oestrogenic activity, should be of great concern. Our screening results can inform environmental protection goals and play an important role in environmental protection and early warnings to human health.

Classification of estrogenic compounds by coupling high content analysis and machine learning algorithms

Environmental toxicants affect human health in various ways. Of the thousands of chemicals present in the environment, those with adverse effects on the endocrine system are referred to as endocrine-disrupting chemicals (EDCs). Here, we focused on a subclass of EDCs that impacts the estrogen receptor (ER), a pivotal transcriptional regulator in health and disease. Estrogenic activity of compounds can be measured by many in vitro or cell-based high throughput assays that record various endpoints from large pools of cells, and increasingly at the single-cell level. To simultaneously capture multiple mechanistic ER endpoints in individual cells that are affected by EDCs, we previously developed a sensitive high throughput/high content imaging assay that is based upon a stable cell line harboring a visible multicopy ER responsive transcription unit and expressing a green fluorescent protein (GFP) fusion of ER. High content analysis generates voluminous multiplex data comprised of minable features that describe numerous mechanistic endpoints. In this study, we present a machine learning pipeline for rapid, accurate, and sensitive assessment of the endocrine-disrupting potential of benchmark chemicals based on data generated from high content analysis. The multidimensional imaging data was used to train a classification model to ultimately predict the impact of unknown compounds on the ER, either as agonists or antagonists. To this end, both linear logistic regression and nonlinear Random Forest classifiers were benchmarked and evaluated for predicting the estrogenic activity of unknown compounds. Furthermore, through feature selection, data visualization, and model discrimination, the most informative features were identified for the classification of ER agonists/antagonists. The results of this data-driven study showed that highly accurate and generalized classification models with a minimum number of features can be constructed without loss of generality, where these machine learning models serve as a means for rapid mechanistic/phenotypic evaluation of the estrogenic potential of many chemicals.

Chemical contaminants or toxicants pose environmental and health-related risks for exposure. The ability to rapidly understand their biological impact, specifically on a key modulator of important physiological and pathological states in the human body is essential for diagnosing and avoiding undesirable health outcomes during environmental emergencies. In this study, we use advanced data analytics for creating statistical models that can accurately predict the endocrinological activity of toxic chemicals based on high throughput/high content image analysis data. We focus on a subclass of chemicals that affect the estrogen receptor (ER), which is a pivotal transcriptional regulator in health and disease. The multidimensional imaging data of these benchmark chemicals are used to train a classification model to ultimately predict the impact of unknown compounds on the ER, either as agonists or antagonists. To this end, we evaluate linear and nonlinear classifiers for predicting the estrogenic activity of unknown compounds and use feature selection, data visualization, and model discrimination methodologies to identify the most informative features for the classification of ER agonists/antagonists.

Discriminant models on mitochondrial toxicity improved by consensus modeling and resolving imbalance in training

Humans and animals may be exposed to tens of thousands of natural and synthetic chemicals during their lifespan. It is difficult to assess risk for all the chemicals with experimental toxicity tests. An alternative approach is to use computational toxicology methods such as quantitative structure-activity relationship (QSAR) modeling. Mitochondrial toxicity is involved in many diseases such as cancer, neurodegeneration, type 2 diabetes, cardiovascular diseases and autoimmune diseases. Thus, it is important to rapidly and efficiently identify chemicals with mitochondrial toxicity. In this study, five machine learning algorithms and twelve types of molecular fingerprints were employed to generate QSAR discriminant models for mitochondrial toxicity. A threshold moving method was adopted to resolve the imbalance issue in the training data. Consensus of the models by an averaging probability strategy improved prediction performance. The best model has correct classification rates of 81.8% and 88.3% in ten-fold cross validation and external validation, respectively. Substructures such as phenol, carboxylic acid, nitro and arylchloride were found informative through analysis of information gain and frequency of substructures. The results demonstrate that resolving imbalance in training and building consensus models can improve classification rates for mitochondrial toxicity prediction.

In silico prediction of chemical neurotoxicity using machine learning

Neurotoxicity is one of the main causes of drug withdrawal, and the biological experimental methods of detecting neurotoxic toxicity are time-consuming and laborious. In addition, the existing computational prediction models of neurotoxicity still have some shortcomings. In response to these shortcomings, we collected a large number of data set of neurotoxicity and used PyBioMed molecular descriptors and eight machine learning algorithms to construct regression prediction models of chemical neurotoxicity. Through the cross-validation and test set validation of the models, it was found that the extra-trees regressor model had the best predictive effect on neurotoxicity ([Formula: see text] = 0.784). In addition, we get the applicability domain of the models by calculating the standard deviation distance and the lever distance of the training set. We also found that some molecular descriptors are closely related to neurotoxicity by calculating the contribution of the molecular descriptors to the models. Considering the accuracy of the regression models, we recommend using the extra-trees regressor model to predict the chemical autonomic neurotoxicity.

Quantitative structure-property relationship of the photoelectrochemical oxidation of phenolic pollutants at modified nanoporous titanium oxide using supervised machine learning

Here we report on an advanced photoelectrochemical (PEC) oxidation of 22 phenolic pollutants based on modified nanoporous TiO2, which was directly grown on a titanium substrate electrochemically. Their degradation rate constants were experimentally determined and their physicochemical properties were computaionally calculated. The quantitative structure-property relationship (QSPR) was elucidated by employing multiple linear regression (MLR) method. A supervised machine learning approach was employed to build QSPR models. The high predictive abilities of the QSPR model were validated via leave-one-out (LOO) method and a strict regimen of statistical validation tests. The significant descriptors identified in the QSPR Model for the phenolic compounds were also assessed using a typical dye pollutant Rhodamine B, further confirming the high effectiveness and predictability of the optimized model. Our study has shown that the integrated effect of the structural, hydrophobic and topological properties along with electronic property should be considered in order to design an efficient PEC catalytic approach for environmental applications.

Modeling and insights into molecular basis of low molecular weight respiratory sensitizers

Respiratory sensitization has been considered an important toxicological endpoint, because of the severe risk to human health. A great part of sensitization events were caused by low molecular weight (< 1000) respiratory sensitizers in the past decades. However, there is currently no widely accepted test method that can identify prospective low molecular weight respiratory sensitisers. Herein, we performed the study of modeling and insights into molecular basis of low molecular weight respiratory sensitizers with a high-quality data set containing 136 respiratory sensitizers and 518 nonsensitizers. We built a number of classification models by using OCHEM tools, and a consensus model was developed based on the ten best individual models. The consensus model showed good predictive ability with a balanced accuracy of 0.78 and 0.85 on fivefold cross-validation and external validation, respectively. The readers can predict the respiratory sensitization of organic compounds via https://ochem.eu/article/114857 . The effect of several molecular properties on respiratory sensitization was also evaluated. The results indicated that these properties differ significantly between respiratory sensitizers and nonsensitizers. Furthermore, 14 privileged substructures responsible for respiratory sensitization were identified. We hope the models and the findings could provide useful help for environmental risk assessment.

Prediction on the mutagenicity of nitroaromatic compounds using quantum chemistry descriptors based QSAR and machine learning derived classification methods

Nitroaromatic compounds (NACs) are an important type of environmental organic pollutants. However, it is lack of sufficient information relating to their potential adverse effects on human health and the environment due to the limited resources. Thus, using in silico technologies to assess their potential hazardous effects is urgent and promising. In this study, quantitative structure activity relationship (QSAR) and classification models were constructed using a set of NACs based on their mutagenicity against Salmonella typhimurium TA100 strain. For QSAR studies, DRAGON descriptors together with quantum chemistry descriptors were calculated for characterizing the detailed molecular information. Based on genetic algorithm (GA) and multiple linear regression (MLR) analyses, we screened descriptors and developed QSAR models. For classification studies, seven machine learning methods along with six molecular fingerprints were applied to develop qualitative classification models. The goodness of fitting, reliability, robustness and predictive performance of all developed models were measured by rigorous statistical validation criteria, then the best QSAR and classification models were chosen. Moreover, the QSAR models with quantum chemistry descriptors were compared to that without quantum chemistry descriptors and previously reported models. Notably, we also obtained some specific molecular properties or privileged substructures responsible for the high mutagenicity of NACs. Overall, the developed QSAR and classification models can be utilized as potential tools for rapidly predicting the mutagenicity of new or untested NACs for environmental hazard assessment and regulatory purposes, and may provide insights into the in vivo toxicity mechanisms of NACs and related compounds.

Determination of Endocrine Disruption Potential of Bisphenol A Alternatives in Food Contact Materials Using In Vitro Assays: State of the Art and Future Challenges

Alternatives to bisphenol A (BPA) are developed for food contact materials as a result of increasing evidence of exposure-correlated harmful effects of BPA. In vitro assays provide the fast, affordable, and mechanism insightful ways to screen endocrine disruption (ED), which is a major concern of new BPA alternatives. In this review, we summarize the safety and regulation information on the alternatives to BPA, review the state of the art of in vitro assays for ED evaluation, highlight their advantages and limitations, and discuss the challenges and future research needs. Our review shows that ligand binding, reporter gene, cell proliferation, and steroidogenesis are four commonly used in vitro assays to determine the ED at the response of receptor, gene transcription, and whole cell level. Major challenges are found from in vitro-in vivo translation and identification of ED chemicals in polymers. More studies on these areas are needed in the future.

Machine Learning Consensus To Predict the Binding to the Androgen Receptor within the CoMPARA Project

The nuclear androgen receptor (AR) is one of the most relevant biological targets of Endocrine Disrupting Chemicals (EDCs), which produce adverse effects by interfering with hormonal regulation and endocrine system functioning. This paper describes novel in silico models to identify organic AR modulators in the context of the Collaborative Modeling Project of Androgen Receptor Activity (CoMPARA), coordinated by the National Center of Computational Toxicology (U.S. Environmental Protection Agency). The collaborative project involved 35 international research groups to prioritize the experimental tests of approximatively 40k compounds, based on the predictions provided by each participant. In this paper, we describe our machine learning approach to predict the binding to AR, which is based on a consensus of a multivariate Bernoulli Naive Bayes, a Random Forest, and N-Nearest Neighbor classification models. The approach was developed in compliance with the Organization of Economic Cooperation and Development (OECD) principles, trained on 1687 ToxCast molecules classified according to 11 in vitro assays, and further validated on a set of 3,882 external compounds. The models provided robust and reliable predictions and were used to gather novel data-driven insights on the structural features related to AR binding, agonism, and antagonism.

Computational QSAR model combined molecular descriptors and fingerprints to predict HDAC1 inhibitors

The dynamic balance between acetylation and deacetylation of histones plays a crucial role in the epigenetic regulation of gene expression. It is equilibrated by two families of enzymes: histone acetyltransferases and histone deacetylases (HDACs). HDACs repress transcription by regulating the conformation of the higher-order chromatin structure. HDAC inhibitors have recently become a class of chemical agents for potential treatment of the abnormal chromatin remodeling process involved in certain cancers. In this study, we constructed a large dataset to predict the activity value of HDAC1 inhibitors. Each compound was represented with seven fingerprints, and computational models were subsequently developed to predict HDAC1 inhibitors via five machine learning methods. These methods include naïve Bayes, κ-nearest neighbor, C4.5 decision tree, random forest, and support vector machine (SVM) algorithms. The best predicting model was CDK fingerprint with SVM, which exhibited an accuracy of 0.89. This model also performed best in five-fold cross-validation. Some representative substructure alerts responsible for HDAC1 inhibitors were identified by using MoSS in KNIME, which could facilitate the identification of HDAC1 inhibitors.

In Silico Prediction of O⁶-Methylguanine-DNA Methyltransferase Inhibitory Potency of Base Analogs with QSAR and Machine Learning Methods

O⁶-methylguanine-DNA methyltransferase (MGMT), a unique DNA repair enzyme, can confer resistance to DNA anticancer alkylating agents that modify the O⁶-position of guanine. Thus, inhibition of MGMT activity in tumors has a great interest for cancer researchers because it can significantly improve the anticancer efficacy of such alkylating agents. In this study, we performed a quantitative structure activity relationship (QSAR) and classification study based on a total of 134 base analogs related to their ED50 values (50% inhibitory concentration) against MGMT. Molecular information of all compounds were described by quantum chemical descriptors and Dragon descriptors. Genetic algorithm (GA) and multiple linear regression (MLR) analysis were combined to develop QSAR models. Classification models were generated by seven machine-learning methods based on six types of molecular fingerprints. Performances of all developed models were assessed by internal and external validation techniques. The best QSAR model was obtained with Q²Loo = 0.83, R² = 0.87, Q²ext = 0.67, and R²ext = 0.69 based on 84 compounds. The results from QSAR studies indicated topological charge indices, polarizability, ionization potential (IP), and number of primary aromatic amines are main contributors for MGMT inhibition of base analogs. For classification studies, the accuracies of 10-fold cross-validation ranged from 0.750 to 0.885 for top ten models. The range of accuracy for the external test set ranged from 0.800 to 0.880 except for PubChem-Tree model, suggesting a satisfactory predictive ability. Three models (Ext-SVM, Ext-Tree and Graph-RF) showed high and reliable predictive accuracy for both training and external test sets. In addition, several representative substructures for characterizing MGMT inhibitors were identified by information gain and substructure frequency analysis method. Our studies might be useful for further study to design and rapidly identify potential MGMT inhibitors.

Multiclassification Prediction of Enzymatic Reactions for Oxidoreductases and Hydrolases Using Reaction Fingerprints and Machine Learning Methods

Drug metabolism is a complex procedure in the human body, including a series of enzymatically catalyzed reactions. However, it is costly and time consuming to investigate drug metabolism experimentally; computational methods are hence developed to predict drug metabolism and have shown great advantages. As the first step, classification of metabolic reactions and enzymes is highly desirable for drug metabolism prediction. In this study, we developed multiclassification models for prediction of reaction types catalyzed by oxidoreductases and hydrolases, in which three reaction fingerprints were used to describe the reactions and seven machine learnings algorithms were employed for model building. Data retrieved from KEGG containing 1055 hydrolysis and 2510 redox reactions were used to build the models, respectively. The external validation data consisted of 213 hydrolysis and 512 redox reactions extracted from the Rhea database. The best models were built by neural network or logistic regression with a 2048-bit transformation reaction fingerprint. The predictive accuracies of the main class, subclass, and superclass classification models on external validation sets were all above 90%. This study will be very helpful for enzymatic reaction annotation and further study on metabolism prediction.

In Silico Prediction of Chemical Toxicity for Drug Design Using Machine Learning Methods and Structural Alerts

During drug development, safety is always the most important issue, including a variety of toxicities and adverse drug effects, which should be evaluated in preclinical and clinical trial phases. This review article at first simply introduced the computational methods used in prediction of chemical toxicity for drug design, including machine learning methods and structural alerts. Machine learning methods have been widely applied in qualitative classification and quantitative regression studies, while structural alerts can be regarded as a complementary tool for lead optimization. The emphasis of this article was put on the recent progress of predictive models built for various toxicities. Available databases and web servers were also provided. Though the methods and models are very helpful for drug design, there are still some challenges and limitations to be improved for drug safety assessment in the future.

In silico prediction of pesticide aquatic toxicity with chemical category approaches

Aquatic toxicity is an important issue in pesticide development. In this study, using nine molecular fingerprints to describe pesticides, binary and ternary classification models were constructed to predict aquatic toxicity of pesticides via six machine learning methods: Naïve Bayes (NB), Artificial Neural Network (ANN), k-Nearest Neighbor (kNN), Classification Tree (CT), Random Forest (RF) and Support Vector Machine (SVM). For the binary models, local models were obtained with 829 pesticides on rainbow trout (RT) and 151 pesticides on lepomis (LP), and global models were constructed on the basis of 1258 diverse pesticides on RT and LP and 278 on other fish species. After analyzing the local binary models, we found that fish species caused influence in terms of accuracy. Considering the data size and predictive range, the 1258 pesticides were also used to build global ternary models. The best local binary models were Maccs_ANN for RT and Maccs_SVM for LP, which exhibited accuracies of 0.90 and 0.90, respectively. For global binary models, the best model was Graph_SVM with an accuracy of 0.89. Accuracy of the best global ternary model Graph_SVM was 0.81, which was a little lower than that of the best global binary model. In addition, several substructural alerts were identified including nitrobenzene, chloroalkene and nitrile, which could significantly correlate with pesticide aquatic toxicity. This study provides a useful tool for an early evaluation of pesticide aquatic toxicity in environmental risk assessment.

Development of estrogen receptor beta binding prediction model using large sets of chemicals

We developed an ER_β binding prediction model to facilitate identification of chemicals specifically bind ER_β or ER_α together with our previously developed ER_α binding model. Decision Forest was used to train ER_β binding prediction model based on a large set of compounds obtained from EADB. Model performance was estimated through 1000 iterations of 5-fold cross validations. Prediction confidence was analyzed using predictions from the cross validations. Informative chemical features for ER_β binding were identified through analysis of the frequency data of chemical descriptors used in the models in the 5-fold cross validations. 1000 permutations were conducted to assess the chance correlation. The average accuracy of 5-fold cross validations was 93.14% with a standard deviation of 0.64%. Prediction confidence analysis indicated that the higher the prediction confidence the more accurate the predictions. Permutation testing results revealed that the prediction model is unlikely generated by chance. Eighteen informative descriptors were identified to be important to ER_β binding prediction. Application of the prediction model to the data from ToxCast project yielded very high sensitivity of 90-92%. Our results demonstrated ER_β binding of chemicals could be accurately predicted using the developed model. Coupling with our previously developed ER_α prediction model, this model could be expected to facilitate drug development through identification of chemicals that specifically bind ER_β or ER_α.

Integration of in silico methods and computational systems biology to explore endocrine-disrupting chemical binding with nuclear hormone receptors

Thousands of potential endocrine-disrupting chemicals present difficult regulatory challenges. Endocrine-disrupting chemicals can interfere with several nuclear hormone receptors associated with a variety of adverse health effects. The U.S. Environmental Protection Agency (U.S. EPA) has released its reviews of Tier 1 screening assay results for a set of pesticides in the Endocrine Disruptor Screening Program (EDSP), and recently, the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP) data. In this study, the predictive ability of QSAR and docking approaches is evaluated using these data sets. This study also presents a computational systems biology approach using carbaryl (1-naphthyl methylcarbamate) as a case study. For estrogen receptor and androgen receptor binding predictions, two commercial and two open source QSAR tools were used, as was the publicly available docking tool Endocrine Disruptome. For estrogen receptor binding predictions, the ADMET Predictor, VEGA, and OCHEM models (specificity: 0.88, 0.88, and 0.86, and accuracy: 0.81, 0.84, and 0.88, respectively) were each more reliable than the MetaDrug™ model (specificity 0.81 and accuracy 0.77). For androgen receptor binding predictions, the Endocrine Disruptome and ADMET Predictor models (specificity: 0.94 and 0.8, and accuracy: 0.78 and 0.71, respectively) were more reliable than the MetaDrug™ model (specificity 0.33 and accuracy 0.4). A consensus approach is proposed that reaches general agreement among the models (specificity 0.94 and accuracy 0.89). This study integrates QSAR, docking, and systems biology approaches as a virtual screening tool for use in risk assessment. As such, this systems biology pathways and network analysis approach provides a means to more critically assess the potential effects of endocrine-disrupting chemicals.

In Silico Prediction of Chemicals Binding to Aromatase with Machine Learning Methods

Environmental chemicals may affect endocrine systems through multiple mechanisms, one of which is via effects on aromatase (also known as CYP19A1), an enzyme critical for maintaining the normal balance of estrogens and androgens in the body. Therefore, rapid and efficient identification of aromatase-related endocrine disrupting chemicals (EDCs) is important for toxicology and environment risk assessment. In this study, on the basis of the Tox21 10K compound library, in silico classification models for predicting aromatase binders/nonbinders were constructed by machine learning methods. To improve the prediction ability of the models, a combined classifier (CC) strategy that combines different independent machine learning methods was adopted. Performances of the models were measured by test and external validation sets containing 1336 and 216 chemicals, respectively. The best model was obtained with the MACCS (Molecular Access System) fingerprint and CC method, which exhibited an accuracy of 0.84 for the test set and 0.91 for the external validation set. Additionally, several representative substructures for characterizing aromatase binders, such as ketone, lactone, and nitrogen-containing derivatives, were identified using information gain and substructure frequency analysis. Our study provided a systematic assessment of chemicals binding to aromatase. The built models can be helpful to rapidly identify potential EDCs targeting aromatase.

A ternary classification using machine learning methods of distinct estrogen receptor activities within a large collection of environmental chemicals

Endocrine-disrupting chemicals (EDCs), which can threaten ecological safety and be harmful to human beings, have been cause for wide concern. There is a high demand for efficient methodologies for evaluating potential EDCs in the environment. Herein an evaluation platform was developed using novel and statistically robust ternary models via different machine learning models (i.e., linear discriminant analysis, classification and regression tree, and support vector machines). The platform is aimed at effectively classifying chemicals with agonistic, antagonistic, or no estrogen receptor (ER) activities. A total of 440 chemicals from the literature were selected to derive and optimize the three-class model. One hundred and nine new chemicals appeared on the 2014 EPA list for EDC screening, which were used to assess the predictive performances by comparing the E-screen results with the predicted results of the classification models. The best model was obtained using support vector machines (SVM) which recognized agonists and antagonists with accuracies of 76.6% and 75.0%, respectively, on the test set (with an overall predictive accuracy of 75.2%), and achieved a 10-fold cross-validation (CV) of 73.4%. The external predicted accuracy validated by the E-screen assay was 87.5%, which demonstrated the application value for a virtual alert for EDCs with ER agonistic or antagonistic activities. It was demonstrated that the ternary computational model could be used as a faster and less expensive method to identify EDCs that act through nuclear receptors, and to classify these chemicals into different mechanism groups.

In silico prediction of serious eye irritation or corrosion potential of chemicals

Chemical fingerprints combined with machine learning methods were used to build binary classification models for predicting the potential EC/EI of compounds.

Integrated in silico strategy for PBT assessment and prioritization under REACH

Chemicals may persist in the environment, bioaccumulate and be toxic for humans and wildlife, posing great concern. These three properties, persistence (P), bioaccumulation (B), and toxicity (T) are the key targets of the PBT-hazard assessment. The European regulation for the Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) requires assessment of PBT-properties for all chemicals that are produced or imported in Europe in amounts exceeding 10 tonnes per year, checking whether the criteria set out in REACH Annex XIII are met, so the substance should therefore be considered to have properties of very high concern. Considering how many substances can fall under the REACH regulation, there is a pressing need for new strategies to identify and screen large numbers fast and inexpensively. An efficient non-testing screening approach to identify PBT candidates is necessary, as a valuable alternative to money- and time-consuming laboratory tests and a good start for prioritization since few tools exist (e.g. the PBT profiler developed by US EPA). The aim of this work was to offer a conceptual scheme for identifying and prioritizing chemicals for further assessment and if appropriate further testing, based on their PBT-potential, using a non-testing screening approach. We integrated in silico models (using existing and developing new ones) in a final algorithm for screening and ranking PBT-potential, which uses experimental and predicted values as well as associated uncertainties. The Multi-Criteria Decision-Making (MCDM) theory was used to integrate the different values. Then we compiled a new set of data containing known PBT and non-PBT substances, in order to check how well our approach clearly differentiated compounds labeled as PBT from those labeled as non-PBT. This indicated that the integrated model distinguished between PBT from non-PBT compounds.

Prediction of Estrogenic Bioactivity of Environmental Chemical Metabolites

The US Environmental Protection Agency's (EPA) Endocrine Disruptor Screening Program (EDSP) is using in vitro data generated from ToxCast/Tox21 high-throughput screening assays to assess the endocrine activity of environmental chemicals. Considering that in vitro assays may have limited metabolic capacity, inactive chemicals that are biotransformed into metabolites with endocrine bioactivity may be missed for further screening and testing. Therefore, there is a value in developing novel approaches to account for metabolism and endocrine activity of both parent chemicals and their associated metabolites. We used commercially available software to predict metabolites of 50 parent compounds, out of which 38 chemicals are known to have estrogenic metabolites, and 12 compounds and their metabolites are negative for estrogenic activity. Three ER QSAR models were used to determine potential estrogen bioactivity of the parent compounds and predicted metabolites, the outputs of the models were averaged, and the chemicals were then ranked based on the total estrogenicity of the parent chemical and metabolites. The metabolite prediction software correctly identified known estrogenic metabolites for 26 out of 27 parent chemicals with associated metabolite data, and 39 out of 46 estrogenic metabolites were predicted as potential biotransformation products derived from the parent chemical. The QSAR models estimated stronger estrogenic activity for the majority of the known estrogenic metabolites compared to their parent chemicals. Finally, the three models identified a similar set of parent compounds as top ranked chemicals based on the estrogenicity of putative metabolites. This proposed in silico approach is an inexpensive and rapid strategy for the detection of chemicals with estrogenic metabolites and may reduce potential false negative results from in vitro assays.

Prediction of in vitro and in vivo oestrogen receptor activity using hierarchical clustering

In this study, hierarchical clustering classification models were developed to predict in vitro and in vivo oestrogen receptor (ER) activity. Classification models were developed for binding, agonist, and antagonist in vitro ER activity and for mouse in vivo uterotrophic ER binding. In vitro classification models yielded balanced accuracies ranging from 0.65 to 0.85 for the external prediction set. In vivo ER classification models yielded balanced accuracies ranging from 0.72 to 0.83. If used as additional biological descriptors for in vivo models, in vitro scores were found to increase the prediction accuracy of in vivo ER models. If in vitro activity was used directly as a surrogate for in vivo activity, the results were poor (balanced accuracy ranged from 0.49 to 0.72). Under-sampling negative compounds in the training set was found to increase the coverage (fraction of chemicals which can be predicted) and increase prediction sensitivity.

Biomarker-based drug safety assessment in the age of systems pharmacology: from foundational to regulatory science

Improved biomarker-based assessment of drug safety is needed in drug discovery and development as well as regulatory evaluation. However, identifying drug safety-related biomarkers such as genes, proteins, miRNA and single-nucleotide polymorphisms remains a big challenge. The advances of 'omics' and computational technologies such as genomics, transcriptomics, metabolomics, proteomics, systems biology, network biology and systems pharmacology enable us to explore drug actions at the organ and organismal levels. Computational and experimental systems pharmacology approaches could be utilized to facilitate biomarker-based drug safety assessment for drug discovery and development and to inform better regulatory decisions. In this article, we review the current status and advances of systems pharmacology approaches for the development of predictive models to identify biomarkers for drug safety assessment.