tcpl: the ToxCast pipeline for high-throughput screening data.

Motivation: Large high-throughput screening (HTS) efforts are widely used in drug development and chemical toxicity screening. Wide use and integration of these data can benefit from an efficient, transparent and reproducible data pipeline. Summary: The tcpl R package and its associated MySQL database provide a generalized platform for efficiently storing, normalizing and dose-response modeling of large high-throughput and high-content chemical screening data. The novel dose-response modeling algorithm has been tested against millions of diverse dose-response series, and robustly fits data with outliers and cytotoxicity-related signal loss. Availability and Implementation: tcpl is freely available on the Comprehensive R Archive Network under the GPL-2 license. Contact: martin.matt@epa.gov.

Screening the ToxCast phase II libraries for alterations in network function using cortical neurons grown on multi-well microelectrode array (mwMEA) plates.

Methods are needed for rapid screening of environmental compounds for neurotoxicity, particularly ones that assess function. To demonstrate the utility of microelectrode array (MEA)-based approaches as a rapid neurotoxicity screening tool, 1055 chemicals from EPA's phase II ToxCast library were evaluated for effects on neural function and cell health. Primary cortical networks were grown on multi-well microelectrode array (mwMEA) plates. On day in vitro 13, baseline activity (40 min) was recorded prior to exposure to each compound (40 µM). Changes in spontaneous network activity [mean firing rate (MFR)] and cell viability (lactate dehydrogenase and CellTiter Blue) were assessed within the same well following compound exposure. Following exposure, 326 compounds altered (increased or decreased) normalized MFR beyond hit thresholds based on 2× the median absolute deviation of DMSO-treated wells. Pharmaceuticals, pesticides, fungicides, chemical intermediates, and herbicides accounted for 86% of the hits. Further, changes in MFR occurred in the absence of cytotoxicity, as only eight compounds decreased cell viability. ToxPrint chemotype analysis identified several structural domains (e.g., biphenyls and alkyl phenols) significantly enriched with MEA actives relative to the total test set. The top 5 enriched ToxPrint chemotypes were represented in 26% of the MEA hits, whereas the top 11 ToxPrints were represented in 34% of MEA hits. These results demonstrate that large-scale functional screening using neural networks on MEAs can fill a critical gap in assessment of neurotoxicity potential in ToxCast assay results. Further, a data-mining approach identified ToxPrint chemotypes enriched in the MEA-hit subset, which define initial structure-activity relationship inferences, establish potential mechanistic associations to other ToxCast assay endpoints, and provide working hypotheses for future studies.

An "EAR" on Environmental Surveillance and Monitoring: A Case Study on the Use of Exposure-Activity Ratios (EARs) to Prioritize Sites, Chemicals, and Bioactivities of Concern in Great Lakes Waters.

Current environmental monitoring approaches focus primarily on chemical occurrence. However, based on concentration alone, it can be difficult to identify which compounds may be of toxicological concern and should be prioritized for further monitoring, in-depth testing, or management. This can be problematic because toxicological characterization is lacking for many emerging contaminants. New sources of high-throughput screening (HTS) data, such as the ToxCast database, which contains information for over 9000 compounds screened through up to 1100 bioassays, are now available. Integrated analysis of chemical occurrence data with HTS data offers new opportunities to prioritize chemicals, sites, or biological effects for further investigation based on concentrations detected in the environment linked to relative potencies in pathway-based bioassays. As a case study, chemical occurrence data from a 2012 study in the Great Lakes Basin along with the ToxCast effects database were used to calculate exposure-activity ratios (EARs) as a prioritization tool. Technical considerations of data processing and use of the ToxCast database are presented and discussed. EAR prioritization identified multiple sites, biological pathways, and chemicals that warrant further investigation. Prioritized bioactivities from the EAR analysis were linked to discrete adverse outcome pathways to identify potential adverse outcomes and biomarkers for use in subsequent monitoring efforts.

Retrospective mining of toxicology data to discover multispecies and chemical class effects: Anemia as a case study.

