Towards reproducible structure-based chemical categories for PFAS to inform and evaluate toxicity and toxicokinetic testing

Using Zebrafish to Screen Developmental Toxicity of Per- and Polyfluoroalkyl Substances (PFAS)

Per- and polyfluoroalkyl substances (PFAS) are found in many consumer and industrial products. While some PFAS, notably perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), are developmentally toxic in mammals, the vast majority of PFAS have not been evaluated for developmental toxicity potential. A concentration-response study of 182 unique PFAS chemicals using the zebrafish medium-throughput, developmental vertebrate toxicity assay was conducted to investigate chemical structural identifiers for toxicity. Embryos were exposed to each PFAS compound (≤100 μM) beginning on the day of fertilization. At 6 days post-fertilization (dpf), two independent observers graded developmental landmarks for each larva (e.g., mortality, hatching, swim bladder inflation, edema, abnormal spine/tail, or craniofacial structure). Thirty percent of the PFAS were developmentally toxic, but there was no enrichment of any OECD structural category. PFOS was developmentally toxic (benchmark concentration [BMC] = 7.48 μM); however, other chemicals were more potent: perfluorooctanesulfonamide (PFOSA), N-methylperfluorooctane sulfonamide (N-MeFOSA), ((perfluorooctyl)ethyl)phosphonic acid, perfluoro-3,6,9-trioxatridecanoic acid, and perfluorohexane sulfonamide. The developmental toxicity profile for these more potent PFAS is largely unexplored in mammals and other species. Based on these zebrafish developmental toxicity results, additional screening may be warranted to understand the toxicity profile of these chemicals in other species.

An artificial intelligence platform for automated PFAS subgroup classification: A discovery tool for PFAS screening

Since structural analyses and toxicity assessments have not been able to keep up with the discovery of unknown per- and polyfluoroalkyl substances (PFAS), there is an urgent need for effective categorization and grouping of PFAS. In this study, we presented PFAS-Atlas, an artificial intelligence-based platform containing a rule-based automatic classification system and a machine learning-based grouping model. Compared with previously developed classification software, the platform's classification system follows the latest Organization for Economic Co-operation and Development (OECD) definition of PFAS and reduces the number of uncategorized PFAS. In addition, the platform incorporates deep unsupervised learning models to visualize the chemical space of PFAS by clustering similar structures and linking related classes. Through real-world use cases, we demonstrate that PFAS-Atlas can rapidly screen for relationships between chemical structure and persistence, bioaccumulation, or toxicity data for PFAS. The platform can also guide the planning of the PFAS testing strategy by showing which PFAS classes urgently require further attention. Ultimately, the release of PFAS-Atlas will benefit both the PFAS research and regulation communities.

Transcriptomics-Based Points of Departure for Daphnia magna Exposed to 18 Per- and Polyfluoroalkyl Substances

Per- and polyfluoroalkyl substances (PFAS) represent a large group of contaminants of concern based on their widespread use, environmental persistence, and potential toxicity. Many traditional models for estimating toxicity, bioaccumulation, and other toxicological properties are not well suited for PFAS. Consequently, there is a need to generate hazard information for PFAS in an efficient and cost-effective manner. In the present study, Daphnia magna were exposed to multiple concentrations of 22 different PFAS for 24 h in a 96-well plate format. Following exposure, whole-body RNA was extracted and extracts, each representing five exposed individuals, were subjected to RNA sequencing. Following analytical measurements to verify PFAS exposure concentrations and quality control on processed cDNA libraries for sequencing, concentration-response modeling was applied to the data sets for 18 of the tested compounds, and the concentration at which a concerted molecular response occurred (transcriptomic point of departure; tPOD) was calculated. The tPODs, based on measured concentrations of PFAS, generally ranged from 0.03 to 0.58 µM (9.9-350 µg/L; interquartile range). In most cases, these concentrations were two orders of magnitude lower than similarly calculated tPODs for human cell lines exposed to PFAS. They were also lower than apical effect concentrations reported for seven PFAS for which some crustacean or invertebrate toxicity data were available, although there were a few exceptions. Despite being lower than most other available hazard benchmarks, D. magna tPODs were, on average, four orders of magnitude greater than the maximum aqueous concentrations of PFAS measured in Great Lakes tributaries. Overall, this high-throughput transcriptomics assay with D. magna holds promise as a component of a tiered hazard evaluation strategy employing new approach methodologies. Environ Toxicol Chem 2024;00:1-16. © 2024 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

