Improving the human hazard characterization of chemicals: a Tox21 update

The role of analytical chemistry in exposure science: Focus on the aquatic environment

Exposure science, in its broadest sense, studies the interactions between stressors (chemical, biological, and physical agents) and receptors (e.g. humans and other living organisms, and non-living items like buildings), together with the associated pathways and processes potentially leading to negative effects on human health and the environment. The aquatic environment may contain thousands of compounds, many of them still unknown, that can pose a risk to ecosystems and human health. Due to the unquestionable importance of the aquatic environment, one of the main challenges in the field of exposure science is the comprehensive characterization and evaluation of complex environmental mixtures beyond the classical/priority contaminants to new emerging contaminants. The role of advanced analytical chemistry to identify and quantify potential chemical risks, that might cause adverse effects to the aquatic environment, is essential. In this paper, we present the strategies and tools that analytical chemistry has nowadays, focused on chromatography hyphenated to (high-resolution) mass spectrometry because of its relevance in this field. Key issues, such as the application of effect direct analysis to reduce the complexity of the sample, the investigation of the huge number of transformation/degradation products that may be present in the aquatic environment, the analysis of urban wastewater as a source of valuable information on our lifestyle and substances we consumed and/or are exposed to, or the monitoring of drinking water, are discussed in this article. The trends and perspectives for the next few years are also highlighted, when it is expected that new developments and tools will allow a better knowledge of chemical composition in the aquatic environment. This will help regulatory authorities to protect water bodies and to advance towards improved regulations that enable practical and efficient abatements for environmental and public health protection.

The NCATS BioPlanet - An Integrated Platform for Exploring the Universe of Cellular Signaling Pathways for Toxicology, Systems Biology, and Chemical Genomics

Chemical genomics aims to comprehensively define, and ultimately predict, the effects of small molecule compounds on biological systems. Chemical activity profiling approaches must consider chemical effects on all pathways operative in mammalian cells. To enable a strategic and maximally efficient chemical profiling of pathway space, we have created the NCATS BioPlanet, a comprehensive integrated pathway resource that incorporates the universe of 1,658 human pathways sourced from publicly available, manually curated sources, which have been subjected to thorough redundancy and consistency cross-evaluation. BioPlanet supports interactive browsing, retrieval, and analysis of pathways, exploration of pathway connections, and pathway search by gene targets, category, and availability of corresponding bioactivity assay, as well as visualization of pathways on a 3-dimensional globe, in which the distance between any two pathways is proportional to their degree of gene component overlap. Using this resource, we propose a strategy to identify a minimal set of 362 biological assays that can interrogate the universe of human pathways. The NCATS BioPlanet is a public resource, which will be continually expanded and updated, for systems biology, toxicology, and chemical genomics, available at http://tripod.nih.gov/bioplanet/.

Improving the Utility of the Tox21 Dataset by Deep Metadata Annotations and Constructing Reusable Benchmarked Chemical Reference Signatures

The Toxicology in the 21st Century (Tox21) project seeks to develop and test methods for high-throughput examination of the effect certain chemical compounds have on biological systems. Although primary and toxicity assay data were readily available for multiple reporter gene modified cell lines, extensive annotation and curation was required to improve these datasets with respect to how FAIR (Findable, Accessible, Interoperable, and Reusable) they are. In this study, we fully annotated the Tox21 published data with relevant and accepted controlled vocabularies. After removing unreliable data points, we aggregated the results and created three sets of signatures reflecting activity in the reporter gene assays, cytotoxicity, and selective reporter gene activity, respectively. We benchmarked these signatures using the chemical structures of the tested compounds and obtained generally high receiver operating characteristic (ROC) scores, suggesting good quality and utility of these signatures and the underlying data. We analyzed the results to identify promiscuous individual compounds and chemotypes for the three signature categories and interpreted the results to illustrate the utility and re-usability of the datasets. With this study, we aimed to demonstrate the importance of data standards in reporting screening results and high-quality annotations to enable re-use and interpretation of these data. To improve the data with respect to all FAIR criteria, all assay annotations, cleaned and aggregate datasets, and signatures were made available as standardized dataset packages (Aggregated Tox21 bioactivity data, 2019).

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.

The Carcinogenome Project: In Vitro Gene Expression Profiling of Chemical Perturbations to Predict Long-Term Carcinogenicity

BACKGROUND

Most chemicals in commerce have not been evaluated for their carcinogenic potential. The de facto gold-standard approach to carcinogen testing adopts the 2-y rodent bioassay, a time-consuming and costly procedure. High-throughput in vitro assays are a promising alternative for addressing the limitations in carcinogen screening.

OBJECTIVES

We developed a screening process for predicting chemical carcinogenicity and genotoxicity and characterizing modes of actions (MoAs) using in vitro gene expression assays.

METHODS

We generated a large toxicogenomics resource comprising [Formula: see text] expression profiles corresponding to 330 chemicals profiled in HepG2 (human hepatocellular carcinoma cell line) at multiple doses and replicates. Predictive models of carcinogenicity and genotoxicity were built using a random forest classifier. Differential pathway enrichment analysis was performed to identify pathways associated with carcinogen exposure. Signatures of carcinogenicity and genotoxicity were compared with external sources, including Drugmatrix and the Connectivity Map.

RESULTS

Among profiles with sufficient bioactivity, our classifiers achieved 72.2% Area Under the ROC Curve (AUC) for predicting carcinogenicity and 82.3% AUC for predicting genotoxicity. Chemical bioactivity, as measured by the strength and reproducibility of the transcriptional response, was not significantly associated with long-term carcinogenicity in doses up to [Formula: see text]. However, sufficient bioactivity was necessary for a chemical to be used for prediction of carcinogenicity. Pathway enrichment analysis revealed pathways consistent with known pathways that drive cancer, including DNA damage and repair. The data is available at https://clue.io/CRCGN_ABC , and a portal for query and visualization of the results is accessible at https://carcinogenome.org .

