CERAPP: Collaborative Estrogen Receptor Activity Prediction Project

Comparing Machine Learning Models for Aromatase (P450 19A1)

Aromatase, or cytochrome P450 19A1, catalyzes the aromatization of androgens to estrogens within the body. Changes in the activity of this enzyme can produce hormonal imbalances that can be detrimental to sexual and skeletal development. Inhibition of this enzyme can occur with drugs and natural products as well as environmental chemicals. Therefore, predicting potential endocrine disruption via exogenous chemicals requires that aromatase inhibition be considered in addition to androgen and estrogen pathway interference. Bayesian machine learning methods can be used for prospective prediction from the molecular structure without the need for experimental data. Herein, the generation and evaluation of multiple machine learning models utilizing different sources of aromatase inhibition data are described. These models are applied to two test sets for external validation with molecules relevant to drug discovery from the public domain. In addition, the performance of multiple machine learning algorithms was evaluated by comparing internal five-fold cross-validation statistics of the training data. These methods to predict aromatase inhibition from molecular structure, when used in concert with estrogen and androgen machine learning models, allow for a more holistic assessment of endocrine-disrupting potential of chemicals with limited empirical data and enable the reduction of the use of hazardous substances.

Visualization and Analysis of the REACH-chemical Space with Generative Topographic Mapping

In the framework of REACH (Registration Evaluation Authorization and restriction of Chemicals) regulation, industries have generated and reported a huge amount of (eco)toxicological data on substance produced or imported in Europe. The registration procedure initiated the creation of a large REACH database of well defined (eco)toxicological properties. Here, the data distribution in the REACH chemical space was analyzed with the help of the Generative Topographic Mapping (GTM) approach. GTM generates 2-dimensional maps on which each compound is represented as a data point. The 3^rd dimension can be used in order to display a distribution of the given (eco)toxicological property, which can further be used for property assessment of new compounds projected on the map. We report the "Universal REACH map" which accommodates 11 endpoints, covering environmental fate and (eco)toxicological properties. This map demonstrates acceptable predictive performance: in cross-validation, balanced accuracy ranges from 0.60 to 0.78. The 11 endpoints profile has been computed for each REACH-registered substance. Some concerns related to acute aquatic toxicity have been identified, whereas for environmental fate and human health endpoints the amount of compounds predicted as of concern was much smaller. It has been demonstrated that superposition of several class landscapes allows to select the zones in the chemical space populated by compounds with a given (eco)toxicological profile.

Comparison of Machine Learning Models for the Androgen Receptor

The androgen receptor (AR) is a target of interest for endocrine disruption research, as altered signaling can affect normal reproductive and neurological development for generations. In an effort to prioritize compounds with alternative methodologies, the U.S. Environmental Protection Agency (EPA) used in vitro data from 11 assays to construct models of AR agonist and antagonist signaling pathways. While these EPA ToxCast AR models require in vitro data to assign a bioactivity score, Bayesian machine learning methods can be used for prospective prediction from molecule structure alone. This approach was applied to multiple types of data corresponding to the EPA's AR signaling pathway with proprietary software, Assay Central. The training performance of all machine learning models, including six other algorithms, was evaluated by internal 5-fold cross-validation statistics. Bayesian machine learning models were also evaluated with external predictions of reference chemicals to compare prediction accuracies to published results from the EPA. The machine learning model group selected for further studies of endocrine disruption consisted of continuous AC50 data from the February 2019 release of ToxCast/Tox21. These efforts demonstrate how machine learning can be used to predict AR-mediated bioactivity and can also be applied to other targets of endocrine disruption.

Development of classification models for predicting inhibition of mitochondrial fusion and fission using machine learning methods

Mitochondrial fusion and fission are processes to maintain mitochondrial function when cells respond to environment stresses. Disruption of mitochondrial fusion and fission influences cell health and can cause adverse events such as neurodegenerative disorders. It is critical to identify environmental chemicals that can disrupt mitochondrial fusion and fission. However, experimentally testing all the chemicals is not practical because experimental methods are time-consuming and costly. Quantitative structure-activity relationship (QSAR) modeling is an attractive approach for evaluation of chemicals disrupting potential on mitochondrial fusion and fission. In this study, QSAR models were developed for differentiating chemicals capable of inhibition of mitochondrial fusion and fission using machine learning algorithms (i.e. random forest, logistic regression, Bernoulli naive Bayes, and deep neural network). One hundred iterations of five-fold cross validations and external validations showed that the best model on mitochondrial fusion had area under the receiver operating characteristic curve (AUC) of 82.8% and 78.1%, respectively; and the best model for mitochondrial fission yielded AUC of 84.3% and 97.5%, respectively. Furthermore, 45 and 56 structural alerts were identified for inhibition of mitochondrial fusion and fission, respectively. The results demonstrated that the models and the structural alerts could be useful for screening chemicals that inhibit mitochondrial fusion and fission.

Machine Learning Models for Estrogen Receptor Bioactivity and Endocrine Disruption Prediction

The U.S. Environmental Protection Agency (EPA) periodically releases in vitro data across a variety of targets, including the estrogen receptor (ER). In 2015, the EPA used these data to construct mathematical models of ER agonist and antagonist pathways to prioritize chemicals for endocrine disruption testing. However, mathematical models require in vitro data prior to predicting estrogenic activity, but machine learning methods are capable of prospective prediction from the molecular structure alone. The current study describes the generation and evaluation of Bayesian machine learning models grouped by the EPA's ER agonist pathway model using multiple data types with proprietary software, Assay Central. External predictions with three test sets of in vitro and in vivo reference chemicals with agonist activity classifications were compared to previous mathematical model publications. Training data sets were subjected to additional machine learning algorithms and compared with rank normalized scores of internal five-fold cross-validation statistics. External predictions were found to be comparable or superior to previous studies published by the EPA. When assessing six additional algorithms for the training data sets, Assay Central performed similarly at a reduced computational cost. This study demonstrates that machine learning can prioritize chemicals for future in vitro and in vivo testing of ER agonism.

NURA: A curated dataset of nuclear receptor modulators

Nuclear receptors (NRs) are key regulators of human health and constitute a relevant target for medicinal chemistry applications as well as for toxicological risk assessment. Several open databases dedicated to small molecules that modulate NRs exist; however, depending on their final aim (i.e., adverse effect assessment or drug design), these databases contain a different amount and type of annotated molecules, along with a different distribution of experimental bioactivity values. Stemming from these considerations, in this work we aim to provide a unified dataset, NURA (NUclear Receptor Activity) dataset, collecting curated information on small molecules that modulate NRs, to be intended for both pharmacological and toxicological applications. NURA contains bioactivity annotations for 15,247 molecules and 11 selected NRs, and it was obtained by integrating and curating data from toxicological and pharmacological databases (i.e., Tox21, ChEMBL, NR-DBIND and BindingDB). Our results show that NURA dataset is a useful tool to bridge the gap between toxicology- and medicinal-chemistry-related databases, as it is enriched in terms of number of molecules, structural diversity and covered atomic scaffolds compared to the single sources. To the best of our knowledge, NURA dataset is the most exhaustive collection of small molecules annotated for their modulation of the chosen nuclear receptors. NURA dataset is intended to support decision-making in pharmacology and toxicology, as well as to contribute to data-driven applications, such as machine learning. The dataset and the data curation pipeline can be downloaded free of charge on Zenodo at the following DOI: https://doi.org/10.5281/zenodo.3991561.

