CERAPP: Collaborative Estrogen Receptor Activity Prediction Project

Predictive Models for Acute Oral Systemic Toxicity: A Workshop to Bridge the Gap from Research to Regulation

In early 2018, the Interagency Coordinating Committee for the Validation of Alternative Methods (ICCVAM) published the "Strategic Roadmap for Establishing New Approaches to Evaluate the Safety of Chemicals and Medical Products in the United States" (ICCVAM 2018). Cross-agency federal workgroups have been established to implement this roadmap for various toxicological testing endpoints, with an initial focus on acute toxicity testing. The ICCVAM acute toxicity workgroup (ATWG) helped organize a global collaboration to build predictive in silico models for acute oral systemic toxicity, based on a large dataset of rodent studies and targeted towards regulatory needs identified across federal agencies. Thirty-two international groups across government, industry, and academia participated in the project, culminating in a workshop in April 2018 held at the National Institutes of Health (NIH). At the workshop, computational modelers and regulatory decision makers met to discuss the feasibility of using predictive model outputs for regulatory use in lieu of acute oral systemic toxicity testing. The models were combined to yield consensus predictions which demonstrated excellent performance when compared to the animal data, and workshop outcomes and follow-up activities to make these tools available and put them into practice are discussed here.

A combination of ternary classification models and reporter gene assays for the comprehensive thyroid hormone disruption profiles of 209 polychlorinated biphenyls

Computational toxicology is widely applied to screen tens and thousands of potential environmental endocrine disruptors (EDCs) for its great efficiency and rapid evaluation in recent years. Polychlorinated biphenyls (PCBs) with 209 congeners have not been comprehensively tested for their ability to interact with the thyroid receptor (TR), which is one of the most extensively studied targets related to the effects of EDCs. In this study, we aimed to determine the thyroid-disrupting mechanisms of 209 PCBs through the combination of a novel computational ternary classification model and luciferase reporter gene assay. The reporter gene assay was performed to examine the hormone activities of 22 PCBs via TR to classify their profiles as agonistic, antagonistic or inactive. Thus, six agonists, eleven antagonists and seven inactive chemicals against TR were identified in in vitro assays. Then, six relevant variables, including molecular structural descriptors and molecular docking scores, were selected for describing compounds. Machine learning methods (i.e., linear discriminant analysis (LDA) and support vector machines (SVM)) were used to build classification models. The optimal model was obtained by using SVM, with an accuracy of 88.24% in the training set, 80.0% in the test set and 75.0% in 10-fold cross-validation. The remaining 187 PCB congeners' hormone activities were predicted using the obtained models. Out of these PCBs, the SVM model predicted 38 agonists and 81 antagonists. The findings revealed that higher chlorinated PCBs tended to have TR-antagonistic activities, whereas lower chlorinated PCBs were agonists. This study provided a reference for further exploring PCBs' thyroid effect.

Detecting the bioaccumulation patterns of chemicals through data-driven approaches

This work investigates the bioaccumulation patterns of 168 organic chemicals in fish, by comparing their bioconcentration factor (BCF), biomagnification factor (BMF) and octanol-water partitioning coefficient (K_OW). It aims to gain insights on the relationships between dietary and non-dietary bioaccumulation in aquatic environment, on the effectiveness of K_OW and BCF to detect compounds that bioaccumulate through diet, as well as to detect the presence of structure-related bioaccumulation patterns. A linear relationship between logBMF and logK_OW was observed (logBMF = 1.14·logBCF - 6.20) up to logK_OW ≈ 4, as well as between logBMF and logBCF (logBMF = 0.96·logBCF - 4.06) up to a logBCF ≈ 5. 10% of compounds do not satisfy the linear BCF-BMF relationship. The deviations from such linear relationships were further investigated with the aid of a self-organizing map and canonical correlation analysis, which allowed us to shed light on some structure-related patterns. Finally, the usage of K_OW- and BCF-based thresholds to detect compounds that accumulate through diet led to many false positives (47%-91% for K_OW), and a moderate number of false negatives (up to 5% for BCF). These results corroborate the need of using the experimental BMF for hazard assessment practices, as well as of developing computational tools for BMF prediction.

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

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

Exploring non-linear distance metrics in the structure-activity space: QSAR models for human estrogen receptor

Background

Quantitative structure-activity relationship (QSAR) models are important tools used in discovering new drug candidates and identifying potentially harmful environmental chemicals. These models often face two fundamental challenges: limited amount of available biological activity data and noise or uncertainty in the activity data themselves. To address these challenges, we introduce and explore a QSAR model based on custom distance metrics in the structure-activity space.

Methods

The model is built on top of the k-nearest neighbor model, incorporating non-linearity not only in the chemical structure space, but also in the biological activity space. The model is tuned and evaluated using activity data for human estrogen receptor from the US EPA ToxCast and Tox21 databases.

Results

The model closely trails the CERAPP consensus model (built on top of 48 individual human estrogen receptor activity models) in agonist activity predictions and consistently outperforms the CERAPP consensus model in antagonist activity predictions.

Discussion

We suggest that incorporating non-linear distance metrics may significantly improve QSAR model performance when the available biological activity data are limited.

Electronic supplementary material

The online version of this article (10.1186/s13321-018-0300-0) contains supplementary material, which is available to authorized users.

Rapid experimental measurements of physicochemical properties to inform models and testing

The structures and physicochemical properties of chemicals are important for determining their potential toxicological effects, toxicokinetics, and route(s) of exposure. These data are needed to prioritize the risk for thousands of environmental chemicals, but experimental values are often lacking. In an attempt to efficiently fill data gaps in physicochemical property information, we generated new data for 200 structurally diverse compounds, which were rigorously selected from the USEPA ToxCast chemical library, and whose structures are available within the Distributed Structure-Searchable Toxicity Database (DSSTox). This pilot study evaluated rapid experimental methods to determine five physicochemical properties, including the log of the octanol:water partition coefficient (known as log(Kow) or logP), vapor pressure, water solubility, Henry's law constant, and the acid dissociation constant (pKa). For most compounds, experiments were successful for at least one property; log(Kow) yielded the largest return (176 values). It was determined that 77 ToxPrint structural features were enriched in chemicals with at least one measurement failure, indicating which features may have played a role in rapid method failures. To gauge consistency with traditional measurement methods, the new measurements were compared with previous measurements (where available). Since quantitative structure-activity/property relationship (QSAR/QSPR) models are used to fill gaps in physicochemical property information, 5 suites of QSPRs were evaluated for their predictive ability and chemical coverage or applicability domain of new experimental measurements. The ability to have accurate measurements of these properties will facilitate better exposure predictions in two ways: 1) direct input of these experimental measurements into exposure models; and 2) construction of QSPRs with a wider applicability domain, as their predicted physicochemical values can be used to parameterize exposure models in the absence of experimental data.

