Exploring sex differences in human health risk assessment for PFNA and PFDA using a PBPK model

Estimation of species- and sex-specific PFAS pharmacokinetics in mice, rats, and non-human primates using a Bayesian hierarchical methodology

The carbon chain length, degree of fluorination, and functional group of per- and polyfluoroalkyl substances (PFAS) influences the bioaccumulation and half-lives of these substances in humans and laboratory animals. Pharmacokinetic (PK) studies using laboratory animals characterize the absorption, distribution, metabolism, and excretion (ADME) of a PFAS and can provide the underlying data for inter-species extrapolation to inform human pharmacokinetics. However, variations in ADME arise due to differences in protein binding and renal and hepatobiliary clearance mechanisms. In particular, sex- and species-specific differences in active transporter abundance and PFAS binding affinity challenge body weight-based extrapolation assumptions from animal models to human PK parameters. Because these protein-dependent changes in ADME do not always scale with species body weight, classic allometric scaling assumptions can fail to account for species-specific transporter-mediated clearance. In addition, study-dependent differences in pharmacokinetic modeling approaches and parameterization techniques can result in large differences among the PK parameters reported in the literature. To better quantify PFAS pharmacokinetics and characterize the underlying uncertainty, we implemented a Bayesian inference hierarchical model to estimate PFAS PK parameters for multiple species (mice, rats, and non-human primates) using numerous single-dose animal studies. Through an alternative parameterization of the one- and two-compartment models, this method improved parameter identifiability and allowed for the use of prior information on PFAS absorption rate, clearance, and volume of distribution. Using reported time-course concentration data, we estimated sex-specific clearance, volume of distribution, and half-life across mice, rats, and non-human primates using a consistent modeling methodology for eight PFAS: PFHxA, PFHxS, PFNA, PFDA, PFBS, PFBA, PFOA, and PFOS. The resulting comparison to available human data demonstrated that standard volume of distribution body-mass scaling (BW1) for PFAS generally agrees with reported human values while standard assumptions for allometric scaling of clearance (BW3/4) are not appropriate for most of the PFAS investigated in this study. In addition, we demonstrated that there may be considerable differences in clearance for PFAS in some species when comparing across different sexes and routes of exposure.

Cross-cutting studies of per- and polyfluorinated alkyl substances (PFAS) in Arctic wildlife and humans

This cross-cutting review focuses on the presence and impacts of per- and polyfluoroalkyl substances (PFAS) in the Arctic. Several PFAS undergo long-range transport via atmospheric (volatile polyfluorinated compounds) and oceanic pathways (perfluorinated alkyl acids, PFAAs), causing widespread contamination of the Arctic. Beyond targeting a few well-known PFAS, applying sum parameters, suspect and non-targeted screening are promising approaches to elucidate predominant sources, transport, and pathways of PFAS in the Arctic environment, wildlife, and humans, and establish their time-trends. Across wildlife species, concentrations were dominated by perfluorooctane sulfonic acid (PFOS), followed by perfluorononanoic acid (PFNA); highest concentrations were present in mammalian livers and bird eggs. Time trends were similar for East Greenland ringed seals (Pusa hispida) and polar bears (Ursus maritimus). In polar bears, PFOS concentrations increased from the 1980s to 2006, with a secondary peak in 2014-2021, while PFNA increased regularly in the Canadian and Greenlandic ringed seals and polar bear livers. Human time trends vary regionally (though lacking for the Russian Arctic), and to the extent local Arctic human populations rely on traditional wildlife diets, such as marine mammals. Arctic human cohort studies implied that several PFAAs are immunotoxic, carcinogenic or contribute to carcinogenicity, and affect the reproductive, endocrine and cardiometabolic systems. Physiological, endocrine, and reproductive effects linked to PFAS exposure were largely similar among humans, polar bears, and Arctic seabirds. For most polar bear subpopulations across the Arctic, modeled serum concentrations exceeded PFOS levels in human populations, several of which already exceeded the established immunotoxic thresholds for the most severe risk category. Data is typically limited to the western Arctic region and populations. Monitoring of legacy and novel PFAS across the entire Arctic region, combined with proactive community engagement and international restrictions on PFAS production remain critical to mitigate PFAS exposure and its health impacts in the Arctic.

Exposure to perfluorodecanoic acid impairs follicular development via inducing granulosa cell necroptosis