Structures of Endocrine-Disrupting Chemicals Determine Binding to and Activation of the Estrogen Receptor α and Androgen Receptor

Endocrine-disrupting chemicals (EDCs) can interact with nuclear receptors, including estrogen receptor α (ERα) and androgen receptor (AR), to affect the normal endocrine system function, causing severe symptoms. Limited studies queried the EDC mechanisms, focusing on limited chemicals or a set of structurally similar compounds. It remained uncertain how hundreds of diverse EDCs could bind to ERα and AR and cause distinct functional consequences. Here, we employed a series of computational methodologies to investigate the structural features of EDCs that bind to and activate ERα and AR based on more than 4000 compounds. We used molecular docking and molecular dynamics simulations to elucidate the functional consequences and validated structure-function correlations experimentally using a time-resolved fluorescence resonance energy-transfer assay. We found that EDCs share three levels of key fragments. Primary (20 for ERα and 18 for AR) and secondary fragments (38 for ERα and 29 for AR) are responsible for the binding to receptors, and tertiary fragments determine the activity type (agonist, antagonist, or mixed). In summary, our study provides a general mechanism for the EDC function. Discovering the three levels of key fragments may drive fast screening and evaluation of potential EDCs from large sets of commercially used synthetic compounds.

An integrated chemical environment with tools for chemical safety testing

Moving toward species-relevant chemical safety assessments and away from animal testing requires access to reliable data to develop and build confidence in new approaches. The Integrated Chemical Environment (ICE) provides tools and curated data centered around chemical safety assessment. This article describes updates to ICE, including improved accessibility and interpretability of in vitro data via mechanistic target mapping and enhanced interactive tools for in vitro to in vivo extrapolation (IVIVE). Mapping of in vitro assay targets to toxicity endpoints of regulatory importance uses literature-based mode-of-action information and controlled terminology from existing knowledge organization systems to support data interoperability with external resources. The most recent ICE update includes Tox21 high-throughput screening data curated using analytical chemistry data and assay-specific parameters to eliminate potential artifacts or unreliable activity. Also included are physicochemical/ADME parameters for over 800,000 chemicals predicted by quantitative structure-activity relationship models. These parameters are used by the new ICE IVIVE tool in combination with the U.S. Environmental Protection Agency's httk R package to estimate in vivo exposures corresponding to in vitro bioactivity concentrations from stored or user-defined assay data. These new ICE features allow users to explore the applications of an expanded data space and facilitate building confidence in non-animal approaches.

Review and analysis of the potential for glyphosate to interact with the estrogen, androgen and thyroid pathways

Glyphosate was recently evaluated for its potential to interact with the estrogen, androgen and thyroid (EAT) hormone pathways, including steroidogenesis, under the United States Environmental Protection Agency's (USEPA) Endocrine Disruptor Screening Program (EDSP), then by Germany, the rapporteur Member State who led the European Annex 1 renewal for glyphosate, and then by the European Food Protection Agency (EFSA) also as part of the Annex 1 renewal for glyphosate. Under the EDSP, 11 Tier 1 assays were run following the USEPA's validated 890-series test guidelines and included five in vitro and six in vivo assays to evaluate the EAT pathways. Steroidogenesis was evaluated as part of the estrogen and androgen pathways. An up-to-date critical review has been conducted that considered results from the EDSP Tier 1 battery, guideline regulatory studies and an in-depth analysis of the literature studies that informed an endocrine assessment. A strength of this evaluation was that it included data across multiple levels of biological organization, and mammalian and nonmammalian test systems. There was strong agreement across the in vitro and in vivo Tier 1 battery, guideline studies and relevant literature studies, demonstrating that glyphosate does not interact with EAT pathways including steroidogenesis. Based on an analysis of the comprehensive toxicology database for glyphosate and the literature, this review has concluded that glyphosate does not have endocrine-disrupting properties through estrogen, androgen, thyroid and steroidogenic modes of action. © 2020 Society of Chemical Industry.

Tautomer Standardization in Chemical Databases: Deriving Business Rules from Quantum Chemistry

Databases of small, potentially bioactive molecules are ubiquitous across the industry and academia. Designed such that each unique compound should appear only once, the multiplicity of ways in which many compounds can be represented means that these databases require methods for standardizing the representation of chemistry. This is commonly achieved through the use of "Chemistry Business Rules", sets of predefined rules that describe the "house style" of the database in question. At Syngenta, the historical approach to the design of chemistry business rules has been to focus on consistency of representation, with chemical relevance given secondary consideration. In this work, we overturn that convention. Through the use of quantum chemistry calculations, we define a set of chemistry business rules for tautomer standardization that reproduces gas-phase energetic preferences. We go on to show that, compared to our historic approach, this method yields tautomers that are in better agreement with those observed experimentally in condensed phases and that are better suited for use in predictive models.

Toward a Computational Ecotoxicity Assay

Thousands of anthropogenic chemicals are released into the environment each year, posing potential hazards to human and environmental health. Toxic chemicals may cause a variety of adverse health effects, triggering immediate symptoms or delayed effects over longer periods of time. It is thus crucial to develop methods that can rapidly screen and predict the toxicity of chemicals to limit the potential harmful impacts of chemical pollutants. Computational methods are being increasingly used in toxicity predictions. Here, the method of molecular docking is assessed for screening potential toxicity of a variety of xenobiotic compounds, including pesticides, pharmaceuticals, pollutants, and toxins derived from the chemical industry. The method predicts the binding energy of pollutants to a set of carefully selected receptors under the assumption that toxicity in many cases is related to interference with biochemical pathways. The strength of the applied method lies in its rapid generation of interaction maps between potential toxins and the targeted enzymes, which could quickly yield molecular-level information and insight into potential perturbation pathways, aiding in the prioritization of chemicals for further tests. Two scoring functions are compared: Autodock Vina and the machine-learning scoring function RF-Score-VS. The results are promising, although hampered by the accuracy of the scoring functions. The strengths and weaknesses of the docking protocol are discussed, as well as future directions for improving the accuracy for the purpose of toxicity predictions.

Selecting a minimal set of androgen receptor assays for screening chemicals

Screening certain environmental chemicals for their ability to interact with endocrine targets, including the androgen receptor (AR), is an important global concern. We previously developed a model using a battery of eleven in vitro AR assays to predict in vivo AR activity. Here we describe a revised mathematical modeling approach that also incorporates data from newly available assays and demonstrate that subsets of assays can provide close to the same level of predictivity. These subset models are evaluated against the full model using 1820 chemicals, as well as in vitro and in vivo reference chemicals from the literature. Agonist batteries of as few as six assays and antagonist batteries of as few as five assays can yield balanced accuracies of 95% or better relative to the full model. Balanced accuracy for predicting reference chemicals is 100%. An approach is outlined for researchers to develop their own subset batteries to accurately detect AR activity using assays that map to the pathway of key molecular and cellular events involved in chemical-mediated AR activation and transcriptional activity. This work indicates in vitro bioactivity and in silico predictions that map to the AR pathway could be used in an integrated approach to testing and assessment for identifying chemicals that interact directly with the mammalian AR.

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.

