Paraquat pharmacokinetics in primates and extrapolation to humans.

By extending our Paraquat (PQ) work to include primates we have implemented a modelling and simulation strategy that has enabled PQ pharmacokinetic data to be integrated into a single physiologically based pharmacokinetic (PBPK) model that enables more confident extrapolation to humans. Because available data suggested there might be differences in PQ kinetics between primates and non-primates, a radiolabelled study was conducted to characterize pharmacokinetics and excretion in Cynomolgus monkeys. Following single intravenous doses of 0.01 or 0.1 mg paraquat dichloride/kg bw, plasma PQ concentration-time profiles were dose-proportional. Excretion up to 48 h (predominantly urinary) was 82.9%, with ca. 10% remaining unexcreted. In vitro blood binding was similar across Cynomolgus monkeys, humans and rat. Our PBPK model for the rat, mouse and dog, employing a single set of PQ-specific parameters, was scaled to Cynomolgus monkeys and well represented the measured plasma concentration-time profiles over 14 days. Addition of a cartilage compartment to the model better captured the percent remaining in the monkeys at 48 h, whilst having negligible effect on model predictions for the other species. The PBPK model performed well for all four species, demonstrating there is little difference in PQ kinetics between non-primates and primates enabling a more confident extrapolation to humans. Scaling of the PBPK model to humans, with addition of a human-specific dermal submodel based on in vitro human dermal absorption data, provides a valuable tool that could be employed in defining internal dosimetry to complement human health risk assessments.

Integration of paraquat pharmacokinetic data across species using PBPK modelling.

Paraquat dichloride (PQ) is a non-selective herbicide which has been the subject of numerous toxicology studies over more than 50 years. This paper describes the development of a physiologically-based pharmacokinetic (PBPK) model of PQ kinetics for the rat, mouse and dog, firstly to aid the interpretation of studies in which no kinetic measurements were made, and secondly to enable the future extension of the model to humans. Existing pharmacokinetic data were used to develop a model for the rat and mouse. Simulations with this preliminary model were then used to identify key data gaps and to design a new blood binding study to reduce uncertainty in critical aspects of the model. The new data provided evidence to support the model structure, and its predictive performance was then assessed against dog and rat datasets not used in model development. The PQ-specific model parameters are the same for all three species, with only the physiological parameters varying between species. This consistency across species provides a strong basis for extrapolation to other species, as demonstrated here for the dog. The model enables a wide range of PQ data to be linked together to provide a broad understanding of PQ pharmacokinetics in rodents and the dog, showing that the key aspects of PQ kinetics in these species are understood and adequately encapsulated within the model.

Development of a rat and human PBPK model for bromate and estimation of human equivalent concentrations in drinking water.

A physiologically based pharmacokinetic (PBPK) model was developed to described uptake, disposition and clearance of bromate in the rat using published experimental data in rat. The rodent bromate model was extrapolated to human using species-specific physiological parameters and standard interspecies scaling of rate constants. The bromate model is kinetically linear (i.e. AUC and Cmax) across the range of drinking water concentrations used in the cancer bioassays (15 to 500 ppm). This is likely the result of the poor oral bioavailability of bromate due to high reduction rates in the intestinal tract. The bromate PBPK model was used to assess the human equivalent drinking water concentration (HEC) consistent with average plasma concentrations in the rodent bioassays. At drinking water concentrations <500 mg/L, the predicted HEC was two to three fold lower than the bioassay concentration and was dependent on the reported drinking water intake reported in the bioassay.

Quantitative bias analysis of the association between subclinical thyroid disease and two perfluoroalkyl substances in a single study.

Exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) has been associated with the occurrence of thyroid disease in some epidemiologic studies. We hypothesized that in a specific epidemiologic study based on the National Health and Nutrition Examination Survey, the association of subclinical thyroid disease with serum concentration of PFOA and PFOS was due to reverse causality. Thyroid hormone affects glomerular filtration, which in turn affects excretion of PFOA and PFOS. We evaluated this by linking a model of thyroid disease status over the lifetime to physiologically based pharmacokinetic models of PFOA and PFOS. Using Monte Carlo methods, we simulated the target study population and analyzed the data using multivariable logistic regression. The target and simulated populations were similar with respect to age, estimated glomerular filtration rate, serum concentrations of PFOA and PFOS, and prevalence of subclinical thyroid disease. Our findings suggest that in the target study the associations with subclinical hypothyroidism were overstated and the results for subclinical hyperthyroidism were, in general, understated. For example, for subclinical hypothyroidism in men, the reported odds ratio per ln(PFOS) increase was 1.98 (95% CI 1.19-3.28), whereas in the simulated data the bias due to reverse causality gave an odds ratio of 1.19 (1.16-1.23). Our results provide evidence of bias due to reverse causality in a specific cross-sectional study of subclinical thyroid disease with exposure to PFOA and PFOS among adults.

Development of a physiologically based pharmacokinetic model of diisononyl phthalate (DiNP) in pregnant rat and human.

A physiologically based pharmacokinetic (PBPK) model for di-isononyl phthalate (DiNP) was developed by adapting the existing models for di(2-ethylhexyl) phthalate (DEHP) and di-butylphthalate (DBP). Both pregnant rat and human time-course plasma and urine data were used to address the hydrolysis of DiNP in intestinal tract, plasma, and liver as well as hepatic oxidative metabolism and conjugation of the monoester and primary oxidative metabolites. Data in both rats and humans were available to inform the uptake and disposition of mono-isononyl phthalate (MiNP) as well as the three primary oxidative metabolites including hydroxy (7-OH)-, oxo (7-OXO)-, and carboxy (7-COX)-monoisononyl phthalate in plasma and urine. The DiNP model was reliable over a wide range of exposure levels in the pregnant rat as well as the two low exposure levels in humans including capturing the nonlinear behavior in the pregnant rat after repeated 750 mg/kg/day dosing. The presented DiNP PBPK model in pregnant rat and human, based upon an extensive kinetic dataset in both species, may provide a basis for assessing human equivalent exposures based upon either rodent or in vitro points of departure.

Incorporation of in vitro metabolism data and physiologically based pharmacokinetic modeling in a risk assessment for chloroprene.

Objective: To develop a physiologically based pharmacokinetic (PBPK) model for chloroprene in the mouse, rat and human, relying only on in vitro data to estimate tissue metabolism rates and partitioning, and to apply the model to calculate an inhalation unit risk (IUR) for chloroprene.Materials and methods: Female B6C3F1 mice were the most sensitive species/gender for lung tumors in the 2-year bioassay conducted with chloroprene. The PBPK model included tissue metabolism rate constants for chloroprene estimated from results of in vitro gas uptake studies using liver and lung microsomes. To assess the validity of the PBPK model, a 6-hr, nose-only chloroprene inhalation study was conducted with female B6C3F1 mice in which both chloroprene blood concentrations and ventilation rates were measured. The PBPK model was then used to predict dose measures - amounts of chloroprene metabolized in lungs per unit time - in mice and humans.Results: The mouse PBPK model accurately predicted in vivo pharmacokinetic data from the 6-hr, nose-only chloroprene inhalation study. The PBPK model was used to conduct a cancer risk assessment based on metabolism of chloroprene to reactive epoxides in the lung, the target tissue in mice. The IUR was over100-fold lower than the IUR from the EPA Integrated Risk Information System (IRIS), which was based on inhaled chloroprene concentration. The different result from the PBPK model risk assessment arises from use of the more relevant tissue dose metric, amount metabolized, rather than inhaled concentrationDiscussion and conclusions: The revised chloroprene PBPK model is based on the best available science, including new test animal in vivo validation, updated literature review and a Markov-Chain Monte Carlo analysis to assess parameter uncertainty. Relying on both mouse and human metabolism data also provides an important advancement in the use of quantitative in vitro to in vivo extrapolation (QIVIVE). Inclusion of the best available science is especially important when deriving a toxicity value based on species extrapolation for the potential carcinogenicity of a reactive metabolite.

Combining transcriptomics and PBPK modeling indicates a primary role of hypoxia and altered circadian signaling in dichloromethane carcinogenicity in mouse lung and liver.