Per- and polyfluoroalkyl substances (PFAS) have attracted significant attention due to their environmental toxicity. However, the detrimental impact of PFAS on the development of the female reproductive system remains controversial. In this study, we investigated the effects of three specific PFAS compounds perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), and perfluorodecanoic acid (PFDA) on ovarian development. Among these compounds, PFDA demonstrated the most pronounced cytotoxic effect on ovarian granulosa cells. The results showed that a 200 μM concentration of PFDA induced cell apoptosis via the intrinsic pathway by elevating reactive oxygen species (ROS) levels and activating Caspase-9 and Caspase-3. Furthermore, 200 μM PFDA triggered necroptosis, a form of regulated cell death (RCD), through the receptor-interacting serine/threonine kinase 1 (RIPK1), receptor interacting protein kinase 3 (RIPK3), and mixed-lineage kinase domain-like protein (MLKL) axis, mediated by inhibition of the canonical apoptosis proteolytic enzyme Caspase-8. In vivo experiments confirmed that mice exposed to PFDA displayed a significantly reduced ovarian index compared to the control group, accompanied by evident follicular atresia. Ovarian tissues from the PFDA-exposed group showed upregulated necroptosis markers, which were effectively mitigated by inhibiting the phosphorylation of RIPK1 at Ser166. Importantly, this study provides the first evidence that PFDA disrupts ovarian development through a novel mechanism involving the RIPK1-mediated necroptosis pathway, alongside the detection of the intrinsic apoptosis pathway. This greatly expands our insight into the effects of PFDA on cell death. This finding highlights the potential public health hazards associated with PFDA exposure and emphasizes the need for further research to fully understand its broader implications.

Hunting Metabolic Biomarkers for Exposure to Per- and Polyfluoroalkyl Substances: A Review

Per- and polyfluoroalkyl substances (PFAS) represent a class of persistent synthetic chemicals extensively utilized across industrial and consumer sectors, raising substantial environmental and human health concerns. Epidemiological investigations have robustly linked PFAS exposure to a spectrum of adverse health outcomes. Altered metabolites stand as promising biomarkers, offering insights into the identification of specific environmental pollutants and their deleterious impacts on human health. However, elucidating metabolic alterations attributable to PFAS exposure and their ensuing health effects has remained challenging. In light of this, this review aims to elucidate potential biomarkers of PFAS exposure by presenting a comprehensive overview of recent metabolomics-based studies exploring PFAS toxicity. Details of PFAS types, sources, and human exposure patterns are provided. Furthermore, insights into PFAS-induced liver toxicity, reproductive and developmental toxicity, cardiovascular toxicity, glucose homeostasis disruption, kidney toxicity, and carcinogenesis are synthesized. Additionally, a thorough examination of studies utilizing metabolomics to delineate PFAS exposure and toxicity biomarkers across blood, liver, and urine specimens is presented. This review endeavors to advance our understanding of PFAS biomarkers regarding exposure and associated toxicological effects.

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

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.

Dietary exposure to per- and polyfluoroalkyl substances: Potential health impacts on human liver

Per- and polyfluoroalkyl substances (PFAS), dubbed "forever chemicals", are widely present in the environment. Environmental contamination and food contact substances are the main sources of PFAS in food, increasing the risk of human dietary exposure. Numerous epidemiological studies have established the link between dietary exposure to PFAS and liver disease. Correspondingly, PFAS induced-hepatotoxicity (e.g., hepatomegaly, cell viability, inflammation, oxidative stress, bile acid metabolism dysregulation and glycolipid metabolism disorder) observed from in vitro models and in vivo rodent studies have been extensively reported. In this review, the pertinent literature of the last 5 years from the Web of Science database was researched. This study summarized the source and fate of PFAS, and reviewed the occurrence of PFAS in food system (natural and processed food). Subsequently, the characteristics of human dietary exposure PFAS (population characteristics, distribution trend, absorption and distribution) were mentioned. Additionally, epidemiologic evidence linking PFAS exposure and liver disease was alluded, and the PFAS-induced hepatotoxicity observed from in vitro models and in vivo rodent studies was comprehensively reviewed. Lastly, we highlighted several critical knowledge gaps and proposed future research directions. This review aims to raise public awareness about food PFAS contamination and its potential risks to human liver health.

Low dose PFDA induces DNA damage and DNA repair inhibition by promoting nuclear cGAS accumulation in ovarian epithelial cells

