Verification of In Vivo Estrogenic Activity for Four Per- and Polyfluoroalkyl Substances (PFAS) Identified as Estrogen Receptor Agonists via New Approach Methodologies.

Given concerns about potential toxicological hazards of the thousands of data-poor per- and polyfluorinated alkyl substances (PFAS) currently in commerce and detected in the environment, tiered testing strategies that employ high-throughput in vitro screening as an initial testing tier have been implemented. The present study evaluated the effectiveness of previous in vitro screening for identifying PFAS capable, or incapable, of inducing estrogenic responses in fish exposed in vivo. Fathead minnows (Pimephales promelas) were exposed for 96 h to five PFAS (perfluorooctanoic acid [PFOA]; 1H,1H,8H,8H-perfluorooctane-1,8-diol [FC8-diol]; 1H,1H,10H,10H-perfluorodecane-1,10-diol [FC10-diol]; 1H,1H,8H,8H-perfluoro-3,6-dioxaoctane-1,8-diol [FC8-DOD]; and perfluoro-2-methyl-3-oxahexanoic acid [HFPO-DA]) that showed varying levels of in vitro estrogenic potency. In agreement with in vitro screening results, exposure to FC8-diol, FC10-diol, and FC8-DOD caused concentration-dependent increases in the expression of transcript coding for vitellogenin and estrogen receptor alpha and reduced expression of insulin-like growth factor and apolipoprotein eb. Once differences in bioconcentration were accounted for, the rank order of potency in vivo matched that determined in vitro. These results provide a screening level benchmark for worst-case estimates of potential estrogenic hazards of PFAS and a basis for identifying structurally similar PFAS to scrutinize for putative estrogenic activity.

Tox21BodyMap: a webtool to map chemical effects on the human body.

To support rapid chemical toxicity assessment and mechanistic hypothesis generation, here we present an intuitive webtool allowing a user to identify target organs in the human body where a substance is estimated to be more likely to produce effects. This tool, called Tox21BodyMap, incorporates results of 9,270 chemicals tested in the United States federal Tox21 research consortium in 971 high-throughput screening (HTS) assays whose targets were mapped onto human organs using organ-specific gene expression data. Via Tox21BodyMap's interactive tools, users can visualize chemical target specificity by organ system, and implement different filtering criteria by changing gene expression thresholds and activity concentration parameters. Dynamic network representations, data tables, and plots with comprehensive activity summaries across all Tox21 HTS assay targets provide an overall picture of chemical bioactivity. Tox21BodyMap webserver is available at https://sandbox.ntp.niehs.nih.gov/bodymap/.

The Tox21 10K Compound Library: Collaborative Chemistry Advancing Toxicology.

Since 2009, the Tox21 project has screened ∼8500 chemicals in more than 70 high-throughput assays, generating upward of 100 million data points, with all data publicly available through partner websites at the United States Environmental Protection Agency (EPA), National Center for Advancing Translational Sciences (NCATS), and National Toxicology Program (NTP). Underpinning this public effort is the largest compound library ever constructed specifically for improving understanding of the chemical basis of toxicity across research and regulatory domains. Each Tox21 federal partner brought specialized resources and capabilities to the partnership, including three approximately equal-sized compound libraries. All Tox21 data generated to date have resulted from a confluence of ideas, technologies, and expertise used to design, screen, and analyze the Tox21 10K library. The different programmatic objectives of the partners led to three distinct, overlapping compound libraries that, when combined, not only covered a diversity of chemical structures, use-categories, and properties but also incorporated many types of compound replicates. The history of development of the Tox21 "10K" chemical library and data workflows implemented to ensure quality chemical annotations and allow for various reproducibility assessments are described. Cheminformatics profiling demonstrates how the three partner libraries complement one another to expand the reach of each individual library, as reflected in coverage of regulatory lists, predicted toxicity end points, and physicochemical properties. ToxPrint chemotypes (CTs) and enrichment approaches further demonstrate how the combined partner libraries amplify structure-activity patterns that would otherwise not be detected. Finally, CT enrichments are used to probe global patterns of activity in combined ToxCast and Tox21 activity data sets relative to test-set size and chemical versus biological end point diversity, illustrating the power of CT approaches to discern patterns in chemical-activity data sets. These results support a central premise of the Tox21 program: A collaborative merging of programmatically distinct compound libraries would yield greater rewards than could be achieved separately.

Harmonized Cross-Species Assessment of Endocrine and Metabolic Disruptors by Ecotox FACTORIAL Assay.

