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

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.

Threats of per- and poly-fluoroalkyl pollutants to susceptible populations

Environmental exposure to per- and poly-fluoroalkyl substances (PFAS) has raised significant global health concerns due to potential hazards in healthy adults. However, the impact of PFAS on susceptible populations, including pregnant individuals, newborns, the older people, and those with underlying health conditions, has been overlooked. These susceptible groups often have physiological changes that make them less resilient to the same exposures. Consequently, there is an urgent need for a comprehensive understanding of the health risks posed by PFAS exposure to these populations. In this review, we delve into the potential health risks of PFAS exposure in these susceptible populations. Equally important, we also examine and discuss the molecular mechanisms that underlie this susceptibility. These mechanisms include the induction of oxidative stress, disruption of the immune system, impairment of cellular metabolism, and alterations in gut microbiota, all of which contribute to the enhanced toxicity of PFAS in susceptible populations. Finally, we address the primary research challenges and unresolved issues that require further investigation. This discussion aims to foster research for a better understanding of how PFAS affect susceptible populations and to pave the way for strategies to minimize their adverse effects.

Estimating the dynamic early life exposure to PFOA and PFOS of the HELIX children: Emerging profiles via prenatal exposure, breastfeeding, and diet

In utero and children's exposure to per- and polyfluoroalkyl substances (PFAS) is a major concern in health risk assessment as early life exposures are suspected to induce adverse health effects. Our work aims to estimate children's exposure (from birth to 12 years old) to PFOA and PFOS, using a Physiologically-Based Pharmacokinetic (PBPK) modelling approach. A model for PFAS was updated to simulate the internal PFAS exposures during the in utero life and childhood, and including individual characteristics and exposure scenarios (e.g., duration of breastfeeding, weight at birth, etc.). Our approach was applied to the HELIX cohort, involving 1,239 mother-child pairs with measured PFOA and PFOS plasma concentrations at two sampling times: maternal and child plasma concentrations (6 to 12 y.o). Our model predicted an increase in plasma concentrations during fetal development and childhood until 2 y.o when the maximum concentrations were reached. Higher plasma concentrations of PFOA than PFOS were predicted until 2 y.o, and then PFOS concentrations gradually became higher than PFOA concentrations. From 2 to 8 y.o, mean concentrations decreased from 3.1 to 1.88 µg/L or ng/mL (PFOA) and from 4.77 to 3.56 µg/L (PFOS). The concentration-time profiles vary with the age and were mostly influenced by in utero exposure (on the first 4 months after birth), breastfeeding (from 5 months to 2 (PFOA) or 5 (PFOS) y.o of the children), and food intake (after 3 (PFOA) or 6 (PFOS) y.o of the children). Similar measured biomarker levels can correspond to large differences in the simulated internal exposures, highlighting the importance to investigate the children's exposure over the early life to improve exposure classification. Our approach demonstrates the possibility to simulate individual internal exposures using PBPK models when measured biomarkers are scarce, helping risk assessors in gaining insight into internal exposure during critical windows, such as early life.

Studying mixture effects on uptake and tissue distribution of PFAS in zebrafish (Danio rerio) using physiologically based kinetic (PBK) modelling

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously distributed in the aquatic environment. They include persistent, mobile, bioaccumulative, and toxic chemicals and it is therefore critical to increase our understanding on their adsorption, distribution, metabolism, excretion (ADME). The current study focused on uptake of seven emerging PFAS in zebrafish (Danio rerio) and their potential maternal transfer. In addition, we aimed at increasing our understanding on mixture effects on ADME by developing a physiologically based kinetic (PBK) model capable of handling co-exposure scenarios of any number of chemicals. All studied chemicals were taken up in the fish to varying degrees, whereas only perfluorononanoate (PFNA) and perfluorooctanoate (PFOA) were quantified in all analysed tissues. Perfluorooctane sulfonamide (FOSA) was measured at concerningly high concentrations in the brain (Cmax over 15 μg/g) but also in the liver and ovaries. All studied PFAS were maternally transferred to the eggs, with FOSA and 6:2 perfluorooctane sulfonate (6,2 FTSA) showing significant (p < 0.02) signs of elimination from the embryos during the first 6 days of development, while perfluorobutane sulfonate (PFBS), PFNA, and perfluorohexane sulfonate (PFHxS) were not eliminated in embryos during this time-frame. The mixture PBK model resulted in >85 % of predictions within a 10-fold error and 60 % of predictions within a 3-fold error. At studied levels of PFAS exposure, competitive binding was not a critical factor for PFAS kinetics. Gill surface pH influenced uptake for some carboxylates but not the sulfonates. The developed PBK model provides an important tool in understanding kinetics under complex mixture scenarios and this use of New Approach Methodologies (NAMs) is critical in future risk assessment of chemicals and early warning systems.

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.