Dichloromethane (DCM) is a lung and liver carcinogen in mice at inhalation exposures≥2000ppm. The modes of action (MOA) of these responses have been attributed to formation of genotoxic, reactive metabolite(s). Here, we examined gene expression in lung and liver from female B6C3F1 mice exposed to 0, 100, 500, 2000, 3000 and 4000ppm DCM for 90days. We also simulated dose measures - rates of DCM oxidation to carbon monoxide (CO) in lung and liver and expected blood carboxyhemoglobin (HbCO) time courses with a PBPK model inclusive of both conjugation and oxidation pathways. Expression of large numbers of genes was altered at 100ppm with maximal changes in the numbers occurring by 500 or 2000ppm. Most changes in genes common to the two tissues were related to cellular metabolism and circadian clock. At the lower concentrations, the changes in metabolism-related genes were discordant - up in liver and down in lung. These processes included organelle biogenesis, TCA cycle, and respiratory electron transport. Changes in circadian cycle genes - primarily transcription factors - showed strong concentration-related response at higher concentrations (Arntl, Npas2, and Clock were down-regulated; Cry2, Wee1, Bhlhe40, Per3, Nr1d1, Nr1d2 and Dbp) were up-regulated with similar directionality in both tissues. Overall, persistently elevated HbCO from DCM oxidation appears to cause extended periods of hypoxia, leading to altered circadian coupling to cellular metabolism. The dose response for altered circadian processes correlates with the cancer outcome. We found no evidence of changes in genes indicative of responses to cytotoxic, DNA-reactive metabolites.

Improving the risk assessment of lipophilic persistent environmental chemicals in breast milk.

Lipophilic persistent environmental chemicals (LPECs) have the potential to accumulate within a woman's body lipids over the course of many years prior to pregnancy, to partition into human milk, and to transfer to infants upon breastfeeding. As a result of this accumulation and partitioning, a breastfeeding infant's intake of these LPECs may be much greater than his/her mother's average daily exposure. Because the developmental period sets the stage for lifelong health, it is important to be able to accurately assess chemical exposures in early life. In many cases, current human health risk assessment methods do not account for differences between maternal and infant exposures to LPECs or for lifestage-specific effects of exposure to these chemicals. Because of their persistence and accumulation in body lipids and partitioning into breast milk, LPECs present unique challenges for each component of the human health risk assessment process, including hazard identification, dose-response assessment, and exposure assessment. Specific biological modeling approaches are available to support both dose-response and exposure assessment for lactational exposures to LPECs. Yet, lack of data limits the application of these approaches. The goal of this review is to outline the available approaches and to identify key issues that, if addressed, could improve efforts to apply these approaches to risk assessment of lactational exposure to these chemicals.

A hybrid CFD-PBPK model for naphthalene in rat and human with IVIVE for nasal tissue metabolism and cross-species dosimetry.

A PBPK model for naphthalene in the rat and human that incorporates a hybrid CFD-PBPK description of the upper respiratory tract was developed to support cross-species dosimetry comparisons of naphthalene concentrations and tissue normalized rate of metabolism in the nasal respiratory and olfactory epithelium, lung and liver. In vitro measurements of metabolic rates from microsomal incubations published for rat and monkey (surrogate for human) were scaled to the specific tissue based on the tissue microsomal content and volume of tissue. The model reproduces time courses for naphthalene blood concentrations from intravenous and inhalation exposures in rats and upper respiratory tract extraction data in both naïve rats and rats pre-treated to inhibit nasal metabolism. This naphthalene model was applied to estimate human equivalent inhalation concentrations (HECs) corresponding to several NOAELs or LOAELs for the non-cancer effects of naphthalene in rats. Two approaches for cross-species extrapolation were compared: (1) equivalence based on tissue naphthalene concentration and (2) equivalence based on amount metabolized per minute (normalized to tissue volume). At the NOAEL of 0.1 ppm, the regional gas dosimetry ratio (RGDR) based on naphthalene concentration was 0.18 for the dorsal olfactory region; however, the RGDR rises to 5.4 when based on the normalized amount metabolized due to the lower of expression of CYP isozymes in the nasal epithelium of primates and humans. The resulting HEC is 0.12 ppm (0.63 mg/m(3)) continuous exposure at the rat NOAEL of 0.1 ppm (6 h/day, 5 days/week).