Predictive toxicity models rely on large amounts of accurate in vivo data. Here, we analyze the quality of in vivo data from the U.S. EPA Toxicity Reference Database (ToxRefDB), using chemical-induced anemia as an example. Considerations include variation in experimental conditions, changes in terminology over time, distinguishing negative from missing results, observer and diagnostic bias, and data transcription errors. Within ToxRefDB, we use hematological data on 658 chemicals tested in one or more of 1738 studies (subchronic rat or chronic rat, mouse, or dog). Anemia was reported most frequently in the rat subchronic studies, followed by chronic studies in dog, rat, and then mouse. Concordance between studies for a positive finding of anemia (same chemical, different laboratories) ranged from 90% (rat subchronic predicting rat chronic) to 40% (mouse chronic predicting rat chronic). Concordance increased with manual curation by 20% on average. We identified 49 chemicals that showed an anemia phenotype in at least two species. These included 14 aniline moiety-containing compounds that were further analyzed for their potential to be metabolically transformed into substituted anilines, which are known anemia-causing chemicals. This analysis should help inform future use of in vivo databases for model development.

ToxCast Chemical Landscape: Paving the Road to 21st Century Toxicology.

The U.S. Environmental Protection Agency's (EPA) ToxCast program is testing a large library of Agency-relevant chemicals using in vitro high-throughput screening (HTS) approaches to support the development of improved toxicity prediction models. Launched in 2007, Phase I of the program screened 310 chemicals, mostly pesticides, across hundreds of ToxCast assay end points. In Phase II, the ToxCast library was expanded to 1878 chemicals, culminating in the public release of screening data at the end of 2013. Subsequent expansion in Phase III has resulted in more than 3800 chemicals actively undergoing ToxCast screening, 96% of which are also being screened in the multi-Agency Tox21 project. The chemical library unpinning these efforts plays a central role in defining the scope and potential application of ToxCast HTS results. The history of the phased construction of EPA's ToxCast library is reviewed, followed by a survey of the library contents from several different vantage points. CAS Registry Numbers are used to assess ToxCast library coverage of important toxicity, regulatory, and exposure inventories. Structure-based representations of ToxCast chemicals are then used to compute physicochemical properties, substructural features, and structural alerts for toxicity and biotransformation. Cheminformatics approaches using these varied representations are applied to defining the boundaries of HTS testability, evaluating chemical diversity, and comparing the ToxCast library to potential target application inventories, such as used in EPA's Endocrine Disruption Screening Program (EDSP). Through several examples, the ToxCast chemical library is demonstrated to provide comprehensive coverage of the knowledge domains and target inventories of potential interest to EPA. Furthermore, the varied representations and approaches presented here define local chemistry domains potentially worthy of further investigation (e.g., not currently covered in the testing library or defined by toxicity "alerts") to strategically support data mining and predictive toxicology modeling moving forward.

Predicting hepatotoxicity using ToxCast in vitro bioactivity and chemical structure.

The U.S. Tox21 and EPA ToxCast program screen thousands of environmental chemicals for bioactivity using hundreds of high-throughput in vitro assays to build predictive models of toxicity. We represented chemicals based on bioactivity and chemical structure descriptors, then used supervised machine learning to predict in vivo hepatotoxic effects. A set of 677 chemicals was represented by 711 in vitro bioactivity descriptors (from ToxCast assays), 4,376 chemical structure descriptors (from QikProp, OpenBabel, PaDEL, and PubChem), and three hepatotoxicity categories (from animal studies). Hepatotoxicants were defined by rat liver histopathology observed after chronic chemical testing and grouped into hypertrophy (161), injury (101) and proliferative lesions (99). Classifiers were built using six machine learning algorithms: linear discriminant analysis (LDA), Naïve Bayes (NB), support vector machines (SVM), classification and regression trees (CART), k-nearest neighbors (KNN), and an ensemble of these classifiers (ENSMB). Classifiers of hepatotoxicity were built using chemical structure descriptors, ToxCast bioactivity descriptors, and hybrid descriptors. Predictive performance was evaluated using 10-fold cross-validation testing and in-loop, filter-based, feature subset selection. Hybrid classifiers had the best balanced accuracy for predicting hypertrophy (0.84 ± 0.08), injury (0.80 ± 0.09), and proliferative lesions (0.80 ± 0.10). Though chemical and bioactivity classifiers had a similar balanced accuracy, the former were more sensitive, and the latter were more specific. CART, ENSMB, and SVM classifiers performed the best, and nuclear receptor activation and mitochondrial functions were frequently found in highly predictive classifiers of hepatotoxicity. ToxCast and ToxRefDB provide the largest and richest publicly available data sets for mining linkages between the in vitro bioactivity of environmental chemicals and their adverse histopathological outcomes. Our findings demonstrate the utility of high-throughput assays for characterizing rodent hepatotoxicants, the benefit of using hybrid representations that integrate bioactivity and chemical structure, and the need for objective evaluation of classification performance.