Development and Evaluation of Amorphous Solid Dispersion of Riluzole with PBPK Model to Simulate the Pharmacokinetic Profile

In the current work, screening of polymers viz. polyacrylic acid (PAA), polyvinyl pyrrolidone vinyl acetate (PVP VA), and hydroxypropyl methyl cellulose acetate succinate (HPMC AS) based on drug-polymer interaction and wetting property was done for the production of a stable amorphous solid dispersion (ASD) of a poorly water-soluble drug Riluzole (RLZ). PAA showed maximum interaction and wetting property hence, was selected for further studies. Solid state characterization studies confirmed the formation of ASD with PAA. Saturation solubility, dissolution profile, and in vivo pharmacokinetic data of the ASD formulation were generated in rats against its marketed tablet Rilutor. The RLZ:PAA ASD showed exponential enhancement in the dissolution of RLZ. Predicted and observed pharmacokinetic data in rats showed enhanced area under curve (AUC) and Cmax in plasma and brain with respect to Rilutor. Furthermore, a physiologically based pharmacokinetic (PBPK) model of rats for Rilutor and RLZ ASD was developed and then extrapolated to humans where physiological parameters were changed along with a biochemical parameter. The partition coefficient was kept similar in both species. The model was used to predict different exposure scenarios, and the simulated data was compared with observed data points. The PBPK model simulated Cmax and AUC was within two times the experimental data for plasma and brain. The Cmax and AUC in the brain increased with ASD compared to Rilutor for humans showing its potential in improving its biopharmaceutical performance and hence enhanced therapeutic efficacy. The model can predict the RLZ concentration in multiple compartments including plasma and liver.

Development of a Physiologically Based Pharmacokinetic Model for Nitrofurantoin in Rabbits, Rats, and Humans

Nitrofurantoin (NFT) is a commonly used antibiotic for the treatment of urinary tract infections that can cause liver toxicity. Despite reports of hepatic adverse events associated with NFT exposure, there is still limited understanding of the interplay between NFT exposure, its disposition, and the risk of developing liver toxicity. In this study, we aim to develop a physiologically based pharmacokinetic (PBPK) model for NFT in three different species (rabbits, rats, and humans) that can be used as a standard tool for predicting drug-induced liver injury (DILI). We created several versions of the PBPK model using previously published kinetics data from rabbits, and integrated enterohepatic recirculation (EHR) using rat data. Our model showed that active tubular secretion and reabsorption in the kidney are critical in explaining the non-linear renal clearance and urine kinetics of NFT. We subsequently extrapolated the PBPK model to humans. Adapting the physiology to humans led to predictions consistent with human kinetics data, considering a low amount of NFT to be excreted into bile. Model simulations predicted that the liver of individuals with a moderate-to-severe glomerular filtration rate (GFR) is exposed to two-to-three-fold higher concentrations of NFT than individuals with a normal GFR, which coincided with a substantial reduction in the NFT urinary concentration. In conclusion, people with renal insufficiency may be at a higher risk of developing DILI due to NFT exposure, while at the same time having a suboptimal therapeutic effect with a high risk of drug resistance. Our PBPK model can in the future be used to predict NFT kinetics in individual patients on the basis of characteristics like age and GFR.

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.

Development of a Rat Physiologically Based Kinetic Model (PBK) for three Organophosphate Flame Retardants (TDCIPP, TCIPP, TCEP)

Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), Tris (1-chloro-2-propyl) phosphate (TCIPP) and tris (2-chloroethyl) phosphate (TCEP) are three widely used organophosphate flame retardants (OPFRs) being frequently detected in human body fluids. Although OPFRs are being detected in human beings, the toxicological effects of their exposure are not clearly understood due to limited data. For this, a physiologically based kinetic model (PBK) was developed in MCSIM integrated with R studio and validated in rats to understand the toxicokinetics of OPFRs for the first time. The model required the enterohepatic recirculation (EHR) mechanism which was included to explain the non-linear data. Model parameters were optimized using the Bayesian framework (Markov Chain Monte Carlo) along with a visual fitting to explain toxicokinetic data. Goodness-of-fit was calculated to evaluate model predictability power in Rstudio. The model can appropriately predict the concentration of OPFRs in several organs like plasma, urine, kidney, etc. within 1-2-fold of experimental data. Slow elimination of OPFRs was observed from adipose tissue and brain at late time points, showing their potential to accumulate upon daily exposure. The use of PBK was demonstrated by reconstructing the oral exposure equivalent to the in-vitro toxic dose to support neurotoxic risk assessment. This version of PBK can be extrapolated to human for toxicological risk assessment. Nonetheless, further investigation is required to understand whether these chemicals follow similar kinetics in humans, which could lead to a greater risk to human health. CODE AVAILABILITY: The model will be available to access through Rshiny using GIThub soon, InSilicoVida/Flame-Retardant-PBPK-Model: It contains organophosphate flame retardant (OPFRs) PBK for TDCIPP, TCIPP and TCEP (github.com).