Single-Cell Distribution Analysis of AR Levels by High-Throughput Microscopy in Cell Models: Application for Testing Endocrine-Disrupting Chemicals

Cell-to-cell variation of protein expression in genetically homogeneous populations is a common biological trait often neglected during analysis of high-throughput (HT) screens and is rarely used as a metric to characterize chemicals. We have captured single-cell distributions of androgen receptor (AR) nuclear levels after perturbations as a means to evaluate assay reproducibility and characterize a small subset of chemicals. AR, a member of the nuclear receptor family of transcription factors, is the central regulator of male reproduction and is involved in many pathophysiological processes. AR protein levels and nuclear localization often increase following ligand binding, with dihydrotestosterone (DHT) being the natural agonist. HT AR immunofluorescence imaging was used in multiple cell lines to define single-cell nuclear values extracted from thousands of cells per condition treated with DHT or DMSO (control). Analysis of numerous biological replicates led to a quality control metric that takes into account the distribution of single-cell data, and how it changes upon treatments. Dose-response experiments across several cell lines showed a large range of sensitivity to DHT, prompting us to treat selected cell lines with 45 Environmental Protection Agency (EPA)-provided chemicals that include many endocrine-disrupting chemicals (EDCs); data from six of the compounds were then integrated with orthogonal assays. Our comprehensive results indicate that quantitative single-cell distribution analysis of AR protein levels is a valid method to detect the potential androgenic and antiandrogenic actions of environmentally relevant chemicals in a sensitive and reproducible manner.

Reaction Library to Predict Direct Photochemical Transformation Products of Environmental Organic Contaminants in Sunlit Aquatic Systems

Cheminformatics-based applications to predict transformation pathways of environmental contaminants are useful to quickly prioritize contaminants with potentially toxic/persistent products. Direct photolysis can be an important degradation pathway for sunlight-absorbing compounds in aquatic systems. In this study, we developed the first freely available direct phototransformation pathway predictive tool, which uses a rule-based reaction library. Journal publications studying diverse contaminants (such as pesticides, pharmaceuticals, and energetic compounds) were systematically compiled to encode 155 reaction schemes into the reaction library. The execution result of this predictive tool was internally evaluated against 390 compounds from the compiled journal publications and externally evaluated against 138 compounds from the regulatory reports. The recall (sensitivity) and precision (selectivity) were 0.62 and 0.35, respectively, for internal evaluation, and 0.56 and 0.20, respectively, for external evaluation, when only the products formed from the first reaction step were counted. This predictive tool could help to narrow the data gaps in chemical registration/evaluation and inform future experimental studies.

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays

Molecular initiating events (MIEs) are key events in adverse outcome pathways that link molecular chemistry to target biology. As they are based on chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modeling. In this work, we aim to link ligand chemical structures to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Computer-Aided Ligand Discovery for Estrogen Receptor Alpha

Breast cancer (BCa) is one of the most predominantly diagnosed cancers in women. Notably, 70% of BCa diagnoses are Estrogen Receptor α positive (ERα+) making it a critical therapeutic target. With that, the two subtypes of ER, ERα and ERβ, have contrasting effects on BCa cells. While ERα promotes cancerous activities, ERβ isoform exhibits inhibitory effects on the same. ER-directed small molecule drug discovery for BCa has provided the FDA approved drugs tamoxifen, toremifene, raloxifene and fulvestrant that all bind to the estrogen binding site of the receptor. These ER-directed inhibitors are non-selective in nature and may eventually induce resistance in BCa cells as well as increase the risk of endometrial cancer development. Thus, there is an urgent need to develop novel drugs with alternative ERα targeting mechanisms that can overcome the limitations of conventional anti-ERα therapies. Several functional sites on ERα, such as Activation Function-2 (AF2), DNA binding domain (DBD), and F-domain, have been recently considered as potential targets in the context of drug research and discovery. In this review, we summarize methods of computer-aided drug design (CADD) that have been employed to analyze and explore potential targetable sites on ERα, discuss recent advancement of ERα inhibitor development, and highlight the potential opportunities and challenges of future ERα-directed drug discovery.

QSAR without borders

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as quantitative structure-activity relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling but it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas outside of traditional QSAR boundaries including synthesis planning, nanotechnology, materials science, biomaterials, and clinical informatics. As modern research methods generate rapidly increasing amounts of data, the knowledge of robust data-driven modelling methods professed within the QSAR field can become essential for scientists working both within and outside of chemical research. We hope that this contribution highlighting the generalizable components of QSAR modeling will serve to address this challenge.

OECD approaches and considerations for regulatory evaluation of endocrine disruptors

Identifying the potential endocrine disruptor hazard of environmental chemicals is a regulatory mandate for many countries. However, due to the adaptive nature of the endocrine system, absence of a single method capable of identifying endocrine disruption, and the latency between exposure to endocrine disrupting chemical during sensitive life stages and the manifestation of adverse responses, satisfying the regulatory requirement needed to identify a chemical as an endocrine disruptor is a challenge. There are now a variety of validated regulatory tests that can be used in combination to provide evidence that a chemical affects the oestrogen, androgen, thyroid, and steroidogenic pathways of vertebrates, but most rely (at least to some extent) on animal testing and require considerable cost and time to produce the necessary data. Emerging research methods are able to evaluate other endocrine pathways, incorporate more sensitive endpoints, and combine multiple alternative methods to predict in vivo outcomes. Some research approaches may also bridge gaps that have been identified in current endocrine regulatory testing. For the near term, considering new endpoints in a regulatory context may require adding them to existing test methods in order to establish relationships between the traditional and the innovative. From the outset, endocrine testing has always required integration of multiple methods that provide data on different levels of biological organisation, thus, the area of endocrine disruption is particularly adaptable to adverse outcome pathway (AOP) frameworks and integrated test methods built around AOPs. Herein, we provide a review of the status of endocrine disruptors in the OECD context, examples where innovation from research is needed to improve or bridge gaps in endocrine testing, and suggestions for regulators and researchers to facilitate uptake of innovate methods for endocrine disruptor regulatory testing. The increase in several human complex human disorders that include an endocrine component and the alarming decrease in wildlife biodiversity are commanding directives to include the best, most informative, innovative approaches to accelerate the rate and throughput of chemical evaluation for endocrine disruption.

Consensus versus Individual QSARs in Classification: Comparison on a Large-Scale Case Study

Consensus strategies have been widely applied in many different scientific fields, based on the assumption that the fusion of several sources of information increases the outcome reliability. Despite the widespread application of consensus approaches, their advantages in quantitative structure–activity relationship (QSAR) modeling have not been thoroughly evaluated, mainly due to the lack of appropriate large-scale data sets. In this study, we evaluated the advantages and drawbacks of consensus approaches compared to single classification QSAR models. To this end, we used a data set of three properties (androgen receptor binding, agonism, and antagonism) for approximately 4000 molecules with predictions performed by more than 20 QSAR models, made available in a large-scale collaborative project. The individual QSAR models were compared with two consensus approaches, majority voting and the Bayes consensus with discrete probability distributions, in both protective and nonprotective forms. Consensus strategies proved to be more accurate and to better cover the analyzed chemical space than individual QSARs on average, thus motivating their widespread application for property prediction. Scripts and data to reproduce the results of this study are available for download.

