Bayesian evaluation of a physiologically based pharmacokinetic (PBPK) model for perfluorooctane sulfonate (PFOS) to characterize the interspecies uncertainty between mice, rats, monkeys, and humans: Development and performance verification

Hexafluoropropylene Oxide Trimer Acid Is an Unsafe Substitute to Perfluorooctanoic Acid Due to Its Remarkable Liver Accumulation in Mice Disclosed by Comprehensive Toxicokinetic Models

Hexafluoropropylene oxide trimer acid (HFPO-TA, C2F5(CF2OCF(CF3))2COOH) is widely used as an alternative to perfluorooctanoic acid (PFOA), but whether it is a safe alternative requires further evaluation. In this study, male mice were exposed to three dosages (0.56, 2.8, and 14 mg/kg) of HFPO-TA via single oral gavage or intravenous injection for 28 days. HFPO-TA was rapidly absorbed into the blood and tissues within 15 min postexposure, with a volume of distribution approximately 3 times higher than PFOA, indicating a greater propensity for tissue distribution. Notably, HFPO-TA was distinctly more accumulated in liver compared to plasma and other tissues and very poorly excreted, with only 2.23% in urine and 7.26% in feces on the 21st day after oral exposure. A physiologically based toxicokinetic model, extrapolated to long-term low-dose exposure, revealed a lower bile clearance rate (8-fold) and higher liver partition coefficient (7-fold) than PFOA, and a higher hepatic first-pass effect of HFPO-TA (5-fold) than PFOA, contributing to its remarkable liver accumulation (5-fold). Molecular docking analysis reveals strong binding affinity of HFPO-TA with typical enterohepatic circulation transport proteins due to its strong hydrophobicity, flexible chain structure, and formation of additional hydrogen bonds, favoring HFPO-TA accumulation in the liver. The results suggest that HFPO-TA may not be a safe substitute for legacy PFAS, and further human exposure risk assessments are warranted.

An in vitro and machine learning framework for quantifying serum albumin binding of per- and polyfluoroalkyl substances

Per- and polyfluoroalkyl substances (PFAS) are a diverse class of anthropogenic chemicals; many are persistent, bioaccumulative, and mobile in the environment. Worldwide, PFAS bioaccumulation causes serious adverse health impacts, yet the physiochemical determinants of bioaccumulation and toxicity for most PFAS are not well understood, largely due to experimental data deficiencies. As most PFAS are proteinophilic, protein binding is a critical parameter for predicting PFAS bioaccumulation and toxicity. Among these proteins, human serum albumin (HSA) is the predominant blood transport protein for many PFAS. We previously demonstrated the utility of an in vitro differential scanning fluorimetry assay for determining relative HSA binding affinities for 24 PFAS. Here, we report HSA affinities for 65 structurally diverse PFAS from 20 chemical classes. We leverage these experimental data, and chemical/molecular descriptors of PFAS, to build 7 machine learning classifier algorithms and 9 regression algorithms, and evaluate their performance to identify the best predictive binding models. Evaluation of model accuracy revealed that the top-performing classifier model, logistic regression, had an AUROC (area under the receiver operating characteristic curve) statistic of 0.936. The top-performing regression model, support vector regression, had an R2 of 0.854. These top-performing models were then used to predict HSA-PFAS binding for chemicals in the EPAPFASINV list of 430 PFAS. These developed in vitro and in silico methodologies represent a high-throughput framework for predicting protein-PFAS binding based on empirical data, and generate directly comparable binding data of potential use in predictive modeling of PFAS bioaccumulation and other toxicokinetic endpoints.

A physiologically based toxicokinetic model for microplastics and nanoplastics in mice after oral exposure and its implications for human dietary exposure assessment

Evidence of microplastics (MPs) and nanoplastics (NPs) in foods and daily-use products, along with their frequent detection in the human body, has raised concerns regarding their potential impact on human health through dietary ingestion. However, there is a lack of quantitative tools to simulate their bioaccumulation and tissue distribution following environmental exposure. To address this gap, we developed the first physiologically based toxicokinetic (PBTK) model for predicting the biodistribution of MPs and NPs in mice following oral exposure under various exposure scenarios. This novel model incorporated key kinetic mass transport processes, such as membrane permeability, albumin binding, and cellular uptake. We identified that the absorption rate in the gastrointestinal tract and fecal excretion rate constant had significant impacts on organ dosimetry. Our regression analysis indicated that the size-dependent dissociation constant and urine clearance rate constant sharply increased by a factor of 3 as NPs particle size increased to 1 µm. Finally, we developed a graphical user interface to enable interactive visualization and analysis for future applications, supporting human dietary exposure and risk assessment using available food consumption data and MPs/NPs residue data. The simulation results offer a mechanistic perspective, enhancing understanding of the internal organ dosimetry burden and health impacts from dietary exposure to MPs and NPs.

Clinical, histological, molecular, and toxicokinetic renal outcomes of per-/polyfluoroalkyl substances (PFAS) exposure: Systematic review and meta-analysis

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are synthetic chemicals present in the environment that can negatively affect health. Kidney is the major target organ of PFAS exposure, yet the renal impact of PFAS is not completely understood. Here we review the effects of PFAS exposure on kidney health to identify gaps in our understanding and mark potential avenues for future research.

METHODS

PubMed and SCOPUS databases were searched for studies that examined the association between PFAS exposure and kidney-related outcomes. We included all epidemiological, animal, and cell studies and categorized outcomes into four categories: clinical, histological, molecular and toxicokinetic.

RESULTS

We identified 169 studies, including 51 on clinical outcomes, 28 on histological changes, 42 on molecular mechanisms, and 68 on toxicokinetics. Perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) exposure were associated with kidney dysfunction, chronic kidney diseases, and increased risk of kidney cancer. Various histological changes were reported, especially in tubular epithelial cells, and the etiology of PFAS-induced kidney injury included various molecular mechanisms. Although PFOA and PFOS are not considered genotoxic, they exhibit several characteristics of carcinogens. Toxicokinetics of PFOA and PFOS differed significantly between species, with renal elimination influenced by various factors such as sex, age, and structure of the compound.

CONCLUSION

Evidence suggests that PFAS, especially PFOA and PFOS, negatively affects kidney health, though gaps in our understanding of such effects call for further research.

Development of a Physiologically Based Pharmacokinetic (PBPK) Model for F-53B in Pregnant Mice and Its Extrapolation to Humans

Chlorinated polyfluorinated ether sulfonic acid (F-53B), a commonly utilized alternative for perfluorooctane sulfonate, was detected in pregnant women and cord blood recently. However, the lack of detailed toxicokinetic information poses a significant challenge in assessing the human risk assessment for F-53B exposure. Our study aimed to develop a physiologically based pharmacokinetic (PBPK) model for pregnant mice, based on toxicokinetic experiments, and extrapolating it to humans. Pregnant mice were administered 80 μg/kg F-53B orally and intravenously on gestational day 13. F-53B concentrations in biological samples were analyzed via ultraperformance liquid chromatography-mass spectrometry. Results showed the highest F-53B accumulation in the brain, followed by the placenta, amniotic fluid, and liver in fetal mice. These toxicokinetic data were applied to F-53B PBPK model development and evaluation, and Monte Carlo simulations were used to characterize the variability and uncertainty in the human population. Most of the predictive values were within a 2-fold range of experimental data (>72%) and had a coefficient of determination (R2) greater than 0.68. The developed mouse model was then extrapolated to the human and evaluated with human biomonitoring data. Our study provides an important step toward improving the understanding of toxicokinetics of F-53B and enhancing the quantitative risk assessments in sensitive populations, particularly in pregnant women and fetuses.

Contribution of Continued Dermal Exposure of PFAS-Containing Sunscreens to Internal Exposure: Extrapolation from In Vitro and In Vivo Tests to Physiologically Based Toxicokinetic Models

Per- and polyfluoroalkyl substances (PFASs) are widely present in sunscreen products as either active ingredients or impurities. They may penetrate the human skin barrier and then pose potential health risks. Herein, we aimed to develop a physiologically based toxicokinetic (PBTK) model capable of predicting the body loading of PFASs after repeated, long-term dermal application of commercial sunscreens. Ten laboratory-prepared sunscreens, generally falling into two categories of water-in-oil (W/O) and oil-in-water (O/W) sunscreens, were subject to in vitro percutaneous penetration test to assess the impacts of four sunscreen ingredients on PFAS penetration. According to the results, two sunscreen formulas representing W/O and O/W types that mostly enhanced PFAS dermal absorption were then selected for a subsequent 30 day in vivo exposure experiment in mice. PBTK models were successfully established based on the time-dependent PFAS concentrations in mouse tissues (R2 = 0.885-0.947) and validated through another 30 day repeated exposure experiment in mice using two commercially available sunscreens containing PFASs (R2 = 0.809-0.835). The PBTK model results suggest that applying sunscreen of the same amount on a larger skin area is more conducive to PFAS permeation, thus enhancing the exposure risk. This emphasizes the need for caution in practical sunscreen application scenarios, particularly during the summer months.

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.