Per- and polyfluoroalkyl substances (PFAS), the versatile anthropogenic chemicals, are popular with the markets and manufactured in large quantities yearly. Accumulation of PFAS has various adverse health effects on human. Albeit certain members of PFAS were found to have genotoxicity in previous studies, the mechanisms underlying their effects on DNA damage repair remain unclear. Here, we investigated the effects of Perfluorodecanoic acid (PFDA) on DNA damage and DNA damage repair in ovarian epithelial cells through a series of in vivo and in vitro experiments. At environmentally relevant concentration, we firstly found that PFDA can cause DNA damage in primary mouse ovarian epithelial cells and IOSE-80 cells. Moreover, nuclear cGAS increased in PFDA-treated cells, which leaded to the efficiency of DNA homologous recombination (HR) decreased and DNA double-strand breaks perpetuated. In vivo experiments also verified that PFDA can induce more DNA double-strand breaks lesions and nuclear cGAS in ovarian tissue. Taken together, our results unveiled that low dose PFDA can cause deleterious effects on DNA and DNA damage repair (DDR) in ovarian epithelial cells and induce genomic instability.

Perfluorodecanoic acid promotes high-fat diet-triggered adiposity and hepatic lipid accumulation by modulating the NLRP3/caspase-1 pathway in male C57BL/6J mice

Perfluorodecanoic acid (PFDA), a chemical contaminant, may casue became obesity, which makes it a public health concern. In this study, we investigated the effects of PFDA on adiposity development and hepatic lipid accumulation in mice fed with a high-fat diet (HFD). Animals were assigned to two diet treatments (low-fat and high-fat); and PFDA was administered through drinking water for 12 weeks. The contaminant promoted body weight gain and adiposity in HFD-fed mice. Moreover, HFD-fed mice exposed to PFDA had impaired glucose metabolism, inflammation and hepatic lipid accumulation compared to mice fed HFD alone. PFDA activated the expression of hepatic NLRP3 and caspase-1, and induced that of SREBP-1c expression in the liver of HFD-fed mice. PFDA exposure in HFD-fed mice significantly inhibited hepatic AMPK expression than animals fed HFD without PFDA exposure. Furthermore, MCC950, an NLRP3 inhibitor, suppressed the upregulation of NLRP3 and caspase-1 expression, and inhibited the expression of SREBP-1c and the accumulation of hepatic lipid in mice exposed to PFDA. Thus, PFDA may enhance HFD-induced adiposity and hepatic lipid accumulation through the NLRP3/caspase-1 pathway. This contaminant may be a key risk factor for obesity development in individuals consuming high-fat foods, particularly Western diet.

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.

Perfluoroalkyl substance exposure is associated with asthma and innate immune cell count in US adolescents stratified by sex

Exposure to perfluoroalkyl substances (PFAS) may be harmful to humans; however, previous studies have been inconsistent regarding the potential for PFAS-induced immunosuppresion. This study explored the relationship between PFAS exposure and risks of asthma, wheezing, and immunosuppression in 12-19 year-olds using the National Health and Nutrition Examination Survey (NHANES) data. Logistic regression models were used to reveal associations between serum PFAS levels and risks of asthma, wheezing, asthma attack, and emergency department visits. Pearson's correlation was used to determine the relationship between serum PFAS levels and leukocyte count. Data were also stratified by sex. We found that medium-low levels of serum perfluorooctane sulfonate (PFOS) (6.90-12.40 ng/mL) and serum perfluorooctanoic acid (PFOA) (2.43-3.60 ng/mL) were negatively related, respectively, to current asthma and wheezing in boys, and to wheezing in girls. Meanwhile, boys with medium-high levels (1.50-3.00 ng/mL) of serum perfluorohexanesulfonate (PFHxS) had a high risk of wheezing. Among asthmatic participants, both medium-high levels (3.75-5.07 ng/mL) of serum PFOA and high levels (> 3.92 ng/mL) of PFHxS correlated with asthma attacks in boys; likewise, medium-low levels (0.70-0.99 ng/mL) of serum PFNA correlated with asthma attacks in girls. Also, PFOA and PFNA levels were weakly positively correlated with basophil count, whereas PFOS levels were weakly negatively correlated with eosinophils in asthmatic boys, indicating that basophils may be important in the immune response to PFAS exposure among asthmatics.

Application of Minimal Physiologically-Based Pharmacokinetic Model to Simulate Lung and Trachea Exposure of Pyronaridine and Artesunate in Hamsters