DISCUSSION

We demonstrated an in vitro screening approach using gene expression profiling to predict carcinogenicity and infer MoAs of chemical perturbations. https://doi.org/10.1289/EHP3986.

Optimization of a Deep-Learning Method Based on the Classification of Images Generated by Parameterized Deep Snap a Novel Molecular-Image-Input Technique for Quantitative Structure-Activity Relationship (QSAR) Analysis

Numerous chemical compounds are distributed around the world and may affect the homeostasis of the endocrine system by disrupting the normal functions of hormone receptors. Although the risks associated with these compounds have been evaluated by acute toxicity testing in mammalian models, the chronic toxicity of many chemicals remains due to high cost of the compounds and the testing, etc. However, computational approaches may be promising alternatives and reduce these evaluations. Recently, deep learning (DL) has been shown to be promising prediction models with high accuracy for recognition of images, speech, signals, and videos since it greatly benefits from large datasets. Recently, a novel DL-based technique called DeepSnap was developed to conduct QSAR analysis using three-dimensional images of chemical structures. It can be used to predict the potential toxicity of many different chemicals to various receptors without extraction of descriptors. DeepSnap has been shown to have a very high capacity in tests using Tox21 quantitative qHTP datasets. Numerous parameters must be adjusted to use the DeepSnap method but they have not been optimized. In this study, the effects of these parameters on the performance of the DL prediction model were evaluated in terms of the loss in validation as an indicator for evaluating the performance of the DL using the toxicity information in the Tox21 qHTP database. The relations of the parameters of DeepSnap such as (1) number of molecules per SDF split into (2) zoom factor percentage, (3) atom size for van der waals percentage, (4) bond radius, (5) minimum bond distance, and (6) bond tolerance, with the validation loss following quadratic function curves, which suggests that optimal thresholds exist to attain the best performance with these prediction models. Using the parameter values set with the best performance, the prediction model of chemical compounds for CAR agonist was built using 64 images, at 105° angle, with AUC of 0.791. Thus, based on these parameters, the proposed DeepSnap-DL approach will be highly reliable and beneficial to establish models to assess the risk associated with various chemicals.

Deep learning for predicting toxicity of chemicals: a mini review

Humans and wildlife inhabit a world with panoply of natural and synthetic chemicals. Alarmingly, only a limited number of chemicals have undergone comprehensive toxicological evaluation due to limitations of traditional toxicity testing. High-throughput screening assays provide a higher-speed alternative for conventional toxicity testing. Advancement of high-throughput bioassay technology has greatly increased chemical toxicity data volumes in the past decade, pushing toxicology research into a "big data" era. However, traditional data analysis methods fail to effectively process large data volumes, presenting both a challenge and an opportunity for toxicologists. Deep learning, a machine learning method leveraging deep neural networks (DNNs), is a proven useful tool for building quantitative structure-activity relationship (QSAR) models for toxicity prediction utilizing these new large datasets. In this mini review, a brief technical background on DNNs is provided, and the current state of chemical toxicity prediction models built with DNNs is reviewed. In addition, relevant toxicity data sources are summarized, possible limitations are discussed, and perspectives on DNN utilization in chemical toxicity prediction are given.

Using Tox21 High-Throughput Screening Assays for the Evaluation of Botanical and Dietary Supplements

Introduction: Recent nationwide surveys found that natural products, including botanical dietary supplements, are used by ∼18% of adults. In many cases, there is a paucity of toxicological data available for these substances to allow for confident evaluations of product safety. The National Toxicology Program (NTP) has received numerous nominations from the public and federal agencies to study the toxicological effects of botanical dietary supplements. The NTP sought to evaluate the utility of in vitro quantitative high-throughput screening (qHTS) assays for toxicological assessment of botanical and dietary supplements. Materials and Methods: In brief, concentration-response assessments of 90 test substances, including 13 distinct botanical species, and individual purported active constituents were evaluated using a subset of the Tox21 qHTS testing panel. The screen included 20 different endpoints that covered a broad range of biologically relevant signaling pathways to detect test article effects upon endocrine activity, nuclear receptor signaling, stress response signaling, genotoxicity, and cell death signaling. Results and Discussion: Botanical dietary supplement extracts induced measurable and diverse activity. Elevated biological activity profiles were observed following treatments with individual chemical constituents relative to their associated botanical extract. The overall distribution of activity was comparable to activities exhibited by compounds present in the Tox21 10K chemical library. Conclusion: Botanical supplements did not exhibit minimal or idiosyncratic activities that would preclude the use of qHTS platforms as a feasible method to screen this class of compounds. However, there are still many considerations and further development required when attempting to use in vitro qHTS methods to characterize the safety profile of botanical/dietary supplements.

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

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

Multi-dimensional in vitro bioactivity profiling for grouping of glycol ethers

High-content screening data derived from physiologically-relevant in vitro models promise to improve confidence in data-integrative groupings for read-across in human health safety assessments. The biological data-based read-across concept is especially applicable to bioactive chemicals with defined mechanisms of toxicity; however, the challenge of data-derived groupings for chemicals that are associated with little or no bioactivity has not been explored. In this study, we apply a suite of organotypic and population-based in vitro models for comprehensive bioactivity profiling of twenty E-Series and P-Series glycol ethers, solvents with a broad variation in toxicity ranging from relatively non-toxic to reproductive and hematopoetic system toxicants. Both E-Series and P-Series glycol ethers elicited cytotoxicity only at high concentrations (mM range) in induced pluripotent stem cell-derived hepatocytes and cardiomyocytes. Population-variability assessment comprised a study of cytotoxicity in 94 human lymphoblast cell lines from 9 populations and revealed differences in inter-individual variability across glycol ethers, but did not indicate population-specific effects. Data derived from various phenotypic and transcriptomic assays revealed consistent bioactivity trends between both cardiomyocytes and hepatocytes, indicating a more universal, rather than cell-type specific mode-of-action for the tested glycol ethers in vitro. In vitro bioactivity-based similarity assessment using Toxicological Priority Index (ToxPi) showed that glycol ethers group according to their alcohol chain length, longer chains were associated with increased bioactivity. While overall in vitro bioactivity profiles did not correlate with in vivo toxicity data on glycol ethers, in vitro bioactivity of E-series glycol ethers were indicative of and correlated with in vivo irritation scores.