Bayesian analysis of physiologically based toxicokinetic (PBTK) modeling for pentachlorophenol exposure in pregnant women

Pentachlorophenol (PCP) is a persistent organic compound that is widely present in the environment. The estimation of internal exposure levels for a given external exposure using toxicokinetic models is key to the human health risk assessment of PCP. The present study developed a physiologically based multicompartmental pharmacokinetic (PBTK) model to describe and predict the behavior of pentachlorophenol (PCP) in an organism. The model consists of stomach, intestines, adipose tissue, kidneys and fast- and poorly perfused tissues that are interconnected via blood circulation. We constructed a PBTK model of PCP in rats and extrapolated it to human dietary PCP exposure. The toxicokinetic data of PCP in human tissues and excreta were obtained from the published literature. Based on the collected PCP dietary survey and internal exposure data of pregnant women in Shanghai, Bayesian statistical analysis was performed for the model using Markov chain Monte Carlo (MCMC) simulation. The posterior distributions of the sensitive parameters were estimated, and the model was parameter optimized and validated using the pregnant women's test dataset. The results showed that the root mean square error (RMSE) improved 37.3% compared to the original model, and a systematic literature search revealed that the optimized model achieved acceptable prediction results for other datasets in China. A PCP metabolism model based on the exposure characteristics of pregnant women in China was constructed in the present study. The model provides a theoretical basis for the study of PCP toxicity and risk assessment.

Utility of life stage-specific chemical risk assessments based on New Approach Methodologies (NAMs)

The safety assessments for chemicals targeted for use or expected to be exposed to specific life stages, including infancy, childhood, pregnancy and lactation, and geriatrics, need to account for extrapolation of data from healthy adults to these populations to assess their human health risk. However, often adequate and relevant toxicity or pharmacokinetic (PK) data of chemicals in specific life stages are not available. For such chemicals, New Approach Methodologies (NAMs), such as physiologically based pharmacokinetic (PBPK) modeling, biologically based dose response (BBDR) modeling, in vitro to in vivo extrapolation (IVIVE), etc. can be used to understand the variability of exposure and effects of chemicals in specific life stages and assess their associated risk. A life stage specific PBPK model incorporates the physiological and biochemical changes associated with each life stage and simulates their impact on the absorption, distribution, metabolism, and elimination (ADME) of these chemicals. In our review, we summarize the parameterization of life stage models based on New Approach Methodologies (NAMs) and discuss case studies that highlight the utility of a life stage based PBPK modeling for risk assessment. In addition, we discuss the utility of artificial intelligence (AI)/machine learning (ML) and other computational models, such as those based on in vitro data, as tools for estimation of relevant physiological or physicochemical parameters and selection of model. We also discuss existing gaps in the available toxicological datasets and current challenges that need to be overcome to expand the utility of NAMs for life stage-specific chemical risk assessment.

Physiologically based kinetic (PBK) modeling of propiconazole using a machine learning-enhanced read-across approach for interspecies extrapolation

A significant challenge in the traditional human health risk assessment of agrochemicals is the uncertainty in quantifying the interspecies differences between animal models and humans. To work toward a more accurate and animal-free risk determination, new approaches such as physiologically based kinetic (PBK) modeling have been used to perform dosimetry extrapolation from animals to humans. However, the regulatory use and acceptance of PBK modeling is limited for chemicals that lack in vivo animal pharmacokinetic (PK) data, given the inability to evaluate models. To address these challenges, this study developed PBK models in the absence of in vivo PK data for the fungicide propiconazole, an activator of constitutive androstane receptor (CAR)/pregnane X receptor (PXR). A fit-for-purpose read-across approach was integrated with hierarchical clustering - an unsupervised machine learning algorithm, to bridge the knowledge gap. The integration allowed the incorporation of a broad spectrum of attributes for analog consideration, and enabled the analog selection in a simple, reproducible, and objective manner. The applicability was evaluated and demonstrated using penconazole (source) and three pseudo-unknown target chemicals (epoxiconazole, tebuconazole and triadimefon). Applying this machine learning-enhanced read-across approach, difenoconazole was selected as the most appropriate analog for propiconazole. A mouse PBK model was developed and evaluated for difenoconazole (source), with the mode of action of CAR/PXR activation incorporated to simulate the in vivo autoinduction of metabolism. The difenoconazole mouse model then served as a template for constructing the propiconazole mouse model. A parallelogram approach was subsequently applied to develop the propiconazole rat and human models, enabling a quantitative assessment of interspecies differences in dosimetry. This integrated approach represents a substantial advancement toward refining risk assessment of propiconazole within the framework of animal alternative safety assessment strategies.

Species Difference? Bovine, Trout, and Human Plasma Protein Binding of Per- and Polyfluoroalkyl Substances

Per- and polyfluoroalkyl substances (PFAS) strongly bind to proteins and lipids in blood, which govern their accumulation and distribution in organisms. Understanding the plasma binding mechanism and species differences will facilitate the quantitative in vitro-to-in vivo extrapolation and improve risk assessment of PFAS. We studied the binding mechanism of 16 PFAS to bovine serum albumin (BSA), trout, and human plasma using solid-phase microextraction. Binding of anionic PFAS to BSA and human plasma was found to be highly concentration-dependent, while trout plasma binding was linear for the majority of the tested PFAS. At a molar ratio of PFAS to protein ν < 0.1 molPFAS/molprotein, the specific protein binding of anionic PFAS dominated their human plasma binding. This would be the scenario for physiological conditions (ν < 0.01), whereas in in vitro assays, PFAS are often dosed in excess (ν > 1) and nonspecific binding becomes dominant. BSA was shown to serve as a good surrogate for human plasma. As trout plasma contains more lipids, the nonspecific binding to lipids affected the affinities of PFAS for trout plasma. Mass balance models that are parameterized with the protein-water and lipid-water partitioning constants (chemical characteristics), as well as the protein and lipid contents of the plasma (species characteristics), were successfully used to predict the binding to human and trout plasma.

Construction of a physiologically based pharmacokinetic model of paclobutrazol and exposure estimation in the human body

Paclobutrazol (PBZ) is a plant growth regulator that can delay plant growth and improve plant resistance and yield. Although it has been widely used in the growth of medicinal plants, human beings may take it by taking traditional Chinese medicine. There are no published studies on PBZ exposure in humans or standardized limits for PBZ in medicinal plants. We measured the solubility, oil-water partition coefficient (logP), and pharmacokinetics of PBZ in rats and established a physiologically based pharmacokinetic (PBPK) model of PBZ in rats. This was followed by extrapolation to healthy Chinese adult males as a theoretical foundation for future risk assessment of PBZ. The results showed that PBZ had low solubility and high fat solubility. Pharmacokinetic experiments showed that PBZ was absorbed rapidly but eliminated slowly in rats. On this basis, the rat PBPK model was successfully constructed and extrapolated to healthy Chinese adult males to predict the plasma concentration-time curve and exposure of PBZ in humans. The construction of the PBPK model of PBZ in this study facilitates the determination of the standard formulation limits and risk assessment of PBZ residues in medicinal plants.

Hepatic injury and ileitis associated with gut microbiota dysbiosis in mice upon F-53B exposure

Chlorinated polyfluorinated ether sulfonate (F-53B), a substitute of perfluorooctane sulfonic acid (PFOS), has attracted significant attention for its link to hepatotoxicity and enterotoxicity. Nevertheless, the underlying mechanisms of F-53B-induced enterohepatic toxicity remain incompletely understood. This study aimed to explore the role of F-53B exposure on enterohepatic injury based on the gut microbiota, pathological and molecular analysis in mice. Here, we exposed C57BL/6 mice to F-53B (0, 4, 40, and 400 μg/L) for 28 days. Our findings revealed a significant accumulation of F-53B in the liver, followed by small intestines, and feces. In addition, F-53B induced pathological collagen fiber deposition and lipoid degeneration, up-regulated the expression of fatty acid β-oxidation-related genes (PPARα and PPARγ, etc), while simultaneously down-regulating pro-inflammatory genes (Nlrp3, IL-1β, and Mcp1) in the liver. Meanwhile, F-53B induced ileal mucosal barrier damage, and an up-regulation of pro-inflammatory genes and mucosal barrier-related genes (Muc1, Muc2, Claudin1, Occludin, Mct1, and ZO-1) in the ileum. Importantly, F-53B distinctly altered gut microbiota compositions by increasing the abundance of Akkermansia and decreasing the abundance of Prevotellaceae_NK3B31_group in the feces. F-53B-altered microbiota compositions were significantly associated with genes related to fatty acid β-oxidation, inflammation, and mucosal barrier. In summary, our results demonstrate that F-53B is capable of inducing hepatic injury, ileitis, and gut microbiota dysbiosis in mice, and the gut microbiota dysbiosis may play an important role in the F-53B-induced enterohepatic toxicity.