Environmental pollution is a threat to humans and wildlife species. Of particular concern are endocrine disrupting chemicals (EDCs). An important target of EDCs is nuclear receptors (NRs) that control endocrine and metabolic responses through transcriptional regulation. Owing in part to structural differences of NRs, adverse effects of EDCs vary significantly among species. Here, we describe a multiplexed reporter assay (the Ecotox FACTORIAL) enabling parallel assessment of compounds' effects on estrogen, androgen, thyroid, and PPARγ receptors of representative mammals, birds, reptiles, amphibians, and fish. The Ecotox FACTORIAL is a single-well assay comprising a set of species-specific, one-hybrid GAL4-NR reporter constructs transiently transfected into test cells. To harmonize cross-species assessments, we used a combination of two approaches. First, we used the same type of test cells for all reporters; second, we implemented a parallel detection of reporter RNAs. The assay demonstrated excellent quality, reproducibility, and insignificant intra-assay variability. Importantly, the EC50 values for NR ligands were consistent with those reported for conventional assays. Using the assay allowed ranking the hazard potential of environmental pollutants (e.g., bisphenols, polycyclic aromatic hydrocarbons, and synthetic progestins) across species. Furthermore, the assay permitted detecting taxa-specific effects of surface water samples. Therefore, the Ecotox FACTORIAL enables harmonized assessment of the endocrine and metabolic disrupting activity of chemicals and surface water in humans as well as in wildlife species.

Nontarget Screening of Per- and Polyfluoroalkyl Substances Binding to Human Liver Fatty Acid Binding Protein.

More than 1000 per- and polyfluoroalkyl substances (PFASs) have been discovered by nontarget analysis (NTA), but their prioritization for health concerns is challenging. We developed a method by incorporating size-exclusion column co-elution (SECC) and NTA, to screen PFASs binding to human liver fatty acid binding protein (hL-FABP). Of 74 PFASs assessed, 20 were identified as hL-FABP ligands in which eight of them have high binding affinities. Increased PFAS binding affinities correlate with stronger responses in electrospray ionization (ESI-) and longer retention times on a C18 column. This is well explained by a mechanistic model, which revealed that both polar and hydrophobic interactions are crucial for binding affinities. Encouraged by this, we then developed an SECC method to identify hL-FABP ligands, and all eight high-affinity ligands were selectively captured from 74 PFASs. The method was further applied to an aqueous film-forming foam (AFFF) product in which 31 new hL-FABP ligands were identified. Suspect and nontargeted screening revealed these ligands as analogues of perfluorosulfonic acids and homologues of alkyl ether sulfates (C8- and C10/EOn, C8H17(C2H4O)nSO4-, and C10H21(C2H4O)nSO4-). The SECC method was then applied to AFFF-contaminated surface waters. In addition to perfluorooctanesulfonic acid and perfluorohexanesulfonic acid, eight other AFFF chemicals were discovered as novel ligands, including four C14- and C15/EOn. This study implemented a high-throughput method to prioritize PFASs and revealed the existence of many previously unknown hL-FABP ligands.

Potential Toxicity of Complex Mixtures in Surface Waters from a Nationwide Survey of United States Streams: Identifying in Vitro Bioactivities and Causative Chemicals.

While chemical analysis of contaminant mixtures remains an essential component of environmental monitoring, bioactivity-based assessments using in vitro systems increasingly are used in the detection of biological effects. Historically, in vitro assessments focused on a few biological pathways, for example, aryl hydrocarbon receptor (AhR) or estrogen receptor (ER) activities. High-throughput screening (HTS) technologies have greatly increased the number of biological targets and processes that can be rapidly assessed. Here we screened extracts of surface waters from a nationwide survey of United States streams for bioactivities associated with 69 different end points using two multiplexed HTS assays. Bioactivity of extracts from 38 streams was evaluated and compared with concentrations of over 700 analytes to identify chemicals contributing to observed effects. Eleven primary biological end points were detected. Pregnane X receptor (PXR) and AhR-mediated activities were the most commonly detected. Measured chemicals did not completely account for AhR and PXR responses. Surface waters with AhR and PXR effects were associated with low intensity, developed land cover. Likewise, elevated bioactivities frequently associated with wastewater discharges included endocrine-related end points ER and glucocorticoid receptor. These results underscore the value of bioassay-based monitoring of environmental mixtures for detecting biological effects that could not be ascertained solely through chemical analyses.

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.