In vitro screening of per- and polyfluorinated substances (PFAS) for interference with seven thyroid hormone system targets across nine assays

The US Environmental Protection Agency is evaluating the ecological and toxicological effects of per- and polyfluorinated chemicals. A number of perfluorinated chemicals have been shown to impact the thyroid axis in vivo suggesting that the thyroid hormone system is a target of these chemicals. The objective of this study was to evaluate the activity of 136 perfluorinated chemicals at seven key molecular initiating events (MIE) within the thyroid axis using nine in vitro assays. The potential MIE targets investigated are Human Iodothyronine Deiodinase 1 (hDIO1), Human Iodothyronine Deiodinase 2 (hDIO2), Human Iodothyronine Deiodinase 3 (hDIO3), Xenopus Iodothyronine Deiodinase (xDIO3); Human Iodotyrosine Deiodinase (hIYD), Xenopus Iodotyrosine Deiodinase (xIYD), Human Thyroid Peroxidase (hTPO); and the serum binding proteins Human Transthyretin (hTTR) and Human Thyroxine Binding Globulin (hTBG). Of the 136 PFAS chemicals tested, 85 had sufficient activity to produce a half-maximal effect concentration (EC50) in at least one of the nine assays. In general, most of these PFAS chemicals did not have strong potency towards the seven MIEs examined, apart from transthyretin binding, for which several PFAS had potency similar to the respective model inhibitor. These data sets identify potentially active PFAS chemicals to prioritize for further testing in orthogonal in vitro assays and at higher levels of biological organization to evaluate their capacity for altering the thyroid hormone system and causing potential adverse health and ecological effects.

Hazard and risk characterization of 56 structurally diverse PFAS using a targeted battery of broad coverage assays using six human cell types

Per- and poly-fluoroalkyl substances (PFAS) are extensively used in commerce leading to their prevalence in the environment. Due to their chemical stability, PFAS are considered to be persistent and bioaccumulative; they are frequently detected in both the environment and humans. Because of this, PFAS as a class (composed of hundreds to thousands of chemicals) are contaminants of very high concern. Little information is available for the vast majority of PFAS, and regulatory agencies lack safety data to determine whether exposure limits or restrictions are needed. Cell-based assays are a pragmatic approach to inform decision-makers on potential health hazards; therefore, we hypothesized that a targeted battery of human in vitro assays can be used to determine whether there are structure-bioactivity relationships for PFAS, and to characterize potential risks by comparing bioactivity (points of departure) to exposure estimates. We tested 56 PFAS from 8 structure-based subclasses in concentration response (0.1-100 μM) using six human cell types selected from target organs with suggested adverse effects of PFAS - human induced pluripotent stem cell (iPSC)-derived hepatocytes, neurons, and cardiomyocytes, primary human hepatocytes, endothelial and HepG2 cells. While many compounds were without effect; certain PFAS demonstrated cell-specific activity highlighting the necessity of using a compendium of in vitro models to identify potential hazards. No class-specific groupings were evident except for some chain length- and structure-related trends. In addition, margins of exposure (MOE) were derived using empirical and predicted exposure data. Conservative MOE calculations showed that most tested PFAS had a MOE in the 1-100 range; ∼20% of PFAS had MOE<1, providing tiered priorities for further studies. Overall, we show that a compendium of human cell-based models can be used to derive bioactivity estimates for a range of PFAS, enabling comparisons with human biomonitoring data. Furthermore, we emphasize that establishing structure-bioactivity relationships may be challenging for the tested PFAS.

Expanded Systematic Evidence Map for Hundreds of Per- and Polyfluoroalkyl Substances (PFAS) and Comprehensive PFAS Human Health Dashboard

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) encompass a class of chemically and structurally diverse compounds that are extensively used in industry and detected in the environment. The US Environmental Protection Agency (US EPA) 2021 PFAS Strategic Roadmap describes national research plans to address the challenge of PFAS.