The Impact of Novel Assessment Methodologies in Toxicology on Green Chemistry and Chemical Alternatives

The field of experimental toxicology is rapidly advancing by incorporating novel techniques and methods that provide a much more granular view into the mechanisms of potential adverse effects of chemical exposures on human health. The data from various in vitro assays and computational models are useful not only for increasing confidence in hazard and risk decisions, but also are enabling better, faster and cheaper assessment of a greater number of compounds, mixtures, and complex products. This is of special value to the field of green chemistry where design of new materials or alternative uses of existing ones is driven, at least in part, by considerations of safety. This article reviews the state of the science and decision-making in scenarios when little to no data may be available to draw conclusions about which choice in green chemistry is "safer." It is clear that there is no "one size fits all" solution and multiple data streams need to be weighed in making a decision. Moreover, the overall level of familiarity of the decision-makers and scientists alike with new assessment methodologies, their validity, value and limitations is evolving. Thus, while the "impact" of the new developments in toxicology on the field of green chemistry is great already, it is premature to conclude that the data from new assessment methodologies have been widely accepted yet.

Computational in Vitro Toxicology Uncovers Chemical Structures Impairing Mitochondrial Membrane Potential

Technological advances in molecular biology have enabled high-throughput screening (HTS) of large chemical libraries. These approaches have provided valuable toxicity data for many physiological responses, including mitochondrial dysfunction. While several quantitative structure-activity relationship (QSAR) models have been developed for mitochondrial dysfunction, there remains a need to identify specific chemical features associated with this response. Thus, the objective of this study was to identify chemical structures associated with altered mitochondrial membrane potential (MMP). To achieve this, we developed computational models to examine the relationship between specific chemotypes (e.g., ToxPrints) and bioactivity in ToxCast/Tox21 HTS assays for altered MMP. The analysis revealed that the "bond:COH_alcohol_aromatic", "bond:COH_alcohol_aromatic_phenol", and "ring:aromatic_benzene" ToxPrints had the highest average correlations (phi coefficient) with ToxCast/Tox21 assay component endpoints for decreased MMP. These structures also constituted a "core" group of ToxPrints for decreased MMP in a force-directed network model and were the most important chemotypes in a random forest (RF) classification model for the "TOX21_MMP_ratio_down" assay component endpoint. Based on multiple lines of evidence, these structures, which are present in numerous chemicals (e.g., aromatic hydrocarbons, pesticides, and industrial chemicals) are likely involved in mitochondrial dysfunction. Because of the hierarchical structure of ToxPrints, these chemotypes were highly convergent and, when excluded from training data, had limited effects on the classification performance as related structures compensated for predictor loss. These results highlight the flexibility of the RF algorithm and ToxPrints for QSAR modeling, which is useful to identify chemicals affecting mitochondrial function.

Application of Benchmark Concentration (BMC) Analysis on Zebrafish Data: A New Perspective for Quantifying Toxicity in Alternative Animal Models

Over the past decade, the zebrafish is increasingly being used as a model to screen for chemical-mediated toxicities including developmental toxicity (DT) and neurotoxicity (NT). One of the major challenges is lack of harmonization in data analysis approaches, thereby posing difficulty in comparing findings across laboratories. To address this, we sought to establish a unified data analysis strategy for both DT and NT data, by adopting the benchmark concentration (BMC) analysis. There are two critical aspects in the BMC analysis: having a toxicity endpoint amenable for BMC and selecting a proper benchmark response (BMR) for the endpoint. For the former, in addition to the typical endpoints in NT assay (eg, hyper/hypo- response quantified by distance moved), we also used endpoints that assess the differences in movement patterns between chemical-treated embryos and control embryos. For the latter, we standardized the selection of BMR, which is analogous to minimum activity threshold, based on intrinsic response variations in the endpoint. When comparing our BMC results with a traditionally used LOAEL method (lowest-observed-adverse-effect level), we found high active compound concordance (100% for DT vs 74% for NT); generally, the BMC was more sensitive than LOAEL (no. of BMC more sensitive/no. of concordant active compounds, 43/50 for DT vs 16/26 for NT). Using the BMC with standardized toxicity endpoints and an appropriate BMR, we may now have a unified data-analysis approach to comparing results across different zebrafish datasets, for a better understanding of strengths and challenges when using the zebrafish as a screening tool.

Advancing internal exposure and physiologically-based toxicokinetic modeling for 21st-century risk assessments

Scientifically sound, risk-informed evaluation of chemicals is essential to protecting public health. Systematically leveraging information from exposure, toxicology, and epidemiology studies can provide a holistic understanding of how real-world exposure to chemicals may impact the health of populations, including sensitive and vulnerable individuals and life-stages. Increasingly, public health policy makers are employing toxicokinetic (TK) modeling tools to integrate these data streams and predict potential human health impact. Development of a suite of tools for predicting internal exposure, including physiologically-based toxicokinetic (PBTK) models, is being driven by needs to address large numbers of data-poor chemicals efficiently, translate bioactivity, and mechanistic information from new in vitro test systems, and integrate multiple lines of evidence to enable scientifically sound, risk-informed decisions. New modeling approaches are being designed "fit for purpose" to inform specific decision contexts, with applications ranging from rapid screening of hundreds of chemicals, to improved prediction of risks during sensitive stages of development. New data are being generated experimentally and computationally to support these models. Progress to meet the demand for internal exposure and PBTK modeling tools will require transparent publication of models and data to build credibility in results, as well as opportunities to partner with decision makers to evaluate and build confidence in use of these for improved decisions that promote safe use of chemicals.