A fixed-dose combination of pyronaridine and artesunate, one of the artemisinin-based combination therapies, has been used as a potent antimalarial treatment regimen. Recently, several studies have reported the antiviral effects of both drugs against severe acute respiratory syndrome coronavirus two (SARS-CoV-2). However, there are limited data on the pharmacokinetics (PKs), lung, and trachea exposures that could be correlated with the antiviral effects of pyronaridine and artesunate. The purpose of this study was to evaluate the pharmacokinetics, lung, and trachea distribution of pyronaridine, artesunate, and dihydroartemisinin (an active metabolite of artesunate) using a minimal physiologically-based pharmacokinetic (PBPK) model. The major target tissues for evaluating dose metrics are blood, lung, and trachea, and the nontarget tissues were lumped together into the rest of the body. The predictive performance of the minimal PBPK model was evaluated using visual inspection between observations and model predictions, (average) fold error, and sensitivity analysis. The developed PBPK models were applied for the multiple-dosing simulation of daily oral pyronaridine and artesunate. A steady state was reached about three to four days after the first dosing of pyronaridine and an accumulation ratio was calculated to be 1.8. However, the accumulation ratio of artesunate and dihydroartemisinin could not be calculated since the steady state of both compounds was not achieved by daily multiple dosing. The elimination half-life of pyronaridine and artesunate was estimated to be 19.8 and 0.4 h, respectively. Pyronaridine was extensively distributed to the lung and trachea with the lung-to-blood and trachea-to-blood concentration ratios (=Cavg,tissue/Cavg,blood) of 25.83 and 12.41 at the steady state, respectively. Also, the lung-to-blood and trachea-to-blood AUC ratios for artesunate (dihydroartemisinin) were calculated to be 3.34 (1.51) and 0.34 (0.15). The results of this study could provide a scientific basis for interpreting the dose–exposure–response relationship of pyronaridine and artesunate for COVID-19 drug repurposing.

A Scoping Assessment of Implemented Toxicokinetic Models of Per- and Polyfluoro-Alkyl Substances, with a Focus on One-Compartment Models

Toxicokinetic (TK) models have been used for decades to estimate concentrations of per-and polyfluoroalkyl substances (PFAS) in serum. However, model complexity has varied across studies depending on the application and the state of the science. This scoping effort seeks to systematically map the current landscape of PFAS TK models by categorizing different trends and similarities across model type, PFAS, and use scenario. A literature review using Web of Science and SWIFT-Review was used to identify TK models used for PFAS. The assessment covered publications from 2005-2020. PFOA, the PFAS for which most models were designed, was included in 69 of the 92 papers, followed by PFOS with 60, PFHxS with 22, and PFNA with 15. Only 4 of the 92 papers did not include analysis of PFOA, PFOS, PFNA, or PFHxS. Within the corpus, 50 papers contained a one-compartment model, 17 two-compartment models were found, and 33 used physiologically based pharmacokinetic (PBTK) models. The scoping assessment suggests that scientific interest has centered around two chemicals-PFOA and PFOS-and most analyses use one-compartment models in human exposure scenarios.

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.

Perfluoroalkyl substances (PFASs) are substrates of the renal human organic anion transporter 4 (OAT4)

Poly- and perfluoroalkyl substances (PFASs) are omnipresent in the environment and have been shown to accumulate in humans. Most PFASs are not biotransformed in animals and humans, so that elimination is largely dependent on non-metabolic clearance via bile and urine. Accumulation of certain PFASs in humans may relate to their reabsorption from the pre-urine by transporter proteins in the proximal tubules of the kidney, such as URAT1 and OAT4. The present study assessed the in vitro transport of 7 PFASs (PFHpA, PFOA, PFNA, PFDA, PFBS, PFHxS and PFOS) applying URAT1- or OAT4-transfected human embryonic kidney (HEK) cells. Virtually no transport of PFASs could be measured in URAT1-transfected HEK cells. All PFASs, except PFBS, showed clear uptake in OAT4-transfected HEK cells. In addition, these in vitro results were further supported by in silico docking and molecular dynamic simulation studies assessing transporter-ligand interactions. Information on OAT4-mediated transport may provide insight into the accumulation potential of PFASs in humans, but other kinetic aspects may play a role and should also be taken into account. Quantitative information on all relevant kinetic processes should be integrated in physiologically based kinetic (PBK) models, to predict congener-specific accumulation of PFASs in humans in a more accurate manner.

Perfluorodecanoic acid promotes adipogenesis via NLRP3 inflammasome-mediated pathway in HepG2 and 3T3-L1 cells