OBJECTIVES

Systematic Evidence Map (SEM) methods were used to survey and summarize available epidemiological and mammalian bioassay evidence that could inform human health hazard identification for a set of 345 PFAS that were identified by the US EPA's Center for Computational Toxicology and Exposure (CCTE) for in vitro toxicity and toxicokinetic assay testing and through interagency discussions on PFAS of interest. This work builds from the 2022 evidence map that collated evidence on a separate set of ∼ 150 PFAS. Like our previous work, this SEM does not include PFAS that are the subject of ongoing or completed assessments at the US EPA.

METHODS

SEM methods were used to search, screen, and inventory mammalian bioassay and epidemiological literature from peer-reviewed and gray literature sources using manual review and machine-learning software. For each included study, study design details and health end points examined were summarized in interactive web-based literature inventories. Some included studies also underwent study evaluation and detailed extraction of health end point data. All underlying data is publicly available online as interactive visuals with downloadable metadata.

RESULTS

More than 13,000 studies were identified from scientific databases. Screening processes identified 121 mammalian bioassay and 111 epidemiological studies that met screening criteria. Epidemiological evidence (available for 12 PFAS) mostly assessed the reproductive, endocrine, developmental, metabolic, cardiovascular, and immune systems. Mammalian bioassay evidence (available for 30 PFAS) commonly assessed effects in the reproductive, whole-body, nervous, and hepatic systems. Overall, 41 PFAS had evidence across mammalian bioassay and epidemiology data streams (roughly 11% of searched chemicals).

DISCUSSION

No epidemiological and/or mammalian bioassay evidence were identified for most of the PFAS included in our search. Results from this SEM, our 2022 SEM on ∼ 150 PFAS, and other PFAS assessment products from the US EPA are compiled into a comprehensive PFAS dashboard that provides researchers and regulators an overview of the current PFAS human health landscape including data gaps and can serve as a scoping tool to facilitate prioritization of PFAS-related research and/or risk assessment activities. https://doi.org/10.1289/EHP13423.

Polanyi adsorption potential theory for estimating PFAS treatment with granular activated carbon

Per- and polyfluoroalkyl substances (PFAS) are a group of chemicals that have gained interest because some PFAS have been shown to have negative health effects and prolonged environmental and biological persistence. Chemicals classified as PFAS have a wide range of chemical moieties that impart widely variable properties, leading to a range of water treatment process efficacies. The Polanyi Potential Adsorption Theory was used to estimate Freundlich isotherm parameters to predict the efficacy of granular activated carbon (GAC) treatment for 428 PFAS chemicals for which the vast majority had no previously published treatment data. This method accounts for the physical/chemical characteristics of the individual PFAS beyond molecular weight or chain length that have previously been employed. From a statistical analysis of available data and model results, many of the 428 PFAS were predicted to be effectively treatable by GAC. Although not directly applicable to full-scale design, the approach demonstrates a systematic method for predicting the effectiveness of GAC where isotherm or column data are not available. This then can be used to prioritize future research.

Plasma Protein Binding Evaluations of Per- and Polyfluoroalkyl Substances for Category-Based Toxicokinetic Assessment