Binding interactions of halo-benzoic acids, halo-benzenesulfonic acids and halo-phenylboronic acids with human transthyretin

The anionic form-dependent binding interaction of halo-phenolic substances with human transthyretin (hTTR) has been observed previously. This indicates that ionizable compounds should be the primary focus in screening potential hTTR disruptors. Here, the potential binding potency of halo-benzoic acids, halo-benzenesulfonic acids/sulfates and halo-phenylboronic acids with hTTR was determined and analyzed by competitive fluorescence displacement assay integrated with computational methods. The laboratorial results indicated that the three test groups of model compounds exhibited a distinct binding affinity to hTTR. All the tested halo-phenylboronic acids, some of the tested halo-benzoic acids and halo-benzenesulfonic acids/sulfates were shown to be inactive with hTTR. Other halo-benzoic acids and halo-benzenesulfonic acids/sulfates were moderate and/or weak hTTR binders. The binding affinity of halo-benzoic acids and halo-benzenesulfonic acids/sulfates with hTTR was similar. The low distribution ability of the model compounds from water to hTTR may be the reason why they exhibited the binding potency observed with hTTR. By introducing other highly hydrophobic compounds, we observed that the binding affinity between compounds and hTTR increased with increasing molecular hydrophobicity. Those results indicated that the highly hydrophobic halo-benzoic acids and halo-benzenesulfonic acids/sulfates may be high-priority hTTR disruptors. Finally, a binary classification model was constructed employing three predictive variables. The sensitivity (S_n), specificity (S_p), predictive accuracy (Q) values of the training set and validation set were >0.83, indicating that the model had good classification performance. Thus, the binary classification model developed here could be used to distinguish whether a given ionizable compound is a potential hTTR binder or not.

Value and limitation of structure-based profilers to characterize developmental and reproductive toxicity potential

The uncertainty regarding the safety of chemicals leaching from food packaging triggers attention. In silico models provide solutions for screening of these chemicals, since many are toxicologically uncharacterized. For hazard assessment, information on developmental and reproductive toxicity (DART) is needed. The possibility to apply in silico toxicology to identify and quantify DART alerts was investigated. Open-source models and profilers were applied to 195 packaging chemicals and analogues. An approach based on DART and estrogen receptor (ER) binding profilers and molecular docking was able to identify all except for one chemical with documented DART properties. Twenty percent of the chemicals in the database known to be negative in experimental studies were classified as positive. The scheme was then applied to 121 untested chemicals. Alerts were identified for sixteen of them, five being packaging substances, the others structural analogues. Read-across was then developed to translate alerts into quantitative toxicological values. They can be used to calculate margins of exposure (MoE), the size of which reflects safety concern. The application of this approach appears valuable for hazard characterization of toxicologically untested packaging migrants. It is an alternative to the use of default uncertainty factor (UF) applied to animal chronic toxicity value to handle absence of DART data in hazard characterization.

A deep neural network combined with molecular fingerprints (DNN-MF) to develop predictive models for hydroxyl radical rate constants of water contaminants

This work combined a Deep Neural Network (DNN) with molecular fingerprints (MF) to develop models to predict the OH radical rate constants of 593 organic contaminants. Molecular descriptors, most often used in establishing quantitative structural-activity relationships (QSARs), were not used here because of their complicated generation processes that rely on advanced physicochemical and computational knowledge. Instead, we only fed the most basic information of the contaminant structures, i.e., MF encoding the types of atoms and how they are connected, to DNN and DNN then developed predictive models automatically. Here, a dataset containing 457 contaminants and their OH rate constants was first used to develop predictive models by DNN-MF. The hence developed models showed comparable accuracy to the traditional QSARs. The root mean square error (RMSE) values of the test sets were 0.358-0.384. The length of 2048 bits for the MF and 3 hidden layers (each with 1024 neurons) were found to be the optimal parameters for DNN. The model containing additional 89 micorpollutants in the training set was then successfully applied to predict the OH rate constants of 17 organophosphorus flame retardants and 29 additional micropollutants, with comparable accuracy to the reported molecular descriptors-based QSARs.

CoMPARA: Collaborative Modeling Project for Androgen Receptor Activity

BACKGROUND

Endocrine disrupting chemicals (EDCs) are xenobiotics that mimic the interaction of natural hormones and alter synthesis, transport, or metabolic pathways. The prospect of EDCs causing adverse health effects in humans and wildlife has led to the development of scientific and regulatory approaches for evaluating bioactivity. This need is being addressed using high-throughput screening (HTS) in vitro approaches and computational modeling.

OBJECTIVES

In support of the Endocrine Disruptor Screening Program, the U.S. Environmental Protection Agency (EPA) led two worldwide consortiums to virtually screen chemicals for their potential estrogenic and androgenic activities. Here, we describe the Collaborative Modeling Project for Androgen Receptor Activity (CoMPARA) efforts, which follows the steps of the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP).

METHODS

The CoMPARA list of screened chemicals built on CERAPP's list of 32,464 chemicals to include additional chemicals of interest, as well as simulated ToxCast™ metabolites, totaling 55,450 chemical structures. Computational toxicology scientists from 25 international groups contributed 91 predictive models for binding, agonist, and antagonist activity predictions. Models were underpinned by a common training set of 1,746 chemicals compiled from a combined data set of 11 ToxCast™/Tox21 HTS in vitro assays.

RESULTS

The resulting models were evaluated using curated literature data extracted from different sources. To overcome the limitations of single-model approaches, CoMPARA predictions were combined into consensus models that provided averaged predictive accuracy of approximately 80% for the evaluation set.

DISCUSSION

The strengths and limitations of the consensus predictions were discussed with example chemicals; then, the models were implemented into the free and open-source OPERA application to enable screening of new chemicals with a defined applicability domain and accuracy assessment. This implementation was used to screen the entire EPA DSSTox database of ∼875,000 chemicals, and their predicted AR activities have been made available on the EPA CompTox Chemicals dashboard and National Toxicology Program's Integrated Chemical Environment. https://doi.org/10.1289/EHP5580.

CoMPARA: Collaborative Modeling Project for Androgen Receptor Activity

BACKGROUND

Endocrine disrupting chemicals (EDCs) are xenobiotics that mimic the interaction of natural hormones and alter synthesis, transport, or metabolic pathways. The prospect of EDCs causing adverse health effects in humans and wildlife has led to the development of scientific and regulatory approaches for evaluating bioactivity. This need is being addressed using high-throughput screening (HTS) in vitro approaches and computational modeling.

OBJECTIVES

In support of the Endocrine Disruptor Screening Program, the U.S. Environmental Protection Agency (EPA) led two worldwide consortiums to virtually screen chemicals for their potential estrogenic and androgenic activities. Here, we describe the Collaborative Modeling Project for Androgen Receptor Activity (CoMPARA) efforts, which follows the steps of the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP).

METHODS