Development and application of the physiologically-based toxicokinetic (PBTK) model for ochratoxin A (OTA) in rats and humans

Ochratoxin A (OTA) is a common fungal toxin frequently detected in food and human plasma samples. Currently, the physiologically based toxicokinetic (PBTK) model plays an active role in dose translation and can improve and enhance the risk assessment of toxins. In this study, the PBTK model of OTA in rats and humans was established based on knowledge of OTA-specific absorption, distribution, metabolism, and excretion (ADME) in order to better explain the disposition of OTA in humans and the discrepancies with other species. The models were calibrated and optimized using the available kinetic and toxicokinetic (TK) data, and independent test datasets were used for model evaluation. Subsequently, sensitivity analyses and population simulations were performed to characterize the extent to which variations in physiological and specific chemical parameters affected the model output. Finally, the constructed models were used for dose extrapolation of OTA, including the rat-to-human dose adjustment factor (DAF) and the human exposure conversion factor (ECF). The results showed that the unbound fraction (Fup) of OTA in plasma of rat and human was 0.02-0.04% and 0.13-4.21%, respectively. In vitro experiments, the maximum enzyme velocity (Vmax) and Michaelis-Menten constant (Km) of OTA in rat and human liver microsomes were 3.86 and 78.17 μg/g min-1, 0.46 and 4.108 μg/mL, respectively. The predicted results of the model were in good agreement with the observed data, and the models in rats and humans were verified. The PBTK model derived a DAF of 0.1081 between rats and humans, whereas the ECF was 2.03. The established PBTK model can be used to estimate short- or long-term OTA exposure levels in rats and humans, with the capacity for dose translation of OTA to provide the underlying data for risk assessment of OTA.

Quantitative cross-species comparison of serum albumin binding of per- and polyfluoroalkyl substances from five structural classes

Per- and polyfluoroalkyl substances (PFAS) are a class of over 8000 chemicals, many of which are persistent, bioaccumulative, and toxic to humans, livestock, and wildlife. Serum protein binding affinity is instrumental in understanding PFAS toxicity, yet experimental binding data is limited to only a few PFAS congeners. Previously, we demonstrated the usefulness of a high-throughput, in vitro differential scanning fluorimetry assay for determination of relative binding affinities of human serum albumin for 24 PFAS congeners from 6 chemical classes. In the current study, we used this assay to comparatively examine differences in human, bovine, porcine, and rat serum albumin binding of 8 structurally informative PFAS congeners from 5 chemical classes. With the exception of the fluorotelomer alcohol 1H, 1H, 2H, 2H-perfluorooctanol (6:2 FTOH), each PFAS congener bound by human serum albumin was also bound by bovine, porcine, and rat serum albumin. The critical role of the charged functional headgroup in albumin binding was supported by the inability of albumin of each species tested to bind 6:2 FTOH. Significant interspecies differences in serum albumin binding affinities were identified for each of the bound PFAS congeners. Relative to human albumin, perfluoroalkyl carboxylic and sulfonic acids were bound with greater affinity by porcine and rat serum albumin, and the perfluoroalkyl ether acid congener bound with lower affinity to porcine and bovine serum albumin. These comparative affinity data for PFAS binding by serum albumin from human, experimental model, and livestock species reduce critical interspecies uncertainty and improve accuracy of predictive bioaccumulation and toxicity assessments for PFAS.

A physiologically based pharmacokinetic/pharmacodynamic model to determine dosage regimens and withdrawal intervals of aditoprim against Streptococcus suis

Introduction: Streptococcus suis (S. suis) is a zoonotic pathogen threatening public health. Aditoprim (ADP), a novel veterinary medicine, exhibits an antibacterial effect against S. suis. In this study, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model was used to determine the dosage regimens of ADP against S. suis and withdrawal intervals. Methods: The PBPK model of ADP injection can predict drug concentrations in plasma, liver, kidney, muscle, and fat. A semi-mechanistic pharmacodynamic (PD) model, including susceptible subpopulation and resistant subpopulation, is successfully developed by a nonlinear mixed-effect model to evaluate antibacterial effects. An integrated PBPK/PD model is conducted to predict the time-course of bacterial count change and resistance development under different ADP dosages. Results: ADP injection, administrated at 20 mg/kg with 12 intervals for 3 consecutive days, can exert an excellent antibacterial effect while avoiding resistance emergence. The withdrawal interval at the recommended dosage regimen is determined as 18 days to ensure food safety. Discussion: This study suggests that the PBPK/PD model can be applied as an effective tool for the antibacterial effect and safety evaluation of novel veterinary drugs.

Estimation of per- and polyfluoroalkyl substances (PFAS) half-lives in human studies: a systematic review and meta-analysis

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) constitute a heterogeneous group of synthetic compounds widely used in industrial applications. The estimation of PFAS half-life (t1/2) is essential to quantify their persistence, their toxicity and mechanism of action in humans.

OBJECTIVES

The purpose of this review is to summarize the evidence on PFAS half-lives in humans from the available literature, and to investigate the limitations and uncertainties characterizing half-life estimation.

METHODS

The search was conducted on PubMed, Scopus, and Embase databases up to July 03, 2023 and was aimed at identifying all papers that estimated PFAS half-life in human populations. We excluded studies on temporal trends or providing estimates of half-life based solely on renal clearance. As persistent and ongoing exposures can influence half-life estimation, we decided to include only studies that were conducted after the main source of exposure to PFAS had ceased. A random-effects meta-analysis was conducted on studies that reported perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid (PFOS) or perfluorohexanesulfonic acid (PFHxS) half-life estimation. Risk of bias was evaluated using the OHAT tool.

RESULTS

A total of 13 articles were included in the review, with 5 studies conducted in exposed general populations and 8 studies conducted in exposed workers; the estimated mean half-life ranged from 1.48 to 5.1 years for PFOA, from 3.4 to 5.7 years for total PFOS, and from 2.84 to 8.5 years for PFHxS. High heterogeneity among studies was observed; potential reasons include the variability among the investigated populations, discrepancies in considering ongoing exposures, variability in PFAS isomeric compositions, accounting for background exposure, time since exposure stopped and methods used for half-life estimation.

DISCUSSION

Despite the efforts made to better understand PFAS toxicokinetics, further studies are needed to identify important characteristics of these persistent chemicals. Biomonitoring studies should focus on persistent and unaccounted sources of exposure to PFAS and on individual characteristics potentially determining half-life, to ensure accurate estimates.

PBPK Modeling as an Alternative Method of Interspecies Extrapolation that Reduces the Use of Animals: A Systematic Review

INTRODUCTION

Physiologically based pharmacokinetic (PBPK) modeling is a computational approach that simulates the anatomical structure of the studied species and presents the organs and tissues as compartments interconnected by arterial and venous blood flows.

AIM

The aim of this systematic review was to analyze the published articles focused on the development of PBPK models for interspecies extrapolation in the disposition of drugs and health risk assessment, presenting to this modeling an alternative to reduce the use of animals.

METHODS

For this purpose, a systematic search was performed in PubMed using the following search terms: "PBPK" and "Interspecies extrapolation". The revision was performed according to PRISMA guidelines.

RESULTS

In the analysis of the articles, it was found that rats and mice are the most commonly used animal models in the PBPK models; however, most of the physiological and physicochemical information used in the reviewed studies were obtained from previous publications. Additionally, most of the PBPK models were developed to extrapolate pharmacokinetic parameters to humans and the main application of the models was for toxicity testing.

CONCLUSION

PBPK modeling is an alternative that allows the integration of in vitro and in silico data as well as parameters reported in the literature to predict the pharmacokinetics of chemical substances, reducing in large quantity the use of animals that are required in traditional studies.

A physiologically based pharmacokinetic model of diethyl phthalates in humans