Profiling 976 ToxCast chemicals across 331 enzymatic and receptor signaling assays.

Understanding potential health risks is a significant challenge due to the large numbers of diverse chemicals with poorly characterized exposures and mechanisms of toxicities. The present study analyzes 976 chemicals (including failed pharmaceuticals, alternative plasticizers, food additives, and pesticides) in Phases I and II of the U.S. EPA's ToxCast project across 331 cell-free enzymatic and ligand-binding high-throughput screening (HTS) assays. Half-maximal activity concentrations (AC50) were identified for 729 chemicals in 256 assays (7,135 chemical-assay pairs). Some of the most commonly affected assays were CYPs (CYP2C9 and CYP2C19), transporters (mitochondrial TSPO, norepinephrine, and dopaminergic), and GPCRs (aminergic). Heavy metals, surfactants, and dithiocarbamate fungicides showed promiscuous but distinctly different patterns of activity, whereas many of the pharmaceutical compounds showed promiscuous activity across GPCRs. Literature analysis confirmed >50% of the activities for the most potent chemical-assay pairs (54) but also revealed 10 missed interactions. Twenty-two chemicals with known estrogenic activity were correctly identified for the majority (77%), missing only the weaker interactions. In many cases, novel findings for previously unreported chemical-target combinations clustered with known chemical-target interactions. Results from this large inventory of chemical-biological interactions can inform read-across methods as well as link potential targets to molecular initiating events in adverse outcome pathways for diverse toxicities.

Real-time growth kinetics measuring hormone mimicry for ToxCast chemicals in T-47D human ductal carcinoma cells.

High-throughput screening (HTS) assays capable of profiling thousands of environmentally relevant chemicals for in vitro biological activity provide useful information on the potential for disrupting endocrine pathways. Disruption of the estrogen signaling pathway has been implicated in a variety of adverse health effects including impaired development, reproduction, and carcinogenesis. The estrogen-responsive human mammary ductal carcinoma cell line T-47D was exposed to 1815 ToxCast chemicals comprising pesticides, industrial chemicals, pharmaceuticals, personal care products, cosmetics, food ingredients, and other chemicals with known or suspected human exposure potential. Cell growth kinetics were evaluated using real-time cell electronic sensing. T-47D cells were exposed to eight concentrations (0.006-100 µM), and measurements of cellular impedance were repeatedly recorded for 105 h. Chemical effects were evaluated based on potency (concentration at which response occurs) and efficacy (extent of response). A linear growth response was observed in response to prototypical estrogen receptor agonists (17ß-estradiol, genistein, bisphenol A, nonylphenol, and 4-tert-octylphenol). Several compounds, including bisphenol A and genistein, induced cell growth comparable in efficacy to that of 17ß-estradiol, but with decreased potency. Progestins, androgens, and corticosteroids invoked a biphasic growth response indicative of changes in cell number or cell morphology. Results from this cell growth assay were compared with results from additional estrogen receptor (ER) binding and transactivation assays. Chemicals detected as active in both the cell growth and ER receptor binding assays demonstrated potencies highly correlated with two ER transactivation assays (r = 0.72; r = 0.70). While ER binding assays detected chemicals that were highly potent or efficacious in the T-47D cell growth and transactivation assays, the binding assays lacked sensitivity in detecting weakly active compounds. In conclusion, this cell-based assay rapidly detects chemical effects on T-47D growth and shows potential, in combination with other HTS assays, to detect environmentally relevant chemicals with potential estrogenic activity.