The CoMPARA list of screened chemicals built on CERAPP's list of 32,464 chemicals to include additional chemicals of interest, as well as simulated ToxCast™ metabolites, totaling 55,450 chemical structures. Computational toxicology scientists from 25 international groups contributed 91 predictive models for binding, agonist, and antagonist activity predictions. Models were underpinned by a common training set of 1,746 chemicals compiled from a combined data set of 11 ToxCast™/Tox21 HTS in vitro assays.

RESULTS

The resulting models were evaluated using curated literature data extracted from different sources. To overcome the limitations of single-model approaches, CoMPARA predictions were combined into consensus models that provided averaged predictive accuracy of approximately 80% for the evaluation set.

DISCUSSION

The strengths and limitations of the consensus predictions were discussed with example chemicals; then, the models were implemented into the free and open-source OPERA application to enable screening of new chemicals with a defined applicability domain and accuracy assessment. This implementation was used to screen the entire EPA DSSTox database of ∼ 875,000 chemicals, and their predicted AR activities have been made available on the EPA CompTox Chemicals dashboard and National Toxicology Program's Integrated Chemical Environment. https://doi.org/10.1289/EHP5580.

The role of fit-for-purpose assays within tiered testing approaches: A case study evaluating prioritized estrogen-active compounds in an in vitro human uterotrophic assay

Chemical risk assessment relies on toxicity tests that require significant numbers of animals, time and costs. For the >30,000 chemicals in commerce, the current scale of animal testing is insufficient to address chemical safety concerns as regulatory and product stewardship considerations evolve to require more comprehensive understanding of potential biological effects, conditions of use, and associated exposures. We demonstrate the use of a multi-level new approach methodology (NAMs) strategy for hazard- and risk-based prioritization to reduce animal testing. A Level 1/2 chemical prioritization based on estrogen receptor (ER) activity and metabolic activation using ToxCast data was used to select 112 chemicals for testing in a Level 3 human uterine cell estrogen response assay (IKA assay). The Level 3 data were coupled with quantitative in vitro to in vivo extrapolation (Q-IVIVE) to support bioactivity determination (as a surrogate for hazard) in a tissue-specific context. Assay AC₅₀s and Q-IVIVE were used to estimate human equivalent doses (HEDs), and HEDs were compared to rodent uterotrophic assay in vivo-derived points of departure (PODs). For substances active both in vitro and in vivo, IKA assay-derived HEDs were lower or equivalent to in vivo PODs for 19/23 compounds (83%). Activity exposure relationships were calculated, and the IKA assay was as or more protective of human health than the rodent uterotrophic assay for all IKA-positive compounds. This study demonstrates the utility of biologically relevant fit-for-purpose assays and supports the use of a multi-level strategy for chemical risk assessment.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

Ion Channels as Therapeutic Targets for Type 1 Diabetes Mellitus

Ion channels are integral proteins expressed in almost all living cells and are involved in muscle contraction and nutrient transport. They play a critical role in the normal functioning of the excitable tissues of the nervous system and regulate the action potential and contraction events. Dysfunction of genes encodes ion channel proteins, which disrupt the channel function and lead to a number of diseases, among which is type 1 diabetes mellitus (T1DM). Therefore, understanding the complex mechanism of ion channel receptors is necessary to facilitate the diagnosis and management of treatment. In this review, we summarize the mechanism of important ion channels and their potential role in the regulation of insulin secretion along with the limitations of ion channels as therapeutic targets. Furthermore, we discuss the recent investigations of the mechanism regulating the ion channels in pancreatic beta cells, which suggest that ion channels are active participants in the regulation of insulin secretion.

Towards accurate high-throughput ligand affinity prediction by exploiting structural ensembles, docking metrics and ligand similarity

MOTIVATION

Nowadays, virtual screening (VS) plays a major role in the process of drug development. Nonetheless, an accurate estimation of binding affinities, which is crucial at all stages, is not trivial and may require target-specific fine-tuning. Furthermore, drug design also requires improved predictions for putative secondary targets among which is Estrogen Receptor alpha (ERα).

RESULTS

VS based on combinations of Structure-Based VS (SBVS) and Ligand-Based VS (LBVS) is gaining momentum to improve VS performances. In this study, we propose an integrated approach using ligand docking on multiple structural ensembles to reflect receptor flexibility. Then, we investigate the impact of the two different types of features (structure-based and ligand molecular descriptors) on affinity predictions using a random forest algorithm. We find that ligand-based features have lower predictive power (rP = 0.69, R2 = 0.47) than structure-based features (rP = 0.78, R2 = 0.60). Their combination maintains high accuracy (rP = 0.73, R2 = 0.50) on the internal test set, but it shows superior robustness on external datasets. Further improvement and extending the training dataset to include xenobiotics, leads to a novel high-throughput affinity prediction method for ERα ligands (rP = 0.85, R2 = 0.71). The presented prediction tool is provided to the community as a dedicated satellite of the @TOME server in which one can upload a ligand dataset in mol2 format and get ligand docked and affinity predicted.

AVAILABILITY AND IMPLEMENTATION

http://edmon.cbs.cnrs.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ensemble QSAR Modeling to Predict Multispecies Fish Toxicity Lethal Concentrations and Points of Departure

QSAR modeling can be used to aid testing prioritization of the thousands of chemical substances for which no ecological toxicity data are available. We drew on the U.S. Environmental Protection Agency's ECOTOX database with additional data from ECHA to build a large data set containing in vivo test data on fish for thousands of chemical substances. This was used to create QSAR models to predict two types of end points: acute LC50 (median lethal concentration) and points of departure similar to the NOEC (no observed effect concentration) for any duration (named the "LC50" and "NOEC" models, respectively). These models used study covariates, such as species and exposure route, as features to facilitate the simultaneous use of varied data types. A novel method of substituting taxonomy groups for species dummy variables was introduced to maximize generalizability to different species. A stacked ensemble of three machine learning methods-random forest, gradient boosted trees, and support vector regression-was implemented to best make use of a large data set with many descriptors. The LC50 and NOEC models predicted end points within 1 order of magnitude 81% and 76% of the time, respectively, and had RMSEs of roughly 0.83 and 0.98 log10(mg/L), respectively. Benchmarks against the existing TEST and ECOSAR tools suggest improved prediction accuracy.

Toxicity profiling and prioritization of plant-derived antimalarial agents

Human malaria is the most widespread mosquito-borne life-threatening disease worldwide. In the absence of effective vaccines, prevention and treatment of malaria only depend on prophylaxis and drug-based therapy either in monotherapy or in combination. Unfortunately, the number of available antimalarial drugs presenting different mechanisms of action is rather limited. In addition, the appearance of drug-resistance in the parasite strains impacts the efficacy of the treatments. As a result, there is a crucial need to find new drugs to circumvent resistance problems. In the quest to identify new antimalarial agents a huge number of plant-derived compounds (PDCs) have been investigated. Surprisingly in the in silico PDC screening programs, toxicity filters are either never used or so simple that their interest is limited. In this context, the goal of this study was to show how to take advantage of validated toxicity QSAR models for refining the selection of PDCs. From an original data set of 507 PDCs collected from the literature, the use of toxicity filters for endocrine disruption, developmental toxicity, and hepatotoxicity in conjunction with classical pharmacokinetic filters allowed us to obtain a list of 31 compounds of potential interest. The pros and cons of such a strategy have been discussed.

In Silico Predictions of Endocrine Disruptors Properties