On selecting a minimal set of in vitro assays to reliably determine estrogen agonist activity.

The US EPA is charged with screening chemicals for their ability to be endocrine disruptors through interaction with the estrogen, androgen and thyroid axes. The agency is exploring the use of high-throughput in vitro assays to use in the Endocrine Disruptor Screening Program (EDSP), potentially as replacements for lower-throughput in vitro and in vivo tests. The first replacement is an integrated computational and experimental model for estrogen receptor (ER) activity, to be used as an alternative to the EDSP Tier 1 in vitro ER binding and transactivation assays and the in vivo uterotrophic bioassay. The ER agonist model uses a set of 16 in vitro assays that incorporate multiple technologies and cell lines and probe multiple points in the ER pathway. Here, we demonstrate that subsets of assays with as few as 4 assays can predict the activity of all 1811 chemicals tested with accuracy equivalent to that of the full 16-assay model. The prediction accuracy against reference chemicals is higher than that of the full chemical set, partly because the larger set contains many chemicals that can cause a variety of types of assay interference There are multiple accurate assay subsets, allowing flexibility in the construction of a multiplexed assay battery. We also discuss the issue of challenging chemicals, i.e. those that can give false positive results in certain assays, and could hence be more problematic when only a few assays are used.

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.

Predictive endocrine testing in the 21st century using in vitro assays of estrogen receptor signaling responses.

Thousands of environmental chemicals are subject to regulatory review for their potential to be endocrine disruptors (ED). In vitro high-throughput screening (HTS) assays have emerged as a potential tool for prioritizing chemicals for ED-related whole-animal tests. In this study, 1814 chemicals including pesticide active and inert ingredients, industrial chemicals, food additives, and pharmaceuticals were evaluated in a panel of 13 in vitro HTS assays. The panel of in vitro assays interrogated multiple end points related to estrogen receptor (ER) signaling, namely binding, agonist, antagonist, and cell growth responses. The results from the in vitro assays were used to create an ER Interaction Score. For 36 reference chemicals, an ER Interaction Score >0 showed 100% sensitivity and 87.5% specificity for classifying potential ER activity. The magnitude of the ER Interaction Score was significantly related to the potency classification of the reference chemicals (p < 0.0001). ERα/ERß selectivity was also evaluated, but relatively few chemicals showed significant selectivity for a specific isoform. When applied to a broader set of chemicals with in vivo uterotrophic data, the ER Interaction Scores showed 91% sensitivity and 65% specificity. Overall, this study provides a novel method for combining in vitro concentration response data from multiple assays and, when applied to a large set of ER data, accurately predicted estrogenic responses and demonstrated its utility for chemical prioritization.

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.

Evaluation of 147 perfluoroalkyl substances for immunotoxic and other (patho)physiological activities through phenotypic screening of human primary cells.

A structurally diverse set of 147 per- and polyfluoroalkyl substances (PFAS) was screened in a panel of 12 human primary cell systems by measuring 148 biomarkers relevant to (patho)physiological pathways to inform hypotheses about potential mechanistic effects of data-poor PFAS in human model systems. This analysis focused on immunosuppressive activity, which was previously reported as an in vivo effect of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), by comparing PFAS responses to four pharmacological immunosuppressants. The PFOS response profile had little correlation with reference immunosuppressants, suggesting in vivo activity does not occur by similar mechanisms. The PFOA response profile did share features with the profile of dexamethasone, although some distinct features were lacking. Other PFAS, including 2,2,3,3-tetrafluoropropyl acrylate, demonstrated more similarity to the reference immunosuppressants but with additional activities not found in the reference immunosuppressive drugs. Correlation of PFAS profiles with a database of environmental chemical responses and pharmacological probes identified potential mechanisms of bioactivity for some PFAS, including responses similar to ubiquitin ligase inhibitors, deubiquitylating enzyme (DUB) inhibitors, and thioredoxin reductase inhibitors. Approximately 21% of the 147 PFAS with confirmed sample quality were bioactive at nominal testing concentrations in the 1-60 micromolar range in these human primary cell systems. These data provide new hypotheses for mechanisms of action for a subset of PFAS and may further aid in development of a PFAS categorization strategy useful in safety assessment.

A gene expression biomarker for predictive toxicology to identify chemical modulators of NF-κB.