High-throughput screening and environmental risk assessment: State of the science and emerging applications

In 2007 the United States National Research Council (NRC) published a vision for toxicity testing in the 21st century that emphasized the use of in vitro high-throughput screening (HTS) methods and predictive models as an alternative to in vivo animal testing. In the present study we examine the state of the science of HTS and the progress that has been made in implementing and expanding on the NRC vision, as well as challenges to implementation that remain. Overall, significant progress has been made with regard to the availability of HTS data, aggregation of chemical property and toxicity information into online databases, and the development of various models and frameworks to support extrapolation of HTS data. However, HTS data and associated predictive models have not yet been widely applied in risk assessment. Major barriers include the disconnect between the endpoints measured in HTS assays and the assessment endpoints considered in risk assessments as well as the rapid pace at which new tools and models are evolving in contrast with the slow pace at which regulatory structures change. Nonetheless, there are opportunities for environmental scientists and policymakers alike to take an impactful role in the ongoing development and implementation of the NRC vision. Six specific areas for scientific coordination and/or policy engagement are identified. Environ Toxicol Chem 2019;38:12-26. Published 2018 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.

VIVD: Virtual in vitro distribution model for the mechanistic prediction of intracellular concentrations of chemicals in in vitro toxicity assays

In vitro toxicity testing routinely uses nominal treatment concentrations as the driver for measured toxicity endpoints. However, test compounds can bind to the plastic of culture vessels or interact with culture media components, such as lipids and albumin. Additionally, volatile compounds may partition into the air above culture media. These processes reduce the free concentrations of compound to which cells are exposed. Models predicting the freely dissolved concentrations by accounting for these interactions have been published. However, these have only been applied to neutral compounds or assume no differential ionisation of test compounds between the media and cell cytoplasm. Herein, we describe an in vitro distribution model, based on the Fick-Nernst Planck equation accounting for differential compound ionisation in culture medium and intracellular water. The model considers permeability of ionised and unionised species and accounts for membrane potential in the partitioning of ionised moieties. By accounting for lipid and protein binding in culture medium, binding to cell culture plastic, air-partitioning, and lipid binding in the cell, the model can predict chemical concentrations (free and total) in medium and cells. The model can improve in vitro in vivo extrapolation of toxicity endpoint by determining intracellular concentrations for translation to in vivo.

Workflow for defining reference chemicals for assessing performance of in vitro assays

Instilling confidence in use of in vitro assays for predictive toxicology requires evaluation of assay performance. Performance is typically assessed using reference chemicals--compounds with defined activity against the test system target. However, developing reference chemical lists has historically been very resource-intensive. We developed a semi-automated process for selecting and annotating reference chemicals across many targets in a standardized format and demonstrate the workflow here. A series of required fields defines the potential reference chemical: the in vitro molecular target, pathway, or phenotype affected; and the chemical's mode (e.g. agonist, antagonist, inhibitor). Activity information was computationally extracted into a database from multiple public sources including non-curated scientific literature and curated chemical-biological databases, resulting in the identification of chemical activity in 2995 biological targets. Sample data from literature sources covering 54 molecular targets ranging from data-poor to data-rich was manually checked for accuracy. Precision rates were 82.7% from curated data sources and 39.5% from automated literature extraction. We applied the final reference chemical lists to evaluating performance of EPA's ToxCast program in vitro bioassays. The level of support, i.e. the number of independent reports in the database linking a chemical to a target, was found to strongly correlate with likelihood of positive results in the ToxCast assays, although individual assay performance had considerable variation. This overall approach allows rapid development of candidate reference chemical lists for a wide variety of targets that can facilitate performance evaluation of in vitro assays as a critical step in imparting confidence in alternative approaches.

Evidence for Cross Species Extrapolation of Mammalian-Based High-Throughput Screening Assay Results

High-throughput screening (HTS) and computational technologies have emerged as important tools for chemical hazard identification. The US Environmental Protection Agency (EPA) launched the Toxicity ForeCaster (ToxCast) Program, which has screened thousands of chemicals in hundreds of mammalian-based HTS assays for biological activity. The data are being used to prioritize toxicity testing on those chemicals likely to lead to adverse effects. To use HTS assays in predicting hazard to both humans and wildlife, it is necessary to understand how broadly these data may be extrapolated across species. The US EPA Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS; https://seqapass.epa.gov/seqapass/ ) tool was used to assess conservation of the 484 protein targets represented in the suite of ToxCast assays and other HTS assays. To demonstrate the utility of the SeqAPASS data for guiding extrapolation, case studies were developed which focused on targets of interest to the US Endocrine Disruptor Screening Program and the Organisation for Economic Cooperation and Development. These case studies provide a line of evidence for conservation of endocrine targets across vertebrate species, with few exceptions, and demonstrate the utility of SeqAPASS for defining the taxonomic domain of applicability for HTS results and identifying organisms for suitable follow-up toxicity tests.

PotentialNet for Molecular Property Prediction

The arc of drug discovery entails a multiparameter optimization problem spanning vast length scales. The key parameters range from solubility (angstroms) to protein-ligand binding (nanometers) to in vivo toxicity (meters). Through feature learning-instead of feature engineering-deep neural networks promise to outperform both traditional physics-based and knowledge-based machine learning models for predicting molecular properties pertinent to drug discovery. To this end, we present the PotentialNet family of graph convolutions. These models are specifically designed for and achieve state-of-the-art performance for protein-ligand binding affinity. We further validate these deep neural networks by setting new standards of performance in several ligand-based tasks. In parallel, we introduce a new metric, the Regression Enrichment Factor EF_χ ^(R), to measure the early enrichment of computational models for chemical data. Finally, we introduce a cross-validation strategy based on structural homology clustering that can more accurately measure model generalizability, which crucially distinguishes the aims of machine learning for drug discovery from standard machine learning tasks.