Development of PBPK models for PFOA and PFOS for human pregnancy and lactation life stages.

Perfluoroalkyl acid carboxylates and sulfonates (PFAA) have many consumer and industrial applications. Developmental toxicity studies in animals have raised concern about potential reproductive/developmental effects of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS); however, in humans conflicting results have been reported for associations between maternal PFAA levels and these outcomes. Risk assessments and interpretation of available human data during gestation and lactation are hindered due to lack of a framework for understanding and estimating maternal, fetal, and neonatal pharmacokinetics (PK). Physiologically based pharmacokinetic (PBPK) models were developed for PFOA and PFOS for the gestation and lactation life stages in humans to understand how the physiological changes associated with development affect pharmacokinetics of these compounds in the mother, fetus, and infant. These models were derived from PBPK models for PFOA/PFOS that were previously developed for adult humans and rats during gestation and lactation and from existing human pregnancy and lactation models developed for other chemicals. The models simulated PFOA and PFOS concentrations in fetal, infant, and maternal plasma and milk, were compared to available data in humans, and also were used to estimate maternal exposure. The models reported here identified several research needs, which include (1) the identification of transporters involved in renal resorption to explain the multiyear half-lives of these compounds in humans, (2) factors affecting clearance of PFOA/PFOS during gestation and lactation, and (3) data to estimate clearance of PFOA/PFOS in infants. These models may help address concerns regarding possible adverse health effects due to PFOA/PFOS exposure in the fetus and infant and may be useful in comparing pharmacokinetics across life stages.

Incorporation of rapid association/dissociation processes in tissues into the monkey and human physiologically based pharmacokinetic models for manganese.

In earlier physiologically based pharmacokinetic (PBPK) models for manganese (Mn), the kinetics of transport of Mn into and out of tissues were primarily driven by slow rates of association and dissociation of Mn with tissue binding sites. However, Mn is known to show rapidly reversible binding in tissues. An updated Mn model for primates, following similar work with rats, was developed that included rapid association/dissociation processes with tissue Mn-binding sites, accumulation of free Mn in tissues after saturation of these Mn-binding sites and rapid rates of entry into tissues. This alternative structure successfully described Mn kinetics in tissues in monkeys exposed to Mn via various routes including oral, inhalation, and intraperitoneal, subcutaneous, or intravenous injection and whole-body kinetics and tissue levels in humans. An important contribution of this effort is showing that the extension of the rate constants for binding and cellular uptake established in the monkey were also able to describe kinetic data from humans. With a consistent model structure for monkeys and humans, there is less need to rely on cadaver data and whole-body tracer studies alone to calibrate a human model. The increased biological relevance of the Mn model structure and parameters provides greater confidence in applying the Mn PBPK models to risk assessment. This model is also well-suited to explicitly incorporate emerging information on the role of transporters in tissue disposition, intestinal uptake, and hepatobiliary excretion of Mn.

Relative contributions of endogenous and exogenous formaldehyde to formation of deoxyguanosine monoadducts and DNA-protein crosslink adducts of DNA in rat nasal mucosa.

Understanding the dose-response for formaldehyde-induced nasal cancer in rats is complicated by (1) the uneven distribution of inhaled formaldehyde across the interior surface of the nasal cavity and, (2) the presence of endogenous formaldehyde (endoF) in the nasal mucosa. In this work, we used computational fluid dynamics (CFD) modeling to predict flux of inhaled (exogenous) formaldehyde (exogF) from air into tissue at the specific locations where DNA adducts were measured. Experimental work has identified DNA-protein crosslink (DPX) adducts due to exogF and deoxyguanosine (DG) adducts due to both exogF and endoF. These adducts can be considered biomarkers of exposure for effects of endoF and exogF on DNA that may be part of the mechanism of tumor formation. We describe a computational model linking CFD-predicted flux of formaldehyde from air into tissue, and the intracellular production of endoF, with the formation of DPX and DG adducts. We assumed that, like exogF, endoF can produce DPX. The model accurately reproduces exogDPX, exogDG, and endoDG data after inhalation from 0.7 to 15 ppm. The dose-dependent concentrations of exogDPX and exogDG are predicted to exceed the concentrations of their endogenous counterparts at about 2 and 6 ppm exogF, respectively. At all concentrations examined, the concentrations of endoDPX and exogDPX were predicted to be at least 10-fold higher than that of their DG counterparts. The modeled dose-dependent concentrations of these adducts are suitable to be used together with data on the dose-dependence of cell proliferation to conduct quantitative modeling of formaldehyde-induced rat nasal carcinogenicity.