Perfluorodecanoic acid (PFDA) is a toxic persistent pollutant that is extensively used in food applications, such as food packaging and cookware. Emerging evidence indicates that PFDA exposure were associated with higher plasma triglyceride concentration in human. In contrast, it is unknown how PFDA might affect adipogenesis. To explore the effects and underlying mechanisms of PFDA on lipid metabolism in this study, both HepG2 cells and 3T3-L1 differentiation model were used. The results showed that PFDA promoted the cellular triglyceride accumulation and triglyceride content in concentration-dependent manners. Furthermore, PFDA activated the NLRP3 inflammasome, which is crucial for the induction of lipogenic genes expression including fatty acid synthase (FAS), hydroxymethylglutaryl coenzyme A synthase (HMGCS), and stearoyl-CoA desaturase 1 (SCD1). Additionally, PFDA-induced adipogenesis was abolished by caspase-1 inhibitor and siNLRP3 in HepG2 cells. Moreover, after PFDA treatment, the expression of SREBP1, an important regulator of lipid metabolism, was increased, as well as its target genes, and PFDA-induced SREBP1 enhanced expression can be abolished by caspase-1 inhibitor and siNLRP3 as well. Together, these results provide to understanding of the potential health implications of exposure to PFDA on lipid accumulation, and suggest that PFDA can promote adipogenesis via an NLRP3 inflammasome-mediated SREBP1 pathway.

Simulating toxicokinetic variability to identify susceptible and highly exposed populations

BACKGROUND

Toxicokinetic (TK) data needed for chemical risk assessment are not available for most chemicals. To support a greater number of chemicals, the U.S. Environmental Protection Agency (EPA) created the open-source R package "httk" (High Throughput ToxicoKinetics). The "httk" package provides functions and data tables for simulation and statistical analysis of chemical TK, including a population variability simulator that uses biometrics data from the National Health and Nutrition Examination Survey (NHANES).

OBJECTIVE

Here we modernize the "HTTK-Pop" population variability simulator based on the currently available data and literature. We provide explanations of the algorithms used by "httk" for variability simulation and uncertainty propagation.

METHODS

We updated and revised the population variability simulator in the "httk" package with the most recent NHANES biometrics (up to the 2017-18 NHANES cohort). Model equations describing glomerular filtration rate (GFR) were revised to more accurately represent physiology and population variability. The model output from the updated "httk" package was compared with the current version.

RESULTS

The revised population variability simulator in the "httk" package now provides refined, more relevant, and better justified estimations.

SIGNIFICANCE

Fulfilling the U.S. EPA's mission to provide open-source data and models for evaluations and applications by the broader scientific community, and continuously improving the accuracy of the "httk" package based on the currently available data and literature.

Per- and polyfluoroalkyl substances (PFASs) in the blood of police and Beagle dogs from Harbin, China: Concentrations and associations with hematological parameters

Per-and polyfluoroalkyl substances (PFASs) have been omnipresent in the environment and marine organisms. However, little is known about these compounds and their associations with hematological parameters in dogs. In this research, we investigated the concentrations and distributions of PFASs in the blood of dogs and explored the associations between PFASs concentrations in blood and hematological parameters. Perfluorooctanoate (PFOA) was found to be the dominant PFAS in the blood (54.23%), followed by perfluorobutyric acid (PFBA) (16.05%) and perfluorooctanesulfonate (PFOS) (12.05%). On average, PFASs concentration in dogs was 3.553 ng/mL (SD: 2.146). Moreover, age is a key factor influencing the levels of PFBA, PFOA, and PFBS in males, as well as seven PFASs (6:2 Cl-PFESA, PFBA, PFOA, PFOS, PFHxS, PFDA, and PFNA) in female dogs' blood. The results revealed that PFHxS in dietary food accounted for most of the total daily PFASs consumption. We also discovered that greater PFASs exposure (including PFOA and PFOS) could significantly increase amylase (AMY) and decrease cholesterol (CHOL) levels. Furthermore, there are linear relationships between PFDA, PFNA and many biochemical parameters (AMY, CHOL, albumin/globulin (A/G), blood urea nitrogen (BUN), alkaline phosphatase (ALP), creatinine (CREA)). Thus, PFAS accumulation has a certain influence on dogs' health, and we must pay attention to the potential threat posed by these elements to dogs.

Human risk assessment of 4-n-nonylphenol (4-n-NP) using physiologically based pharmacokinetic (PBPK) modeling: analysis of gender exposure differences and application to exposure analysis related to large exposure variability in population