Phthalates are a family of industrial and consumer product chemicals, among which diethyl phthalate (DEP) has been widely used. DEP is metabolized into the active metabolite monoethyl phthalate (MEP) and exposure to DEP may induce male reproductive toxicity, developmental toxicity and hepatotoxicity. To better assess the toxicity of DEP and MEP, it is important to understand and predict their internal concentrations, especially in reproductive organs. Here we present a human physiologically based pharmacokinetic (PBPK) model of DEP. Implemented in R, the PBPK model consists of seven tissue compartments, including blood, gut, liver, fat, skin, gonad, and rest of body (RB). In the blood both DEP and MEP partition into free and bound forms, and tissue distribution is considered as blood flow-limited. DEP is metabolized in the gut and liver into MEP which is further glucuronidated and cleared through the urine. The chemical-specific parameters of the model were predicted in silico or estimated based on published human urinary MEP data after exposure to DEP in the air at 250 or 300 μg/m3 for 3 or 6 h through inhalation and dermal absorption. Sensitivity analysis identified important parameters including partition coefficients of DEP for fat, RB, and skin compartments, and the rate constants for glucuronidation of MEP and urinary excretion, with regard to Cmax, area under the curve (AUC), and clearance half-lives of DEP and MEP. A subset of the sensitive parameters was then included in hierarchical population Bayesian Markov chain Monte Carlo (MCMC) simulations to characterize the uncertainty and variability of these parameters. The model is consistent with the notion that dermal absorption represents a significant route of exposure to DEP in ambient air and clothing can be an effective barrier. The developed human PBPK model can be utilized upon further refinement as a quantitative tool for DEP risk assessment.

Exploring Route-Specific Pharmacokinetics of PFAS in Mice by Coupling in Vivo Tests and Physiologically Based Toxicokinetic Models

BACKGROUND

Oral ingestion, inhalation, and skin contact are important exposure routes for humans to uptake per- and polyfluoroalkyl substances (PFAS). However, nasal and dermal exposure to PFAS remains unclear, and accurately predicting internal body burden of PFAS in humans via multiple exposure pathways is urgently required.

OBJECTIVES

We aimed to develop multiple physiologically based toxicokinetic (PBTK) models to unveil the route-specific pharmacokinetics and bioavailability of PFAS via respective oral, nasal, and dermal exposure pathways using a mouse model and sought to predict the internal concentrations in various tissues through multiple exposure routes and extrapolate it to humans.

METHODS

Mice were administered the mixed solution of perfluorohexane sulfonate, perfluorooctane sulfonate, and perfluorooctanoic acid through oral, nasal, and dermal exposure separately or jointly. The time-dependent concentrations of PFAS in plasma and tissues were determined to calibrate and validate the individual and combined PBTK models, which were applied in single- and repeated-dose scenarios.

RESULTS

The developed route-specific PBTK models successfully simulated the tissue concentrations of PFAS in mice following single or joint exposure routes as well as long-term repeated dose scenarios. The time to peak concentration of PFAS in plasma via dermal exposure was much longer (34.1-83.0 h) than that via nasal exposure (0.960 h). The bioavailability of PFAS via oral exposure was the highest (73.2%-98.0%), followed by nasal (33.9%-66.8%) and dermal exposure (4.59%-7.80%). This model was extrapolated to predict internal levels in human under real environment.

DISCUSSION

Based on these data, we predict the following: PFAS were absorbed quickly via nasal exposure, whereas a distinct hysteresis effect was observed for dermal exposure. Almost all the PFAS to which mice were exposed via gastrointestinal route were absorbed into plasma, which exhibited the highest bioavailability. Exhalation clearance greatly depressed the bioavailability of PFAS via nasal exposure, whereas the lowest bioavailability in dermal exposure was because of the interception of PFAS within the skin layers. https://doi.org/10.1289/EHP11969.

Quantitative Cross-Species Comparison of Serum Albumin Binding of Per- and Polyfluoroalkyl Substances from Five Structural Classes

Per- and polyfluoroalkyl substances (PFAS) are a class of over 8,000 chemicals that are persistent, bioaccumulative, and toxic to humans, livestock, and wildlife. Serum protein binding affinity is instrumental in understanding PFAS toxicity, yet experimental binding data is limited to only a few PFAS congeners. Previously, we demonstrated the usefulness of a high-throughput, in vitro differential scanning fluorimetry assay for determination of relative binding affinities of human serum albumin for 24 PFAS congeners from 6 chemical classes. In the current study, we used this differential scanning fluorimetry assay to comparatively examine differences in human, bovine, porcine, and rat serum albumin binding of 8 structurally informative PFAS congeners from 5 chemical classes. With the exception of the fluorotelomer alcohol 1H,1H,2H,2H-perfluorooctanol (6:2 FTOH), each PFAS congener bound by human serum albumin was also bound by bovine, porcine, and rat serum albumin. The critical role of the charged functional headgroup in albumin binding was supported by the inability of serum albumin of each species tested to bind 6:2 FTOH. Significant interspecies differences in serum albumin binding affinities were identified for each of the bound PFAS congeners. Relative to human albumin, perfluoroalkyl carboxylic and sulfonic acids were bound with greater affinity by porcine and rat serum albumin, and perfluoroalkyl ether congeners bound with lower affinity to porcine and bovine serum albumin. These comparative affinity data for PFAS binding by serum albumin from human, experimental model and livestock species reduce critical interspecies uncertainty and improve accuracy of predictive toxicity assessments for PFAS.

Development and application of an interactive generic physiologically based pharmacokinetic (igPBPK) model for adult beef cattle and lactating dairy cows to estimate tissue distribution and edible tissue and milk withdrawal intervals for per- and polyfluoroalkyl substances (PFAS)

Humans can be exposed to per- and polyfluoroalkyl substances (PFAS) through dietary intake from milk and edible tissues from food animals. This study developed a physiologically based pharmacokinetic (PBPK) model to predict tissue and milk residues and estimate withdrawal intervals (WDIs) for multiple PFAS including PFOA, PFOS and PFHxS in beef cattle and lactating dairy cows. Results showed that model predictions were mostly within a two-fold factor of experimental data for plasma, tissues, and milk with an estimated coefficient of determination (R2) of >0.95. The predicted muscle WDIs for beef cattle were <1 day for PFOA, 449 days for PFOS, and 69 days for PFHxS, while the predicted milk WDIs in dairy cows were <1 day for PFOA, 1345 days for PFOS, and zero day for PFHxS following a high environmental exposure scenario (e.g., 49.3, 193, and 161 ng/kg/day for PFOA, PFOS, and PFHxS, respectively, for beef cattle for 2 years). The model was converted to a web-based interactive generic PBPK (igPBPK) platform to provide a user-friendly dashboard for predictions of tissue and milk WDIs for PFAS in cattle. This model serves as a foundation for extrapolation to other PFAS compounds to improve safety assessment of cattle-derived food products.

Dynamic exposure and body burden models for per- and polyfluoroalkyl substances (PFAS) enable management of food safety risks in cattle

With increasing global focus on planetary boundaries, food safety and quality, the presence of per- and polyfluoroalkyl substances (PFAS) in the food chain presents a challenge for the sustainable production and supply of quality assured food. Consumption of food is the primary PFAS exposure route for the general population. At contaminated sites, PFAS have been reported in a range of agricultural commodities including cattle. Consumer exposure assessments are complicated by the lack of validated modelling approaches to estimate PFAS bioaccumulation in cattle. Previous studies have shown that PFAS bioaccumulation in livestock is influenced by environmental, spatial and temporal factors that necessitate a dynamic modelling approach. This work presents an integrated exposure and population toxicokinetic (PopTK) model for cattle that estimates serum and tissue concentrations of PFAS over time. Daily exposures were estimated from intakes of water, pasture, and soil, and considered animal growth, seasonal variability (pasture moisture content and temperature) and variable PFAS concentrations across paddocks. Modelled serum and tissue estimates were validated against monitoring data from Australian and Swedish cattle farms. The models were also used to develop and test practical management options for reducing PFAS exposure and to prioritise remediation for farms. Model outputs for exposure management scenarios (testing cattle rotation and targeted supplementation of feed and water) showed potential for marked reductions in consumer exposures from cattle produce.

PPARα/ACOX1 as a novel target for hepatic lipid metabolism disorders induced by per- and polyfluoroalkyl substances: An integrated approach

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are persistent and ubiquitous environmental contaminants with well-documented hepatotoxicity. However, the mechanistic linkage between PFAS exposure and non-alcoholic fatty liver disease (NAFLD) remains largely elusive.

OBJECTIVES

This study aimed to explore PFAS-to-NAFLD link and the relevant molecular mechanisms.

METHODS

The cross-sectional analyses using National Health and Nutrition Examination Survey (NHANES) data were conducted to investigate the association between PFAS exposure and NAFLD. A combination of in silico toxicological analyses, bioinformatics approaches, animal experiments, and in vitro assays was used to explore the molecular initiating events (MIEs) and key events (KEs) in PFAS-induced hepatic lipid metabolism disorders.

RESULTS

