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

An in vitro-in silico workflow for predicting renal clearance of PFAS

Per- and poly-fluoroalkyl substances (PFAS) have a wide range of elimination half-lives (days to years) in humans, thought to be in part due to variation in proximal tubule reabsorption. While human biomonitoring studies provide important data for some PFAS, renal clearance (CLrenal) predictions for hundreds of PFAS in commerce requires experimental studies with in vitro models and physiologically-based in vitro-to-in vivo extrapolation (IVIVE). Options for studying renal proximal tubule pharmacokinetics include cultures of renal proximal tubule epithelial cells (RPTECs) and/or microphysiological systems. This study aimed to compare CLrenal predictions for PFAS using in vitro models of varying complexity (96-well plates, static 24-well Transwells and a fluidic microphysiological model, all using human telomerase reverse transcriptase-immortalized and OAT1-overexpressing RPTECs combined with in silico physiologically-based IVIVE. Three PFAS were tested: one with a long half-life (PFOS) and two with shorter half-lives (PFHxA and PFBS). PFAS were added either individually (5 μM) or as a mixture (2 μM of each substance) for 48 h. Bayesian methods were used to fit concentrations measured in media and cells to a three-compartmental model to obtain the in vitro permeability rates, which were then used as inputs for a physiologically-based IVIVE model to estimate in vivo CLrenal. Our predictions for human CLrenal of PFAS were highly concordant with available values from in vivo human studies. The relative values of CLrenal between slow- and faster-clearance PFAS were most highly concordant between predictions from 2D culture and corresponding in vivo values. However, the predictions from the more complex model (with or without flow) exhibited greater concordance with absolute CLrenal. Overall, we conclude that a combined in vitro-in silico workflow can predict absolute CLrenal values, and effectively distinguish between PFAS with slow and faster clearance, thereby allowing prioritization of PFAS with a greater potential for bioaccumulation in humans.

Emerging organic contaminants in drinking water systems: Human intake, emerging health risks, and future research directions

Few earlier reviews on emerging organic contaminants (EOCs) in drinking water systems (DWS) focused on their detection, behaviour, removal and fate. Reviews on multiple exposure pathways, human intake estimates, and health risks including toxicokinetics, and toxicodynamics of EOCs in DWS are scarce. This review presents recent advances in human intake and health risks of EOCs in DWS. First, an overview of the evidence showing that DWS harbours a wide range of EOCs is presented. Multiple human exposure to EOCs occurs via ingestion of drinking water and beverages, inhalation and dermal pathways are discussed. A potential novel exposure may occur via the intravenous route in dialysis fluids. Analysis of global data on pharmaceutical pollution in rivers showed that the cumulative concentrations (μg L-1) of pharmaceuticals (mean ± standard error of the mean) were statistically more than two times significantly higher (p = 0.011) in South America (11.68 ± 5.29), Asia (9.97 ± 3.33), Africa (9.48 ± 2.81) and East Europe (8.09 ± 4.35) than in high-income regions (2.58 ± 0.48). Maximum cumulative concentrations of pharmaceuticals (μg L-1) decreased in the order; Asia (70.7) had the highest value followed by South America (68.8), Africa (51.3), East Europe (32.0) and high-income regions (17.1) had the least concentration. The corresponding human intake via ingestion of untreated river water was also significantly higher in low- and middle-income regions than in their high-income counterparts. For each region, the daily intake of pharmaceuticals was highest in infants, followed by children and then adults. A critique of the human health hazards, including toxicokinetics and toxicodynamics of EOCs is presented. Emerging health hazards of EOCs in DWS include; (1) long-term latent and intergenerational effects, (2) the interactive health effects of EOC mixtures, (3) the challenges of multifinality and equifinality, and (4) the Developmental Origins of Health and Disease hypothesis. Finally, research needs on human health hazards of EOCs in DWS are presented.

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.

A Comparison of In Vitro Points of Departure with Human Blood Levels for Per- and Polyfluoroalkyl Substances (PFAS)

Per- and polyfluoroalkyl substances (PFAS) are widely used, and their fluorinated state contributes to unique uses and stability but also long half-lives in the environment and humans. PFAS have been shown to be toxic, leading to immunosuppression, cancer, and other adverse health outcomes. Only a small fraction of the PFAS in commerce have been evaluated for toxicity using in vivo tests, which leads to a need to prioritize which compounds to examine further. Here, we demonstrate a prioritization approach that combines human biomonitoring data (blood concentrations) with bioactivity data (concentrations at which bioactivity is observed in vitro) for 31 PFAS. The in vitro data are taken from a battery of cell-based assays, mostly run on human cells. The result is a Bioactive Concentration to Blood Concentration Ratio (BCBCR), similar to a margin of exposure (MoE). Chemicals with low BCBCR values could then be prioritized for further risk assessment. Using this method, two of the PFAS, PFOA (Perfluorooctanoic Acid) and PFOS (Perfluorooctane Sulfonic Acid), have BCBCR values < 1 for some populations. An additional 9 PFAS have BCBCR values < 100 for some populations. This study shows a promising approach to screening level risk assessments of compounds such as PFAS that are long-lived in humans and other species.

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.

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.

Risk-based prioritization of PFAS using phenotypic and transcriptomic data from human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes

Per- and polyfluoroalkyl substances (PFAS) are chemicals with important applications; they are persistent in the environment and may pose human health hazards. Regulatory agencies are con­sidering restrictions and bans of PFAS; however, little data exists for informed decisions. Several prioritization strategies were proposed for evaluation of potential hazards of PFAS. Structure-based grouping could expedite the selection of PFAS for testing; still, the hypothesis that structure-effect relationships exist for PFAS requires confirmation. We tested 26 structurally diverse PFAS from 8 groups using human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes, and tested concentration-response effects on cell function and gene expression. Few phenotypic effects were observed in hepatocytes, but negative chronotropy was observed in cardiomyocytes for 8 PFAS. Substance- and cell type-dependent transcriptomic changes were more prominent but lacked substantial group-specific effects. In hepatocytes, we found upregulation of stress-related and extracellular matrix organization pathways, and down-regulation of fat metabolism. In car­diomyocytes, contractility-related pathways were most affected. We derived phenotypic and transcriptomic points of departure and compared them to predicted PFAS exposures. Conservative estimates for bioactivity and exposure were used to derive a bioactivity-to-exposure ratio (BER) for each PFAS; 23 of 26 PFAS had BER > 1. Overall, these data suggest that structure-based PFAS grouping may not be sufficient to predict their biological effects. Testing of individual PFAS may be needed for scientifically-supported decision-making. Our proposed strategy of using two human cell types and considering phenotypic and transcriptomic effects, combined with dose-response analysis and calculation of BER, may be used for PFAS prioritization.

Binding of Per- and Polyfluoroalkyl Substances (PFAS) to Serum Proteins: Implications for Toxicokinetics in Humans

Per- and polyfluoroalkyl substances (PFAS) are a diverse class of highly persistent anthropogenic chemicals that are detectable in the serum of most humans. PFAS exposure has been associated with many adverse effects on human health including immunotoxicity, increased risk of certain cancers, and metabolic disruption. PFAS binding to the most abundant blood serum proteins (human serum albumin [HSA] and globulins) is thought to affect transport to active sites, toxicity, and elimination half-lives. However, few studies have investigated the competitive binding of PFAS to these proteins in human serum. Here, we use C18 solid-phase microextraction fibers to measure HSA-water and globulin-water distribution coefficients (DHSA/w, Dglob/w) for PFAS with carbon chains containing 4 to 13 perfluorinated carbons (ηpfc = 4-13) and several functional head-groups. PFAS with ηpfc < 7 were highly bound to HSA relative to globulins, whereas PFAS with ηpfc ≥ 7 showed a greater propensity for binding to globulins. Experimentally measured DHSA/w and Dglob/w and concentrations of serum proteins successfully predicted the variability in PFAS binding in human serum. We estimated that the unbound fraction of serum PFAS varied by up to a factor of 2.5 among individuals participating in the 2017-2018 U.S. National Health and Nutrition Examination Survey. These results suggest that serum HSA and globulins are important covariates for epidemiological studies aimed at understanding the effects of PFAS exposure.

In Vitro and In Silico Analysis of the Bindings between Legacy and Novel Per- and Polyfluoroalkyl Substances and Human Serum Albumin

Per- and polyfluoroalkyl substances (PFASs) are emerging contaminants of concern that can enter the human body through a variety of pathways and thereby cause harmful effects. Exposure of pregnant women to PFASs could even affect both the mother and the child. Human serum albumin (HSA) is considered to be the primary transport protein for a variety of substances in body fluids. It can be bound to different contaminants and might result in possible effects on human health. Yet, few studies are available on the binding affinity of legacy PFASs and their novel alternatives to HSA. In this study, the binding mechanisms of HSA to both legacy PFASs and their novel alternatives were investigated using fluorescence spectroscopy, together with further molecular docking. The results show that all the target PFASs were statically quenched against HSA with binding ratios of 1:1. The binding constants of long-chain PFASs and novel alternatives of perfluoroalkanesulfonic acids (PFSAs) were greater than 102, whereas those of short-chain PFASs alternatives and novel alternatives of perfluorocarboxylic acids (PFCAs) were less than 102. In general, the binding affinities of PFCAs on HSA were less than that of PFSAs, while the binding affinities of short-chain PFASs alternatives on HSA were smaller than those of long-chain PFASs and their novel alternatives. Therefore, bindings to HSA could be considered as an important influencing factor for the bioaccumulation of legacy and novel PFASs in the human body.

Binding Affinity and Mechanism of Six PFAS with Human Serum Albumin: Insights from Multi-Spectroscopy, DFT and Molecular Dynamics Approaches

Per- and Polyfluoroalkyl Substances (PFAS) bioaccumulate in the human body, presenting potential health risks and cellular toxicity. Their transport mechanisms and interactions with tissues and the circulatory system require further investigation. This study investigates the interaction mechanisms of six PFAS with Human Serum Albumin (HSA) using multi-spectroscopy, DFT and a molecular dynamics approach. Multi-spectral analysis shows that perfluorononanoic acid (PFNA) has the best binding capabilities with HSA. The order of binding constants (298 K) is as follows: "Perfluorononanoic Acid (PFNA, 7.81 × 106 L·mol-1) > Perfluoro-2,5-dimethyl-3,6-dioxanonanoic Acid (HFPO-TA, 3.70 × 106 L·mol-1) > Perfluorooctanoic Acid (PFOA, 2.27 × 105 L·mol-1) > Perfluoro-3,6,9-trioxadecanoic Acid (PFO3DA, 1.59 × 105 L·mol-1) > Perfluoroheptanoic Acid (PFHpA, 4.53 × 103 L·mol-1) > Dodecafluorosuberic Acid (DFSA, 1.52 × 103 L·mol-1)". Thermodynamic analysis suggests that PFNA and PFO3DA's interactions with HSA are exothermic, driven primarily by hydrogen bonds or van der Waals interactions. PFHpA, DFSA, PFOA, and HFPO-TA's interactions with HSA, on the other hand, are endothermic processes primarily driven by hydrophobic interactions. Competitive probe results show that the main HSA-PFAS binding site is in the HSA structure's subdomain IIA. These findings are also consistent with the findings of molecular docking. Molecular dynamics simulation (MD) analysis further shows that the lowest binding energy (-38.83 kcal/mol) is fund in the HSA-PFNA complex, indicating that PFNA binds more readily with HSA. Energy decomposition analysis also indicates that van der Waals and electrostatic interactions are the main forces for the HSA-PFAS complexes. Correlation analysis reveals that DFT quantum chemical descriptors related to electrostatic distribution and characteristics like ESP and ALIE are more representative in characterizing HSA-PFAS binding. This study sheds light on the interactions between HSA and PFAS. It guides health risk assessments and control strategies against PFAS, serving as a critical starting point for further public health research.

Dose Response, Dosimetric, and Metabolic Evaluations of Replacement PFAS Perfluoro-(2,5,8-trimethyl-3,6,9-trioxadodecanoic) Acid (HFPO-TeA)

Few studies are available on the environmental and toxicological effects of perfluoroalkyl ether carboxylic acids (PFECAs), such as GenX, which are replacing legacy PFAS in manufacturing processes. To collect initial data on the toxicity and toxicokinetics of a longer-chain PFECA, male and female Sprague Dawley rats were exposed to perfluoro-(2,5,8-trimethyl-3,6,9-trioxadodecanoic) acid (HFPO-TeA) by oral gavage for five days over multiple dose levels (0.3-335.2 mg/kg/day). Clinically, we observed mortality at doses >17 mg/kg/day and body weight changes at doses ≤17 mg/kg/day. For the 17 mg/kg/day dose level, T3 and T4 thyroid hormone concentrations were significantly decreased (p < 0.05) from controls and HFPO-TeA plasma concentrations were significantly different between sexes. Non-targeted analysis of plasma and in vitro hepatocyte assay extractions revealed the presence of another GenX oligomer, perfluoro-(2,5-dimethyl-3,6-dioxanonanoic) acid (HFPO-TA). In vitro to in vivo extrapolation (IVIVE) parameterized with in vitro toxicokinetic data predicted steady-state blood concentrations that were within seven-fold of those observed in the in vivo study, demonstrating reasonable predictivity. The evidence of thyroid hormone dysregulation, sex-based differences in clinical results and dosimetry, and IVIVE predictions presented here suggest that the replacement PFECA HFPO-TeA induces a complex and toxic exposure response in rodents.

The contribution of PFAS to thyroid hormone-displacing activity in Dutch waters: A comparison between two in vitro bioassays with chemical analysis

Per- and polyfluoroalkyl substances (PFAS) are a group of xenobiotics that are widely distributed throughout the aquatic environment. Many PFAS are possible thyroid hormone (TH) system disrupting compounds, because they have the capacity to -amongst other- inhibit the TH thyroxine (T4) from binding to its transport protein transthyretin (TTR). This study investigated the occurrence of TH-displacing activity in the Dutch water cycle, and more specifically, the contribution of PFAS to this effect. Over one year of monitoring data of 29 PFAS (linear and branched) showed the continuous presence of PFAS in drinking waters and their surface water sources. Secondly, the FITC-T4 and TTR-TRβ-CALUX bioassays were mutually compared using positive (HPLC-grade water spiked with PFOA) and negative control samples (HPLC-grade water), as well as relative potency factors (RPFs) of up to 20 PFAS congeners. Both assays were found to be suitable for measuring TH-displacing activity in water samples. As a third aim, a field study was performed in the Dutch water cycle that was comprised of samples from drinking water, surface water, PFAS contaminated sites, and 2 wastewater treatment plants. All samples were analyzed with 1. chemical analysis for 29 PFAS, 2. the FITC-T4 bioassay, and 3. the TTR-TRβ-CALUX bioassay. The bioassays mutually showed good correlation (R2 0.85). Bioanalytical equivalent concentrations (BEQ) based on chemically-determined concentrations and RPFs (BEQchem) revealed that analyzed PFAS only explained ≤4.1 % of their activity in water extracts measured by both bioassays (BEQbio). This indicated that as yet unknown compounds contribute to the majority of the measured TH-displacing activity. Moreover, water treatment processes (e.g. DW production from SW) showed a larger contribution of target PFAS to the BEQbio. This could be a first lead to identify unknown compounds that contribute to this activity, and as such, enable the assessment of possible risks associated by the occurrence of TH-displacing activity in water.

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.