As a toxic substance, 4-n-nonylphenol (4-n-NP) or 4-nonylphenol (4-NP) is widely present in the environment. 4-n-NP is a single substance with a linear-alkyl side chain, but 4-NP usually refers to a random mixture containing various branched types. Unfortunately, human risk assessment and/or exposure level analysis for 4-n-NP (or 4-NP) were almost nonexistent, and related research was urgently needed. This study aimed to analyze the various exposures of 4-n-NP (or 4-NP) through development of a physiologically based-pharmacokinetic (PBPK) model considering gender difference in pharmacokinetics of 4-n-NP and its application to human risk assessment studies. A PBPK model was newly developed considering gender differences in 4-n-NP pharmacokinetics and applied to a human risk assessment for each gender. Exposure analysis was performed using a PBPK model that considered gender differences in 4-n-NP (or 4-NP) exposure and high variabilities in several countries. Furthermore, an extended application was attempted as a human risk assessment for random mixture 4-NP, which is difficult to accurately evaluate in reality. External-exposure and margin-of-safety estimated with the same internal exposure amount differed between genders, meaning the need for a differentiated risk assessment considering gender. Exposure analysis based on biomonitoring data confirmed large variability in exposure to 4-n-NP (or 4-NP) by country, group, and period. External-exposures estimated using PBPK model varied widely, ranging from 0.039 to 63.875 mg/kg/day (for 4-n-NP or 4-NP). By country, 4-n-NP (or 4-NP) exposure was higher in females than in males and the margin-of-safety tended to be low. Overall, exposure to 4-n-NP (or 4-NP) in populations was largely not safe, suggesting need for ongoing management and monitoring. Considering low in vivo accumulation confirmed by PBPK model, risk reduction of 4-n-NP is possible by reducing its use.

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.

Half-lives of several polyfluoroalkyl substances (PFAS) in cattle serum and tissues

Cattle that were at steady-state serum polyfluoroalkyl substances (PFAS) concentrations due to several years of exposure to water contaminated by residues of Aqueous Film-Forming (AFFF) firefighting foam had perfluorooctane sulphonate (PFOS) isomers, perfluoroheptane sulphonate (PFHpS), perfluorohexane sulphonate (PFHxS), perfluorononanoic acid (PFNA) and perfluorodecanoic acid (PFDA) in serum. Elimination serum half-lives were determined in five heifers from serial blood sampling over 215 days. Eleven additional animals that had blood sampled on day 19 (d19) were euthanised on d63. PFAS half-life estimates from the serial blood sampling and from d19/d63 data were not significantly different. The combined (n = 16) serum half-lives (in days) were: total PFOS (tPFOS, 74.1 ± 13.4), PFHpS (45.7 ± 9.4), PFHxS (9.3 ± 1.3), PFNA (12.3 ± 3.2) and PFDA (60.4 ± 10.4). The half-lives of linear PFOS (L-PFOS, 69.4 ± 11.6) and mono branched PFOS isomers (m-PFOS, 83.6 ± 19) were not significantly different from tPFOS, but for the di-branched isomers (di-PFOS), the serum half-life was significantly lower (29.9 ± 5.8). Animal age (1.4-12.3 years old) and serum concentration at the start of depuration did not influence half-lives, and there was no difference between steers and heifers. Consideration of serum and tissue PFAS concentrations at d63 and d215 indicated there was no difference in tPFOS depuration from serum or muscle, but elimination from liver and kidney may be slightly longer. Depuration of PFHpS is essentially the same in serum, kidney and liver, and it is expected depletion from muscle would be comparable. The short half-life of di-PFOS, PFHxS and PFNA did not allow an assessment of clearance from tissues because they were not measurable at d215 but based on the results for PFOS and PFHpS, elimination of PFHxS from tissues is expected to mirror that from serum. Human health risk assessment implications are discussed.

A Model Template Approach for Rapid Evaluation and Application of Physiologically Based Pharmacokinetic Models for Use in Human Health Risk Assessments: A Case Study on Per- and Polyfluoroalkyl Substances

Physiologically based pharmacokinetic (PBPK) models are commonly used in risk assessments to perform inter- and intraspecies extrapolations as well as to extrapolate between different dosing scenarios; however, they must first undergo quality assurance review, which can be a time-consuming process, especially when model code is not readily available. We developed and implemented (using R and MCSim) a PBPK model template capable of replicating published model results for several chemical-specific PBPK models. This model template allows for faster quality assurance review because the general model equations only need to be reviewed once, and application to a specific chemical then only requires reviewing input parameters. The model template can implement PBPK models with oral and intravenous exposure routes, varying numbers of tissue compartments, renal reabsorption, and multiple elimination pathways, including fecal, urinary, and biliary. Using the model template, we reproduced published model simulation results for perfluorohexanesulfonic acid, perfluorononanoic acid, perfluorodecanoic acid, perfluorooctanoate, and perflouorooctane sulfonate. We also show that the template can be a useful tool for identifying potential model errors. Thus, the model template allows for faster evaluation and review of published PBPK models and provides a proof of concept for using this approach with broader classes of chemical-specific PBPK models.