Target-specific toxicity knowledgebase (TsTKb): a novel toolkit for in silico predictive toxicology

As the number of man-made chemicals increases at an unprecedented pace, efforts of quickly screening and accurately evaluating their potential adverse biological effects have been hampered by prohibitively high costs of in vivo/vitro toxicity testing. While it is unrealistic and unnecessary to test every uncharacterized chemical, it remains a major challenge to develop alternative in silico tools with high reliability and precision in toxicity prediction. To address this urgent need, we have developed a novel mode-of-action-guided, molecular modeling-based, and machine learning-enabled modeling approach for in silico chemical toxicity prediction. Here we introduce the core element of this approach, Target-specific Toxicity Knowledgebase (TsTKb), which consists of two main components: Chemical Mode of Action (ChemMoA) database and a suite of prediction model libraries.

Characterizing sources of variability in zebrafish embryo screening protocols

There is a need for fast, efficient, and cost-effective hazard identification and characterization of chemical hazards. This need is generating increased interest in the use of zebrafish embryos as both a screening tool and an alternative to mammalian test methods. A Collaborative Workshop on Aquatic Models and 21st Century Toxicology identified the lack of appropriate and consistent testing protocols as a challenge to the broader application of the zebrafish embryo model. The National Toxicology Program established the Systematic Evaluation of the Application of Zebrafish in Toxicology (SEAZIT) initiative to address the lack of consistent testing guidelines and identify sources of variability for zebrafish-based assays. This report summarizes initial SEAZIT information-gathering efforts. Investigators in academic, government, and industry laboratories that routinely use zebrafish embryos for chemical toxicity testing were asked about their husbandry practices and standard protocols. Information was collected about protocol components including zebrafish strains, feed, system water, disease surveillance, embryo exposure conditions, and endpoints. Literature was reviewed to assess issues raised by the investigators. Interviews revealed substantial variability across design parameters, data collected, and analysis procedures. The presence of the chorion and renewal of exposure media (static versus static-renewal) were identified as design parameters that could potentially influence study outcomes and should be investigated further with studies to determine chemical uptake from treatment solution into embryos. The information gathered in this effort provides a basis for future SEAZIT activities to promote more consistent practices among researchers using zebrafish embryos for toxicity evaluation.

ChemDIS 2: an update of chemical-disease inference system

Computational inference of affected functions, pathways and diseases for chemicals could largely accelerate the evaluation of potential effects of chemical exposure on human beings. Previously, we have developed a ChemDIS system utilizing information of interacting targets for chemical-disease inference. With the target information, testable hypotheses can be generated for experimental validation. In this work, we present an update of ChemDIS 2 system featured with more updated datasets and several new functions, including (i) custom enrichment analysis function for single omics data; (ii) multi-omics analysis function for joint analysis of multi-omics data; (iii) mixture analysis function for the identification of interaction and overall effects; (iv) web application programming interface (API) for programmed access to ChemDIS 2. The updated ChemDIS 2 system capable of analyzing more than 430 000 chemicals is expected to be useful for both drug development and risk assessment of environmental chemicals.

Database URL: ChemDIS 2 is freely accessible via https://cwtung.kmu.edu.tw/chemdis

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.

A Mechanistic Framework for Integrating Chemical Structure and High-Throughput Screening Results to Improve Toxicity Predictions

Adverse Outcome Pathways (AOPs) establish a connection between a molecular initiating event (MIE) and an adverse outcome. Detailed understanding of the MIE provides the ideal data for determining chemical properties required to elicit the MIE. This study utilized high-throughput screening data from the ToxCast program, coupled with chemical structural information, to generate chemical clusters using three similarity methods pertaining to nine MIEs within an AOP network for hepatic steatosis. Three case studies demonstrate the utility of the mechanistic information held by the MIE for integrating biological and chemical data. Evaluation of the chemical clusters activating the glucocorticoid receptor identified activity differences in chemicals within a cluster. Comparison of the estrogen receptor results with previous work showed that bioactivity data and structural alerts can be combined to improve predictions in a customizable way where bioactivity data are limited. The aryl hydrocarbon receptor (AHR) highlighted that while structural data can be used to offset limited data for new screening efforts, not all ToxCast targets have sufficient data to define robust chemical clusters. In this context, an alternative to additional receptor assays is proposed where assays for proximal key events downstream of AHR activation could be used to enhance confidence in active calls. These case studies illustrate how the AOP framework can support an iterative process whereby in vitro toxicity testing and chemical structure can be combined to improve toxicity predictions. In vitro assays can inform the development of structural alerts linking chemical structure to toxicity. Consequently, structurally related chemical groups can facilitate identification of assays that would be informative for a specific MIE. Together, these activities form a virtuous cycle where the mechanistic basis for the in vitro results and the breadth of the structural alerts continually improve over time to better predict activity of chemicals for which limited toxicity data exist.

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.

A systems biology approach to predictive developmental neurotoxicity of a larvicide used in the prevention of Zika virus transmission

The need to prevent developmental brain disorders has led to an increased interest in efficient neurotoxicity testing. When an epidemic of microcephaly occurred in Brazil, Zika virus infection was soon identified as the likely culprit. However, the pathogenesis appeared to be complex, and a larvicide used to control mosquitoes responsible for transmission of the virus was soon suggested as an important causative factor. Yet, it is challenging to identify relevant and efficient tests that are also in line with ethical research defined by the 3Rs rule (Replacement, Reduction and Refinement). Especially in an acute situation like the microcephaly epidemic, where little toxicity documentation is available, new and innovative alternative methods, whether in vitro or in silico, must be considered. We have developed a network-based model using an integrative systems biology approach to explore the potential developmental neurotoxicity, and we applied this method to examine the larvicide pyriproxyfen widely used in the prevention of Zika virus transmission. Our computational model covered a wide range of possible pathways providing mechanistic hypotheses between pyriproxyfen and neurological disorders via protein complexes, thus adding to the plausibility of pyriproxyfen neurotoxicity. Although providing only tentative evidence and comparisons with retinoic acid, our computational systems biology approach is rapid and inexpensive. The case study of pyriproxyfen illustrates its usefulness as an initial or screening step in the assessment of toxicity potentials of chemicals with incompletely known toxic properties.