Predictive Models in Ecotoxicology: Bridging the Gap Between Scientific Progress and Regulatory Applicability

This special series is the outcome of the session "Predictive models in ecotoxicology: Bridging the gap between scientific progress and regulatory applicability," held at the 27th SETAC Europe annual meeting (Brussels, May 2017). In this foreword the rationale behind the special series, the reasons for proposing it, and its objectives are described briefly. Integr Environ Assess Manag 2018;14:601-603. © 2018 SETAC.

Integrating tools for non-targeted analysis research and chemical safety evaluations at the US EPA

Tens-of-thousands of chemicals are registered in the U.S. for use in countless processes and products. Recent evidence suggests that many of these chemicals are measureable in environmental and/or biological systems, indicating the potential for widespread exposures. Traditional public health research tools, including in vivo studies and targeted analytical chemistry methods, have been unable to meet the needs of screening programs designed to evaluate chemical safety. As such, new tools have been developed to enable rapid assessment of potentially harmful chemical exposures and their attendant biological responses. One group of tools, known as "non-targeted analysis" (NTA) methods, allows the rapid characterization of thousands of never-before-studied compounds in a wide variety of environmental, residential, and biological media. This article discusses current applications of NTA methods, challenges to their effective use in chemical screening studies, and ways in which shared resources (e.g., chemical standards, databases, model predictions, and media measurements) can advance their use in risk-based chemical prioritization. A brief review is provided of resources and projects within EPA's Office of Research and Development (ORD) that provide benefit to, and receive benefits from, NTA research endeavors. A summary of EPA's Non-Targeted Analysis Collaborative Trial (ENTACT) is also given, which makes direct use of ORD resources to benefit the global NTA research community. Finally, a research framework is described that shows how NTA methods will bridge chemical prioritization efforts within ORD. This framework exists as a guide for institutions seeking to understand the complexity of chemical exposures, and the impact of these exposures on living systems.

"MS-Ready" structures for non-targeted high-resolution mass spectrometry screening studies

Chemical database searching has become a fixture in many non-targeted identification workflows based on high-resolution mass spectrometry (HRMS). However, the form of a chemical structure observed in HRMS does not always match the form stored in a database (e.g., the neutral form versus a salt; one component of a mixture rather than the mixture form used in a consumer product). Linking the form of a structure observed via HRMS to its related form(s) within a database will enable the return of all relevant variants of a structure, as well as the related metadata, in a single query. A Konstanz Information Miner (KNIME) workflow has been developed to produce structural representations observed using HRMS ("MS-Ready structures") and links them to those stored in a database. These MS-Ready structures, and associated mappings to the full chemical representations, are surfaced via the US EPA's Chemistry Dashboard ( https://comptox.epa.gov/dashboard/ ). This article describes the workflow for the generation and linking of ~ 700,000 MS-Ready structures (derived from ~ 760,000 original structures) as well as download, search and export capabilities to serve structure identification using HRMS. The importance of this form of structural representation for HRMS is demonstrated with several examples, including integration with the in silico fragmentation software application MetFrag. The structures, search, download and export functionality are all available through the CompTox Chemistry Dashboard, while the MetFrag implementation can be viewed at https://msbi.ipb-halle.de/MetFragBeta/ .

Novel application of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene sets associated with disease: use case in breast carcinogenesis

Advances in technology within biomedical sciences have led to an inundation of data across many fields, raising new challenges in how best to integrate and analyze these resources. For example, rapid chemical screening programs like the US Environmental Protection Agency's ToxCast and the collaborative effort, Tox21, have produced massive amounts of information on putative chemical mechanisms where assay targets are identified as genes; however, systematically linking these hypothesized mechanisms with in vivo toxicity endpoints like disease outcomes remains problematic. Herein we present a novel use of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene associations with biological concepts as represented by Medical Subject Headings (MeSH terms) in PubMed. Resources that tag genes to articles were integrated, then cross-species orthologs were identified using UniRef50 clusters. MeSH term frequency was normalized to reflect the MeSH tree structure, and then the resulting GeneID-MeSH associations were ranked using NPMI. The resulting network, called Entity MeSH Co-occurrence Network (EMCON), is a scalable resource for the identification and ranking of genes for a given topic of interest. The utility of EMCON was evaluated with the use case of breast carcinogenesis. Topics relevant to breast carcinogenesis were used to query EMCON and retrieve genes important to each topic. A breast cancer gene set was compiled through expert literature review (ELR) to assess performance of the search results. We found that the results from EMCON ranked the breast cancer genes from ELR higher than randomly selected genes with a recall of 0.98. Precision of the top five genes for selected topics was calculated as 0.87. This work demonstrates that EMCON can be used to link in vitro results to possible biological outcomes, thus aiding in generation of testable hypotheses for furthering understanding of biological function and the contribution of chemical exposures to disease.

QSAR Modeling of ToxCast Assays Relevant to the Molecular Initiating Events of AOPs Leading to Hepatic Steatosis

Nonalcoholic hepatic steatosis is a worldwide epidemiological concern since it is among the most prominent hepatic diseases. Indeed, research in toxicology and epidemiology has gathered evidence that exposure to endocrine disruptors can perturb cellular homeostasis and cause this disease. Therefore, assessing the likelihood of a chemical to trigger hepatic steatosis is a matter of the utmost importance. However, systematic in vivo testing of all the chemicals humans are exposed to is not feasible for ethical and economical reasons. In this context, predicting the molecular initiating events (MIE) leading to hepatic steatosis by QSAR modeling is an issue of practical relevance in modern toxicology. In this article, we present QSAR models based on random forest classifiers and DRAGON molecular descriptors for the prediction of in vitro assays that are relevant to MIEs leading to hepatic steatosis. These assays were provided by the ToxCast program and proved to be predictive for the detection of chemical-induced steatosis. During the modeling process, special attention was paid to chemical and toxicological data curation. We adopted two modeling strategies (undersampling and balanced random forests) to develop robust QSAR models from unbalanced data sets. The two modeling approaches gave similar results in terms of predictivity, and most of the models satisfy a minimum percentage of correctly predicted chemicals equal to 75%. Finally, and most importantly, the developed models proved to be useful as an effective in silico screening test for hepatic steatosis.