Biomarkers of poly- and perfluoroalkyl substances (PFAS) in Sub-Arctic and Arctic communities in Canada

Polyfluoroalkyl substances and perfluoroalkyl substances (PFAS) are a family of anthropogenic chemicals that are used in food packaging, waterproof clothing, and firefighting foams for their water and oil resistant properties. Though levels of some PFAS appear to be decreasing in Canada's south, environmental levels have been increasing in the Arctic due to long-range transport. However, the implications of this on human exposures in sub-Arctic and Arctic populations in Canada have yet to be established. To address this data gap, human biomonitoring research was completed in Old Crow, Yukon, and the Dehcho region, Northwest Territories. Blood samples were collected from adults residing in seven northern First Nations and were analyzed by liquid chromatography mass spectrometry. A total of nine PFAS were quantified: perfluorooctanoic acid (PFOA), perfluorooctane sulphonic acid (PFOS), perfluorohexane sulphonic acid (PFHxS), perfluorononanoic acid (PFNA), perfluorodecanoic acid (PFDA), and perfluoroundecanoic acid (PFUdA), perfluorobutanoic acid (PFBA), perfluorohexanoic acid (PFHxA), and perfluorobutane sulphonic acid (PFBS). In the Dehcho (n = 124), five PFAS had a detection rate greater than 50% including PFOS, PFOA, PFHxS, PFNA, and PFDA. In addition to these PFAS, PFUdA was also detected in at least half of the samples collected in Old Crow (n = 54). Generally, male participants had higher concentrations of PFAS compared to female participants, and PFAS concentrations tended to increase with age. For most PFAS, Old Crow and Dehcho levels were similar or lower to those measured in the general Canadian population (as measured through the Canadian Health Measures Survey or CHMS) and other First Nations populations in Canada (as measured through the First Nations Biomonitoring Initiative or FNBI). The key exception to this was for PFNA which, relative to the CHMS (0.51 μg/L), was approximately 1.8 times higher in Old Crow (0.94 μg/L) and 2.8 times higher in Dehcho (1.42 μg/L) than observed in the general Canadian population. This project provides baseline PFAS levels for participating communities, improving understanding of human exposures to PFAS in Canada. Future research should investigate site-specific PFNA exposure sources and monitor temporal trends in these regions.

Human risk assessment of di-isobutyl phthalate through the application of a developed physiologically based pharmacokinetic model of di-isobutyl phthalate and its major metabolite mono-isobutyl phthalate

Di-isobutyl phthalate (DiBP) is a substance used in the production of objects frequently used in human life. Mono-isobutyl phthalate (MiBP), a major in vivo metabolite of DiBP, is a biomarker for DiBP exposure assessment. Therefore, risk assessment studies on DiBP and MiBP, which have not yet been reported in detail, are needed. The aim of this study was to develop and evaluate a physiologically based pharmacokinetic (PBPK) model for DiBP and MiBP in rats and extend this to human risk assessment based on human exposure. Pharmacokinetic studies were performed in male rats following the administration of 5-100 mg/kg DiBP, and these results were used for the development and validation of the PBPK model. In addition, the previous pharmacokinetic results in female rats following DiBP administration and the pharmacokinetic results in both males and females according to multiple exposures to DiBP were used to develop and validate the PBPK model. The metabolism of DiBP to MiBP in the body was very significant and rapid, and the biodistribution of MiBP was broad and major. Furthermore, the amount of MiBP in the body showed a correlation with DiBP exposure, and from this, a PBPK model was developed to evaluate the external exposure of DiBP from the internal exposure of MiBP. The predicted rat plasma, urine, fecal, and tissue concentrations using the developed PBPK model fitted well with the observed values. The established PBPK model for rats was extrapolated to a human PBPK model of DiBP and MiBP based on human physiological parameters and allometric scaling. The reference dose of 0.512 mg/kg/day of DiBP and external doses of 6.14-280.90 μg/kg/day DiBP for human risk assessment were estimated using Korean biomonitoring values. Valuable insight and approaches to assessing human health risks associated with DiBP exposure were provided by this study.