Aggregating data for computational toxicology applications: The U.S. Environmental Protection Agency (EPA) Aggregated Computational Toxicology Resource (ACToR) System.

Computational toxicology combines data from high-throughput test methods, chemical structure analyses and other biological domains (e.g., genes, proteins, cells, tissues) with the goals of predicting and understanding the underlying mechanistic causes of chemical toxicity and for predicting toxicity of new chemicals and products. A key feature of such approaches is their reliance on knowledge extracted from large collections of data and data sets in computable formats. The U.S. Environmental Protection Agency (EPA) has developed a large data resource called ACToR (Aggregated Computational Toxicology Resource) to support these data-intensive efforts. ACToR comprises four main repositories: core ACToR (chemical identifiers and structures, and summary data on hazard, exposure, use, and other domains), ToxRefDB (Toxicity Reference Database, a compilation of detailed in vivo toxicity data from guideline studies), ExpoCastDB (detailed human exposure data from observational studies of selected chemicals), and ToxCastDB (data from high-throughput screening programs, including links to underlying biological information related to genes and pathways). The EPA DSSTox (Distributed Structure-Searchable Toxicity) program provides expert-reviewed chemical structures and associated information for these and other high-interest public inventories. Overall, the ACToR system contains information on about 400,000 chemicals from 1100 different sources. The entire system is built using open source tools and is freely available to download. This review describes the organization of the data repository and provides selected examples of use cases.

Early Drug-Induced Liver Injury Risk Screening: "Free," as Good as It Gets.

For all the promise of and need for clinical drug-induced liver injury (DILI) risk screening systems, demonstrating the predictive value of these systems versus readily available physicochemical properties and inherent dosing information has not been thoroughly evaluated. Therefore, we utilized a systematic approach to evaluate the predictive value of in vitro safety assays including bile salt export pump transporter inhibition and cytotoxicity in HepG2 and transformed human liver epithelial along with physicochemical properties. We also evaluated the predictive value of in vitro ADME assays including hepatic partition coefficient (Kp) and its unbound counterpart because they provide insight on hepatic accumulation potential. The datasets comprised of 569 marketed drugs with FDA DILIrank annotation (most vs less/none), dose and physicochemical information, 384 drugs with Kp and plasma protein binding data, and 279 drugs with safety assay data. For each dataset and combination of input parameters, we developed random forest machine learning models and measured model performance using the receiver operator characteristic area under the curve (ROC AUC). The median ROC AUC across the various data and parameters sets ranged from 0.67 to 0.77 with little evidence of additive predictivity when including safety or ADME assay data. Subsequent machine learning models consistently demonstrated daily dose, fraction sp3 or ionization, and cLogP/D inputs produced the best, simplest model for predicting clinical DILI risk with an ROC AUC of 0.75. This systematic framework should be used for future assay predictive value assessments and highlights the need for continued improvements to clinical DILI risk annotation.

Predictive model of rat reproductive toxicity from ToxCast high throughput screening.

The U.S. Environmental Protection Agency's ToxCast research program uses high throughput screening (HTS) for profiling bioactivity and predicting the toxicity of large numbers of chemicals. ToxCast Phase I tested 309 well-characterized chemicals in more than 500 assays for a wide range of molecular targets and cellular responses. Of the 309 environmental chemicals in Phase I, 256 were linked to high-quality rat multigeneration reproductive toxicity studies in the relational Toxicity Reference Database. Reproductive toxicants were defined here as having achieved a reproductive lowest-observed-adverse-effect level of less than 500 mg kg(-1) day(-1). Eight-six chemicals were identified as reproductive toxicants in the rat, and 68 of those had sufficient in vitro bioactivity to model. Each assay was assessed for univariate association with the identified reproductive toxicants. Significantly associated assays were linked to gene sets and used for the subsequent predictive modeling. Using linear discriminant analysis and fivefold cross-validation, a robust and stable predictive model was produced capable of identifying rodent reproductive toxicants with 77% ± 2% and 74% ± 5% (mean ± SEM) training and test cross-validation balanced accuracies, respectively. With a 21-chemical external validation set, the model was 76% accurate, further indicating the model's potential for prioritizing the many thousands of environmental chemicals with little to no hazard information. The biological features of the model include steroidal and nonsteroidal nuclear receptors, cytochrome P450 enzyme inhibition, G protein-coupled receptors, and cell signaling pathway readouts-mechanistic information suggesting additional targeted, integrated testing strategies and potential applications of in vitro HTS to risk assessment.