Uncertainty quantification in ToxCast high throughput screening

High throughput screening (HTS) projects like the U.S. Environmental Protection Agency's ToxCast program are required to address the large and rapidly increasing number of chemicals for which we have little to no toxicity measurements. Concentration-response parameters such as potency and efficacy are extracted from HTS data using nonlinear regression, and models and analyses built from these parameters are used to predict in vivo and in vitro toxicity of thousands of chemicals. How these predictions are impacted by uncertainties that stem from parameter estimation and propagated through the models and analyses has not been well explored. While data size and complexity makes uncertainty quantification computationally expensive for HTS datasets, continued advancements in computational resources have allowed these computational challenges to be met. This study uses nonparametric bootstrap resampling to calculate uncertainties in concentration-response parameters from a variety of HTS assays. Using the ToxCast estrogen receptor model for bioactivity as a case study, we highlight how these uncertainties can be propagated through models to quantify the uncertainty in model outputs. Uncertainty quantification in model outputs is used to identify potential false positives and false negatives and to determine the distribution of model values around semi-arbitrary activity cutoffs, increasing confidence in model predictions. At the individual chemical-assay level, curves with high variability are flagged for manual inspection or retesting, focusing subject-matter-expert time on results that need further input. This work improves the confidence of predictions made using HTS data, increasing the ability to use this data in risk assessment.

Multi-Descriptor Read Across (MuDRA): A Simple and Transparent Approach for Developing Accurate Quantitative Structure-Activity Relationship Models