The cross-sectional analyses with NHANES data revealed the significant association between PFAS exposure and hepatic steatosis/NAFLD. The in silico toxicological analyses showed that PPARα activation induced by perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS), prototypical representatives of PFAS, is the critical MIE associated with NAFLD-predominant liver diseases. Transcriptome-based bioinformatic annotation and analyses identified that transcriptional upregulation of hepatic acyl-CoA oxidase 1 (ACOX1) in PPARα-regulated peroxisomal β-oxidation pathway was the KE involved with PFOA/PFOS-perturbed hepatic lipid metabolic pathways in humans, mice and rats. The in vivo and in vitro assays further verified that ACOX1-mediated oxidative stress contributed to mitochondrial compromise and lipid accumulation in PFOA/PFOS-exposed mouse hepatocytes, which could be mitigated by co-treatment with ACOX1 inhibitor and mitochondria ROS scavenger. Additionally, we observed that besides PFOA and PFOS, hepatic ACOX1 exhibited good-fit response to short-term exposures of long-chain (C7-C10) perfluoroalkyl carboxylic acids (PFHpA, PFNA, PFDA) and perfluoroalkyl sulfonic acids (PFHpS, PFDS) in human hepatocyte spheroids through benchmark dose (BMD) modeling.

CONCLUSION

Our study unveils a novel molecular target for PFAS-induced hepatic lipid metabolic disorders, shedding new light on prediction, assessment, and mitigation of PFAS hepatotoxicity.

A State-of-the-Science Review of Interactions of Per- and Polyfluoroalkyl Substances (PFAS) with Renal Transporters in Health and Disease: Implications for Population Variability in PFAS Toxicokinetics

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous in the environment and have been shown to cause various adverse health impacts. In animals, sex- and species-specific differences in PFAS elimination half-lives have been linked to the activity of kidney transporters. However, PFAS molecular interactions with kidney transporters are still not fully understood. Moreover, the impact of kidney disease on PFAS elimination remains unclear.

OBJECTIVES

This state-of-the-science review integrated current knowledge to assess how changes in kidney function and transporter expression from health to disease could affect PFAS toxicokinetics and identified priority research gaps that should be addressed to advance knowledge.

METHODS

We searched for studies that measured PFAS uptake by kidney transporters, quantified transporter-level changes associated with kidney disease status, and developed PFAS pharmacokinetic models. We then used two databases to identify untested kidney transporters that have the potential for PFAS transport based on their endogenous substrates. Finally, we used an existing pharmacokinetic model for perfluorooctanoic acid (PFOA) in male rats to explore the influence of transporter expression levels, glomerular filtration rate (GFR), and serum albumin on serum half-lives.

RESULTS

The literature search identified nine human and eight rat kidney transporters that were previously investigated for their ability to transport PFAS, as well as seven human and three rat transporters that were confirmed to transport specific PFAS. We proposed a candidate list of seven untested kidney transporters with the potential for PFAS transport. Model results indicated PFOA toxicokinetics were more influenced by changes in GFR than in transporter expression.

DISCUSSION

Studies on additional transporters, particularly efflux transporters, and on more PFAS, especially current-use PFAS, are needed to better cover the role of transporters across the PFAS class. Remaining research gaps in transporter expression changes in specific kidney disease states could limit the effectiveness of risk assessment and prevent identification of vulnerable populations. https://doi.org/10.1289/EHP11885.

Effects of perfluorododecanoic acid on testicular function in mice

Perfluorodecanoic acid (PFDoA) is a widely distributed environmental pollutant that can affect the functions of many organs. However, systematic evaluations of the effects of PFDoA on testicular functions are lacking. The aim of this study was to investigate the effects of PFDoA on mouse testicular functions, including spermatogenesis, testosterone synthesis, and stem Leydig cells (SLCs) in the interstitial tissue of the testis. PFDoA (0, 2, 5, 10 mg/kg/d) was administered via gavage to 2-month-old mice for 4 weeks. Serum hormone levels and sperm quality were assayed. Furthermore, to investigate the mechanisms by which PFDoA affects testosterone synthesis and spermatogenesis in vivo, the expression of StAR and P450scc in testicular tissue was measured by immunofluorescence staining and quantitative real-time PCR. In addition, the levels of SLC markers, including nestin and CD51, were studied. PFDoA decreased the luteinizing hormone concentration and sperm quality. Although the difference was not statistically significant, mean testosterone levels showed a downward trend. The expression of StAR, P450scc, CD51, and nestin was also suppressed in the PFDoA-treated groups compared with the control group. Our study suggested that PFDoA exposure can decrease testosterone biosynthesis, and even reduce the number of SLCs. These results indicated that PFDoA suppressed the main functions of testis, and further researches are required to identify strategies for preventing or reducing the effect of PFDoA on testicular function.

Human health risk assessment of 6:2 Cl-PFESA through quantitative in vitro to in vivo extrapolation by integrating cell-based assays, an epigenetic key event, and physiologically based pharmacokinetic modeling

Human health risk assessment of chemicals is essential but often relies on time-consuming and animal and labor-extensive procedures. Here, we develop a population-based, quantitative in vitro to in vivo extrapolation (QIVIVE) approach which depended on cellular effects monitored by in vitro assays, considered chemical internal concentration determined by LC-MS/MS, extrapolated into in vivo target tissue concentration through physiologically based pharmacokinetic (PBPK) modelling, and assessed populational health risk using in silico modelling. By applying this QIVIVE approach to 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA), as a representative of the emerging pollutants, we find that 6:2 Cl-PFESA disturbed lipid homeostasis in HepG2 cells through enhancement of lipid accumulation and fatty acid β-oxidation, during which miR-93-5p served as a key event towards toxicity and thus, could serve as an efficient toxicity marker for risk assessment; further, the disruption potency of lipid homeostasis of 6:2 Cl-PFESA for the most of studied populations in China might be of moderate concern. Together, our approach improved the reliability of QIVIVE during human health risk assessment, which can readily be used for other chemicals.

Facile synthesis of Cu2+-immobilized magnetic covalent organic frameworks for highly efficient enrichment and sensitive determination of five phthalate monoesters from mouse plasma with HPLC-MS/MS

Development of a simple, highly selective, and sensitive analytical method for phthalate monoesters (mPAEs) remains a challenge due to the complexity of biological samples. To address this issue, Cu²⁺ immobilized magnetic covalent organic frameworks (Fe₃O₄@TtDt@Cu²⁺ composites) with core-shell structures were prepared to enhance the enrichment efficiency of mPAEs by a facile approach synthesis of COFs shells with inherent bifunctional groups on Fe₃O₄ NPs and further Cu²⁺ immobilization. The composites exhibit high specific surface area (348.1 m² g^-1), outstanding saturation magnetization (34.94 emu g^-1), ordered mesoporous structure, Cu²⁺ immobilization, and excellent thermal stability. Accordingly, a magnetic solid-phase extraction (MSPE) pretreatment technique based on Cu²⁺ immobilized COF composites combined with high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) was established, and key parameters including the adsorbent amount, adsorption time, elution solvent, etc. were examined in detail. The developed analytical method showed wide linear ranges (10-8000 ng L^-1), low limit of detections (LODs, 2-10 ng L^-1), and good correlation coefficients (R² ≥ 0.9904) for the five mPAEs. Furthermore, the analytical method was also successfully applied to the highly sensitive detection of metabolite mPAEs in mouse plasma samples, indicating the promising application of the Fe₃O₄@TtDt@Cu²⁺ composites as a quick and efficient adsorbent in the sample pretreatment.

Machine learning and artificial intelligence in physiologically based pharmacokinetic modeling

Physiologically based pharmacokinetic (PBPK) models are useful tools in drug development and risk assessment of environmental chemicals. PBPK model development requires the collection of species-specific physiological, and chemical-specific absorption, distribution, metabolism, and excretion (ADME) parameters, which can be a time-consuming and expensive process. This raises a need to create computational models capable of predicting input parameter values for PBPK models, especially for new compounds. In this review, we summarize an emerging paradigm for integrating PBPK modeling with machine learning (ML) or artificial intelligence (AI)-based computational methods. This paradigm includes 3 steps (1) obtain time-concentration PK data and/or ADME parameters from publicly available databases, (2) develop ML/AI-based approaches to predict ADME parameters, and (3) incorporate the ML/AI models into PBPK models to predict PK summary statistics (eg, area under the curve and maximum plasma concentration). We also discuss a neural network architecture "neural ordinary differential equation (Neural-ODE)" that is capable of providing better predictive capabilities than other ML methods when used to directly predict time-series PK profiles. In order to support applications of ML/AI methods for PBPK model development, several challenges should be addressed (1) as more data become available, it is important to expand the training set by including the structural diversity of compounds to improve the prediction accuracy of ML/AI models; (2) due to the black box nature of many ML models, lack of sufficient interpretability is a limitation; (3) Neural-ODE has great potential to be used to generate time-series PK profiles for new compounds with limited ADME information, but its application remains to be explored. Despite existing challenges, ML/AI approaches will continue to facilitate the efficient development of robust PBPK models for a large number of chemicals.