Evaluation and Optimization of Pharmacokinetic Models for in Vitro to in Vivo Extrapolation of Estrogenic Activity for Environmental Chemicals

BACKGROUND

To effectively incorporate in vitro data into regulatory use, confidence must be established in the quantitative extrapolation of in vitro activity to relevant end points in animals or humans.

OBJECTIVE

Our goal was to evaluate and optimize in vitro to in vivo extrapolation (IVIVE) approaches using in vitro estrogen receptor (ER) activity to predict estrogenic effects measured in rodent uterotrophic studies.

METHODS

We evaluated three pharmacokinetic (PK) models with varying complexities to extrapolate in vitro to in vivo dosimetry for a group of 29 ER agonists, using data from validated in vitro [U.S. Environmental Protection Agency (U.S. EPA) ToxCast™ ER model] and in vivo (uterotrophic) methods. In vitro activity values were adjusted using mass-balance equations to estimate intracellular exposure via an enrichment factor (EF), and steady-state model calculations were adjusted using fraction of unbound chemical in the plasma ([Formula: see text]) to approximate bioavailability. Accuracy of each model-adjustment combination was assessed by comparing model predictions with lowest effect levels (LELs) from guideline uterotrophic studies.

RESULTS

We found little difference in model predictive performance based on complexity or route-specific modifications. Simple adjustments, applied to account for in vitro intracellular exposure (EF) or chemical bioavailability ([Formula: see text]), resulted in significant improvements in the predictive performance of all models.

CONCLUSION

Computational IVIVE approaches accurately estimate chemical exposure levels that elicit positive responses in the rodent uterotrophic bioassay. The simplest model had the best overall performance for predicting both oral (PPK_EF) and injection (PPK_[Formula: see text]) LELs from guideline uterotrophic studies, is freely available, and can be parameterized entirely using freely available in silico tools. https://doi.org/10.1289/EHP1655.

Machine Learning Based Toxicity Prediction: From Chemical Structural Description to Transcriptome Analysis

Toxicity prediction is very important to public health. Among its many applications, toxicity prediction is essential to reduce the cost and labor of a drug's preclinical and clinical trials, because a lot of drug evaluations (cellular, animal, and clinical) can be spared due to the predicted toxicity. In the era of Big Data and artificial intelligence, toxicity prediction can benefit from machine learning, which has been widely used in many fields such as natural language processing, speech recognition, image recognition, computational chemistry, and bioinformatics, with excellent performance. In this article, we review machine learning methods that have been applied to toxicity prediction, including deep learning, random forests, k-nearest neighbors, and support vector machines. We also discuss the input parameter to the machine learning algorithm, especially its shift from chemical structural description only to that combined with human transcriptome data analysis, which can greatly enhance prediction accuracy.

High-Content Assay Multiplexing for Muscle Toxicity Screening in Human-Induced Pluripotent Stem Cell-Derived Skeletal Myoblasts

Skeletal muscle-associated toxicity is an underresearched area in the field of high-throughput toxicity screening; hence, the potential adverse effects of drugs and chemicals on skeletal muscle are largely unknown. Novel organotypic microphysiological in vitro models are being developed to replicate the contractile function of skeletal muscle; however, the throughput and a need for specialized equipment may limit the utility of these tissue chip models for screening. In addition, recent developments in stem cell biology have resulted in the generation of induced pluripotent stem cell (iPSC)-derived skeletal myoblasts that enable high-throughput in vitro screening. This study set out to develop a high-throughput multiplexed assay using iPSC-derived skeletal myoblasts that can be used as a first-pass screen to assess the potential for chemicals to affect skeletal muscle. We found that cytotoxicity and cytoskeletal integrity are most useful and reproducible assays for the skeletal myoblasts when evaluating overall cellular health or gauging disruptions in actin polymerization following 24 h of exposure. Both assays are based on high-content imaging and quantitative image processing to derive quantitative phenotypes. Both assays showed good to excellent assay robustness and reproducibility measured by interplate and interday replicability, coefficients of variation of negative controls, and Z'-factors for positive control chemicals. Concentration response assessment of muscle-related toxicants showed specificity of the observed effects compared to the general cytotoxicity. Overall, this study establishes a high-throughput multiplexed assay using skeletal myoblasts that may be used for screening and prioritization of chemicals for more complex tissue chip-based and in vivo evaluation.

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.

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.

Infer the in vivo point of departure with ToxCast in vitro assay data using a robust learning approach

The development and application of high throughput in vitro assays is an important development for risk assessment in the twenty-first century. However, there are still significant challenges to incorporate in vitro assays into routine toxicity testing practices. In this paper, a robust learning approach was developed to infer the in vivo point of departure (POD) with in vitro assay data from ToxCast and Tox21 projects. Assay data from ToxCast and Tox21 projects were utilized to derive the in vitro PODs for several hundred chemicals. These were combined with in vivo PODs from ToxRefDB regarding the rat and mouse liver to build a high-dimensional robust regression model. This approach separates the chemicals into a majority, well-predicted set; and a minority, outlier set. Salient relationships can then be learned from the data. For both mouse and rat liver PODs, over 93% of chemicals have inferred values from in vitro PODs that are within ± 1 of the in vivo PODs on the log₁₀ scale (the target learning region, or TLR) and R² of 0.80 (rats) and 0.78 (mice) for these chemicals. This is comparable with extrapolation between related species (mouse and rat), which has 93% chemicals within the TLR and the R² being 0.78. Chemicals in the outlier set tend to also have more biologically variable characteristics. With the continued accumulation of high throughput data for a wide range of chemicals, predictive modeling can provide a valuable complement for adverse outcome pathway based approach in risk assessment.