Multiple approaches to quantitative structure-activity relationship (QSAR) modeling using various statistical or machine learning techniques and different types of chemical descriptors have been developed over the years. Oftentimes models are used in consensus to make more accurate predictions at the expense of model interpretation. We propose a simple, fast, and reliable method termed Multi-Descriptor Read Across (MuDRA) for developing both accurate and interpretable models. The method is conceptually related to the well-known kNN approach but uses different types of chemical descriptors simultaneously for similarity assessment. To benchmark the new method, we have built MuDRA models for six different end points (Ames mutagenicity, aquatic toxicity, hepatotoxicity, hERG liability, skin sensitization, and endocrine disruption) and compared the results with those generated with conventional consensus QSAR modeling. We find that models built with MuDRA show consistently high external accuracy similar to that of conventional QSAR models. However, MuDRA models excel in terms of transparency, interpretability, and computational efficiency. We posit that due to its methodological simplicity and reliable predictive accuracy, MuDRA provides a powerful alternative to a much more complex consensus QSAR modeling. MuDRA is implemented and freely available at the Chembench web portal ( https://chembench.mml.unc.edu/mudra ).

Beyond the RfD: Broad Application of a Probabilistic Approach to Improve Chemical Dose-Response Assessments for Noncancer Effects

BACKGROUND

The National Academies recommended risk assessments redefine the traditional noncancer Reference Dose (RfD) as a probabilistically derived risk-specific dose, a framework for which was recently developed by the World Health Organization (WHO).

OBJECTIVES

Our aim was to assess the feasibility and implications of replacing traditional RfDs with probabilistic estimates of the human dose associated with an effect magnitude M and population incidence I (HD_MI).

METHODS

We created a comprehensive, curated database of RfDs derived from animal data and developed a standardized, automated, web-accessible probabilistic dose-response workflow implementing the WHO framework.

RESULTS

We identified 1,464 RfDs and associated endpoints, representing 608 chemicals across many types of effects. Applying our standardized workflow resulted in 1,522 HD_MI values. Traditional RfDs are generally within an order of magnitude of the HD_MI lower confidence bound for I=1% and M values commonly used for benchmark doses. The greatest contributor to uncertainty was lack of benchmark dose estimates, followed by uncertainty in the extent of human variability. Exposure at the traditional RfD frequently implies an upper 95% confidence bound of several percent of the population affected. Whether such incidences are considered acceptable is likely to vary by chemical and risk context, especially given the wide range of severity of the associated effects, from clinical chemistry to mortality.

CONCLUSIONS

Overall, replacing RfDs with HD_MI estimates can provide a more consistent, scientifically rigorous, and transparent basis for risk management decisions, as well as support additional decision contexts such as economic benefit-cost analysis, risk-risk tradeoffs, life-cycle impact analysis, and emergency response. https://doi.org/10.1289/EHP3368.

Conditional Toxicity Value (CTV) Predictor: An In Silico Approach for Generating Quantitative Risk Estimates for Chemicals

BACKGROUND

Human health assessments synthesize human, animal, and mechanistic data to produce toxicity values that are key inputs to risk-based decision making. Traditional assessments are data-, time-, and resource-intensive, and they cannot be developed for most environmental chemicals owing to a lack of appropriate data.

OBJECTIVES

As recommended by the National Research Council, we propose a solution for predicting toxicity values for data-poor chemicals through development of quantitative structure-activity relationship (QSAR) models.

METHODS

We used a comprehensive database of chemicals with existing regulatory toxicity values from U.S. federal and state agencies to develop quantitative QSAR models. We compared QSAR-based model predictions to those based on high-throughput screening (HTS) assays.

RESULTS

QSAR models for noncancer threshold-based values and cancer slope factors had cross-validation-based Q² of 0.25-0.45, mean model errors of 0.70-1.11 log₁₀ units, and applicability domains covering >80% of environmental chemicals. Toxicity values predicted from QSAR models developed in this study were more accurate and precise than those based on HTS assays or mean-based predictions. A publicly accessible web interface to make predictions for any chemical of interest is available at http://toxvalue.org.

CONCLUSIONS

An in silico tool that can predict toxicity values with an uncertainty of an order of magnitude or less can be used to quickly and quantitatively assess risks of environmental chemicals when traditional toxicity data or human health assessments are unavailable. This tool can fill a critical gap in the risk assessment and management of data-poor chemicals. https://doi.org/10.1289/EHP2998.

A Workflow for Identifying Metabolically Active Chemicals to Complement in vitro Toxicity Screening

The new paradigm of toxicity testing approaches involves rapid screening of thousands of chemicals across hundreds of biological targets through use of in vitro assays. Such assays may lead to false negatives when the complex metabolic processes that render a chemical bioactive in a living system are unable to be replicated in an in vitro environment. In the current study, a workflow is presented for complementing in vitro testing results with in silico and in vitro techniques to identify inactive parents that may produce active metabolites. A case study applying this workflow involved investigating the influence of metabolism for over 1,400 chemicals considered inactive across18 in vitro assays related to the estrogen receptor (ER) pathway. Over 7,500 first-generation and second-generation metabolites were generated for these in vitro inactive chemicals using an in silico software program. Next, a consensus model comprised of four individual quantitative structure activity relationship (QSAR) models was used to predict ER-binding activity for each of the metabolites. Binding activity was predicted for ~8-10% of metabolites in each generation, with these metabolites linked to 259 in vitro inactive parent chemicals. Metabolites were enriched in substructures consisting of alcohol, aromatic, and phenol bonds relative to their inactive parent chemicals, suggesting these features are potentially favorable for ER-binding. The workflow presented here can be used to identify parent chemicals that can be potentially bioactive, to aid confidence in high throughput risk screening.

In silico identification of protein targets for chemical neurotoxins using ToxCast in vitro data and read-across within the QSAR toolbox

This study evaluates the application of QSAR Toolbox and ToxCast screening data to identify neurological targets for pyrethroids.

There are many mechanisms of neurotoxicity that are initiated by the interaction of chemicals with different neurological targets. Under the U.S. Environmental Protection Agency's ToxCast program, the biological activity of thousands of chemicals was screened in biochemical and cell-based assays in a high-throughput manner. Two hundred sixteen assays in the ToxCast screening database were identified as targeting a total of 123 proteins having neurological functions according to the Gene Ontology database. Data from these assays were imported into the Organization for Economic Co-operation and Development QSAR Toolbox and used to predict neurological targets for chemical neurotoxins. Two sets of data were generated: one set was used to classify compounds as active or inactive and another set, composed of AC₅₀s for only active compounds, was used to predict AC₅₀ values for unknown chemicals. Chemical grouping and read-across within the QSAR Toolbox were used to identify neurologic targets and predict interactions for pyrethroids, a class of compounds known to elicit neurotoxic effects in humans. The classification prediction results showed 79% accuracy while AC₅₀ predictions demonstrated mixed accuracy compared with the ToxCast screening data.

OPERA models for predicting physicochemical properties and environmental fate endpoints

The collection of chemical structure information and associated experimental data for quantitative structure–activity/property relationship (QSAR/QSPR) modeling is facilitated by an increasing number of public databases containing large amounts of useful data. However, the performance of QSAR models highly depends on the quality of the data and modeling methodology used. This study aims to develop robust QSAR/QSPR models for chemical properties of environmental interest that can be used for regulatory purposes. This study primarily uses data from the publicly available PHYSPROP database consisting of a set of 13 common physicochemical and environmental fate properties. These datasets have undergone extensive curation using an automated workflow to select only high-quality data, and the chemical structures were standardized prior to calculation of the molecular descriptors. The modeling procedure was developed based on the five Organization for Economic Cooperation and Development (OECD) principles for QSAR models. A weighted k-nearest neighbor approach was adopted using a minimum number of required descriptors calculated using PaDEL, an open-source software. The genetic algorithms selected only the most pertinent and mechanistically interpretable descriptors (2–15, with an average of 11 descriptors). The sizes of the modeled datasets varied from 150 chemicals for biodegradability half-life to 14,050 chemicals for logP, with an average of 3222 chemicals across all endpoints. The optimal models were built on randomly selected training sets (75%) and validated using fivefold cross-validation (CV) and test sets (25%). The CV Q² of the models varied from 0.72 to 0.95, with an average of 0.86 and an R² test value from 0.71 to 0.96, with an average of 0.82. Modeling and performance details are described in QSAR model reporting format and were validated by the European Commission’s Joint Research Center to be OECD compliant. All models are freely available as an open-source, command-line application called OPEn structure–activity/property Relationship App (OPERA). OPERA models were applied to more than 750,000 chemicals to produce freely available predicted data on the U.S. Environmental Protection Agency’s CompTox Chemistry Dashboard.

High-throughput in-silico prediction of ionization equilibria for pharmacokinetic modeling

Chemical ionization plays an important role in many aspects of pharmacokinetic (PK) processes such as protein binding, tissue partitioning, and apparent volume of distribution at steady state (Vdss). Here, estimates of ionization equilibrium constants (i.e., pKa) were analyzed for 8132 pharmaceuticals and 24,281 other compounds to which humans might be exposed in the environment. Results revealed broad differences in the ionization of pharmaceutical chemicals and chemicals with either near-field (in the home) or far-field sources. The utility of these high-throughput ionization predictions was evaluated via a case-study of predicted PK Vdss for 22 compounds monitored in the blood and serum of the U.S. population by the U.S. Centers for Disease Control and Prevention National Health and Nutrition Examination Survey (NHANES). The chemical distribution ratio between water and tissue was estimated using predicted ionization states characterized by pKa. Probability distributions corresponding to ionizable atom types (IATs) were then used to analyze the sensitivity of predicted Vdss on predicted pKa using Monte Carlo methods. 8 of the 22 compounds were predicted to be ionizable. For 5 of the 8 the predictions based upon ionization are significantly different from what would be predicted for a neutral compound. For all but one (foramsulfuron), the probability distribution of predicted Vdss generated by IAT sensitivity analysis spans both the neutral prediction and the prediction using ionization. As new data sets of chemical-specific information on metabolism and excretion for hundreds of chemicals are being made available (e.g., Wetmore et al., 2015), high-throughput methods for calculating Vdss and tissue-specific PK distribution coefficients will allow the rapid construction of PK models to provide context for both biomonitoring data and high-throughput toxicity screening studies such as Tox21 and ToxCast.

Assessment of the Endocrine-Disrupting Effects of Trichloroethylene and Its Metabolites Using in Vitro and in Silico Approaches

Trichloroethylene (TCE) is a ubiquitous environmental contaminant, which may have effects on both ecosystem and human health. TCE has been reported to cause several toxic effects, but little effort has been made to assess the ecological risks of TCE or its major metabolites: trichloroethanol (TCOH), trichloroacetic acid, and oxalic acid (OA). In this study, the endocrine-disrupting potential of TCE and its metabolites were investigated using in vitro and in silico approaches. We examined alterations in the steroidogenesis pathway using the NCI-H295R cell line and utilized receptor-mediated luciferase reporter cell lines to identify effects on estrogen and androgen receptors. Molecular docking was also used to explore chemical interactions with these receptors. All test chemicals except OA significantly increased 17β-estradiol production which can be attributed to an up-regulation of 17β-hydroxysteroid dehydrogenase. Moreover, TCOH exhibited significant antiestrogenic activity with a RIC20 (20% relative inhibitory concentration) of 3.7 × 10-7 M. Molecular docking simulation supported this finding with lower docking scores for TCOH, indicating that hydrogen bonds may stabilize the interaction between TCOH and the estrogen receptor binding pocket. These findings suggest that TCE contamination poses an endocrine-disrupting threat, which has implications for both ecological and human health.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Molecular Descriptors for Structure-Activity Applications: A Hands-On Approach

Molecular descriptors capture diverse parts of the structural information of molecules and they are the support of many contemporary computer-assisted toxicological and chemical applications. After briefly introducing some fundamental concepts of structure-activity applications (e.g., molecular descriptor dimensionality, classical vs. fingerprint description, and activity landscapes), this chapter guides the readers through a step-by-step explanation of molecular descriptors rationale and application. To this end, the chapter illustrates a case study of a recently published application of molecular descriptors for modeling the activity on cytochrome P450.

Molecular Docking for Predictive Toxicology

Molecular docking is an in silico method widely applied in drug discovery programs to predict the binding mode of a given molecule interacting with a specific biological target. This computational technique is today emerging also in the field of predictive toxicology for regulatory purposes, being for instance successfully applied to develop classification models for the prediction of the endocrine disruptor potential of chemicals. Herein, we describe the protocol for adapting molecular docking to the purposes of predictive toxicology.

Computer-aided prediction of biological activity spectra for chemical compounds: opportunities and limitation

An essential characteristic of chemical compounds is their biological activity since its presence can become the basis for the use of the substance for therapeutic purposes, or, on the contrary, limit the possibilities of its practical application due to the manifestation of side action and toxic effects. Computer assessment of the biological activity spectra makes it possible to determine the most promising directions for the study of the pharmacological action of particular substances, and to filter out potentially dangerous molecules at the early stages of research. For more than 25 years, we have been developing and improving the computer program PASS (Prediction of Activity Spectra for Substances), designed to predict the biological activity spectrum of substance based on the structural formula of its molecules. The prediction is carried out by the analysis of structure-activity relationships for the training set, which currently contains information on structures and known biological activities for more than one million molecules. The structure of the organic compound is represented in PASS using Multilevel Neighborhoods of Atoms descriptors; the activity prediction for new compounds is performed by the naive Bayes classifier and the structure-activity relationships determined by the analysis of the training set. We have created and improved both local versions of the PASS program and freely available web resources based on PASS (http://www.way2drug.com). They predict several thousand biological activities (pharmacological effects, molecular mechanisms of action, specific toxicity and adverse effects, interaction with the unwanted targets, metabolism and action on molecular transport), cytotoxicity for tumor and non-tumor cell lines, carcinogenicity, induced changes of gene expression profiles, metabolic sites of the major enzymes of the first and second phases of xenobiotics biotransformation, and belonging to substrates and/or metabolites of metabolic enzymes. The web resource Way2Drug is used by over 18,000 researchers from more than 90 countries around the world, which allowed them to obtain over 600,000 predictions and publish about 500 papers describing the obtained results. The analysis of the published works shows that in some cases the interpretation of the prediction results presented by the authors of these publications requires an adjustment. In this work, we provide the theoretical basis and consider, on particular examples, the opportunities and limitations of computer-aided prediction of biological activity spectra.

Impact of Molecular Descriptors on Computational Models

Molecular descriptors encode a wide variety of molecular information and have become the support of many contemporary chemoinformatic and bioinformatic applications. They grasp specific molecular features (e.g., geometry, shape, pharmacophores, or atomic properties) and directly affect computational models, in terms of outcome, performance, and applicability. This chapter aims to illustrate the impact of different molecular descriptors on the structural information captured and on the perceived chemical similarity among molecules. After introducing the fundamental concepts of molecular descriptor theory and application, a step-by-step retrospective virtual screening procedure guides users through the fundamental processing steps and discusses the impact of different types of molecular descriptors.

Dark matter in host-microbiome metabolomics: Tackling the unknowns-A review

The "dark matter" in metabolomics (unknowns) represents an exciting frontier with significant potential for discovery in relation to biochemistry, yet it also presents one of the largest challenges to overcome. This focussed review takes a close look at the current state-of-the-art and future challenges in tackling the unknowns with specific focus on the human gut microbiome and host-microbe interactions. Metabolomics, like metabolism itself, is a very dynamic discipline, with many workflows and methods under development, both in terms of chemical analysis and post-analysis data processing. Here, we look at developments in the mutli-omic analyses and the use of mass spectrometry to investigate the exchange of metabolites between the host and the microbiome as well as the environment within the microbiome. A case study using HuMiX, a microfluidics-based human-microbial co-culture system that enables the co-culture of human and microbial cells under controlled conditions, is used to highlight opportunities and current limitations. Common definitions, approaches, databases and elucidation techniques from both the environmental and metabolomics fields are covered, with perspectives on how to merge these, as the boundaries blur between the fields. While reflecting on the number of unknowns remaining to be conquered in typical complex samples measured with mass spectrometry (often orders of magnitude above the "knowns"), we provide an outlook on future perspectives and challenges in elucidating the relevant "dark matter".

Development of Decision Forest Models for Prediction of Drug-Induced Liver Injury in Humans Using A Large Set of FDA-approved Drugs

Drug-induced liver injury (DILI) presents a significant challenge to drug development and regulatory science. The FDA’s Liver Toxicity Knowledge Base (LTKB) evaluated >1000 drugs for their likelihood of causing DILI in humans, of which >700 drugs were classified into three categories (most-DILI, less-DILI, and no-DILI). Based on this dataset, we developed and compared 2-class and 3-class DILI prediction models using the machine learning algorithm of Decision Forest (DF) with Mold2 structural descriptors. The models were evaluated through 1000 iterations of 5-fold cross-validations, 1000 bootstrapping validations and 1000 permutation tests (that assessed the chance correlation). Furthermore, prediction confidence analysis was conducted, which provides an additional parameter for proper interpretation of prediction results. We revealed that the 3-class model not only had a higher resolution to estimate DILI risk but also showed an improved capability to differentiate most-DILI drugs from no-DILI drugs in comparison with the 2-class DILI model. We demonstrated the utility of the models for drug ingredients with warnings very recently issued by the FDA. Moreover, we identified informative molecular features important for assessing DILI risk. Our results suggested that the 3-class model presents a better option than the binary model (which most publications are focused on) for drug safety evaluation.

A systematic evaluation of analogs and automated read-across prediction of estrogenicity: A case study using hindered phenols

Read-across is an important data gap filling technique used within category and analog approaches for regulatory hazard identification and risk assessment. Although much technical guidance is available that describes how to develop category/analog approaches, practical principles to evaluate and substantiate analog validity (suitability) are still lacking. This case study uses hindered phenols as an example chemical class to determine: (1) the capability of three structure fingerprint/descriptor methods (PubChem, ToxPrints and MoSS MCSS) to identify analogs for read-across to predict Estrogen Receptor (ER) binding activity and, (2) the utility of data confidence measures, physicochemical properties, and chemical R-group properties as filters to improve ER binding predictions. The training dataset comprised 462 hindered phenols and 257 non- hindered phenols. For each chemical of interest (target), source analogs were identified from two datasets (hindered and non-hindered phenols) that had been characterized by a fingerprint/descriptor method and by two cut-offs: (1) minimum similarity distance (range: 0.1 - 0.9) and, (2) N closest analogs (range: 1 - 10). Analogs were then filtered using: (1) physicochemical properties of the phenol (termed global filtering) and, (2) physicochemical properties of the R-groups neighboring the active hydroxyl group (termed local filtering). A read-across prediction was made for each target chemical on the basis of a majority vote of the N closest analogs. The results demonstrate that: (1) concordance in ER activity increases with structural similarity, regardless of the structure fingerprint/descriptor method, (2) increased data confidence significantly improves read-across predictions, and (3) filtering analogs using global and local properties can help identify more suitable analogs. This case study illustrates that the quality of the underlying experimental data and use of endpoint relevant chemical descriptors to evaluate source analogs are critical to achieving robust read-across predictions.

Predictive Structure-Based Toxicology Approaches To Assess the Androgenic Potential of Chemicals

We present a practical and easy-to-run in silico workflow exploiting a structure-based strategy making use of docking simulations to derive highly predictive classification models of the androgenic potential of chemicals. Models were trained on a high-quality chemical collection comprising 1689 curated compounds made available within the CoMPARA consortium from the US Environmental Protection Agency and were integrated with a two-step applicability domain whose implementation had the effect of improving both the confidence in prediction and statistics by reducing the number of false negatives. Among the nine androgen receptor X-ray solved structures, the crystal 2PNU (entry code from the Protein Data Bank) was associated with the best performing structure-based classification model. Three validation sets comprising each 2590 compounds extracted by the DUD-E collection were used to challenge model performance and the effectiveness of Applicability Domain implementation. Next, the 2PNU model was applied to screen and prioritize two collections of chemicals. The first is a small pool of 12 representative androgenic compounds that were accurately classified based on outstanding rationale at the molecular level. The second is a large external blind set of 55450 chemicals with potential for human exposure. We show how the use of molecular docking provides highly interpretable models and can represent a real-life option as an alternative nontesting method for predictive toxicology.

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_α.

Predicting Organ Toxicity Using in Vitro Bioactivity Data and Chemical Structure

Animal testing alone cannot practically evaluate the health hazard posed by tens of thousands of environmental chemicals. Computational approaches making use of high-throughput experimental data may provide more efficient means to predict chemical toxicity. Here, we use a supervised machine learning strategy to systematically investigate the relative importance of study type, machine learning algorithm, and type of descriptor on predicting in vivo repeat-dose toxicity at the organ-level. A total of 985 compounds were represented using chemical structural descriptors, ToxPrint chemotype descriptors, and bioactivity descriptors from ToxCast in vitro high-throughput screening assays. Using ToxRefDB, a total of 35 target organ outcomes were identified that contained at least 100 chemicals (50 positive and 50 negative). Supervised machine learning was performed using Naïve Bayes, k-nearest neighbor, random forest, classification and regression trees, and support vector classification approaches. Model performance was assessed based on F1 scores using 5-fold cross-validation with balanced bootstrap replicates. Fixed effects modeling showed the variance in F1 scores was explained mostly by target organ outcome, followed by descriptor type, machine learning algorithm, and interactions between these three factors. A combination of bioactivity and chemical structure or chemotype descriptors were the most predictive. Model performance improved with more chemicals (up to a maximum of 24%), and these gains were correlated (ρ = 0.92) with the number of chemicals. Overall, the results demonstrate that a combination of bioactivity and chemical descriptors can accurately predict a range of target organ toxicity outcomes in repeat-dose studies, but specific experimental and methodologic improvements may increase predictivity.

Evaluation of estrogen receptor alpha activation by glyphosate-based herbicide constituents

The safety, including the endocrine disruptive capability, of glyphosate-based herbicides (GBHs) is a matter of intense debate. We evaluated the estrogenic potential of glyphosate, commercial GBHs and polyethoxylated tallowamine adjuvants present as co-formulants in GBHs. Glyphosate (≥10,000 μg/L or 59 μM) promoted proliferation of estrogen-dependent MCF-7 human breast cancer cells. Glyphosate also increased the expression of an estrogen response element-luciferase reporter gene (ERE-luc) in T47D-KBluc cells, which was blocked by the estrogen antagonist ICI 182,780. Commercial GBH formulations or their adjuvants alone did not exhibit estrogenic effects in either assay. Transcriptomics analysis of MCF-7 cells treated with glyphosate revealed changes in gene expression reflective of hormone-induced cell proliferation but did not overlap with an ERα gene expression biomarker. Calculation of glyphosate binding energy to ERα predicts a weak and unstable interaction (-4.10 kcal mol-1) compared to estradiol (-25.79 kcal mol-1), which suggests that activation of this receptor by glyphosate is via a ligand-independent mechanism. Induction of ERE-luc expression by the PKA signalling activator IBMX shows that ERE-luc is responsive to ligand-independent activation, suggesting a possible mechanism of glyphosate-mediated activation. Our study reveals that glyphosate, but not other components present in GBHs, can activate ERα in vitro, albeit at relatively high concentrations.

Grouping of CoMPARA data with respect to compounds from the carcinogenic potency database

Methods for clustering and measures of similarity of chemical structures have become an important supporting tool in chemoinformatics. They represent the basis for categorization of chemicals and read-across, where a molecular property is estimated from 'similar molecules'. This study proposes a clustering scheme within the given dataset with respect to a reference dataset. The scheme was applied on two datasets ToxCast_AR_Agonist and ToxCast_AR_Antagonists with 1654 and 1522 compounds, respectively. The compounds are tested to androgen receptor activity (AR) in 11 high throughput screening assays. The carcinogenic dataset was used as the reference set.

Molecular dynamics simulations and applications in computational toxicology and nanotoxicology

Nanotoxicology studies toxicity of nanomaterials and has been widely applied in biomedical researches to explore toxicity of various biological systems. Investigating biological systems through in vivo and in vitro methods is expensive and time taking. Therefore, computational toxicology, a multi-discipline field that utilizes computational power and algorithms to examine toxicology of biological systems, has gained attractions to scientists. Molecular dynamics (MD) simulations of biomolecules such as proteins and DNA are popular for understanding of interactions between biological systems and chemicals in computational toxicology. In this paper, we review MD simulation methods, protocol for running MD simulations and their applications in studies of toxicity and nanotechnology. We also briefly summarize some popular software tools for execution of MD simulations.

Ligand efficiency metrics in drug discovery: the pros and cons from a practical perspective

INTRODUCTION

Ligand efficiency metrics are almost universally accepted as a valuable indicator of compound quality and an aid to reduce attrition. Areas covered: In this review, the authors describe ligand efficiency metrics giving a balanced overview on their merits and points of weakness in order to enable the readers to gain an informed opinion. Relevant theoretical breakthroughs and drug-like properties are also illustrated. Several recent exemplary case studies are discussed in order to illustrate the main fields of application of ligand efficiency metrics. Expert opinion: As a medicinal chemist guide, ligand efficiency metrics perform in a context- and chemotype-dependent manner; thus, they should not be used as a magic box. Since the 'big bang' of efficiency metrics occurred more or less ten years ago and the average time to develop a new drug is over the same period, the next few years will give a clearer outlook on the increased rate of success, if any, gained by means of these new intriguing tools.

Navigating through the minefield of read-across tools: A review of in silico tools for grouping

Read-across is a popular data gap filling technique used within analogue and category approaches for regulatory purposes. In recent years there have been many efforts focused on the challenges involved in read-across development, its scientific justification and documentation. Tools have also been developed to facilitate read-across development and application. Here, we describe a number of publicly available read-across tools in the context of the category/analogue workflow and review their respective capabilities, strengths and weaknesses. No single tool addresses all aspects of the workflow. We highlight how the different tools complement each other and some of the opportunities for their further development to address the continued evolution of read-across.

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.

Chemical Risk Assessment: Traditional vs Public Health Perspectives

Preventing adverse health effects of environmental chemical exposure is fundamental to protecting individual and public health. When done efficiently and properly, chemical risk assessment enables risk management actions that minimize the incidence and effects of environmentally induced diseases related to chemical exposure. However, traditional chemical risk assessment is faced with multiple challenges with respect to predicting and preventing disease in human populations, and epidemiological studies increasingly report observations of adverse health effects at exposure levels predicted from animal studies to be safe for humans. This discordance reinforces concerns about the adequacy of contemporary risk assessment practices for protecting public health. It is becoming clear that to protect public health more effectively, future risk assessments will need to use the full range of available data, draw on innovative methods to integrate diverse data streams, and consider health endpoints that also reflect the range of subtle effects and morbidities observed in human populations. Considering these factors, there is a need to reframe chemical risk assessment to be more clearly aligned with the public health goal of minimizing environmental exposures associated with disease.

Prediction of Hydrolysis Products of Organic Chemicals under Environmental pH Conditions

Cheminformatics-based software tools can predict the molecular structure of transformation products using a library of transformation reaction schemes. This paper presents the development of such a library for abiotic hydrolysis of organic chemicals under environmentally relevant conditions. The hydrolysis reaction schemes in the library encode the process science gathered from peer-reviewed literature and regulatory reports. Each scheme has been ranked on a scale of one to six based on the median half-life in a data set compiled from literature-reported hydrolysis rates. These ranks are used to predict the most likely transformation route when more than one structural fragment susceptible to hydrolysis is present in a molecule of interest. Separate rank assignments are established for pH 5, 7, and 9 to represent standard conditions in hydrolysis studies required for registration of pesticides in Organisation for Economic Co-operation and Development (OECD) member countries. The library is applied to predict the likely hydrolytic transformation products for two lists of chemicals, one representative of chemicals used in commerce and the other specific to pesticides, to evaluate which hydrolysis reaction pathways are most likely to be relevant for organic chemicals found in the natural environment.

High-Performance Prediction of Human Estrogen Receptor Agonists Based on Chemical Structures

Many agonists for the estrogen receptor are known to disrupt endocrine functioning. We have developed a computational model that predicts agonists for the estrogen receptor ligand-binding domain in an assay system. Our model was entered into the Tox21 Data Challenge 2014, a computational toxicology competition organized by the National Center for Advancing Translational Sciences. This competition aims to find high-performance predictive models for various adverse-outcome pathways, including the estrogen receptor. Our predictive model, which is based on the random forest method, delivered the best performance in its competition category. In the current study, the predictive performance of the random forest models was improved by strictly adjusting the hyperparameters to avoid overfitting. The random forest models were optimized from 4000 descriptors simultaneously applied to 10,000 activity assay results for the estrogen receptor ligand-binding domain, which have been measured and compiled by Tox21. Owing to the correlation between our model's and the challenge's results, we consider that our model currently possesses the highest predictive power on agonist activity of the estrogen receptor ligand-binding domain. Furthermore, analysis of the optimized model revealed some important features of the agonists, such as the number of hydroxyl groups in the molecules.

Critical Assessment of Small Molecule Identification 2016: automated methods

BACKGROUND

The fourth round of the Critical Assessment of Small Molecule Identification (CASMI) Contest ( www.casmi-contest.org ) was held in 2016, with two new categories for automated methods. This article covers the 208 challenges in Categories 2 and 3, without and with metadata, from organization, participation, results and post-contest evaluation of CASMI 2016 through to perspectives for future contests and small molecule annotation/identification.

RESULTS

The Input Output Kernel Regression (CSI:IOKR) machine learning approach performed best in "Category 2: Best Automatic Structural Identification-In Silico Fragmentation Only", won by Team Brouard with 41% challenge wins. The winner of "Category 3: Best Automatic Structural Identification-Full Information" was Team Kind (MS-FINDER), with 76% challenge wins. The best methods were able to achieve over 30% Top 1 ranks in Category 2, with all methods ranking the correct candidate in the Top 10 in around 50% of challenges. This success rate rose to 70% Top 1 ranks in Category 3, with candidates in the Top 10 in over 80% of the challenges. The machine learning and chemistry-based approaches are shown to perform in complementary ways.

CONCLUSIONS

The improvement in (semi-)automated fragmentation methods for small molecule identification has been substantial. The achieved high rates of correct candidates in the Top 1 and Top 10, despite large candidate numbers, open up great possibilities for high-throughput annotation of untargeted analysis for "known unknowns". As more high quality training data becomes available, the improvements in machine learning methods will likely continue, but the alternative approaches still provide valuable complementary information. Improved integration of experimental context will also improve identification success further for "real life" annotations. The true "unknown unknowns" remain to be evaluated in future CASMI contests. Graphical abstract .

Rigorous 3-dimensional spectral data activity relationship approach modeling strategy for ToxCast estrogen receptor data classification, validation, and feature extraction

The estrogenic potential (expressed as a score composite of 18 high throughput screening bioassays) of 1528 compounds from the ToxCast database was modeled by a 3-dimensional spectral data activity relationship approach (3D-SDAR). Due to a lack of ¹⁷ O nuclear magnetic resonance (NMR) simulation software, the most informative carbon-carbon 3D-SDAR fingerprints were augmented with indicator variables representing oxygen atoms from carbonyl and carboxamide, ester, sulfonyl, nitro, aliphatic hydroxyl, and phenolic hydroxyl groups. To evaluate the true predictive performance of the authors' model the United States Environmental Protection Agency provided them with a blind test set consisting of 2008 compounds. Of these, 543 had available literature data-their binding affinity served to estimate the external classification accuracy of the developed model: predictive accuracy of 0.62, sensitivity of 0.71, and specificity of 0.53 were obtained. Compared with alternative modeling techniques, the authors' model displayed very little reduction in performance between the modeling and the prediction set. A 3D-SDAR mapping technique allowed identification of structural features essential for estrogenicity: 1) the presence of a phenolic OH group or cyclohexenone, 2) a second aromatic or phenolic ring at a distance of 6 Å to 8 Å from the oxygen of the first phenol ring, 3) the presence of a methyl group approximately 6 Å away from the centroid of a phenol ring, and 4) a carbonyl group in close proximity (∼4 Å measured to the centroid) to 1 of the phenol rings. Environ Toxicol Chem 2017;36:823-830. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