Endocrine-disrupting chemicals (EDCs) are a broad class of molecules present in our environment that are suspected to cause adverse effects in the endocrine system by interfering with the synthesis, transport, degradation, or action of endogenous ligands. The characterization of the harmful interaction between environmental compounds and their potential cellular targets and the development of robust in vivo, in vitro, and in silico screening methods are important for assessment of the toxic potential of large numbers of chemicals. In this context, computer-aided technologies that will allow for activity prediction of endocrine disruptors and environmental risk assessments are being developed. These technologies must be able to cope with diverse data and connect chemistry at the atomic level with the biological activity at the cellular, organ, and organism levels. Quantitative structure-activity relationship methods became popular for toxicity issues. They correlate the chemical structure of compounds with biological activity through a number of molecular descriptors (e.g., molecular weight and parameters to account for hydrophobicity, topology, or electronic properties). Chemical structure analysis is a first step; however, modeling intermolecular interactions and cellular behavior will also be essential. The increasing number of three-dimensional crystal structures of EDCs' targets has provided a wealth of structural information that can be used to predict their interactions with EDCs using docking and scoring procedures. In the present review, we have described the various computer-assisted approaches that use ligands and targets properties to predict endocrine disruptor activities.

Novel Consensus Architecture To Improve Performance of Large-Scale Multitask Deep Learning QSAR Models

Advances in the development of high-throughput screening and automated chemistry have rapidly accelerated the production of chemical and biological data, much of them freely accessible through literature aggregator services such as ChEMBL and PubChem. Here, we explore how to use this comprehensive mapping of chemical biology space to support the development of large-scale quantitative structure-activity relationship (QSAR) models. We propose a new deep learning consensus architecture (DLCA) that combines consensus and multitask deep learning approaches together to generate large-scale QSAR models. This method improves knowledge transfer across different target/assays while also integrating contributions from models based on different descriptors. The proposed approach was validated and compared with proteochemometrics, multitask deep learning, and Random Forest methods paired with various descriptors types. DLCA models demonstrated improved prediction accuracy for both regression and classification tasks. The best models together with their modeling sets are provided through publicly available web services at https://predictor.ncats.io .

Integrating in silico models and read-across methods for predicting toxicity of chemicals: A step-wise strategy

In silico methods and models are increasingly used for predicting properties of chemicals for hazard identification and hazard characterisation in the absence of experimental toxicity data. Many in silico models are available and can be used individually or in an integrated fashion. Whilst such models offer major benefits to toxicologists, risk assessors and the global scientific community, the lack of a consistent framework for the integration of in silico results can lead to uncertainty and even contradictions across models and users, even for the same chemicals. In this context, a range of methods for integrating in silico results have been proposed on a statistical or case-specific basis. Read-across constitutes another strategy for deriving reference points or points of departure for hazard characterisation of untested chemicals, from the available experimental data for structurally-similar compounds, mostly using expert judgment. Recently a number of software systems have been developed to support experts in this task providing a formalised and structured procedure. Such a procedure could also facilitate further integration of the results generated from in silico models and read-across. This article discusses a framework on weight of evidence published by EFSA to identify the stepwise approach for systematic integration of results or values obtained from these "non-testing methods". Key criteria and best practices for selecting and evaluating individual in silico models are also described, together with the means to combining the results, taking into account any limitations, and identifying strategies that are likely to provide consistent results.

An Approach for Using In Vitro and In Silico Data to Identify Pharmaceuticals with Potential (Anti-)Estrogenic Activity in Aquatic Vertebrates at Environmentally Relevant Concentrations

Endocrine-active pharmaceuticals can cause adverse reproductive and developmental effects in nontarget organisms. Aquatic vertebrates may be susceptible to the effects of such pharmaceuticals given that the structure of hormone receptors and the physiology of the endocrine system are highly conserved across vertebrates. To aid in the regulatory review of the environmental impact of drugs, we demonstrate an approach to screen and support the prioritization of pharmaceuticals based on their ability to interact with estrogen receptors (ERs) at environmentally relevant concentrations. Tox21 in vitro results from ER agonist and antagonist assays were retrieved for 1123 pharmaceuticals. In silico predictions from the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP) models were used to estimate ER agonist and antagonist activity for an additional 170 pharmaceuticals not tested in the Tox21 assay platform. The estrogenic effect ratio (EER) and anti-estrogenic effect ratio (AEER) were calculated by comparing the activity concentration at half-maximal response (AC50) for ER agonism and antagonism, respectively, with estimated pharmaceutical concentrations in fish tissue based on estimates of environmental exposures. A total of 73 and 127 pharmaceuticals were identified as ER agonists and antagonists, respectively. As expected, 17β-estradiol and 17α-ethinylestradiol displayed EERs  > 1, and raloxifene and bazedoxifene acetate displayed AEERs > 1, thus indicating that these pharmaceuticals have the potential to reach fish tissue levels that exceed concentrations estimated to interact with ERs. Four pharmaceuticals displayed EERs between 0.1 and 1, and 6 displayed AEERs between 0.1 and 1. This approach may help determine the need for submission of environmental assessment data for new drug applications and support prioritization of pharmaceuticals with the potential to disrupt endocrine signaling in vertebrates. Environ Toxicol Chem 2019;38:2154-2168. © 2019 SETAC.

Determining the endocrine disruption potential of industrial chemicals using an integrative approach: Public databases, in vitro exposure, and modeling receptor interactions

Environmental and occupational exposure to industrial chemicals has been linked to toxic and carcinogenic effects in animal models and human studies. However, current toxicology testing does not thoroughly explore the endocrine disrupting effects of industrial chemicals, which may have low dose effects not predicted when determining the limit of toxicity. The objective of this study was to evaluate the endocrine disrupting potential of a broad range of chemicals used in the petrochemical sector. Therefore, 139 chemicals were classified for reproductive toxicity based on the United Nations Globally Harmonized System for hazard classification. These chemicals were evaluated in PubMed for reported endocrine disrupting activity, and their endocrine disrupting potential was estimated by identifying chemicals with active nuclear receptor endpoints publicly available databases. Evaluation of ToxCast data suggested that these chemicals preferentially alter the activity of the estrogen receptor (ER). Four chemicals were prioritized for in vitro testing using the ER-positive, immortalized human uterine Ishikawa cell line and a range of concentrations below the reported limit of toxicity in humans. We found that 2,6-di-tert-butyl-p-cresol (BHT) and diethanolamine (DEA) repressed the basal expression of estrogen-responsive genes PGR, NPPC, and GREB1 in Ishikawa cells, while tetrachloroethylene (PCE) and 2,2'-methyliminodiethanol (MDEA) induced the expression of these genes. Furthermore, low-dose combinations of PCE and MDEA produced additive effects. All four chemicals interfered with estradiol-mediated induction of PGR, NPPC, and GREB1. Molecular docking demonstrated that these chemicals could bind to the ligand binding site of ERα, suggesting the potential for direct stimulatory or inhibitory effects. We found that these chemicals altered rates of proliferation and regulated the expression of cell proliferation associated genes. These findings demonstrate previously unappreciated endocrine disrupting effects and underscore the importance of testing the endocrine disrupting potential of chemicals in the future to better understand their potential to impact public health.