PBPK modeling to evaluate maximum tolerated doses: A case study with 3-chloroallyl alcohol.

Introduction: A physiologically based pharmacokinetic (PBPK) model for 3-chloroallyl alcohol (3-CAA) was developed and used to evaluate the design of assays for the in vivo genotoxicity of 3-CAA. Methods: Model development was supported by read across from a published PBPK model for ethanol. Read across was motivated by the expectation that 3-CAA, which like ethanol is a primary alcohol, is metabolized largely by hepatic alcohol dehydrogenases. The PBPK model was used to evaluate how two metrics of tissue dosimetry, maximum blood concentration (Cmax; mg/L) and area under the curve (AUC; mg-hr/L) vary with dose of 3-CAA and with dose route (oral gavage, drinking water). Results: The model predicted that oral gavage results in a 6-fold higher Cmax than the same dose administered in drinking water, but in similar AUCs. Predicted Cmax provided the best correlation with severe toxicity (e.g., lethality) from 3-CAA, consistent with the production of a reactive metabolite. Therefore, drinking water administration can achieve higher sustained concentration without severe toxicity in vivo. Discussion: This evaluation is significant because cytotoxicity is a potential confounder of mutagenicity testing. The PBPK model can be used to ensure that studies meet OECD and USEPA test guidelines and that the highest dose used is not associated with severe toxicity. In addition, PBPK modeling provides assurance of target tissue (e.g., bone marrow) exposure even in the absence of laboratory data, by defining the relationship between applied dose and target tissue dose based on accepted principles of pharmacokinetics, relevant physiology and biochemistry of the dosed animals, and chemical-specific information.

Quantitative in vitro to in vivo extrapolation of cell-based toxicity assay results.

The field of toxicology is currently undergoing a global paradigm shift to use of in vitro approaches for assessing the risks of chemicals and drugs in a more mechanistic and high throughput manner than current approaches relying primarily on in vivo testing. However, reliance on in vitro data entails a number of new challenges associated with translating the in vitro results to corresponding in vivo exposures. Physiologically based pharmacokinetic (PBPK) modeling provides an effective framework for conducting quantitative in vitro to in vivo extrapolation (QIVIVE). Their physiological structure facilitates the incorporation of in silico- and in vitro-derived chemical-specific parameters in order to predict in vivo absorption, distribution, metabolism and excretion. In particular, the combination of in silico- and in vitro parameter estimation with PBPK modeling can be used to predict the in vivo exposure conditions that would produce chemical concentrations in the target tissue equivalent to the concentrations at which effects were observed with in vitro assays of tissue/organ toxicity. This review describes the various elements of QIVIVE and highlights key aspects of the process, with an emphasis on extrapolation of in vitro metabolism data to predict in vivo clearance as the key element. Other important elements include characterization of free concentration in the toxicity assay and potential complications associated with intestinal absorption and renal clearance. Examples of successful QIVIVE approaches are described ranging from a simple steady-state approach that is suitable for a high throughput environment to more complicated approaches requiring full PBPK models.

Development of a physiologically based pharmacokinetic model for inhalation of jet fuels in the rat.