The nuclear factor-kappa B (NF-κB) is a transcription factor with important roles in inflammation, immune response, and oncogenesis. Dysregulation of NF-κB signaling is associated with inflammation and certain cancers. We developed a gene expression biomarker predictive of NF-κB modulation and used the biomarker to screen a large compendia of gene expression data. The biomarker consists of 108 genes responsive to tumor necrosis factor α in the absence but not the presence of IκB, an inhibitor of NF-κB. Using a set of 450 profiles from cells treated with immunomodulatory factors with known NF-κB activity, the balanced accuracy for prediction of NF-κB activation was > 90%. The biomarker was used to screen a microarray compendium consisting of 12,061 microarray comparisons from human cells exposed to 2,672 individual chemicals to identify chemicals that could cause toxic effects through NF-κB. There were 215 and 49 chemicals that were identified as putative or known NF-κB activators or suppressors, respectively. NF-κB activators were also identified using two high-throughput screening assays; 165 out of the ~3,800 chemicals (ToxCast assay) and 55 out of ~7,500 unique compounds (Tox21 assay) were identified as potential activators. A set of 32 chemicals not previously associated with NF-κB activation and which partially overlapped between the different screens were selected for validation in wild-type and NFKB1-null HeLa cells. Using RT-qPCR and targeted RNA-Seq, 31 of the 32 chemicals were confirmed to be NF-κB activators. These results comprehensively identify a set of chemicals that could cause toxic effects through NF-κB.

Comprehensive assessment of NR ligand polypharmacology by a multiplex reporter NR assay.

Nuclear receptors (NR) are ligand-modulated transcription factors that regulate multiple cell functions and thus represent excellent drug targets. However, due to a considerable NR structural homology, NR ligands often interact with multiple receptors. Here, we describe a multiplex reporter assay (the FACTORIAL NR) that enables parallel assessment of NR ligand activity across all 48 human NRs. The assay comprises one-hybrid GAL4-NR reporter modules transiently transfected into test cells. To evaluate the reporter activity, we assessed their RNA transcripts. We used a homogeneous RNA detection approach that afforded equal detection efficacy and permitted the multiplex detection in a single-well format. For validation, we examined a panel of selective NR ligands and polypharmacological agonists and antagonists of the progestin, estrogen, PPAR, ERR, and ROR receptors. The assay produced highly reproducible NR activity profiles (r > 0.96) permitting quantitative assessment of individual NR responses. The inferred EC50 values agreed with the published data. The assay showed excellent quality ( = 0.73) and low variability ( = 7.2%). Furthermore, the assay permitted distinguishing direct and non-direct NR responses to ligands. Therefore, the FACTORIAL NR enables comprehensive evaluation of NR ligand polypharmacology.

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.

A Roadmap to the Structure-Related Metabolism Pathways of Per- and Polyfluoroalkyl Substances in the Early Life Stages of Zebrafish (Danio rerio).

BACKGROUND: Thousands of per- and polyfluoroalkyl substances (PFAS) with diverse structures have been detected in the ambient environment. Apart from a few well-studied PFAS, the structure-related toxicokinetics of a broader set of PFAS remain unclear. OBJECTIVES: To understand the toxicokinetics of PFAS, we attempted to characterize the metabolism pathways of 74 structurally diverse PFAS samples from the U.S. Environmental Protection Agency's PFAS screening library. METHODS: Using the early life stages of zebrafish (Danio rerio) as a model, we determined the bioconcentration factors and phenotypic toxicities of 74 PFAS. Then, we applied high-resolution mass spectrometry-based nontargeted analysis to identify metabolites of PFAS in zebrafish larvae after 5 d of exposure by incorporating retention time and mass spectra. In vitro enzymatic activity experiments with human recombinant liver carboxylesterase (hCES1) were employed to validate the structure-related hydrolysis of 11 selected PFAS. RESULTS: Our findings identified five structural categories of PFAS prone to metabolism. The metabolism pathways of PFAS were highly related to their structures as exemplified by fluorotelomer alcohols that the predominance of ß-oxidation or taurine conjugation pathways were primarily determined by the number of hydrocarbons. Hydrolysis was identified as a major metabolism pathway for diverse PFAS, and perfluoroalkyl carboxamides showed the highest in vivo hydrolysis rates, followed by carboxyesters and sulfonamides. The hydrolysis of PFAS was verified with recombinant hCES1, with strong substrate preferences toward perfluoroalkyl carboxamides. CONCLUSIONS: We suggest that the roadmap of the structure-related metabolism pathways of PFAS established in this study would provide a starting point to inform the potential health risks of other PFAS. https://doi.org/10.1289/EHP7169.