New approach methodologies (NAMs) that make use of in vitro screening and in silico approaches to inform chemical evaluations rely on in vitro toxicokinetic (TK) data to translate in vitro bioactive concentrations to exposure metrics reflective of administered dose. With 1364 per- and polyfluoroalkyl substances (PFAS) identified as of interest under Section 8 of the U.S. Toxic Substances Control Act (TSCA) and concern over the lack of knowledge regarding environmental persistence, human health, and ecological effects, the utility of NAMs to understand potential toxicities and toxicokinetics across these data-poor compounds is being evaluated. To address the TK data deficiency, 71 PFAS selected to span a wide range of functional groups and physico-chemical properties were evaluated for in vitro human plasma protein binding (PPB) by ultracentrifugation with liquid chromatography-mass spectrometry analysis. For the 67 PFAS successfully evaluated by ultracentrifugation, fraction unbound in plasma (fup) ranged from less than 0.0001 (pentadecafluorooctanoyl chloride) to 0.7302 (tetrafluorosuccinic acid), with over half of the PFAS showing PPB exceeding 99.5% (fup < 0.005). Category-based evaluations revealed that perfluoroalkanoyl chlorides and perfluorinated carboxylates (PFCAs) with 6-10 carbons were the highest bound, with similar median values for alkyl, ether, and polyether PFCAs. Interestingly, binding was lower for the PFCAs with a carbon chain length of ≥11. Lower binding also was noted for fluorotelomer carboxylic acids when compared to their carbon-equivalent perfluoroalkyl acids. Comparisons of the fup value derived using two PPB methods, ultracentrifugation or rapid equilibrium dialysis (RED), revealed RED failure for a subset of PFAS of high mass and/or predicted octanol-water partition coefficients exceeding 4 due to failure to achieve equilibrium. Bayesian modeling was used to provide uncertainty bounds around fup point estimates for incorporation into TK modeling. This PFAS PPB evaluation and grouping exercise across 67 structures greatly expand our current knowledge and will aid in PFAS NAM development.

Category-Based Toxicokinetic Evaluations of Data-Poor Per- and Polyfluoroalkyl Substances (PFAS) using Gas Chromatography Coupled with Mass Spectrometry

Concern over per- and polyfluoroalkyl substances (PFAS) has increased as more is learned about their environmental presence, persistence, and bioaccumulative potential. The limited monitoring, toxicokinetic (TK), and toxicologic data available are inadequate to inform risk across this diverse domain. Here, 73 PFAS were selected for in vitro TK evaluation to expand knowledge across lesser-studied PFAS alcohols, amides, and acrylates. Targeted methods developed using gas chromatography-tandem mass spectrometry (GC-MS/MS) were used to measure human plasma protein binding and hepatocyte clearance. Forty-three PFAS were successfully evaluated in plasma, with fraction unbound (f_up) values ranging from 0.004 to 1. With a median f_up of 0.09 (i.e., 91% bound), these PFAS are highly bound but exhibit 10-fold lower binding than legacy perfluoroalkyl acids recently evaluated. Thirty PFAS evaluated in the hepatocyte clearance assay showed abiotic loss, with many exceeding 60% loss within 60 min. Metabolic clearance was noted for 11 of the 13 that were successfully evaluated, with rates up to 49.9 μL/(min × million cells). The chemical transformation simulator revealed potential (bio)transformation products to consider. This effort provides critical information to evaluate PFAS for which volatility, metabolism, and other routes of transformation are likely to modulate their environmental fates.

Evaluation of Per- and Polyfluoroalkyl Substances (PFAS) In Vitro Toxicity Testing for Developmental Neurotoxicity

Per- and polyfluoroalkyl substances (PFAS) are a diverse set of commercial chemicals widely detected in humans and the environment. However, only a limited number of PFAS are associated with epidemiological or experimental data for hazard identification. To provide developmental neurotoxicity (DNT) hazard information, the work herein employed DNT new approach methods (NAMs) to generate in vitro screening data for a set of 160 PFAS. The DNT NAMs battery was comprised of the microelectrode array neuronal network formation assay (NFA) and high-content imaging (HCI) assays to evaluate proliferation, apoptosis, and neurite outgrowth. The majority of PFAS (118/160) were inactive or equivocal in the DNT NAMs, leaving 42 active PFAS that decreased measures of neural network connectivity and neurite length. Analytical quality control indicated 43/118 inactive PFAS samples and 10/42 active PFAS samples were degraded; as such, careful interpretation is required as some negatives may have been due to loss of the parent PFAS, and some actives may have resulted from a mixture of parent and/or degradants of PFAS. PFAS containing a perfluorinated carbon (C) chain length ≥8, a high C:fluorine ratio, or a carboxylic acid moiety were more likely to be bioactive in the DNT NAMs. Of the PFAS positives in DNT NAMs, 85% were also active in other EPA ToxCast assays, whereas 79% of PFAS inactives in the DNT NAMs were active in other assays. These data demonstrate that a subset of PFAS perturb neurodevelopmental processes in vitro and suggest focusing future studies of DNT on PFAS with certain structural feature descriptors.