The pharmacokinetic behavior of the majority of jet fuel constituents has not been previously described in the framework of a physiologically based pharmacokinetic (PBPK) model for inhalation exposure. Toxic effects have been reported in multiple organ systems, though exposure methods varied across studies, utilizing either vaporized or aerosolized fuels. The purpose of this work was to assess the pharmacokinetics of aerosolized and vaporized fuels, and develop a PBPK model capable of describing both types of exposures. To support model development, n-tetradecane and n-octane exposures were conducted at 89 mg/m(3) aerosol+vapor and 1000-5000 ppm vapor, respectively. Exposures to JP-8 and S-8 were conducted at ~900-1000 mg/m(3), and ~200 mg/m(3) to a 50:50 blend of both fuels. Sub-models were developed to assess the behavior of representative constituents and grouped unquantified constituents, termed "lumps", accounting for the remaining fuel mass. The sub-models were combined into the first PBPK model for petroleum and synthetic jet fuels. Inhalation of hydrocarbon vapors was described with simple gas-exchange assumptions for uptake and exhalation. For aerosol droplets systemic uptake occurred in the thoracic region. Visceral tissues were described using perfusion and diffusion-limited equations. The model described kinetics at multiple fuel concentrations, utilizing a chemical "lumping" strategy to estimate parameters for fractions of speciated and unspeciated hydrocarbons and gauge metabolic interactions. The model more accurately simulated aromatic and lower molecular weight (MW) n-alkanes than some higher MW chemicals. Metabolic interactions were more pronounced at high (~2700-1000 mg/m(3)) concentrations. This research represents the most detailed assessment of fuel pharmacokinetics to date.

Development of multi-route physiologically-based pharmacokinetic models for ethanol in the adult, pregnant, and neonatal rat.

Biofuel blends of 10% ethanol (EtOH) and gasoline are common in the USA, and higher EtOH concentrations are being considered (15-85%). Currently, no physiologically-based pharmacokinetic (PBPK) models are available to describe the kinetics of EtOH-based biofuels. PBPK models were developed to describe life-stage differences in the kinetics of EtOH alone in adult, pregnant, and neonatal rats for inhalation, oral, and intravenous routes of exposure, using data available in the open literature. Whereas ample data exist from gavage and intravenous routes of exposure, kinetic data from inhalation exposures are limited, particularly at concentrations producing blood and target tissue concentrations associated with developmental neurotoxicity. Compared to available data, the three models reported in this paper accurately predicted the kinetics of EtOH, including the absorption, peak concentration, and clearance across multiple datasets. In general, model predictions for adult and pregnant animals matched inhalation and intravenous datasets better than gavage data. The adult model was initially better able to predict the time-course of blood concentrations than was the neonatal model. However, after accounting for age-related changes in gastric uptake using the calibrated neonate model, simulations consistently reproduced the early kinetic behavior in blood. This work provides comprehensive multi-route life-stage models of EtOH pharmacokinetics and represents a first step in development of models for use with gasoline-EtOH blends, with additional potential applicability in investigation of the pharmacokinetics of EtOH abuse, addiction, and toxicity.

A comprehensive physiologically based pharmacokinetic (PBPK) model for nicotine in humans from using nicotine-containing products with different routes of exposure.

Physiologically based pharmacokinetic (PBPK) modeling can be a useful tool for characterizing nicotine pharmacokinetics (PK) from use of tobacco products. We expand a previously published PBPK model to simulate a nicotine PK profile, following single or multiple use of various tobacco products [cigarettes, smokeless tobacco, and electronic nicotine delivery systems, or a nicotine inhaler (NICOTROL)] The uptake route in the model was designed to allow for three uptake compartments: buccal cavity (BC), upper respiratory tract (URT) (conducting and transitional airways) and lower respiratory tract (alveolar region). Within each region, the model includes product-specific descriptions of the flux of nicotine into plasma, as well as the flux of nicotine from the BC and URT to the gastrointestinal tract. These descriptions are based on regional deposition and diffusion models of nicotine into plasma, which depends on the product type. Regional deposition flux combined with regional differences in physiological parameters (e.g., blood perfusion ratio and tissue thickness) play a key role in the product-specific PK profile of nicotine. The current model describes the slower flux of nicotine into plasma across the BC and URT, as well as the rapid flux known to occur in the alveolar region. Overall, the addition of the BC and respiratory tract compartments to the nicotine model provided simulation results that are comparable to the nicotine time-course plasma concentrations reported from clinical studies for the four product categories simulated.

Evaluation and prediction of pharmacokinetics of PFOA and PFOS in the monkey and human using a PBPK model.