In vitro cardiotoxicity assessment of environmental chemicals using an organotypic human induced pluripotent stem cell-derived model

An important target area for addressing data gaps through in vitro screening is the detection of potential cardiotoxicants. Despite the fact that current conservative estimates relate at least 23% of all cardiovascular disease cases to environmental exposures, the identities of the causative agents remain largely uncharacterized. Here, we evaluate the feasibility of a combinatorial in vitro/in silico screening approach for functional and mechanistic cardiotoxicity profiling of environmental hazards using a library of 69 representative environmental chemicals and drugs. Human induced pluripotent stem cell-derived cardiomyocytes were exposed in concentration-response for 30min or 24h and effects on cardiomyocyte beating and cellular and mitochondrial toxicity were assessed by kinetic measurements of intracellular Ca²⁺ flux and high-content imaging using the nuclear dye Hoechst 33342, the cell viability marker Calcein AM, and the mitochondrial depolarization probe JC-10. More than half of the tested chemicals exhibited effects on cardiomyocyte beating after 30min of exposure. In contrast, after 24h, effects on cell beating without concomitant cytotoxicity were observed in about one third of the compounds. Concentration-response data for in vitro bioactivity phenotypes visualized using the Toxicological Prioritization Index (ToxPi) showed chemical class-specific clustering of environmental chemicals, including pesticides, flame retardants, and polycyclic aromatic hydrocarbons. For environmental chemicals with human exposure predictions, the activity-to-exposure ratios between modeled blood concentrations and in vitro bioactivity were between one and five orders of magnitude. These findings not only demonstrate that some ubiquitous environmental pollutants might have the potential at high exposure levels to alter cardiomyocyte function, but also indicate similarities in the mechanism of these effects both within and among chemicals and classes.