Estimating toxicity-related biological pathway altering doses for high-throughput chemical risk assessment.

We describe a framework for estimating the human dose at which a chemical significantly alters a biological pathway in vivo, making use of in vitro assay data and an in vitro-derived pharmacokinetic model, coupled with estimates of population variability and uncertainty. The quantity we calculate, the biological pathway altering dose (BPAD), is analogous to current risk assessment metrics in that it combines dose-response data with analysis of uncertainty and population variability to arrive at conservative exposure limits. The analogy is closest when perturbation of a pathway is a key event in the mode of action (MOA) leading to a specified adverse outcome. Because BPADs are derived from relatively inexpensive, high-throughput screening (HTS) in vitro data, this approach can be applied to high-throughput risk assessments (HTRA) for thousands of data-poor environmental chemicals. We envisage the first step of HTRA to be an assessment of in vitro concentration-response relationships across biologically important pathways to derive biological pathway altering concentrations (BPAC). Pharmacokinetic (PK) modeling is then used to estimate the in vivo doses required to achieve the BPACs in the blood at steady state. Uncertainty and variability are incorporated in both the BPAC and the PK parameters and then combined to yield a probability distribution for the dose required to perturb the critical pathway. We finally define the BPADL as the lower confidence bound of this pathway-altering dose. This perspective outlines a framework for using HTRA to estimate BPAD values; provides examples of the use of this approach, including a comparison of BPAD values with published dose-response data from in vivo studies; and discusses challenges and alternative formulations.

Impact of environmental chemicals on key transcription regulators and correlation to toxicity end points within EPA's ToxCast program.

Exposure to environmental chemicals adds to the burden of disease in humans and wildlife to a degree that is difficult to estimate and, thus, mitigate. The ability to assess the impact of existing chemicals for which little to no toxicity data are available or to foresee such effects during early stages of chemical development and use, and before potential exposure occurs, is a pressing need. However, the capacity of the current toxicity evaluation approaches to meet this demand is limited by low throughput and high costs. In the context of EPA's ToxCast project, we have evaluated a novel cellular biosensor system (Factorial (1) ) that enables rapid, high-content assessment of a compound's impact on gene regulatory networks. The Factorial biosensors combined libraries of cis- and trans-regulated transcription factor reporter constructs with a highly homogeneous method of detection enabling simultaneous evaluation of multiplexed transcription factor activities. Here, we demonstrate the application of the technology toward determining bioactivity profiles by quantitatively evaluating the effects of 309 environmental chemicals on 25 nuclear receptors and 48 transcription factor response elements. We demonstrate coherent transcription factor activity across nuclear receptors and their response elements and that Nrf2 activity, a marker of oxidative stress, is highly correlated to the overall promiscuity of a chemical. Additionally, as part of the ToxCast program, we identify molecular targets that associate with in vivo end points and represent modes of action that can serve as potential toxicity pathway biomarkers and inputs for predictive modeling of in vivo toxicity.

Xenobiotic-metabolizing enzyme and transporter gene expression in primary cultures of human hepatocytes modulated by ToxCast chemicals.