A New CSRML Structure-Based Fingerprint Method for Profiling and Categorizing Per- and Polyfluoroalkyl Substances (PFAS)

The term PFAS encompasses diverse per- and polyfluorinated alkyl (and increasingly aromatic) chemicals spanning industrial processes, commercial uses, environmental occurrence, and potential concerns. With increased chemical curation, currently exceeding 14,000 structures in the PFASSTRUCTV5 inventory on EPA's CompTox Chemicals Dashboard, has come increased motivation to profile, categorize, and analyze the PFAS structure space using modern cheminformatics approaches. Making use of the publicly available ToxPrint chemotypes and ChemoTyper application, we have developed a new PFAS-specific fingerprint set consisting of 129 TxP_PFAS chemotypes coded in CSRML, a chemical-based XML-query language. These are split into two groups, the first containing 56 mostly bond-type ToxPrints modified to incorporate attachment to either a CF group or F atom to enforce proximity to the fluorinated portion of the chemical. This focus resulted in a dramatic reduction in TxP_PFAS chemotype counts relative to the corresponding ToxPrint counts (averaging 54%). The remaining TxP_PFAS chemotypes consist of various lengths and types of fluorinated chains, rings, and bonding patterns covering indications of branching, alternate halogenation, and fluorotelomers. Both groups of chemotypes are well represented across the PFASSTRUCT inventory. Using the ChemoTyper application, we show how the TxP_PFAS chemotypes can be visualized, filtered, and used to profile the PFASSTRUCT inventory, as well as to construct chemically intuitive, structure-based PFAS categories. Lastly, we used a selection of expert-based PFAS categories from the OECD Global PFAS list to evaluate a small set of analogous structure-based TxP_PFAS categories. TxP_PFAS chemotypes were able to recapitulate the expert-based PFAS category concepts based on clearly defined structure rules that can be computationally implemented and reproducibly applied to process large PFAS inventories without need to consult an expert. The TxP_PFAS chemotypes have the potential to support computational modeling, harmonize PFAS structure-based categories, facilitate communication, and allow for more efficient and chemically informed exploration of PFAS chemicals moving forward.

Targeted Per- and Polyfluoroalkyl substances (PFAS) assessments for high throughput screening: Analytical and testing considerations to inform a PFAS stock quality evaluation framework

Per- and polyfluoroalkyl substances (PFAS) represent a large chemical class lacking hazard, toxicokinetic, and exposure information. To accelerate PFAS hazard evaluation, new approach methodologies (NAMs) comprised of in vitro high-throughput toxicity screening, toxicokinetic data, and computational modeling are being employed in read across strategies to evaluate the larger PFAS landscape. A critical consideration to ensure robust evaluations is a parallel assessment of the quality of the screening stock solutions, where dimethyl sulfoxide (DMSO) is often the diluent of choice. Challenged by the lack of commercially available reference standards for many of the selected PFAS and reliance on mass spectrometry approaches for such an evaluation, we developed a high-throughput framework to evaluate the quality of screening stocks for 205 PFAS selected for these NAM efforts. Using mass spectrometry coupled with either liquid or gas chromatography, a quality scoring system was developed that incorporated observations during mass spectral examination to provide a simple pass or fail notation. Informational flags were used to further describe findings regarding parent analyte presence through accurate mass identification, evidence of contaminants and/or degradation, or further describe characteristics such as isomer presence. Across the PFAS-DMSO stocks tested, 148 unique PFAS received passing quality scores to allow for further in vitro testing whereas 57 received a failing score primarily due to detection issues or confounding effects of DMSO. Principle component analysis indicated vapor pressure and Henry's Law Constant as top indicators for a failed quality score for those analyzed by gas chromatography. Three PFAS in the hexafluoropropylene oxide family failed due to degradation in DMSO. As the PFAS evaluated spanned over 20 different structural categories, additional commentary describes analytical observations across specific groups related to PFAS stock composition, detection, stability, and methodologic considerations that will be useful for informing future analytical assessment and downstream HTS efforts. The high-throughput stock quality scoring workflow presented holds value as a tool to evaluate chemical presence and quality efficiently and for informing data inclusion in PFAS or other NAM screening efforts.