An Overview of Physiologically-Based Pharmacokinetic Models for Forensic Science

A physiologically-based pharmacokinetic (PBPK) model represents the structural components of the body with physiologically relevant compartments connected via blood flow rates described by mathematical equations to determine drug disposition. PBPK models are used in the pharmaceutical sector for drug development, precision medicine, and the chemical industry to predict safe levels of exposure during the registration of chemical substances. However, one area of application where PBPK models have been scarcely used is forensic science. In this review, we give an overview of PBPK models successfully developed for several illicit drugs and environmental chemicals that could be applied for forensic interpretation, highlighting the gaps, uncertainties, and limitations.

A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species

Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t½) have been observed in some cases. Knowledge of chemical-specific t½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.

In vivo relative bioavailability of perfluorooctanoic acid (PFOA) and its alternative hexafluoropropylene oxide trimer acid (HFPO-TA): Influence of food and mechanisms exploration

The extensive use of perfluorooctanoic acid (PFOA), and its substitute hexafluoropropylene oxide trimer acid (HFPO-TA) has resulted in their frequent detection in environmental samples. However, little is known of their bioavailability via oral ingestion and the influence of food co-ingestion on absorption. Here, the relative bioavailability (RBA) of PFOA and HFPO-TA in soil was measured using an in vivo mouse model in the presence of food with different nutritional statuses (n = 11). PFOA and HFPO-TA RBA in soil was variable depending on nutrient co-administration, ranging from 29.8-95.5 % and 43.9-68.0 %, respectively. For both PFOA and HFPO-TA, a significantly negative correlation was observed between RBA and protein content in food (r = 0.57-0.72), while a positive correlation was observed with carbohydrate content (r = 0.51-0.57). Mechanistic studies showed that protein in food decreased PFOA and HFPO-TA RBA by down-regulating the expression of fatty acid binding protein 1 (FABP1) and up-regulating the expression of multidrug resistance associated protein 4 (Mrp4) in the liver, which are responsible for the absorption and efflux of PFOA and HFPO-TA. Dietary carbohydrates promoted albumin synthesis and up-regulated FABP1 expression thereby enhancing absorption and increasing PFOA and HFPO-TA RBA. This study provides an insight into potential dietary strategies for reducing exposure to per- and polyfluoroalkyl substances.

Physiologically based pharmacokinetic (PBPK) modeling of BDE-209 following oral exposure in Chinese population

The wide usage of decabromodiphenyl ether (BDE-209) as additive brominated flame retardant has caused its widespread occurrence in the environment and high exposure risk in humans. Estimating its internal exposure dose and reconstruction of external exposure dose using physiologically based pharmacokinetic (PBPK) modelling approach is a key step in the risk assessment of BDE-209. However, the PBPK model for BDE-209 is currently unavailable. This study has established two oral permeability-limited PBPK models of BDE-209 without enterohepatic recirculation (EHR) (model 1) and with EHR (model 2) for Chinese population. Using the in vitro experiments, the average binding of BDE-209 to human plasma protein (99.64% ± 2.97%) was obtained. Moreover, blood sample analysis and systematic literature review were performed to obtain internal and external exposure data of BDE-209 used for model calibration and validation. The predictions of both models were within 2-fold of the observed, and a longer half-life of serum BDE-209 was observed in model 2 than model 1. Based on the models, a human biomonitoring guidance value (HBM-GV) of 93.61 μg/g lw was derived for BDE-209, and there is no health risk found for Chinese population currently. This study provides new quantitative assessment tools for health risk assessment of BDE-209.

Commentary: cumulative risk assessment of perfluoroalkyl carboxylic acids and perfluoralkyl sulfonic acids: what is the scientific support for deriving tolerable exposures by assembling 27 PFAS into 1 common assessment group?

This commentary proposes an approach to risk assessment of mixtures of per- and polyfluorinated alkyl substances (PFAS) as EFSA was tasked to derive a tolerable intake for a group of 27 PFAS. The 27 PFAS to be considered contain different functional groups and have widely variable physicochemical (PC) properties and toxicokinetics and thus should not treated as one group based on regulatory guidance for risk assessment of mixtures. The proposed approach to grouping is to split the 27 PFAS into two groups, perfluoroalkyl carboxylates and perfluoroalkyl sulfonates, and apply a relative potency factor approach (as proposed by RIVM) to obtain two separate group TDIs based on liver toxicity in rodents since liver toxicity is a sensitive response of rodents to PFAS. Short chain PFAS and other PFAS structures should not be included in the groups due to their low potency and rapid elimination. This approach is in better agreement with scientific and regulatory guidance for mixture risk assessment.

Immunomodulation and exposure to per- and polyfluoroalkyl substances: an overview of the current evidence from animal and human studies

Per- and polyfluoroalkyl substances (PFAS) have been widely used and represent a class of environmental persistent chemicals. An association of a reduction of vaccination efficacy with PFAS serum levels in humans was used by the European Food Safety Authority as a key effect for PFAS risk assessment. The data support for using this association is reviewed by a critical analysis of the respective human epidemiology and the available animal studies on the immunomodulation of PFAS. Based on an analysis of the available human epidemiology, the overall level of evidence regarding associations between PFAS serum levels and reduced antibody response remains weak. Absence of an association between an increase in clinical infections and PFAS serum levels and the limited understanding of the importance of antibody levels as an isolated data point further support this conclusion. Animal toxicity studies with PFAS focusing on immunomodulation also provide only limited support for immunomodulation as an important endpoint in PFAS toxicity. While immunomodulation is observed after PFAS administration, generally at blood concentrations several orders of magnitude above those seen in environmentally exposed humans, the relevance of these observation is hampered by the high doses required to influence immune endpoints, the limited number of endpoints assessed, and inconsistent results. The limitations of the current database on associations of human PFAS exposures outlined here indicate that more evidence is required to select immunomodulation as a critical endpoint for human PFAS risk assessment.

Development of a multi-route physiologically based pharmacokinetic (PBPK) model for nanomaterials: a comparison between a traditional versus a new route-specific approach using gold nanoparticles in rats

BACKGROUND

Physiologically based pharmacokinetic (PBPK) modeling is an important tool in predicting target organ dosimetry and risk assessment of nanoparticles (NPs). The methodology of building a multi-route PBPK model for NPs has not been established, nor systematically evaluated. In this study, we hypothesized that the traditional route-to-route extrapolation approach of PBPK modeling that is typically used for small molecules may not be appropriate for NPs. To test this hypothesis, the objective of this study was to develop a multi-route PBPK model for different sizes (1.4-200 nm) of gold nanoparticles (AuNPs) in adult rats following different routes of administration (i.e., intravenous (IV), oral gavage, intratracheal instillation, and endotracheal inhalation) using two approaches: a traditional route-to-route extrapolation approach for small molecules and a new approach that is based on route-specific data that we propose to be applied generally to NPs.

RESULTS

We found that the PBPK model using this new approach had superior performance than the traditional approach. The final PBPK model was optimized rigorously using a Bayesian hierarchical approach with Markov chain Monte Carlo simulations, and then converted to a web-based interface using R Shiny. In addition, quantitative structure-activity relationships (QSAR) based multivariate linear regressions were established to predict the route-specific key biodistribution parameters (e.g., maximum uptake rate) based on the physicochemical properties of AuNPs (e.g., size, surface area, dose, Zeta potential, and NP numbers). These results showed the size and surface area of AuNPs were the main determinants for endocytic/phagocytic uptake rates regardless of the route of administration, while Zeta potential was an important parameter for the estimation of the exocytic release rates following IV administration.

CONCLUSIONS

This study suggests that traditional route-to-route extrapolation approaches for PBPK modeling of small molecules are not applicable to NPs. Therefore, multi-route PBPK models for NPs should be developed using route-specific data. This novel PBPK-based web interface serves as a foundation for extrapolating to other NPs and to humans to facilitate biodistribution estimation, safety, and risk assessment of NPs.

An Interactive Generic Physiologically Based Pharmacokinetic (igPBPK) Modeling Platform to Predict Drug Withdrawal Intervals in Cattle and Swine: A Case Study on Flunixin, Florfenicol and Penicillin G