Primary human hepatocyte cultures are useful in vitro model systems of human liver because when cultured under appropriate conditions the hepatocytes retain liver-like functionality such as metabolism, transport, and cell signaling. This model system was used to characterize the concentration- and time-response of the 320 ToxCast chemicals for changes in expression of genes regulated by nuclear receptors. Fourteen gene targets were monitored in quantitative nuclease protection assays: six representative cytochromes P-450, four hepatic transporters, three Phase II conjugating enzymes, and one endogenous metabolism gene involved in cholesterol synthesis. These gene targets are sentinels of five major signaling pathways: AhR, CAR, PXR, FXR, and PPARalpha. Besides gene expression, the relative potency and efficacy for these chemicals to modulate cellular health and enzymatic activity were assessed. Results demonstrated that the culture system was an effective model of chemical-induced responses by prototypical inducers such as phenobarbital and rifampicin. Gene expression results identified various ToxCast chemicals that were potent or efficacious inducers of one or more of the 14 genes, and by inference the 5 nuclear receptor signaling pathways. Significant relative risk associations with rodent in vivo chronic toxicity effects are reported for the five major receptor pathways. These gene expression data are being incorporated into the larger ToxCast predictive modeling effort.

Analysis of eight oil spill dispersants using rapid, in vitro tests for endocrine and other biological activity.

The Deepwater Horizon oil spill has led to the use of >1 M gallons of oil spill dispersants, which are mixtures of surfactants and solvents. Because of this large scale use there is a critical need to understand the potential for toxicity of the currently used dispersant and potential alternatives, especially given the limited toxicity testing information that is available. In particular, some dispersants contain nonylphenol ethoxylates (NPEs), which can degrade to nonylphenol (NP), a known endocrine disruptor. Given the urgent need to generate toxicity data, we carried out a series of in vitro high-throughput assays on eight commercial dispersants. These assays focused on the estrogen and androgen receptors (ER and AR), but also included a larger battery of assays probing other biological pathways. Cytotoxicity in mammalian cells was also quantified. No activity was seen in any AR assay. Two dispersants showed a weak ER signal in one assay (EC50 of 16 ppm for Nokomis 3-F4 and 25 ppm for ZI-400). NPs and NPEs also had a weak signal in this same ER assay. Note that Corexit 9500, the currently used product, does not contain NPEs and did not show any ER activity. Cytotoxicity values for six of the dispersants were statistically indistinguishable, with median LC50 values approximately 100 ppm. Two dispersants, JD 2000 and SAF-RON GOLD, were significantly less cytotoxic than the others with LC50 values approaching or exceeding 1000 ppm.

Variability in in vivo studies: Defining the upper limit of performance for predictions of systemic effect levels.

New approach methodologies (NAMs) for chemical hazard assessment are often evaluated via comparison to animal studies; however, variability in animal study data limits NAM accuracy. The US EPA Toxicity Reference Database (ToxRefDB) enables consideration of variability in effect levels, including the lowest effect level (LEL) for a treatment-related effect and the lowest observable adverse effect level (LOAEL) defined by expert review, from subacute, subchronic, chronic, multi-generation reproductive, and developmental toxicity studies. The objectives of this work were to quantify the variance within systemic LEL and LOAEL values, defined as potency values for effects in adult or parental animals only, and to estimate the upper limit of NAM prediction accuracy. Multiple linear regression (MLR) and augmented cell means (ACM) models were used to quantify the total variance, and the fraction of variance in systemic LEL and LOAEL values explained by available study descriptors (e.g., administration route, study type). The MLR approach considered each study descriptor as an independent contributor to variance, whereas the ACM approach combined categorical descriptors into cells to define replicates. Using these approaches, total variance in systemic LEL and LOAEL values (in log10-mg/kg/day units) ranged from 0.74 to 0.92. Unexplained variance in LEL and LOAEL values, approximated by the residual mean square error (MSE), ranged from 0.20-0.39. Considering subchronic, chronic, or developmental study designs separately resulted in similar values. Based on the relationship between MSE and R-squared for goodness-of-fit, the maximal R-squared may approach 55 to 73% for a NAM-based predictive model of systemic toxicity using these data as reference. The root mean square error (RMSE) ranged from 0.47 to 0.63 log10-mg/kg/day, depending on dataset and regression approach, suggesting that a two-sided minimum prediction interval for systemic effect levels may have a width of 58 to 284-fold. These findings suggest quantitative considerations for building scientific confidence in NAM-based systemic toxicity predictions.