Comprehensive Analyses and Prioritization of Tox21 10K Chemicals Affecting Mitochondrial Function by in-Depth Mechanistic Studies

BACKGROUND

A central challenge in toxicity testing is the large number of chemicals in commerce that lack toxicological assessment. In response, the Tox21 program is re-focusing toxicity testing from animal studies to less expensive and higher throughput in vitro methods using target/pathway-specific, mechanism-driven assays.

OBJECTIVES

Our objective was to use an in-depth mechanistic study approach to prioritize and characterize the chemicals affecting mitochondrial function.

METHODS

We used a tiered testing approach to prioritize for more extensive testing 622 compounds identified from a primary, quantitative high-throughput screen of 8,300 unique small molecules, including drugs and industrial chemicals, as potential mitochondrial toxicants by their ability to significantly decrease the mitochondrial membrane potential (MMP). Based on results from secondary MMP assays in HepG2 cells and rat hepatocytes, 34 compounds were selected for testing in tertiary assays that included formation of reactive oxygen species (ROS), upregulation of p53 and nuclear erythroid 2-related factor 2/antioxidant response element (Nrf2/ARE), mitochondrial oxygen consumption, cellular Parkin translocation, and larval development and ATP status in the nematode Caenorhabditis elegans.

RESULTS

A group of known mitochondrial complex inhibitors (e.g., rotenone) and uncouplers (e.g., chlorfenapyr), as well as potential novel complex inhibitors and uncouplers, were detected. From this study, we identified four not well-characterized potential mitochondrial toxicants (lasalocid, picoxystrobin, pinacyanol, and triclocarban) that merit additional in vivo characterization.

CONCLUSIONS

The tier-based approach for identifying and mechanistically characterizing mitochondrial toxicants can potentially reduce animal use in toxicological testing. https://doi.org/10.1289/EHP2589.

Confirmation of high-throughput screening data and novel mechanistic insights into VDR-xenobiotic interactions by orthogonal assays

High throughput screening (HTS) programs have demonstrated that the Vitamin D receptor (VDR) is activated and/or antagonized by a wide range of structurally diverse chemicals. In this study, we examined the Tox21 qHTS data set generated against VDR for reproducibility and concordance and elucidated functional insights into VDR-xenobiotic interactions. Twenty-one potential VDR agonists and 19 VDR antagonists were identified from a subset of >400 compounds with putative VDR activity and examined for VDR functionality utilizing select orthogonal assays. Transient transactivation assay (TT) using a human VDR plasmid and Cyp24 luciferase reporter construct revealed 20/21 active VDR agonists and 18/19 active VDR antagonists. Mammalian-2-hybrid assay (M2H) was then used to evaluate VDR interactions with co-activators and co-regulators. With the exception of a select few compounds, VDR agonists exhibited significant recruitment of co-regulators and co-activators whereas antagonists exhibited considerable attenuation of recruitment by VDR. A unique set of compounds exhibiting synergistic activity in antagonist mode and no activity in agonist mode was identified. Cheminformatics modeling of VDR-ligand interactions were conducted and revealed selective ligand VDR interaction. Overall, data emphasizes the molecular complexity of ligand-mediated interactions with VDR and suggest that VDR transactivation may be a target site of action for diverse xenobiotics.

Biocompatibility and Cytotoxic Evaluation of New Sorbent Cartridges for Blood Hemoperfusion

BACKGROUND/AIMS

The use of adsorption cartridges for hemoperfusion (HP) is rapidly evolving. For these devices, the potential induced cytotoxicity is an important issue. The aim of this study was to investigate potential in vitro cytotoxic effects of different sorbent cartridges, HA130, HA230, HA330, HA380 (Jafron, China), on U937 monocytes.

METHODS

Monocytes were exposed to the sorbent material in static and dynamic manners. In static test, cell medium samples were collected after 24 h of incubation in the cartridges. In dynamic test, HP modality has been carried out and samples at 30, 60, 90, and 120 min were collected.

RESULTS

Compared to control samples, there was no evidence of increased necrosis or apoptosis in monocytes exposed to the cartridges both in the static and dynamic tests.

CONCLUSION

Our in vitro testing suggests that HA cartridges carry an optimal level of biocompatibility and their use in HP is not associated with adverse reactions or signs of cytotoxicity.

Differentiating Pathway-Specific From Nonspecific Effects in High-Throughput Toxicity Data: A Foundation for Prioritizing Adverse Outcome Pathway Development

The U.S. Environmental Protection Agency's ToxCast program has screened thousands of chemicals for biological activity, primarily using high-throughput in vitro bioassays. Adverse outcome pathways (AOPs) offer a means to link pathway-specific biological activities with potential apical effects relevant to risk assessors. Thus, efforts are underway to develop AOPs relevant to pathway-specific perturbations detected in ToxCast assays. Previous work identified a "cytotoxic burst" (CTB) phenomenon wherein large numbers of the ToxCast assays begin to respond at or near test chemical concentrations that elicit cytotoxicity, and a statistical approach to defining the bounds of the CTB was developed. To focus AOP development on the molecular targets corresponding to ToxCast assays indicating pathway-specific effects, we conducted a meta-analysis to identify which assays most frequently respond at concentrations below the CTB. A preliminary list of potentially important, target-specific assays was determined by ranking assays by the fraction of chemical hits below the CTB compared with the number of chemicals tested. Additional priority assays were identified using a diagnostic-odds-ratio approach which gives greater ranking to assays with high specificity but low responsivity. Combined, the two prioritization methods identified several novel targets (e.g., peripheral benzodiazepine and progesterone receptors) to prioritize for AOP development, and affirmed the importance of a number of existing AOPs aligned with ToxCast targets (e.g., thyroperoxidase, estrogen receptor, aromatase). The prioritization approaches did not appear to be influenced by inter-assay differences in chemical bioavailability. Furthermore, the outcomes were robust based on a variety of different parameters used to define the CTB.

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.