Open-source QSAR models for pKa prediction using multiple machine learning approaches

BACKGROUND

The logarithmic acid dissociation constant pKa reflects the ionization of a chemical, which affects lipophilicity, solubility, protein binding, and ability to pass through the plasma membrane. Thus, pKa affects chemical absorption, distribution, metabolism, excretion, and toxicity properties. Multiple proprietary software packages exist for the prediction of pKa, but to the best of our knowledge no free and open-source programs exist for this purpose. Using a freely available data set and three machine learning approaches, we developed open-source models for pKa prediction.

METHODS

The experimental strongest acidic and strongest basic pKa values in water for 7912 chemicals were obtained from DataWarrior, a freely available software package. Chemical structures were curated and standardized for quantitative structure-activity relationship (QSAR) modeling using KNIME, and a subset comprising 79% of the initial set was used for modeling. To evaluate different approaches to modeling, several datasets were constructed based on different processing of chemical structures with acidic and/or basic pKas. Continuous molecular descriptors, binary fingerprints, and fragment counts were generated using PaDEL, and pKa prediction models were created using three machine learning methods, (1) support vector machines (SVM) combined with k-nearest neighbors (kNN), (2) extreme gradient boosting (XGB) and (3) deep neural networks (DNN).

RESULTS

The three methods delivered comparable performances on the training and test sets with a root-mean-squared error (RMSE) around 1.5 and a coefficient of determination (R2) around 0.80. Two commercial pKa predictors from ACD/Labs and ChemAxon were used to benchmark the three best models developed in this work, and performance of our models compared favorably to the commercial products.

CONCLUSIONS

This work provides multiple QSAR models to predict the strongest acidic and strongest basic pKas of chemicals, built using publicly available data, and provided as free and open-source software on GitHub.

Overview of the effects of chemical mixtures with endocrine disrupting activity in the context of real-life risk simulation: An integrative approach (Review)

Research over the past years has indicated that chronic human exposure to very low doses of various chemical species in mixtures and administered via different routes (percutaneous, orally, etc.) should be the main focus of new biochemical and toxicological studies. Humans have daily contact with various chemicals, such as food additives, pesticides from fruits/vegetables, antibiotics (and other veterinary drugs) from meat, different types of preservatives from cosmetics, to name a few. Simultaneous exposure to this wide array of chemicals does not produce immediate effects, but summative effect/s over time that may be clinically manifested several years thereafter. Classical animal studies designed to test the toxic outcome of a single chemical are not suitable to assess, and then extrapolate to humans, the effects of a whole mixture of chemicals. Testing the aftermath of a combination of chemicals, at low doses, around or below the no observed adverse effect is stressed by many toxicologists. Thus, there is a need to reformulate the design of biochemical and toxicological studies in order to perform real-life risk simulation. This review discuss the potential use of computational methods as a complementary tool for in vitro and in vivo toxicity tests with a high predictive potential that could contribute to reduce animal testing, cost and time, when assessing the effects of chemical combinations. This review focused on the use of these methods to predict the potential endocrine disrupting activity of a mixture of chemicals.

Predicting estrogen receptor binding of chemicals using a suite of in silico methods - Complementary approaches of (Q)SAR, molecular docking and molecular dynamics

With the aim of obtaining reliable estimates of Estrogen Receptor (ER) binding for diverse classes of compounds, a weight of evidence approach using estimates from a suite of in silico models was assessed. The predictivity of a simple Majority Consensus of (Q)SAR models was assessed using a test set of compounds with experimental Relative Binding Affinity (RBA) data. Molecular docking was also carried out and the binding energies of these compounds to the ERα receptor were determined. For a few selected compounds, including a known full agonist and antagonist, the intrinsic activity was determined using low-mode molecular dynamics methods. Individual (Q)SAR model predictivity varied, as expected, with some models showing high sensitivity, others higher specificity. However, the Majority Consensus (Q)SAR prediction showed a high accuracy and reasonably balanced sensitivity and specificity. Molecular docking provided quantitative information on strength of binding to the ERα receptor. For the 50 highest binding affinity compounds with positive RBA experimental values, just 5 of them were predicted to be non-binders by the Majority QSAR Consensus. Furthermore, agonist-specific assay experimental values for these 5 compounds were negative, which indicates that they may be ER antagonists. We also showed different scenarios of combining (Q)SAR results with Molecular docking classification of ER binding based on cut-off values of binding energies, providing a rational combined strategy to maximize terms of toxicological interest.

Combining in Silico Tools with Multicriteria Analysis for Alternatives Assessment of Hazardous Chemicals: A Case Study of Decabromodiphenyl Ether Alternatives

Alternatives assessment is applied for minimizing the risk of unintentionally replacing a hazardous chemical with another hazardous chemical. Central challenges are the diversity of properties to consider and the lack of high-quality experimental data. To address this, a novel alternatives assessment procedure was developed based on in silico data and multicriteria decision analysis (MCDA) methods. As a case study, 16 alternatives to the flame retardant decabromodiphenyl ether were considered. The hazard properties included persistence (P), bioaccumulation potential (B), toxicities (T), and mobility in water (M). Databases were consulted and 2866 experimental data points were collected for the target chemicals; however, these were mostly replicate data points for some hazard criteria for a subset of alternatives. Therefore, in silico data and three MCDA strategies were tested including heat mapping, multiattribute utility theory (MAUT), and Elimination Et Choix Traduisant la REalité (ELECTRE III). The heat map clearly showed that none of the target chemicals are hazard-free, whereas MAUT and ELECTRE III agreed on ranking the "least worst" choices. This study identified several challenges and the complexity in the alternatives assessment processes motivating more case studies combining in silico and MCDA approaches.

The Next Generation Blueprint of Computational Toxicology at the U.S. Environmental Protection Agency

The U.S. Environmental Protection Agency (EPA) is faced with the challenge of efficiently and credibly evaluating chemical safety often with limited or no available toxicity data. The expanding number of chemicals found in commerce and the environment, coupled with time and resource requirements for traditional toxicity testing and exposure characterization, continue to underscore the need for new approaches. In 2005, EPA charted a new course to address this challenge by embracing computational toxicology (CompTox) and investing in the technologies and capabilities to push the field forward. The return on this investment has been demonstrated through results and applications across a range of human and environmental health problems, as well as initial application to regulatory decision-making within programs such as the EPA's Endocrine Disruptor Screening Program. The CompTox initiative at EPA is more than a decade old. This manuscript presents a blueprint to guide the strategic and operational direction over the next 5 years. The primary goal is to obtain broader acceptance of the CompTox approaches for application to higher tier regulatory decisions, such as chemical assessments. To achieve this goal, the blueprint expands and refines the use of high-throughput and computational modeling approaches to transform the components in chemical risk assessment, while systematically addressing key challenges that have hindered progress. In addition, the blueprint outlines additional investments in cross-cutting efforts to characterize uncertainty and variability, develop software and information technology tools, provide outreach and training, and establish scientific confidence for application to different public health and environmental regulatory decisions.

Predicting estrogen receptor activation by a group of substituted phenols: An integrated approach to testing and assessment case study