ToxRefDB version 2.0: Improved utility for predictive and retrospective toxicology analyses.

The Toxicity Reference Database (ToxRefDB) structures information from over 5000 in vivo toxicity studies, conducted largely to guidelines or specifications from the US Environmental Protection Agency and the National Toxicology Program, into a public resource for training and validation of predictive models. Herein, ToxRefDB version 2.0 (ToxRefDBv2) development is described. Endpoints were annotated (e.g. required, not required) according to guidelines for subacute, subchronic, chronic, developmental, and multigenerational reproductive designs, distinguishing negative responses from untested. Quantitative data were extracted, and dose-response modeling for nearly 28,000 datasets from nearly 400 endpoints using Benchmark Dose (BMD) Modeling Software were generated and stored. Implementation of controlled vocabulary improved data quality; standardization to guideline requirements and cross-referencing with United Medical Language System (UMLS) connects ToxRefDBv2 observations to vocabularies linked to UMLS, including PubMed medical subject headings. ToxRefDBv2 allows for increased connections to other resources and has greatly enhanced quantitative and qualitative utility for predictive toxicology.

Profiling 58 compounds including cosmetic-relevant chemicals using ToxRefDB and ToxCast.

Safety assessment for cosmetic-relevant chemicals (CRCs) in the European Union has been reshaped by restrictions on animal testing, and new approach methodologies (NAMs) for predicting toxicity are critical to ensure new cosmetic product safety. To demonstrate NAMs for safety assessment, we surveyed in vitro bioactivity and in vivo systemic toxicity data in the US Environmental Protection Agency's (EPA's) Toxicity Forecaster (ToxCast) and Toxicity Reference databases (ToxRefDB), respectively, for 58 chemicals identified as CRCs, including cosmetic ingredients as well as trace contaminants. CRCs were diverse in use types as suggested by broad chemical use categories. In terms of both target organ effects and study type, the median of the lowest effect level (LEL) doses in ToxRefDB for CRCs tended to be slightly higher than the median for the remaining 928 chemicals with study data in ToxRefDB, though the ranges of LELs were similar. For 17 of the 58 CRCs, high-throughput toxicokinetic data were used to calculate administered equivalent doses (AEDs) in mg/kg/day units for the in vitro bioactivity observed, and these AEDs served as conservative estimators of the systemic LELs observed in vivo. This work suggests that NAMs for bioactivity may inform a conservative point-of-departure estimate for diverse CRCs.

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.

Use of high-throughput in vitro toxicity screening data in cancer hazard evaluations by IARC Monograph Working Groups.

Evidence regarding carcinogenic mechanisms serves a critical role in International Agency for Research on Cancer (IARC) Monograph evaluations. Three recent IARC Working Groups pioneered inclusion of the US Environmental Protection Agency (EPA) ToxCast program high-throughput screening (HTS) data to supplement other mechanistic evidence. In Monograph V110, HTS profiles were compared between perfluorooctanoic acid (PFOA) and prototypical activators across multiple nuclear receptors. For Monograph V112-113, HTS assays were mapped to 10 key characteristics of carcinogens identified by an IARC expert group, and systematically considered as an additional mechanistic data stream. Both individual assay results and ToxPi-based rankings informed mechanistic evaluations. Activation of multiple nuclear receptors in HTS assays showed that PFOA targets not only peroxisome proliferator activated receptors, but also other receptors. ToxCast assays substantially covered 5 of 10 key characteristics, corroborating literature evidence of "induces oxidative stress" and "alters cell proliferation, cell death or nutrient supply" and filling gaps for "modulates receptor-mediated effects." Thus, ToxCast HTS data were useful both in evaluating specific mechanistic hypotheses and in contributing to the overall evaluation of mechanistic evidence. However, additional HTS assays are needed to provide more comprehensive coverage of the 10 key characteristics of carcinogens that form the basis of current IARC mechanistic evaluations.