Violative chemical residues in edible tissues from food-producing animals are of global public health concern. Great efforts have been made to develop physiologically based pharmacokinetic (PBPK) models for estimating withdrawal intervals (WDIs) for extralabel prescribed drugs in food animals. Existing models are insufficient to address the food safety concern as these models are either limited to 1 specific drug or difficult to be used by non-modelers. This study aimed to develop a user-friendly generic PBPK platform that can predict tissue residues and estimate WDIs for multiple drugs including flunixin, florfenicol, and penicillin G in cattle and swine. Mechanism-based in silico methods were used to predict tissue/plasma partition coefficients and the models were calibrated and evaluated with pharmacokinetic data from Food Animal Residue Avoidance Databank (FARAD). Results showed that model predictions were, in general, within a 2-fold factor of experimental data for all 3 drugs in both species. Following extralabel administration and respective U.S. FDA-approved tolerances, predicted WDIs for both cattle and swine were close to or slightly longer than FDA-approved label withdrawal times (eg, predicted 8, 28, and 7 days vs labeled 4, 28, and 4 days for flunixin, florfenicol, and penicillin G in cattle, respectively). The final model was converted to a web-based interactive generic PBPK platform. This PBPK platform serves as a user-friendly quantitative tool for real-time predictions of WDIs for flunixin, florfenicol, and penicillin G following FDA-approved label or extralabel use in both cattle and swine, and provides a basis for extrapolating to other drugs and species.

Insights into the Competitive Mechanisms of Per- and Polyfluoroalkyl Substances Partition in Liver and Blood

Some per- and polyfluoroalkyl substances (PFASs) tend to be accumulated in liver and cause hepatotoxicity. However, the difficulty to directly measure liver concentrations of PFASs in humans hampers our understanding of their hepatotoxicity and mechanisms of action. We investigated the partitioning of 11 PFASs between liver and blood in male CD-1 mice. Although accumulation of the perfluoroalkanesulfonic acids (PFSAs) in mice serum was higher than their carboxylic acids (PFCAs) counterparts as expected, the liver-blood partition coefficients (R_L/S) of PFSAs were lower than the PFCAs R_L/S, implying a competition between liver and blood. The in vitro experiments further indicated that the partitioning was dominantly determined by their competitive binding between human liver fatty acid binding protein (hL-FABP) and serum albumin (HSA). The binding affinities (K_d) of PFASs to both proteins were measured. The correlations between the R_L/S and log K_d (hL-FABP)/log K_d (HSA) were stronger than those with log K_d (hL-FABP) alone, magnifying that the partitioning was dominantly controlled by competitive binding between hL-FABP and HSA. Therefore, the liver concentrations of the selected PFASs in humans could be predicted from the available serum concentrations, which is important for assessing their hepatotoxicity.

Comparative Pharmacokinetics of Sulfadiazine and Its Metabolite N4-Acetyl Sulfadiazine in Grass Carp (Ctenopharyngodon idella) at Different Temperatures after Oral Administration

In this study, the plasma pharmacokinetics and tissue disposition of sulfadiazine (SDZ) and its main metabolite, N₄-acetyl sulfadiazine (ACT-SDZ), were compared between 18 and 24 °C following a single oral administration of SDZ at 50 mg/kg in grass carp (Ctenopharyngodon idella). The plasma and tissues were sampled from 0.167 h up to 96 h and analyzed by ultra-performance liquid chromatography with an ultraviolet detector. The pharmacokinetic parameters were estimated using a one-compartmental approach. Results showed that pharmacokinetics of SDZ and ACT-SDZ in plasma and tissues were notably influenced by the increase of temperature. The increased temperature shortened the absorption half-life (K01_HL) of SDZ and ACT-SDZ in gill, kidney, and plasma, but increased in liver and muscle + skin. The elimination half-life (K10_HF) and the area under concentration-time curve (AUC_0–∞) of SDZ and ACT-SDZ all presented a declined trend. The apparent volume of distribution (V_F) of SDZ in plasma was increased from 0.93 to 1.64 L/kg, and the apparent systemic total body clearance (Cl_F) was also increased from 0.01 to 0.05 L/h/kg. Overall, the rise of temperature decreased K10_HF, AUC_0–∞ of SDZ, and ACT-SDZ in plasma and tissues, but increased V_F and Cl_F in the plasma for SDZ.

Integration of Toxicogenomics and Physiologically Based Pharmacokinetic Modeling in Human Health Risk Assessment of Perfluorooctane Sulfonate

Toxicogenomics and physiologically based pharmacokinetic (PBPK) models are useful approaches in chemical risk assessment, but the methodology to incorporate toxicogenomic data into a PBPK model to inform risk assessment remains to be developed. This study aimed to develop a probabilistic human health risk assessment approach by integrating toxicogenomic dose-response data and PBPK modeling using perfluorooctane sulfonate (PFOS) as a case study. Based on the available human in vitro and mouse in vivo toxicogenomic data, we identified the differentially expressed genes (DEGs) at each exposure paradigm/duration. Kyoto Encyclopedia of Genes and Genomes and disease ontology enrichment analyses were conducted on the DEGs to identify significantly enriched pathways and diseases. The dose-response data of DEGs were analyzed using the Bayesian benchmark dose (BMD) method. Using a previously published PBPK model, the gene BMDs were converted to human equivalent doses (HEDs), which were summarized to pathway and disease HEDs and then extrapolated to reference doses (RfDs) by considering an uncertainty factor of 30 for mouse in vivo data and 10 for human in vitro data. The results suggested that the median RfDs at different exposure paradigms were similar to the 2016 U.S. Environmental Protection Agency's recommended RfD, while the RfDs for the most sensitive pathways and diseases were closer to the recent European Food Safety Authority's guidance values. In conclusion, genomic dose-response data and PBPK modeling can be integrated to become a useful alternative approach in risk assessment of environmental chemicals. This approach considers multiple endpoints, provides toxicity mechanistic insights, and does not rely on apical toxicity endpoints.

Effect of Enterohepatic Circulation on the Accumulation of Per- and Polyfluoroalkyl Substances: Evidence from Experimental and Computational Studies

The pharmacokinetic characteristics of per- and polyfluoroalkyl substances (PFAS) affect their distribution and bioaccumulation in biological systems. The enterohepatic circulation leads to reabsorption of certain chemicals from bile back into blood and the liver and thus influences their elimination, yet its influence on PFAS bioaccumulation remains unclear. We explored the role of enterohepatic circulation in PFAS bioaccumulation by examining tissue distribution of various PFAS in wild fish and a rat model. Computational models were used to determine the reabsorbed fractions of PFAS by calculating binding affinities of PFAS for key transporter proteins of enterohepatic circulation. The results indicated that higher concentrations were observed in blood, the liver, and bile compared to other tissues for some PFAS in fish. Furthermore, exposure to a PFAS mixture on the rat model showed that the reabsorption phenomenon appeared during 8-12 h for most long-chain PFAS. Molecular docking calculations suggest that PFAS can bind to key transporter proteins via electrostatic and hydrophobic interactions. Further regression analysis adds support to the hypothesis that binding affinity of the apical sodium-dependent bile acid transporter is the most important variable to predict the human half-lives of PFAS. This study demonstrated the critical role of enterohepatic circulation in reabsorption, distribution, and accumulation of PFAS.

Physiologically based pharmacokinetic (PBPK) modeling of perfluorohexane sulfonate (PFHxS) in humans

Per- and polyfluoroalkyl substances (PFAS) are persistent, man-made compounds prevalent in the environment and consistently identified in human biomonitoring samples. In particular, perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and perfluorohexane sulfonic acid (PFHxS) have been identified at U.S. Air Force installations. The study of human toxicokinetics and physiologically based pharmacokinetic (PBPK) modeling of PFHxS has been less robust and has been limited in scope and application as compared to PFOS and PFOA. The primary goal of the current effort was to develop a PBPK model describing PFHxS disposition in humans that can be applied to retrospective, current, and future human health risk assessment of PFHxS. An existing model developed for PFOS and PFOA was modified and key parameter values for exposure and toxicokinetics were calibrated for PFHxS prediction based on human biomonitoring data, particularly general population serum levels from the U.S. Centers for Disease Prevention and Control (CDC) National Health and Nutrition Examination Survey (NHANES). Agreement between the model and the calibration and evaluation data was excellent and recapitulated observed trends across sex, age, and calendar years. Confidence in the model is greatest for application to adults in the 2000-2018 time frame and for shorter-term future projections.

Bayesian Refinement of the Permeability-Limited Physiologically Based Pharmacokinetic Model for Perfluorooctanoic Acid in Male Rats

