Sex-specific risk assessment of PFHxS using a physiologically based pharmacokinetic model

Towards process-based modelling and parameterisation of bioaccumulation in humans across PFAS congeners

Per- and polyfluoroalkyl substances (PFAS) are a large class of stable toxic chemicals which have ended up in the environment and in organisms in significant concentrations. Toxicokinetic models are needed to facilitate extrapolation of bioaccumulation data across PFAS congeners and species. For the present study, we carried out an inventory of accumulation processes specific for PFAS, deviating from traditional Persistent Organic Pollutants (POPs). In addition, we reviewed toxicokinetic models on PFAS reported in literature, classifying them according to the number of compartments distinguished as a one-compartment model (1-CM), two-compartment model (2- CM) or a multi-compartment model, (multi-CM) as well as the accumulation processes included and the parameters used. As the inventory showed that simple 1-CMs were lacking, we developed a generic 1-CM of ourselves to include PFAS specific processes and validated the model for legacy perfluoroalkyl acids. Predicted summed elimination constants were accurate for long carbon chains (>C6), indicating that the model properly represented toxicokinetic processes for most congeners. Results for urinary elimination rate constants were mixed, which might be caused by the exclusion of reabsorption processes (renal reabsorption, enterohepatic circulation). The 1-CM needs to be improved further in order to better predict individual elimination pathways. Besides that, more data on PFAS-transporter specific processes are needed to extrapolate across PFAS congeners and species.

Evaluation of the Effect of Perfluorohexane Sulfonate on the Proliferation of Human Liver Cells

Perfluorohexane sulfonate (PFHxS) is a widely detected replacement for legacy long-chain perfluoroalkyl substances (PFAS) in the environment and human blood samples. Its potential toxicity led to its recent classification as a globally regulated persistent organic pollutant. Although animal studies have shown a positive association between PFHxS levels and hepatic steatosis and hepatocellular hypertrophy, the link with liver toxicity, including end-stage liver cancer, remains inconclusive. In this study, we examined the effects of PFHxS on the proliferation of Hep3B (human hepatocellular carcinoma) and SK-Hep1 (human liver sinusoidal endothelial cells). Cells were exposed to different PFHxS concentrations for 24-48 h to assess viability and 12-14 days to measure colony formation. The viability of both cell lines increased at PFHxS concentrations <200 μM, decreased at >400 μM, and was highest at 50 μM. Colony formation increased at <300 μM and decreased at 500 μM PFHxS. Consistent with the effect on cell proliferation, PFHxS increased the expression of proliferating cell nuclear antigen (PCNA) and cell-cycle molecules (CDK2, CDK4, cyclin E, and cyclin D1). In summary, PFHxS exhibited a biphasic effect on liver cell proliferation, promoting survival and proliferation at lower concentrations and being cytotoxic at higher concentrations. This suggests that PFHxS, especially at lower concentrations, might be associated with HCC development and progression.

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.

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.

Per- and polyfluoroalkyl substances exposure and its influence on the intestinal barrier: An overview on the advances

Per- and polyfluoroalkyl substances (PFAS) are a class of artificially synthetic organic compounds that are hardly degraded in the natural environment. PFAS have been widely used for many decades, and the persistence and potential toxicity of PFAS are an emerging concern in the world. PFAS exposed via diet can be readily absorbed by the intestine and enter the circulatory system or accumulate directly at intestinal sites, which could interact with the intestine and cause the destruction of intestinal barrier. This review summarizes current relationships between PFAS exposure and intestinal barrier damage with a focus on more recent toxicological studies. Exposure to PFAS could cause inflammation in the gut, destruction of the gut epithelium and tight junction structure, reduction of the mucus layer, and induction of the toxicity of immune cells. PFAS accumulation could also induce microbial disorders and metabolic products changes. In addition, there are limited studies currently, and most available studies converge on the health risk of PFAS exposure for human intestinal disease. Therefore, more efforts are deserved to further understand potential associations between PFAS exposure and intestinal dysfunction and enable better assessment of exposomic toxicology and health risks for humans in the future.

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.

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.

In Vivo Transformation of a Novel Polyfluoroether Surfactant