The Role of "Physiologically Based Pharmacokinetic Model (PBPK)" New Approach Methodology (NAM) in Pharmaceuticals and Environmental Chemical Risk Assessment

Physiologically Based Pharmacokinetic (PBPK) models are mechanistic tools generally employed in the pharmaceutical industry and environmental health risk assessment. These models are recognized by regulatory authorities for predicting organ concentration-time profiles, pharmacokinetics and daily intake dose of xenobiotics. The extension of PBPK models to capture sensitive populations such as pediatric, geriatric, pregnant females, fetus, etc., and diseased populations such as those with renal impairment, liver cirrhosis, etc., is a must. However, the current modelling practices and existing models are not mature enough to confidently predict the risk in these populations. A multidisciplinary collaboration between clinicians, experimental and modeler scientist is vital to improve the physiology and calculation of biochemical parameters for integrating knowledge and refining existing PBPK models. Specific PBPK covering compartments such as cerebrospinal fluid and the hippocampus are required to gain mechanistic understanding about xenobiotic disposition in these sub-parts. The PBPK model assists in building quantitative adverse outcome pathways (qAOPs) for several endpoints such as developmental neurotoxicity (DNT), hepatotoxicity and cardiotoxicity. Machine learning algorithms can predict physicochemical parameters required to develop in silico models where experimental data are unavailable. Integrating machine learning with PBPK carries the potential to revolutionize the field of drug discovery and development and environmental risk. Overall, this review tried to summarize the recent developments in the in-silico models, building of qAOPs and use of machine learning for improving existing models, along with a regulatory perspective. This review can act as a guide for toxicologists who wish to build their careers in kinetic modeling.

Unravelling sex-specific BPA toxicokinetics in children using a pediatric PBPK model

Bisphenol A (BPA) is a widely known endocrine disruptor (ED) found in many children's products such as toys, feeding utensils, and teething rings. Recent epidemiology association studies have shown postnatal BPA exposure resulted in developing various diseases such as diabetes, obesity, and neurodegeneration, etc., later in their lives. However, little is known about its sex-specific metabolism and consequently internal exposure. The aim of this study was to develop a sex-specific pediatric physiologically based pharmacokinetic model (PBPK) for BPA to compare their toxicokinetic differences. First, the published adult PBPK model was re-validated, and then this model was extended by interpolation to incorporate pediatric sex specific physiological and biochemical parameters. We used both the classical body weight and ontogeny-based scaling approach to interpolate the metabolic process. Then, the pharmacokinetic attributes of the models using the two-scaling approach mentioned above were compared with adult model. Further, a sex-specific PBPK model with an ontogeny scaling approach was preferred to evaluate the pharmacokinetic differences. Moreover, this model was used to reconstruct the BPA exposure from two cohorts (Helix and PBAT Cohort) from 7 EU countries. The half-life of BPA was found to be almost the same in boys and girls at the same exposure levels. Our model estimated BPA children's exposure to be about 1500 times higher than the tolerable daily intake (TDI) recently set by European Food Safety Authority (EFSA) i.e., 0.04 ng/kg BW/day. The model demonstrated feasibility of extending the adult PBPK to sex-specific pediatric, thus investigate a gender-specific health risk assessment.

Framework for risk assessment of PFAS utilizing experimental studies and in-silico models

Perfluoroalkyl substances (PFAS), especially PFOS and PFOA, are two widely used synthetic chemicals that can impact human health based on evidence from animal and epidemiologic studies. In this paper, we have reviewed and summarized the influence of PFAS exposure on health, pointing the quality of evidence, and applied translational techniques to integrate evidence for PFAS policy making. This is the first review where highly referred articles on PFAS used for policymaking by several regulatory agencies were collected and evaluated based on the review guidelines developed by the US National Toxicology Program's Office of Health Assessment and Translation (OHAT) review guidelines. Several limitations were observed, including co-exposure to multiple chemicals and limited measurement of primary and secondary outcomes related to specific toxicity. However, data from all the studies provided a moderate to strong level of confidence for link between PFAS exposure and different adverse outcomes. Secondly, for translating the risk to humans, an in-silico model and scaling approach was utilized. Physiologically based pharmacokinetic model (PBPK) was used to calculate the human equivalent dose (HED) from two widely accepted studies and compared with tolerable daily intakes (TDIs) established by various regulatory agencies. Inter-species dose extrapolation was done to compare with human the relevance of dosing scenarios used in animals. Overall, a framework for translation of risk was proposed based on the conclusions of this review with the goal of improving policymaking. The current paper can improve the methodological protocols for PFAS experimental studies and encourage the utilization of in-silico models for translating risk.

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.