Expediting toxicity testing with increased precision, predictive power, and clinical utility

Federal Government management of health risks associated with the use of therapeutics and unintended environmental chemical exposures must be expedited to meet public health needs. Although US agencies initiated the Tox21 strategy over a decade ago to expedite toxicity testing and improve the reliability of risk assessments, recent status reports indicate that achieving its goals is still decades away. Emerging technologies create an opportunity to both expedite toxicity testing and improve its predictive power. The way forward may be an augmentation of the strategy aimed at enhancing the resolution and scope of Tox21 and exploring the adaptability of real-time chemical sensor, digital imaging, and other technologies to toxicity testing. Among the anticipated returns on the associated investment would likely be enhanced accuracy in prediction, reductions in the time needed to conduct hazard identifications and toxicity assessments, and an overall increase in the precision and reliability of the risk assessment process. This in turn expedites risk management decisions and reduces scientific uncertainty and the need to incorporate margins of safety that can add cost without necessarily returning improved health protection.

Application of the key characteristics of carcinogens in cancer hazard identification

Smith et al. (Env. Health Perspect. 124: 713, 2016) identified 10 key characteristics (KCs), one or more of which are commonly exhibited by established human carcinogens. The KCs reflect the properties of a cancer-causing agent, such as 'is genotoxic,' 'is immunosuppressive' or 'modulates receptor-mediated effects,' and are distinct from the hallmarks of cancer, which are the properties of tumors. To assess feasibility and limitations of applying the KCs to diverse agents, methods and results of mechanistic data evaluations were compiled from eight recent IARC Monograph meetings. A systematic search, screening and evaluation procedure identified a broad literature encompassing multiple KCs for most (12/16) IARC Group 1 or 2A carcinogens identified in these meetings. Five carcinogens are genotoxic and induce oxidative stress, of which pentachlorophenol, hydrazine and malathion also showed additional KCs. Four others, including welding fumes, are immunosuppressive. The overall evaluation was upgraded to Group 2A based on mechanistic data for only two agents, tetrabromobisphenol A and tetrachloroazobenzene. Both carcinogens modulate receptor-mediated effects in combination with other KCs. Fewer studies were identified for Group 2B or 3 agents, with the vast majority (17/18) showing only one or no KCs. Thus, an objective approach to identify and evaluate mechanistic studies pertinent to cancer revealed strong evidence for multiple KCs for most Group 1 or 2A carcinogens but also identified opportunities for improvement. Further development and mapping of toxicological and biomarker endpoints and pathways relevant to the KCs can advance the systematic search and evaluation of mechanistic data in carcinogen hazard identification.

High-Throughput H295R Steroidogenesis Assay: Utility as an Alternative and a Statistical Approach to Characterize Effects on Steroidogenesis

The U.S. Environmental Protection Agency Endocrine Disruptor Screening Program and the Organization for Economic Co-operation and Development (OECD) have used the human adrenocarcinoma (H295R) cell-based assay to predict chemical perturbation of androgen and estrogen production. Recently, a high-throughput H295R (HT-H295R) assay was developed as part of the ToxCast program that includes measurement of 11 hormones, including progestagens, corticosteroids, androgens, and estrogens. To date, 2012 chemicals have been screened at 1 concentration; of these, 656 chemicals have been screened in concentration-response. The objectives of this work were to: (1) develop an integrated analysis of chemical-mediated effects on steroidogenesis in the HT-H295R assay and (2) evaluate whether the HT-H295R assay predicts estrogen and androgen production specifically via comparison with the OECD-validated H295R assay. To support application of HT-H295R assay data to weight-of-evidence and prioritization tasks, a single numeric value based on Mahalanobis distances was computed for 654 chemicals to indicate the magnitude of effects on the synthesis of 11 hormones. The maximum mean Mahalanobis distance (maxmMd) values were high for strong modulators (prochloraz, mifepristone) and lower for moderate modulators (atrazine, molinate). Twenty-five of 28 reference chemicals used for OECD validation were screened in the HT-H295R assay, and produced qualitatively similar results, with accuracies of 0.90/0.75 and 0.81/0.91 for increased/decreased testosterone and estradiol production, respectively. The HT-H295R assay provides robust information regarding estrogen and androgen production, as well as additional hormones. The maxmMd from this integrated analysis may provide a data-driven approach to prioritizing lists of chemicals for putative effects on steroidogenesis.

Transcriptomic and phenotypic profiling in developing zebrafish exposed to thyroid hormone receptor agonists

There continues to be a need to develop in vivo high-throughput screening (HTS) and computational methods to screen chemicals for interaction with the estrogen, androgen, and thyroid pathways and as complements to in vitro HTS assays. This study explored the utility of an embryonic zebrafish HTS approach to identify and classify endocrine bioactivity using phenotypically-anchored transcriptome profiling. Transcriptome analysis was conducted on zebrafish embryos exposed to 25 estrogen-, androgen-, or thyroid-active chemicals at concentrations that elicited adverse malformations or mortality at 120 h post-fertilization in 80% of animals exposed. Analysis of the top 1000 significant differentially expressed transcripts and developmental toxicity profiles across all treatments identified a unique transcriptional and phenotypic signature for thyroid hormone receptor agonists. This unique signature has the potential to be used as a tiered in vivo HTS and may aid in identifying chemicals that interact with the thyroid hormone receptor.

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.