Per- and polyfluoroalkyl substances are a class of fluorochemicals that can degrade into perfluoroalkyl acids, which are well known to be persistent in the environment. It is thus important that novel fluorinated surfactants be designed to degrade into small, nonbioaccumulative products. We report the biotransformation and elimination kinetics of one such novel polyfluorinated surfactant, di(polyfluoroether thioether(S)-oate) sulfonate (diFESOS), and its metabolites. Biotransformation was investigated in vitro using S9 liver fractions and in vivo in Sprague-Dawley rats. Rats dosed by oral gavage with diFESOS were found to have relatively fast elimination kinetics, with half-lives on the order of hours, compared with legacy fluorinated surfactants such as the disubstituted polyfluoroalkyl phosphates that have half-lives on the order of days. To interrogate degradation of the polyfluorinated chain, rats were then dosed with a polyfluoroether thioether alcohol (a suspected product of carboxylate cleavage of diFESOS) either orally or intravenously, and the novel metabolite 2H-3:2 polyfluoroether sulfonic acid (2H-3:2 PFESA) was identified. Perfluoropropionic acid was detected in rat urine and is likely a terminal product. The blood of orally dosed rats contained higher levels of metabolites than the blood of intravenously dosed rats, suggesting the importance of metabolism in the gut and liver. Elimination kinetics of all the novel metabolites were faster than their fully fluorinated counterparts. Environ Toxicol Chem 2021;40:3328-3336. © 2021 SETAC.

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.

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.

Endocrine Disruptor Potential of Short- and Long-Chain Perfluoroalkyl Substances (PFASs)-A Synthesis of Current Knowledge with Proposal of Molecular Mechanism

Endocrine disruptors are a group of chemical compounds that, even in low concentrations, cause a hormonal imbalance in the body, contributing to the development of various harmful health disorders. Many industry compounds, due to their important commercial value and numerous applications, are produced on a global scale, while the mechanism of their endocrine action has not been fully understood. In recent years, per- and polyfluoroalkyl substances (PFASs) have gained the interest of major international health organizations, and thus more and more studies have been aimed to explain the toxicity of these compounds. PFASs were firstly synthesized in the 1950s and broadly used in the industry in the production of firefighting agents, cosmetics and herbicides. The numerous industrial applications of PFASs, combined with the exceptionally long half-life of these substances in the human body and extreme environmental persistence, result in a common and chronic exposure of the general population to their action. Available data have suggested that human exposure to PFASs can occur during different stages of development and may cause short- or/and long-term health effects. This paper synthetizes the current literature reports on the presence, bioaccumulation and, particularly, endocrine toxicity of selected long- and short-chain PFASs, with a special emphasis on the mechanisms underlying their endocrine actions.

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.

Occurrence, partition and environmental risk assessment of per- and polyfluoroalkyl substances in water and sediment from the Baiyangdian Lake, China

This work examined the contamination of poly- and perfluorinated compounds (PFASs) in the water and sediment of the Baiyangdian Lake. The total concentration of PFASs in the surface water varied from 140.5 to 1828.5 ng/L, and the highest concentration of PFASs were observed near the entrance of Fuhe river. The topmost contaminant was sodium perfluorohexanesulfonate (PFHxS) and perfluorooctanoic acid (PFOA) in the north and south of the Baiyangdian Lake respectively, which indicated different contamination sources. The total concentration of PFASs in the sediment varied from 0.48 to 30 ng/g, and the distribution of PFASs in the sediment was similar with that in the surface water. The concentrations of polyfluoroalkyl phosphoric diesters (diPAPs) were three to four orders of magnitude lower than those of perfluorocarboxylates (PFCAs) and PFSAs. Although the pore water and the surface water had similar ΣPFASs, the concentration of perfluorodecanoic acid (PFDA) in pore water was 1.4 to 4.4 times higher than that in surface water, and the concentration of perfluoropentanoic acid (PFPeA) in pore water was 20–70% that in surface water. The results of ecological risk assessment showed that the PFASs were currently of no immediate risk to the aquatic life.

Occurrence and source apportionment of novel and legacy poly/perfluoroalkyl substances in Hai River basin in China using receptor models and isomeric fingerprints