Physiologically based pharmacokinetic (PBPK) modeling is a powerful technique to inform risk assessment of xenobiotic substances such as perfluorooctanoic acid (PFOA). In our previous study, a permeability-limited PBPK model was developed to simulate the toxicokinetics and tissue distribution of PFOA in male rats. However, due to limited information on some key model parameters (e.g., protein binding and active transport rates), the uncertainty of the permeability-limited PBPK model was quite high. To address this issue, a hierarchical Bayesian analysis with Markov chain Monte Carlo (MCMC) was applied to reduce the uncertainty of parameters and improve the performance of the PBPK model. With the optimized posterior parameters, the PBPK model was evaluated by comparing its prediction with experimental data from three different studies. The results show that the uncertainties of the posterior model parameters were reduced substantially. In addition, most of the PBPK model predictions were improved: with the posterior parameters, most of the predicted plasma toxicokinetics (e.g., half-life) and tissue distribution fell well within a factor of 2.0 of the experimental data. Finally, the Bayesian framework could provide insights into the molecular mechanisms driving PFOA toxicokinetics: PFOA-protein binding, membrane permeability, and active transport.

An integrated exposure and pharmacokinetic modeling framework for assessing population-scale risks of phthalates and their substitutes

To effectively incorporate in vitro-in silico-based methods into the regulation of consumer product safety, a quantitative connection between product phthalate concentrations and in vitro bioactivity data must be established for the general population. We developed, evaluated, and demonstrated a modeling framework that integrates exposure and pharmacokinetic models to convert product phthalate concentrations into population-scale risks for phthalates and their substitutes. A probabilistic exposure model was developed to generate the distribution of multi-route exposures based on product phthalate concentrations, chemical properties, and human activities. Pharmacokinetic models were developed to simulate population toxicokinetics using Bayesian analysis via the Markov chain Monte Carlo method. Both exposure and pharmacokinetic models demonstrated good predictive capability when compared with worldwide studies. The distributions of exposures and pharmacokinetics were integrated to predict the population distributions of internal dosimetry. The predicted distributions showed reasonable agreement with the U.S. biomonitoring surveys of urinary metabolites. The "source-to-outcome" local sensitivity analysis revealed that food contact materials had the greatest impact on body burden for di(2-ethylhexyl) adipate (DEHA), di-2-ethylhexyl phthalate (DEHP), di(isononyl) cyclohexane-1,2-dicarboxylate (DINCH), and di(2-propylheptyl) phthalate (DPHP), whereas the body burden of diethyl phthalate (DEP) was most sensitive to the concentration in personal care products. The upper bounds of predicted plasma concentrations showed no overlap with ToxCast in vitro bioactivity values. Compared with the in vitro-to-in vivo extrapolation (IVIVE) approach, the integrated modeling framework has significant advantages in mapping product phthalate concentrations to multi-route risks, and thus is of great significance for regulatory use with a relatively low input requirement. Further integration with new approach methodologies will facilitate these in vitro-in silico-based risk assessments for a broad range of products containing an equally broad range of chemicals.

Risk Assessment of Perfluorooctane Sulfonate (PFOS) using Dynamic Age Dependent Physiologically based Pharmacokinetic Model (PBPK) across Human Lifetime

The widespread use of Perfluorooctane sulfonate (PFOS) in everyday life, its long half-life, and the lipophilicity that makes it easily accumulate in the body, raises the question of its safe exposure among different population groups. There are currently enough epidemiological studies showing evidence of PFOS exposure and its associated adverse effects on humans. Moreover, it is already known that physiological changes along with age e.g. organ volume, renal blood flow, cardiac output and albumin concentrations affect chemicals body burden. Human biomonitoring cohort studies have reported PFOS concentrations in blood and autopsy tissue data with PFOS present in sensitive organs across all human lifespan. However, to interpret such biomonitoring data in the context of chemical risk assessment, it is necessary to have a mechanistic framework that explains show the physiological changes across age affects the concentration of chemical inside different tissues of the human body. PBPK model is widely and successfully used in the field of risk assessment. The objective of this manuscript is to develop a dynamic age-dependent PBPK model as an extension of the previously published adult PFOS model and utilize this model to predict and compare the PFOS tissue distribution and plasma concentration across different age groups. Different cohort study data were used for exposure dose reconstruction and evaluation of time-dependent concentration in sensitive organs. Predicted plasma concentration followed trends observed in biomonitoring data and model predictions showed the increased disposition of PFOS in the geriatric population. PFOS model is sensitive to parameters governing renal resorption and elimination across all ages, which is related to PFOS half-life in humans. This model provides an effective framework for improving the quantitative risk assessment of PFOS throughout the human lifetime, particularly in susceptible age groups. The dynamic age-dependent PBPK model provides a step forward for developing such kind of dynamic model for other perfluoroalkyl substances.

Why is elevation of serum cholesterol associated with exposure to perfluoroalkyl substances (PFAS) in humans? A workshop report on potential mechanisms

Serum concentrations of cholesterol are positively correlated with exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) in humans. The associated change in cholesterol is small across a broad range of exposure to PFOA and PFOS. Animal studies generally have not indicated a mechanism that would account for the association in humans. The extent to which the relationship is causal is an open question. Nonetheless, the association is of particular importance because increased serum cholesterol has been considered as an endpoint to derive a point of departure in at least one recent risk assessment. To gain insight into potential mechanisms for the association, both causal and non-causal, an expert workshop was held Oct 31 and Nov 1, 2019 to discuss relevant data and propose new studies. In this report, we summarize the relevant background data, the discussion among the attendees, and their recommendations for further research.

Application of physiologically based pharmacokinetic models to promote the development of veterinary drugs with high efficacy and safety

Physiologically based pharmacokinetic (PBPK) models have become important tools for the development of novel human drugs. Food-producing animals and pets comprise an important part of human life, and the development of veterinary drugs (VDs) has greatly impacted human health. Owing to increased affordability of and demand for drug development, VD manufacturing companies should have more PBPK models required to reduce drug production costs. So far, little attention has been paid on applying PBPK models for the development of VDs. This review begins with the development processes of VDs; then summarizes case studies of PBPK models in human or VD development; and analyzes the application, potential, and advantages of PBPK in VD development, including candidate screening, formulation optimization, food effects, target-species safety, and dosing optimization. Then, the challenges of applying the PBPK model to VD development are discussed. Finally, future opportunities of PBPK models in designing dosing regimens for intracellular pathogenic infections and for efficient oral absorption of VDs are further forecasted. This review will be relevant to readers who are interested in using a PBPK model to develop new VDs.

Bayesian population physiologically-based pharmacokinetic model for robustness evaluation of withdrawal time in tilapia aquaculture administrated to florfenicol

The antimicrobial residues of aquacultural production is a growing public concern, leading to reexamine the method for establishing robust withdrawal time and ensuring food safety. Our study aims to develop the optimizing population physiologically-based pharmacokinetic (PBPK) model for assessing florfenicol residues in the tilapia tissues, and for evaluating the robustness of the withdrawal time (WT). Fitting with published pharmacokinetic profiles that experimented under temperatures of 22 and 28 °C, a PBPK model was constructed by applying with the Bayesian Markov chain Monte Carol (MCMC) algorithm to estimate WTs under different physiological, environmental and dosing scenarios. Results show that the MCMC algorithm improves the estimates of uncertainty and variability of PBPK-related parameters, and optimizes the simulation of the PBPK model. It is noteworthy that posterior sets generated from temperature-associated datasets to be respectively used for simulating residues under corresponding temperature conditions. Simulating the residues under regulated regimen and overdosing scenarios for Taiwan, the estimated WTs were 12-16 days at 22 °C and 9-12 days at 28 °C, while for the USA, the estimated WTs were 14-18 and 11-14 days, respectively. Comparison with the regulated WT of 15 days, results indicate that the current WT has well robustness and resilience in the environment of higher temperatures. The optimal Bayesian population PBPK model provides effective analysis for determining WTs under scenario-specific conditions. It is a new insight into the increasing body of literature on developing the Bayesian-PBPK model and has practical implications for improving the regulation of food safety.

Suggestions for Improving the Characterization of Risk from Exposures to Per and Polyfluorinated Alkyl Substances

Many state and federal environmental and health agencies have developed risk-based criteria for assessing the risk of adverse health effects of per- and polyfluorinated alkyl substances (PFAS) exposure to humans and the environment. However, the criteria that have been developed vary; drinking water criteria developed for perfluorooctanoic acid, for example, can vary by up to 750 fold. This is due to differences and variability in the data and information used, study/endpoint selection, assumptions and magnitude of uncertainty factors used in the absence and extrapolation of critical effect data, differences in underlying approaches to addressing exposure within criteria development, and/or policy decisions on levels of acceptable risk. We have critically evaluated the methods used to develop these criteria while focusing on derivation and application of drinking water criteria and discuss a range of improvements to risk-characterization practice recently presented at a Focused Topic Meeting on PFAS conducted by the Society of Environmental Toxicology and Chemistry in Durham, North Carolina, USA, 12 to 15 August 2019. We propose methods that consider maximizing the use of disparate data streams, seeking patterns, and proposing biologically based approaches to evidence integration toward informed criteria development. Environ Toxicol Chem 2021;40:883-898. © 2020 SETAC.