Perfluoroalkyl acid carboxylates and sulfonates (PFAAs) have many consumer and industrial applications. The persistence and widespread distribution of these compounds in humans have brought them under intense scrutiny. Limited pharmacokinetic data is available in humans; however, human data exists for two communities with drinking water contaminated by PFAAs. Also, there is toxicological and pharmacokinetic data for monkeys, which can be quite useful for cross-species extrapolation to humans. The goal of this research was to develop a physiologically-based pharmacokinetic (PBPK) model for PFOA and PFOS for monkeys and then scale this model to humans in order to describe available human drinking water data. The monkey model simulations were consistent with available PK data for monkeys. The monkey model was then extrapolated to the human and then used to successfully simulate the data collected from residents of two communities exposed to PFOA in drinking water. Human PFOS data is minimal; however, using the half-life estimated from occupational exposure, our model exhibits reasonable agreement with the available human serum PFOS data. It is envisioned that our PBPK model will be useful in supporting human health risk assessments for PFOA and PFOS by aiding in understanding of human pharmacokinetics.

A Kinetic Analysis of DNA-Deoxy Guanine Adducts in the Nasal Epithelium Produced by Inhaled Formaldehyde in Rats-Assessing Contributions to Adduct Production From Both Endogenous and Exogenous Sources of Formaldehyde.

Although formaldehyde is a normal constituent of tissues, lifetime inhalation exposures at 6 h/day, 5 days/week at concentrations ≥6 ppm caused a nonlinear increase in nasal tumors in rats with incidence reaching close to 50% at 15 ppm. Studies with heavy isotope labeled [13CD2]-formaldehyde permit quantification of both the mass-labeled exogenous and endogenous DNA-formaldehyde reaction products. An existing pharmacokinetic model developed initially to describe 14C-DNA-protein crosslinks (DPX) provided a template for describing the time course of mass-labeled adducts. Published datasets included both DPX and N2-HO13CD2-dG adducts measured after a single 6-h exposure to 0.7, 2, 6, 9, 10, or 15 ppm formaldehyde, after multi-day exposures to 2 ppm for 6 h/day, 7 days/week with interim sacrifices up to 28 days, and after 28-day exposures for 6 h/day, 7 days/week to 0.3, 0.03, or 0.001 ppm. The existing kinetic model overpredicted endogenous adducts in the nasal epithelium after 1-day [13CD2]-formaldehyde exposure, requiring adjustment of parameters for rates of tissue metabolism and background formaldehyde. After refining tissue formaldehyde parameters, we fit the model to both forms of adducts by varying key parameters and optimizing against all 3 studies. Fitting to all these studies required 2 nonlinear pathways-one for high-exposure saturation of clearance in the nasal epithelial tissues and another for extracellular clearance that restricts uptake into the epithelial tissue for inhaled concentrations below 0.7 ppm. This refined pharmacokinetic model for endogenous and exogenous formaldehyde acetal adducts can assist in updating biologically based dose-response models for formaldehyde carcinogenicity.

Application of a combined aggregate exposure pathway and adverse outcome pathway (AEP-AOP) approach to inform a cumulative risk assessment: A case study with phthalates.

Advancements in measurement and modeling capabilities are providing unprecedented access to estimates of chemical exposure and bioactivity. With this influx of new data, there is a need for frameworks that help organize and disseminate information on chemical hazard and exposure in a manner that is accessible and transparent. A case study approach was used to demonstrate integration of the Adverse Outcome Pathway (AOP) and Aggregate Exposure Pathway (AEP) frameworks to support cumulative risk assessment of co-exposure to two phthalate esters that are ubiquitous in the environment and that are associated with disruption of male sexual development in the rat: di(2-ethylhexyl) phthalate (DEHP) and di-n-butyl phthalate (DnBP). A putative AOP was developed to guide selection of an in vitro assay for derivation of bioactivity values for DEHP and DnBP and their metabolites. AEPs for DEHP and DnBP were used to extract key exposure data as inputs for a physiologically based pharmacokinetic (PBPK) model to predict internal metabolite concentrations. These metabolite concentrations were then combined using in vitro-based relative potency factors for comparison with an internal dose metric, resulting in an estimated margin of safety of ~13,000. This case study provides an adaptable workflow for integrating exposure and toxicity data by coupling AEP and AOP frameworks and using in vitro and in silico methodologies for cumulative risk assessment.