A variety of fluorinated alternatives are being manufactured and applied as a consequence of stringent regulations on legacy poly/perfluoroalkyl substances (PFASs). In this study, 26 emerging and legacy PFASs were measured in the surface water (including dissolved phase and suspended particulate matter) and sediments taken from Hai River basin, China. The total concentrations of PFASs (∑PFASs) ranged from 1.74 to 172 ng/L, with perfluorooctanonate (PFOA) as the dominant compound (15.2% of the ∑PFASs, median value). Emerging PFASs, such as dimer acid of hexafluoropropylene oxide dimer acid (HFPO-DA) and trimer acid (HFPO-TA), were widely detected in the water samples. Specifically, chlorinated polyfluorinated ether sulfonate (F-53B) was observed to be predominant in some sediment samples. A receptor model, Unmix, was introduced to identify the sources of PFASs in the surface water, and the results indicated that fire-fighting foam/fluoropolymer processing aids (36.6%) were the dominant source. The field-based sediment-water (organic carbon normalized) coefficients, K_oc, were correlated to the carbon chain lengths of the PFASs. A technique coupling one-way analysis of variance with chemical mass balance model was developed to trace the manufacturing sources of PFOA. Electrochemical fluorination (ECF) was the major PFOA manufacturing source with considerable contribution by telomerization. For the first time, the isomers of perfluorooctane sulfonamide (PFOSA) were quantified in the environmental samples. The lower proportion of branched (br-) PFOSA isomers and higher percentage of br-perfluorooctane sulfonate (PFOS) isomers in the water samples relative to their corresponding commercial products, provided more direct evidences that br-PFOSA isomers were biotransformed more easily than n-PFOSA, explaining the observed enrichment of br-PFOS in the aquatic environment.

Derivation of a chronic reference dose for perfluorohexane sulfonate (PFHxS) for reproductive toxicity in mice

Perfluorohexane sulfonate (PFHxS) is a six-carbon perfluoroalkyl sulfonic acid that was used as an industrial surfactant, but is now found as an environmental contaminant worldwide. In addition to its use as an industrial surfactant, it is a legacy contaminant from the use of aqueous film-forming foams. Despite its widespread occurrence in the environment and evidence of biological activity associated with PFHxS and similar perfluoroalkyl sulfonic acids in rodents, there is no oral toxicity value currently available from the IRIS Database. To derive an oral reference dose (RfD) for PFHxS, available toxicity studies were reviewed using a weight-of-evidence approach. A 42-day mouse reproductive study was chosen as the critical study for the derivation of the oral RfD. Benchmark dose modeling was utilized to derive a point of departure (POD) for a reduction in litter size. A 95% lower confidence limit on the benchmark dose (BMDL) of 13,900 ng/mL (serum PFHxS) was modeled for a reduction in litter size. An oral RfD for PFHxS of 4.0 ng/kg/d was calculated by conversion of the BMDL to a human equivalent oral dose using a human half-life adjusted dosimetric conversion factor and the application of a total uncertainty factor of 300. Additional research is needed to better characterize the toxicity associated with oral exposure to PFHxS and refine the development of toxicity values.

Integration of Food Animal Residue Avoidance Databank (FARAD) empirical methods for drug withdrawal interval determination with a mechanistic population-based interactive physiologically based pharmacokinetic (iPBPK) modeling platform: example for flunixin meglumine administration

Violative chemical residues in animal-derived food products affect food safety globally and have impact on the trade of international agricultural products. The Food Animal Residue Avoidance Databank program has been developing scientific tools to provide appropriate withdrawal interval (WDI) estimations after extralabel drug use in food animals for the past three decades. One of the tools is physiologically based pharmacokinetic (PBPK) modeling, which is a mechanistic-based approach that can be used to predict tissue residues and WDIs. However, PBPK models are complicated and difficult to use by non-modelers. Therefore, a user-friendly PBPK modeling framework is needed to move this field forward. Flunixin was one of the top five violative drug residues identified in the United States from 2010 to 2016. The objective of this study was to establish a web-based user-friendly framework for the development of new PBPK models for drugs administered to food animals. Specifically, a new PBPK model for both cattle and swine after administration of flunixin meglumine was developed. Population analysis using Monte Carlo simulations was incorporated into the model to predict WDIs following extralabel administration of flunixin meglumine. The population PBPK model was converted to a web-based interactive PBPK (iPBPK) framework to facilitate its application. This iPBPK framework serves as a proof-of-concept for further improvements in the future and it can be applied to develop new models for other drugs in other food animal species, thereby facilitating the application of PBPK modeling in WDI estimation and food safety assessment.