EFSA Panel on Contaminants in the Food Chain (EFSA CONTAM Panel), Schrenk D, Bignami M, Bodin L, Chipman JK, del Mazo J, Grasl‐Kraupp B, Hogstrand C, Hoogenboom L(R, Leblanc J, Nebbia CS, Nielsen E, Ntzani E, Petersen A, Sand S, Vleminckx C, Wallace H, Barregård L, Ceccatelli S, Cravedi J, Halldorsson TI, Haug LS, Johansson N, Knutsen HK, Rose M, Roudot A, Van Loveren H, Vollmer G, Mackay K, Riolo F, Schwerdtle T. Risk to human health related to the presence of perfluoroalkyl substances in food

The European Commission asked EFSA for a scientific evaluation on the risks to human health related to the presence of perfluoroalkyl substances (PFASs) in food. Based on several similar effects in animals, toxicokinetics and observed concentrations in human blood, the CONTAM Panel decided to perform the assessment for the sum of four PFASs: PFOA, PFNA, PFHxS and PFOS. These made up half of the lower bound (LB) exposure to those PFASs with available occurrence data, the remaining contribution being primarily from PFASs with short half-lives. Equal potencies were assumed for the four PFASs included in the assessment. The mean LB exposure in adolescents and adult age groups ranged from 3 to 22, the 95th percentile from 9 to 70 ng/kg body weight (bw) per week. Toddlers and 'other children' showed a twofold higher exposure. Upper bound exposure was 4- to 49-fold higher than LB levels, but the latter were considered more reliable. 'Fish meat', 'Fruit and fruit products' and 'Eggs and egg products' contributed most to the exposure. Based on available studies in animals and humans, effects on the immune system were considered the most critical for the risk assessment. From a human study, a lowest BMDL 10 of 17.5 ng/mL for the sum of the four PFASs in serum was identified for 1-year-old children. Using PBPK modelling, this serum level of 17.5 ng/mL in children was estimated to correspond to long-term maternal exposure of 0.63 ng/kg bw per day. Since accumulation over time is important, a tolerable weekly intake (TWI) of 4.4 ng/kg bw per week was established. This TWI also protects against other potential adverse effects observed in humans. Based on the estimated LB exposure, but also reported serum levels, the CONTAM Panel concluded that parts of the European population exceed this TWI, which is of concern.

Risk assessment for humans using physiologically based pharmacokinetic model of diethyl phthalate and its major metabolite, monoethyl phthalate

Diethyl phthalate (DEP) belongs to phthalates with short alkyl chains. It is a substance frequently used to make various products. Thus, humans are widely exposed to DEP from the surrounding environment such as food, soil, air, and water. As previously reported in many studies, DEP is an endocrine disruptor with reproductive toxicity. Monoethyl phthalate (MEP), a major metabolite of DEP in vivo, is a biomarker for DEP exposure assessment. It is also an endocrine disruptor with reproductive toxicity, similar to DEP. However, toxicokinetic studies on both MEP and DEP have not been reported in detail yet. Therefore, the objective of this study was to evaluate and develop physiologically based pharmacokinetic (PBPK) model for both DEP and MEP in rats and extend this to human risk assessment based on human exposure. This study was conducted in vivo after intravenous or oral administration of DEP into female (2 mg/kg dose) and male (0.1-10 mg/kg dose) rats. Biological samples consisted of urine, plasma, and 11 different tissues. These samples were analyzed using UPLC-ESI-MS/MS method. For DEP, the tissue to plasma partition coefficient was the highest in the kidney, followed by that in the liver. For MEP, the tissue to plasma partition coefficient was the highest in the liver. It was less than unity in all other tissues. Plasma, urine, and fecal samples were also obtained after IV administration of MEP (10 mg/kg dose) to male rats. All results were reflected in a model developed in this study, including in vivo conversion from DEP to MEP. Predicted concentrations of DEP and MEP in rat urine, plasma, and tissue samples using the developed PBPK model fitted well with observed values. We then extrapolated the PBPK model in rats to a human PBPK model of DEP and MEP based on human physiological parameters. Reference dose of 0.63 mg/kg/day (or 0.18 mg/kg/day) for DEP and external doses of 0.246 μg/kg/day (pregnant), 0.193 μg/kg/day (fetus), 1.005-1.253 μg/kg/day (adults), 0.356-0.376 μg/kg/day (adolescents), and 0.595-0.603 μg/kg/day (children) for DEP for human risk assessment were estimated using Korean biomonitoring values. Our study provides valuable insight into human health risk assessment regarding DEP exposure.