Traditional approaches for chemical risk assessment cannot keep pace with the number of substances requiring assessment. Thus, in a global effort to expedite and modernize chemical risk assessment, New Approach Methodologies (NAMs) are being explored and developed. Included in this effort is the OECD Integrated Approaches for Testing and Assessment (IATA) program, which provides a forum for OECD member countries to develop and present case studies illustrating the application of NAM in various risk assessment contexts. Here, we present an IATA case study for the prediction of estrogenic potential of three target phenols: 4-tert-butylphenol, 2,4-di-tert-butylphenol and octabenzone. Key features of this IATA include the use of two computational approaches for analogue selection for read-across, data collected from traditional and NAM sources, and a workflow to generate predictions regarding the targets' ability to bind the estrogen receptor (ER). Endocrine disruption can occur when a chemical substance mimics the activity of natural estrogen by binding to the ER and, if potency and exposure are sufficient, alters the function of the endocrine system to cause adverse effects. The data indicated that of the three target substances that were considered herein, 4-tert-butylphenol is a potential endocrine disruptor. Further, this IATA illustrates that the NAM approach explored is health protective when compared to in vivo endpoints traditionally used for human health risk assessment.

Predictive models in ecotoxicology: Bridging the gap between scientific progress and regulatory applicability-Remarks and research needs

This paper concludes a special series of 7 articles (4 on toxicokinetic-toxicodynamic [TK-TD] models and 3 on quantitative structure-activity relationship [QSAR] models) published in previous issues of Integrated Environmental Assessment and Management (IEAM). The present paper summarizes the special series articles and highlights their contribution to the topic of increasing the regulatory applicability of effect models. For both TK-TD and QSAR approaches, we then describe the main research needs. The use of TK-TD models for describing sublethal effects must be better developed, particularly through the improvement of the dynamic energy budget (DEBtox) approach. The potential of TK-TD models for moving from lower (molecular) to higher (population) hierarchical levels is highlighted as a promising research line. Some relevant issues to improve the acceptance of QSAR models at the regulatory level are also described, such as increased transparency of the performance assessment and of the modeling algorithms, model documentation, relevance of the chosen target for regulatory needs, and improved mechanistic interpretability. Integr Environ Assess Manag 2019;00:000-000. © 2019 SETAC.

Homology models of mouse and rat estrogen receptor-α ligand-binding domain created by in silico mutagenesis of a human template: molecular docking with 17ß-estradiol, diethylstilbestrol, and paraben analogs

Crystal structures exist for human, but not rodent, estrogen receptor-α ligand-binding domain (ERα-LBD). Consequently, rodent studies involving binding of compounds to ERα-LBD are limited in their molecular-level interpretation and extrapolation to humans. Because the sequences of rodent and human ERα-LBDs are > 95% identical, we expected their 3D structures and ligand binding to be highly similar. To test this hypothesis, we used the human ERα-LBD structure (PDB 3UUD) as a template to produce rat and mouse homology models. Employing the rodent models and human structure, we generated docking poses of 23 Group A ligands (17ß-estradiol, diethylstilbestrol, and 21 paraben analogs) in AutoDock Vina for interspecies comparisons. Ligand RMSDs (Å) (median, 95% CI) were 0.49 (0.21-1.82) (human-mouse) and 1.19 (0.22-1.82) (human-rat), well below the 2.0-2.5 Å range for equivalent docking poses. Numbers of interspecies ligand-receptor residue contacts were highly similar, with Sorensen Sc (%) = 96.8 (90.0-100) (human-mouse) and 97.7 (89.5-100) (human-rat). Likewise, numbers of interspecies ligand-receptor residue contacts were highly correlated: Pearson r = 0.913 (human-mouse) and 0.925 (human-rat). Numbers of interspecies ligand-receptor atom contacts were even more tightly correlated: r = 0.979 (human-mouse) and 0.986 (human-rat). Pyramid plots of numbers of ligand-receptor atom contacts by residue exhibited high interspecies symmetry and had Spearman r s = 0.977 (human-mouse) and 0.966 (human-rat). Group B ligands included 15 ring-substituted parabens recently shown experimentally to exhibit decreased binding to human ERα and to exert increased antimicrobial activity. Ligand efficiencies calculated from docking ligands into human ERα-LBD were well correlated with those derived from published experimental data (Pearson partial r p = 0.894 and 0.918; Groups A and B, respectively). Overall, the results indicate that our constructed rodent ERα-LBDs interact with ligands in like manner to the human receptor, thus providing a high level of confidence in extrapolations of rodent to human ligand-receptor interactions.

Development, validation and integration of in silico models to identify androgen active chemicals

Humans are exposed to large numbers of environmental chemicals, some of which potentially interfere with the endocrine system. The identification of potential endocrine disrupting chemicals (EDCs) has gained increasing priority in the assessment of environmental hazards. The U.S. Environmental Protection Agency (U.S. EPA) has developed the Endocrine Disruptor Screening Program (EDSP) which aims to prioritize and screen potential EDCs. The Toxicity Forecaster (ToxCast) program has generated data using in vitro high-throughput screening (HTS) assays measuring activity of chemicals at multiple points along the androgen receptor (AR) activity pathway. In the present study, using a large and diverse data set of 1667 chemicals provided by the U.S. EPA from the combined ToxCast AR assays in the framework of the Collaborative Modeling Project for Androgen Receptor Activity (CoMPARA). Two models were built using ADMET Predictor™; one is based on Artificial Neural Networks (ANNs) technology and the other uses a Support Vector Machine (SVM) algorithm; one model is a Decision Tree (DT) developed in R; and two models make use of differently combined sets of structural alerts (SAs) automatically extracted by SARpy. We used two strategies to integrate predictions from single models; one is based on a majority vote approach and the other on prediction convergence. These strategies led to enhanced statistical performance in most cases. Moreover, the majority vote approach improved prediction coverage when one or more single models were not able to provide any estimations. This study integrates multiple in silico approaches as a virtual screening tool for use in risk assessment of endocrine disrupting chemicals.

In Silico Prediction of Endocrine Disrupting Chemicals Using Single-Label and Multilabel Models

Endocrine disruption (ED) has become a serious public health issue and also poses a significant threat to the ecosystem. Due to complex mechanisms of ED, traditional in silico models focusing on only one mechanism are insufficient for detection of endocrine disrupting chemicals (EDCs), let alone offering an overview of possible action mechanisms for a known EDC. To remove these limitations, in this study both single-label and multilabel models were constructed across six ED targets, namely, AR (androgen receptor), ER (estrogen receptor alpha), TR (thyroid receptor), GR (glucocorticoid receptor), PPARg (peroxisome proliferator-activated receptor gamma), and aromatase. Two machine learning methods were used to build the single-label models, with multiple random under-sampling combining voting classification to overcome the challenge of data imbalance. Four methods were explored to construct the multilabel models that can predict the interaction of one EDC against multiple targets simultaneously. The single-label models of all the six targets have achieved reasonable performance with balanced accuracy (BA) values from 0.742 to 0.816. Each top single-label model was then joined to predict the multilabel test set with BA values from 0.586 to 0.711. The multilabel models could offer a significant boost over the single-label baselines with BA values for the multilabel test set from 0.659 to 0.832. Therefore, we concluded that single-label models could be employed for identification of potential EDCs, while multilabel ones are preferable for prediction of possible mechanisms of known EDCs.