A comprehensive review of regulatory test methods for endocrine adverse health effects

Development of new endocrine disruption-relevant test methods has been the subject of intensive research efforts for the past several decades, prompted in part by mandates in the 1996 Food Quality Protection Act (FQPA). While scientific understanding and test methods have advanced, questions remain on whether current scientific methods are capable of adequately addressing the complexities of the endocrine system for regulatory health and ecological risk assessments. The specific objective of this article is to perform a comprehensive, detailed evaluation of the adequacy of current test methods to inform regulatory risk assessments of whether a substance has the potential to perturb endocrine-related pathways resulting in human adverse effects. To that end,  approximately 42 existing test guidelines (TGs) were considered in the evaluation of coverage for endocrine-related adverse effects. In addition to evaluations of whether test methods are adequate to capture endocrine-related effects, considerations of further enhancements to current test methods, along with the need to develop novel test methods to address existing test method gaps are described. From this specific evaluation, up to 35 test methods are capable of informing whether a chemical substance perturbs known endocrine related biological pathways. Based on these findings, it can be concluded that current validated test methods are adequate to discern substances that may perturb the endocrine system, resulting in an adverse health effect. Together, these test methods predominantly form the core data requirements of a typical food-use pesticide registration submission. It is recognized, however, that the current state of science is rapidly advancing and there is a need to update current test methods to include added enhancements to ensure continued coverage and public health and environmental protection.

Integrated Framework for Identifying Toxic Transformation Products in Complex Environmental Mixtures

Complex environmental mixtures consist of hundreds to thousands of unknown and unregulated organic compounds that may have toxicological relevance, including transformation products (TPs) of anthropogenic organic pollutants. Non-targeted analysis and suspect screening analysis offer analytical approaches for potentially identifying these toxic transformation products. However, additional tools and strategies are needed in order to reduce the number of chemicals of interest and focus analytical efforts on chemicals that may pose risks to humans and the environment. This brief review highlights recent developments in this field and suggests an integrated framework that incorporates complementary instrumental techniques, computational chemistry, and toxicity analysis, for prioritizing and identifying toxic TPs in the environment.