Search

Fitted algebraic equations with probability distributions are presented to predict the maturation of selected physiological parameters that occur postnatally in preterm and term neonates through early infancy. The neonatal population, in particular preterm neonates, remains a challenge for the development of physiologically based pharmacokinetic (PBPK) models. The need to make available these building blocks for PBPK models for this age group is important for this field to advance. While the equations do not represent an exhaustive evaluation of anatomy, equations with 90% prediction intervals, where data permits, are presented for body weight, organ and tissue growth, blood flows, cardiac output, plasma albumin concentration, and glomerular filtration rate. Ontogeny measurements for blood flows to organs and tissues are limited. Longitudinal anthropomorphic measurements versus intrauterine or newborn measurements in preterm neonates would help reduce the uncertainty in the growth trajectories of neonates.

Arsenic is a ubiquitous contaminant, with typical human dietary intake below 1 μg/kg bw/d and extreme drinking water exposures up to ∼50 μg/kg bw/d. The formation and binding of trivalent metabolites are central to arsenic toxicity and strong human evidence suggests special concern for early life exposures in the etiology of adult diseases, especially cancer. This study measured the metabolism and disposition of arsenite in neonatal mice to understand the role of maturation in metabolic activation and detoxification of arsenic. Many age-related differences were observed after gavage administration of arsenite, with consistent evidence in blood and tissues for higher exposures to trivalent arsenic species in neonatal mice related to the immaturity of metabolic and/or excretory functions. The evidence for greater tissue binding of arsenic species in young mice is consistent with enhanced susceptibility to toxicity based on metabolic and toxicokinetic differences alone. Lactational transfer from arsenite-dosed dams to suckling mice was minimal, based on no dosing-related changes in the levels of arsenic species in pup blood or milk collected from the dams. Animal models evaluating whole-life exposure to inorganic arsenic must use direct dosing in early neonatal life to predict accurately potential toxicity from early life exposures in children.

Chronic exposure to inorganic arsenic (iAs), a contaminant of water and food supplies, is associated with many adverse health effects. A notable feature of iAs metabolism is sequential methylation reactions which produce mono- and di-methylated arsenicals that can contain arsenic in either the trivalent (III) or pentavalent (V) valence states. Because methylated arsenicals containing trivalent arsenic are more potent toxicants than their pentavalent counterparts, the ability to distinguish between the +3 and +5 valence states is a crucial property for physiologically based pharmacokinetic (PBPK) models of arsenicals to possess if they are to be of use in risk assessment. Unfortunately, current analytic techniques for quantifying arsenicals in tissues disrupt the valence state; hence, pharmacokinetic studies in animals, used for model calibration, only reliably provide data on the sum of the +3 and +5 valence forms of a given metabolite. In this paper we show how mathematical modeling can be used to overcome this obstacle and present a PBPK model for the dimethylated metabolite of iAs, which exists as either dimethylarsinous acid, (CH(3))(2)As(III)OH (abbreviated DMA(III)) or dimethylarsinic acid, (CH(3))(2)As(V)(O)OH (abbreviated DMA(V)). The model distinguishes these two forms and sets a lower bound on how much of an organ’s DMA burden is present in the more reactive and toxic trivalent valence state. We conjoin the PBPK model to a simple model for DMA(III)-induced oxidative stress in liver and use this extended model to predict cytotoxicity in liver in response to the high oral dose of DMA(V). The model incorporates mechanistic details derived from in vitro studies and is iteratively calibrated with lumped-valence-state PK data for intravenous or oral dosing with DMA(V). Model formulation leads us to predict that orally administered DMA(V) undergoes extensive reduction in the gastrointestinal (GI) tract to the more toxic trivalent DMA(III).

The estimated glomerular filtration rate (eGFR) equations based on serum creatinine (SCR) have been used for pediatric dose adjustment in drug labeling. This study evaluated the performance of those equations in estimating individual clearance of drugs that are predominantly eliminated by glomerular filtration, using clinical data from the renally eliminated drugs gadobutrol, gadoterate, amikacin, and vancomycin. The eGFR was compared with the observed drug clearance (CL) in 352 pediatric patients from birth to 12 years of age. Multiple eGFR equations overestimated the drug CL on average, including the original and bedside Schwartz equations, which showed an average eGFR/CL ratio between 1 and 3. Further analysis with bedside Schwartz equation showed a higher eGFR/CL ratio in the subjects with a lower SCR or CL. Supraphysiological eGFR as high as 380 mL/min/1.73 m2 was obtained using the bedside Schwartz equation for some of the subjects, most of whom are children < 2 years of age with SCR < 0.2 mg/dL. Excluding the subjects with supraphysiological eGFR from the analysis did not change the overall trend of overestimation. In conclusion, Schwartz equations led to an overestimation of drug clearance for the drugs evaluated. When greater precision is required in predicting eGFR for pediatric patients, such as in drug dosing, revised k constants for the Schwartz equation or new methods of glomerular filtration rate estimation may be necessary.

The objective of this study was to develop pediatric physiologically based pharmacokinetic (PBPK) models for pantoprazole and esomeprazole. Pediatric PBPK models were developed by Simcyp version 15 by incorporating cytochrome P450 (CYP)2C19 maturation and auto‐inhibition. The predicted‐to‐observed pantoprazole clearance (CL) ratio ranged from 0.96–1.35 in children 1–17 years of age and 0.43–0.70 in term infants. The predicted‐to‐observed esomeprazole CL ratio ranged from 1.08–1.50 for children 6–17 years of age, and 0.15–0.33 for infants. The prediction was markedly improved by assuming no auto‐inhibition of esomeprazole in infants in the PBPK model. Our results suggested that the CYP2C19 auto‐inhibition model was appropriate for esomeprazole in adults and older children but could not be directly extended to infants. A better understanding of the complex interplay of enzyme maturation, inhibition, and compensatory mechanisms for CYP2C19 is necessary for PBPK modeling in infants.

Arsenic is a ubiquitous contaminant, with typical dietary intake below 1 μg/kg bw/d and drinking water exposures up to 50 μg/kg bw/d. Arsenic exposures are associated with human diseases and doses of toxicological concern are similar to typical dietary intake. Metabolism of arsenite to dimethylarsinate (DMAV) by arsenite-3-methyltransferase (As3MT) promotes clearance, but also generates reactive trivalent intermediates that bind extensively to cellular thiols. This study measured pentavalent and trivalent arsenic species in blood and tissues after oral and intravenous administration of arsenite (50 μg/kg bw). After oral administration, the intestine and liver contained elevated levels of AsIII and MMAIII, relative to erythrocytes, lung, and kidney, suggesting incomplete conversion to DMA during first-pass metabolism. However, blood concentrations of the predominant species, DMA, were similar for oral and intravenous dosing. While all tissues examined contained DMAIII, muscle, brain, and plasma had undetectable levels of MMAIII. Tissue levels of arsenic species were similar following intravenous vs. oral administration, except lower in the intestine. The results confirm the role of metabolism in producing fluxes of putatively toxic trivalent arsenic intermediates. Tissue dosimetry suggests that the intestine, liver, lung, and kidney could be more susceptible to effects of bound arsenic, relative to muscle and brain.

Physiologically based pharmacokinetic (PBPK) models are developed from compound-independent information to describe important anatomical and physiological characteristics of an individual or population of interest. Modeling pediatric populations is challenging because of the rapid changes that occur during growth, particularly in the first few weeks and months after birth. Neonates who are born premature pose several unique challenges in PBPK model development. To provide appropriate descriptions for body weight (BW) and height (Ht) for age and appropriate incremental gains in PBPK models of the developing preterm and full term neonate, anthropometric measurements collected longitudinally from 1,063 preterm and 158 full term neonates were combined with 2,872 cross-sectional measurements obtained from the NHANES 2007-2010 survey. Age-specific polynomial growth equations for BW and Ht were created for male and female neonates with corresponding gestational birth ages of 25, 28, 31, 34, and 40 weeks. Model-predicted weights at birth were within 20% of published fetal/neonatal reference standards. In comparison to full term neonates, postnatal gains in BW and Ht were slower in preterm subgroups, particularly in those born at earlier gestational ages. Catch up growth for BW in neonates born at 25, 28, 31, and 34 weeks gestational age was complete by 13, 8, 6, and 2 months of life (males) and by 10, 6, 5, and 2 months of life (females), respectively. The polynomial growth equations reported in this paper represent extrauterine growth in full term and preterm neonates and differ from the intrauterine growth standards that were developed for the healthy unborn fetus.

Polyfluorinated compounds (PFCs) are authorized for use as greaseproofing agents in food contact paper. As C8-PFCs (8-carbons) are known to accumulate in tissues, shorter-chain C6-PFCs (6-carbons) have replaced C8-PFCs in many food contact applications. However, the potential of C6-PFCs for human biopersistence has not been fully evaluated. For the first time, we provide internal exposure estimates to key metabolites of 6:2 fluorotelomer alcohol (6:2 FTOH), a monomeric component of C6-PFCs, to extend our understanding of exposure beyond estimates of external exposure. Pharmacokinetic data from published rat and human studies on 6:2 FTOH were used to estimate clearance and area under the curve (AUC) for its metabolites: 5:3 fluorotelomer carboxylic acid (5:3 A), perfluorohexanoic acid (PFHxA) and perfluoroheptanoic acid (PFHpA). Internal exposure to 5:3 A was the highest of evaluated metabolites across species and it had the slowest clearance. Additionally, 5:3 A clearance decreased with increasing 6:2 FTOH exposure. Our analysis provides insight into association of increased internal 5:3 A exposure with high biopersistence potential of 6:2 FTOH. Our results identify 5:3 A as an important biomarker of internal 6:2 FTOH exposure for use in biomonitoring studies, and are potentially useful for toxicological assessment of chronic dietary 6:2 FTOH exposure.

Background: Per- and polyfluoroalkyl substances (PFAS) are considered chemicals of emerging concern, in part due to their environmental and biological persistence and the potential for widespread human exposure. In 2007, a PFAS manufacturer near Decatur, Alabama notified the United States Environmental Protection Agency (EPA) it had discharged PFAS into a wastewater treatment plant, resulting in environmental contamination and potential exposures to the local community. Objectives: To characterize PFAS exposure over time, the Agency for Toxic Substances and Disease Registry (ATSDR) collected blood and urine samples from local residents. Methods: Eight PFAS were measured in serum in 2010 (n=153). Eleven PFAS were measured in serum, and five PFAS were measured in urine (n=45) from some of the same residents in 2016. Serum concentrations were compared to nationally representative data and change in serum concentration over time was evaluated. Biological half-lives were estimated for perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid (PFOS), and perfluorohexane sulfonic acid (PFHxS) using a one-compartment pharmacokinetic model. Results: In 2010 and 2016, geometric mean PFOA and PFOS serum concentrations were elevated in participants compared to the general U.S. population. In 2016, the geometric mean PFHxS serum concentration was elevated compared to the general U.S. population. Geometric mean serum concentrations of PFOA, PFOS, and perfluorononanoic acid (PFNA) were significantly (p≤0.0001) lower (49%, 53%, and 58%, respectively) in 2016 compared to 2010. Half-lives for PFOA, PFOS, and PFHxS were estimated to be 3.9, 3.3, and 15.5years, respectively. Concentrations of PFOA in serum and urine were highly correlated (r=0.75) in males. Conclusions: Serum concentrations of some PFAS are decreasing in this residentially exposed community, but remain elevated compared to the U.S. general population. Keywords: PFAS, PFOA, PFOS, Biomonitoring, Half-life

Physiologically based pharmacokinetic (PBPK) modeling is a computational process that simulates the absorption, distribution, metabolism, and excretion of a substance in the body of an organism based on the interrelationships among key physiological, biochemical, and physicochemical factors using mathematical equations. This chapter provides a summary of the history and recent efforts of PBPK modeling in several applications, including (1) drug discovery and development, (2) toxicology and risk assessment, (3) veterinary pharmacology and animal-derived food safety assessment, and (4) nanomedicine and nanotoxicology. The history of PBPK modeling software programs, books, and the Society of Toxicology Biological Modeling Specialty Section is also briefly introduced. The authors’ perspective on the current and future applications of PBPK modeling in different fields, especially in toxicology and risk assessment, is presented.

This chapter provides the reader new to physiologically based pharmacokinetic (PBPK) modeling a perspective on the evolution of this computational tool for use in chemical toxicology and how to get started learning PBPK modeling with a hands-on exercise. A step-by-step process is presented to learn the components of a PBPK model that represent the body of a given species and how the chemical-specific information is used in a PBPK model. Rate equations for compartments that represent the body are discussed along with equations used to simulate chemical exposures via inhalation, intravenous administration, and oral bolus administration (laboratory animals). The computer software code is provided for a PBPK model written in Magnolia freeware for the chemical trichloroethylene (TCE). Rats are exposed to TCE via inhalation or intravenous administration. Step-by-step exercises are presented for the beginner to experience the wonderful world of simulation.

Accounting for metabolism in a physiologically based pharmacokinetic (PBPK) model is important from a pharmacokinetics perspective. Hepatic and extrahepatic enzymatic biotransformations occur across species and are usually dependent upon age and perhaps sex. Metabolism is either linked causally to toxicity or detoxification. In this chapter, brief descriptions of experimental methods are given for how quantitative metabolic studies are conducted and how the information is used in PBPK models. A challenging simulation exercise is provided for the chemical bisphenol A (BPA), which undergoes extensive hepatic and extrahepatic Phase II conjugation. Phase II conjugation is an inactivation step for BPA.

Oseltamivir is an antiviral drug approved to treat influenza in humans. Although the dosing regimen of this drug is well established for non-pregnant patients, it is not clear if the significant physiological alterations associated with pregnancy affect the pharmacokinetics of oseltamivir and, thus, warrant different dosing regimens to assure efficacy. In this study, we investigated the suitability of rhesus macaques as an animal model for studying oseltamivir pharmacokinetics during all trimesters of pregnancy in comparison to pre-pregnant conditions. Specifically, we compared the pharmacokinetics of oseltamivir and its pharmacologically active metabolite oseltamivir carboxylate in rhesus monkeys after intravenous and nasogastric administration of 2.5 mg oseltamivir phosphate/kg body weight given prior to and during the first, second, and third trimesters of pregnancy. Pregnancy had only a modest effect upon the pharmacokinetic parameters of oseltamivir and oseltamivir carboxylate. Monkeys treated intravenously in the third trimester had a reduction in Vd and CL, compared to non-pregnant monkeys. These changes did not occur in the other two trimesters. Pregnant monkeys treated intravenously had 20–25% decrease in AUC0-∞ of oseltamivir carboxylate and a corresponding increase in Vd and CL. Pregnant monkeys treated nasogastrically with oseltamivir phosphate demonstrated a pattern that recapitulated intravenous dosing. Taken together these data indicate that rhesus monkeys are an acceptable model for studying drug-pregnancy interactions.

Many alloys used in cardiovascular device applications contain high levels of nickel, which if released in sufficient quantities, can lead to adverse health effects. While nickel release from these devices is typically characterized through the use of in-vitro immersion tests, it is unclear if the rate at which nickel is released from a device during in-vitro testing is representative of the release rate following implantation in the body. To address this uncertainty, we have developed a novel biokinetic model that combines a traditional toxicokinetic compartment model with a physics-based model to estimate nickel release from an implanted device. This model links the rate of in-vitro nickel release from a cardiovascular device to serum nickel concentrations, an easily measured endpoint, to estimate the rate and extent of in-vivo nickel release from an implanted device. The model was initially parameterized using data in the literature on in-vitro nickel release from a nickel-containing alloy (nitinol) and baseline serum nickel levels in humans. The results of this first step were then used to validate specific components of the model. The remaining unknown quantities were fit using serum values reported in patients following implantation with nitinol atrial occluder devices. The model is not only consistent with levels of nickel in serum and urine of patients following treatment with the atrial occluders, but also the optimized parameters in the model were all physiologically plausible. The congruity of the model with available data suggests that it can provide a framework to interpret nickel biomonitoring data and use data from in-vitro nickel immersion tests to estimate in-vivo nickel release from implanted cardiovascular devices.

Modeling and simulation approaches are increasingly being utilized in pediatric drug development. Physiologically based pharmacokinetic (PBPK) modeling offers an enhanced ability to predict age-related changes in pharmacokinetics in the pediatric population. In the current study, adult PBPK models were developed for the renally excreted drugs linezolid and emtricitabine. PBPK models were then utilized to predict pharmacokinetics in pediatric patients for various age groups from the oldest to the youngest patients in a stepwise approach. Pharmacokinetic predictions for these two drugs in the pediatric population, including infants and neonates, were within a twofold range of clinical observations. Based on this study, linezolid and emtricitabine pediatric PBPK models incorporating the ontogeny in renal maturation describe the pharmacokinetic differences between adult and pediatric populations, even though the contribution of renal clearance to the total clearance of two drugs was very different (30 % for linezolid vs. 86 % for emtricitabine). These results suggest that PBPK modeling may provide one option to help predict the pharmacokinetics of renally excreted drugs in neonates and infants.

Physiologically-based pharmacokinetic (PBPK) modeling analysis does not stand on its own for regulatory purposes but is a robust tool to support drug/chemical safety assessment. While the development of PBPK models have grown steadily since their emergence, only a handful of models have been accepted to support regulatory purposes due to obstacles such as the lack of a standardized template for reporting PBPK analysis. Here, we expand the existing guidances designed for pharmaceutical applications by recommending additional elements that are relevant to environmental chemicals. This harmonized reporting template can be adopted and customized by public health agencies receiving PBPK model submission, and it can also serve as general guidance for submitting PBPK-related studies for publication in journals or other modeling sharing purposes. The current effort represents one of several ongoing collaborations among the PBPK modeling and risk assessment communities to promote, when appropriate, incorporating PBPK modeling to characterize the influence of pharmacokinetics on safety decisions made by regulatory agencies.

Renal elimination and the resulting clearance of perfluorooctanoic acid (PFOA) from the serum exhibit pronounced sex differences in the adult rat. The literature suggests that this is largely due to hormonally regulated expression of organic anion transporters (OATs) on the apical and basolateral membranes of the proximal tubule cells that facilitate excretion and reabsorption of PFOA from the filtrate into the blood. Previously developed PBPK models of PFOA exposure in the rat have not been parameterized to specifically account for transporter-mediated renal elimination. We developed a PBPK model for PFOA in the male and female rat to explore the role of Oat1, Oat3, and Oatp1a1 in sex-specific renal reabsorption and excretion of PFOA. Descriptions of the kinetic behavior of these transporters were extrapolated from in vitro studies and the model was used to simulate time-course serum, liver, and urine data for intravenous (IV) and oral exposures in both sexes. Model predicted concentrations of PFOA in the liver, serum, and urine showed good agreement with experimental data for both the male and female rat indicating that in vitro derived physiological descriptions of transporter-mediated renal reabsorption can successfully predict sex-dependent excretion of PFOA in the rat. This study supports the hypothesis that sex-specific serum half-lives for PFOA are largely driven by expression of transporters in the kidney and contributes to the development of PBPK modeling as a tool for evaluating the role of transporters in renal clearance.

2-Phenoxyethanol (PhE) has been shown to induce hepatotoxicity, renal toxicity, and hemolysis at dosages ≥ 400 mg/kg/day in subchronic and chronic studies in multiple species. To reduce uncertainty associated with interspecies extrapolations and to evaluate the margin of exposure (MOE) for use of PhE in cosmetics and baby products, a physiologically-based pharmacokinetic (PBPK) model of PhE and its metabolite 2-phenoxyacetic acid (PhAA) was developed. The PBPK model incorporated key kinetic processes describing the absorption, distribution, metabolism and excretion of PhE and PhAA following oral and dermal exposures. Simulations of repeat dose rat studies facilitated the selection of systemic AUC as the appropriate dose metric for evaluating internal exposures to PhE and PhAA in rats and humans. Use of the PBPK model resulted in refinement of the total default UF for extrapolation of the animal data to humans from 100 to 25. Based on very conservative assumptions for product composition and aggregate product use, model-predicted exposures to PhE and PhAA resulting from adult and infant exposures to cosmetic products are significantly below the internal dose of PhE observed at the NOAEL dose in rats. Calculated MOEs for all exposure scenarios were above the PBPK-refined UF of 25.

Extensive first-pass metabolism of ingested bisphenol A (BPA) in the gastro-intestinal tract and liver restricts blood concentrations of bioactive BPA to

Here we present data to evaluate potential absorption of Bisphenol A through non-metabolizing tissues of the upper digestive tract. Concurrent serum and urine concentrations of d6-BPA, and its glucuronide and sulfate conjugates, were measured over a 24h period in 10 adult male volunteers following ingestion of 30μg d6-BPA/kg body weight in soup. The pharmacokinetic behavior of BPA and its metabolites in this cohort (rapid absorption, complete elimination, evidence against sublingual absorption) was reported. This Data in Brief article contains the corresponding individual pharmacokinetic data, reports the demographics of the cohort and provides additional details related to the analytical methods employed and is related to [4].

The development and application of physiologically based pharmacokinetic (PBPK) models in chemical toxicology have grown steadily since their emergence in the 1980s. However, critical evaluation of PBPK models to support public health decision-making across federal agencies has thus far occurred for only a few environmental chemicals. In order to encourage decision-makers to embrace the critical role of PBPK modeling in risk assessment, several important challenges require immediate attention from the modeling community. The objective of this contemporary review is to highlight 3 of these challenges, including: (1) difficulties in recruiting peer reviewers with appropriate modeling expertise and experience; (2) lack of confidence in PBPK models for which no tissue/plasma concentration data exist for model evaluation; and (3) lack of transferability across modeling platforms. Several recommendations for addressing these 3 issues are provided to initiate dialog among members of the PBPK modeling community, as these issues must be overcome for the field of PBPK modeling to advance and for PBPK models to be more routinely applied in support of public health decision-making.

Early-life stages and maternal reproductive states are generally perceived as potentially sensitive to the disrupting effects of drugs or chemicals. Physiological models play an important role in neurodevelopmental toxicology safety assessments by predicting the dosimetry of chemicals or drugs during rapid maturation of biological systems and, in some cases, pharmacodynamic responses. Brief reviews of PBPK models are provided for the pesticides atrazine, chlorpyrifos, and deltamethrin; the essential-nutrient manganese; and the environmental contaminant perchlorate combined with the essential nutrient iodine.

Manufacturing of perfluorooctanoic acid (PFOA), a synthetic chemical with a long half-life in humans, peaked between 1970 and 2002, and has since diminished. In the United States, PFOA is detected in the blood of >99% of people tested, but serum concentrations have decreased since 1999. Much is known about exposure to PFOA in drinking water; however, the impact of non-drinking water PFOA exposure on serum PFOA concentrations is not well characterized. The objective of this research is to apply physiologically based pharmacokinetic (PBPK) modeling and Monte Carlo analysis to evaluate the impact of historic non-drinking water PFOA exposure on serum PFOA concentrations. In vitro to in vivo extrapolation was utilized to inform descriptions of PFOA transport in the kidney. Monte Carlo simulations were incorporated to evaluate factors that account for the large inter-individual variability of serum PFOA concentrations measured in individuals from North Alabama in 2010 and 2016, and the Mid-Ohio River Valley between 2005 and 2008. Predicted serum PFOA concentrations were within two-fold of experimental data. With incorporation of Monte Carlo simulations, the model successfully tracked the large variability of serum PFOA concentrations measured in populations from the Mid-Ohio River Valley. Simulation of exposure in a population of 45 adults from North Alabama successfully predicted 98% of individual serum PFOA concentrations measured in 2010 and 2016, respectively, when non-drinking water ingestion of PFOA exposure was included. Variation in serum PFOA concentrations may be due to inter-individual variability in the disposition of PFOA and potentially elevated historical non-drinking water exposures.

Many physiologically based pharmacokinetic (PBPK) models for environmental chemicals, drugs, and nanomaterials have been developed to aid risk and safety assessments using acslX. However, acslX has been rendered sunset since November 2015. Alternative modeling tools and tutorials are needed for future PBPK applications. This forum article aimed to: (1) demonstrate the performance of 4 PBPK modeling software packages (acslX, Berkeley Madonna, MATLAB, and R language) tested using 2 existing models (oxytetracycline and gold nanoparticles); (2) provide a tutorial of PBPK model code conversion from acslX to Berkeley Madonna, MATLAB, and R language; (3) discuss the advantages and disadvantages of each software package in the implementation of PBPK models in toxicology, and (4) share our perspective about future direction in this field. Simulation results of plasma/tissue concentrations/amounts of oxytetracycline and gold from different models were compared visually and statistically with linear regression analyses. Simulation results from the original models were correlated well with results from the recoded models, with time-concentration/amount curves nearly superimposable and determination coefficients of 0.86-1.00. Step-by-step explanations of the recoding of the models in different software programs are provided in the Supplementary Data. In summary, this article presents a tutorial of PBPK model code conversion for a small molecule and a nanoparticle among 4 software packages, and a performance comparison of these software packages in PBPK model implementation. This tutorial helps beginners learn PBPK modeling, provides suggestions for selecting a suitable tool for future projects, and may lead to the transition from acslX to alternative modeling tools.

Our previous report on pharmacokinetic (PK) evaluation of 6:2 fluorotelomer alcohol (6:2 FTOH) examined the biopersistence potential of its metabolites based on data published from single inhalation and occupational 6:2 FTOH exposure studies. We calculated internal exposure estimates of three key metabolites of 6:2 FTOH, of which 5:3 fluorotelomer carboxylic acid (5:3 acid) had the highest internal exposure and the slowest clearance. No oral repeated 6:2 FTOH exposure data were available at the time to fully characterize the biopersistence potential of the metabolite 5:3 acid. We recently received additional data on 6:2 FTOH and 5:3 acid, which included a 90-day toxicokinetic study report on repeated oral 6:2 FTOH exposure to rats. We reviewed the study and analyzed the reported 5:3 acid concentrations in plasma, liver, and fat using one-compartment PK modeling and calculated elimination rate constants (kel), elimination half-lives (t1/2) and times to steady state (tss) of 5:3 acid at three 6:2 FTOH doses. Our results showed that tss of 5:3 acid in plasma and evaluated tissues were approximately close to 1 year, such that the majority of highest values were observed at the lowest 6:2 FTOH dose, indicating its association with the biopersistence of 6:2 FTOH. The results of our PK analysis are the first to characterize biopersistence potential of the 5:3 acid after repeated oral exposure to the parent compound 6:2 FTOH based on steady state PK parameters, and therefore, may have an impact on future study designs when conducting toxicity assays for such compounds.

The widespread use and high abuse liability of tobacco products has received considerable public health attention, in particular for youth, who are vulnerable to nicotine addiction. In this study, adult and adolescent squirrel monkeys were used to evaluate age-related metabolism and pharmacokinetics of nicotine after intravenous administration. A physiologically-based pharmacokinetic (PBPK) model was created to characterize the pharmacokinetic behaviors of nicotine and its metabolites, cotinine, trans-3'-hydroxycotinine (3'-OH cotinine), and trans-3'-hydroxycotinine glucuronide (3'-OH cotinine glucuronide) for both adult and adolescent squirrel monkeys. The PBPK nicotine model was first calibrated for adult squirrel monkeys utilizing in vitro nicotine metabolic data, plasma concentration-time profiles and cumulative urinary excretion data for nicotine and metabolites. Further model refinement was conducted when the calibrated adult model was scaled to the adolescents, because adolescents appeared to clear nicotine and cotinine more rapidly relative to adults. More specifically, the resultant model parameters representing systemic clearance of nicotine and cotinine for adolescent monkeys were approximately two- to three-fold of the adult values on a per body weight basis. The nonhuman primate PBPK model in general captured experimental observations that were used for both model calibration and evaluation, with acceptable performance metrics for precision and bias. The model also identified differences in nicotine pharmacokinetics between adolescent and adult nonhuman primates which might also be present in humans.

The use of toxicokinetic (TK) data is becoming more prevalent in the evaluation of food ingredient safety as more TK information is being incorporated in safety data packages. Data demonstrating “1) the extent of absorption, 2) tissue distribution, 3) pathways and rates of metabolism, and 4) rate(s) of elimination” of food ingredients and their metabolites of intermediate and high toxicological potential may be useful for planning and designing toxicity studies, selecting doses for toxicity studies, addressing species differences, and understanding the potential modes of action to evaluate their safety. TK data reported in the literature or generated from mechanistic TK studies can be analyzed using mathematical methods, including compartment and noncompartment TK methods, whose predictions can enhance interpretation of observed effects. Because of recent advancements, several approaches have been developed to improve sensitivity of analyses of available TK data and reduce uncertainty for evaluating safety of food ingredients. An example of advanced TK methods is physiologically-based TK (PBTK) modeling that incorporates physiological/biochemical parameters into a TK framework to predict internal exposure. In this review, we discuss the utility of some TK methods and explore their relevance and potential value for food ingredient safety evaluation. We also describe the strengths and limitations of these TK methods and discuss current challenges and opportunities for expanding their application for evaluating safety of food ingredients. This review represents a state of science report, and not a guidance document, on the utility and relevance of TK methods for the safety evaluation of food ingredients.

A probabilistic-based seven-compartment fit for purpose PBPK model for the preterm and term neonates was constructed and used with the combination antimicrobial drugs piperacillin (PIP) and tazobactam (TAZ) as a case study. Anthropomorphic maturation equations with distributions were used to account for variability in growth patterns for neonates born premature (gestation age ≥ 25 weeks). Renal excretion, the primary pathway for elimination of these drugs from the body, was described by passive diffusion (glomerular filtration, GFR) and active tubular secretion. Tubular secretion rates were calculated using in-vitro to in-vivo extrapolation methods for OAT1 and OAT3 protein transporters. Competitive inhibition of each anion (PIP and TAZ) for tubular secretion by the OAT transporters suggested that PIP decreased the systemic clearance rate of TAZ. After accounting for variability in maturation of growth (organs, kidney function and blood flows) model predictions did not account for all the observed variability in the pharmacokinetic plasma data for PIP and TAZ. For the 31 preterm and term neonates, the PBPK model for PIP and TAZ, with no adjustments in model parameters, was judged successful for 48% and 43% of the preterm neonates and 25% of the term neonates. Further refinements in model parameters to increase model performance include experiments to determine plasma protein binding in preterm and term neonates and solubility or partitioning of PIP and TAZ into tissues. A better understanding of the role of illness on modifying pharmacokinetic behavior of PIP and TAZ is needed.

Mammalian estrogen receptors modulate many physiological processes. Chemicals with structural features similar to estrogens can interact with estrogen receptors to produce biological effects similar to those caused by endogenous estrogens in the body. Bisphenol A (BPA) is a structural analogue of estrogen that binds to estrogen receptors. Exposure to BPA in humans is virtually ubiquitous in industrialized societies, but BPA is rapidly detoxified by metabolism and does not accumulate in the body. Whether or not serum concentrations of BPA in humans are sufficiently high to disrupt normal estrogen-related biology is the subject of intense political and scientific debate. Here we show a convergence of robust methods for measuring or calculating serum concentrations of BPA in humans from 93 published studies of more than 30,000 individuals in 19 countries across all life stages. Typical serum BPA concentrations are orders of magnitude lower than levels measurable by modern analytical methods and below concentrations required to occupy more than 0.0009% of Type II Estrogen Binding Sites, GPR30, ERα or ERβ receptors. Occupancies would be higher, but ≤0.04%, for the highest affinity receptor, ERRγ. Our results show limited or no potential for estrogenicity in humans, and question reports of measurable BPA in human serum.

Bisphenol A (BPA) has received considerable attention throughout the last decade due to its widespread use in consumer products. For the first time a physiologically based pharmacokinetic (PBPK) model was developed in neonatal and adult rats to quantitatively evaluate age-dependent pharmacokinetics of BPA and its phase II metabolites. The PBPK model was calibrated in adult rats using studies on BPA metabolism and excretion in the liver and gastrointestinal tract, and pharmacokinetic data with BPA in adult rats. For immature rats the hepatic and gastrointestinal metabolism of BPA was inferred from studies on the maturation of phase II enzymes coupled with serum time course data in pups. The calibrated model predicted the measured serum concentrations of BPA and BPA conjugates after administration of 100μg/kg of d6-BPA in adult rats (oral gavage and intravenous administration) and postnatal days 3, 10, and 21 pups (oral gavage). The observed age-dependent BPA serum concentrations were partially attributed to the immature metabolic capacity of pups. A comparison of the dosimetry of BPA across immature rats and monkeys suggests that dose adjustments would be necessary to extrapolate toxicity studies from neonatal rats to infant humans.

Bromate (BrO3 −), a by-product of ozonation of drinking water, induces nephrotoxicity in male rats at much lower doses than in female rats. This difference appears to be related to the development of α-2u-globulin nephropathy in males. To determine sex-dependent changes in mRNA and protein expression in the renal cortex attributable to α-2u-globulin nephropathy, we performed microarray and immunohistochemical analyses in proximal renal tubules of male and female F344 rats treated with KBrO3 for 28 days. Particular attention was paid to molecular biomarkers of renal tubular injury. Microarray analysis of male and female rats treated with BrO3 − at low doses (125 mg/L KBrO3) displayed marked sex-dependent changes in renal gene expression. The greatest differences were seen in genes encoding for cellular differentiation, apoptosis, ion transport, and cell proliferation. Differences by sex were especially prominent for the cell cycle checkpoint gene p21, the renal injury protein Kim-1, and the kidney injury and cancer biomarker protein osteopontin. Dose-related nephrotoxicity, assessed by hematoxylin and eosin staining, was greater in males compared to female rats, as was cellular proliferation, assessed by bromodeoxyuridine staining. The fraction of proximal renal cells with elevated 8-oxodeoxyguanosine (8-OH-dG) was only increased at the high dose and did not differ by sex. Dose-dependent increases in the expression of osteopontin were detected immunohistochemically only in male rats and were localized in proximal tubule cells. Similarly, BrO3 − treatment increased clusterin and Kim-1 staining in the proximal tubules; however, staining for these proteins did not differ appreciably between males and females. These data demonstrate both qualitative and quantitative differences in the response of male versus female kidneys to BrO3 −-treatment. These sex-dependent effects likely contribute to renal carcinogenesis of BrO3 − in the male rat.

A biologically based dose-response model (BBDR) for the hypothalamic pituitary thyroid (HPT) axis was developed in the near-term pregnant mother and fetus. This model was calibrated to predict serum levels of iodide, total thyroxine (T4), free thyroxine (fT4), and total triiodothyronine (T3) in the mother and fetus for a range of dietary iodide intake. The model was extended to describe perchlorate, an environmental and food contaminant, that competes with the sodium iodide symporter protein for thyroidal uptake of iodide. Using this mode-of-action framework, simulations were performed to determine the daily ingestion rates of perchlorate that would be associated with hypothyroxinemia or onset of hypothyroidism for varying iodide intake. Model simulations suggested that a maternal iodide intake of 75 to 250 µg/day and an environmentally relevant exposure of perchlorate (~0.1 µg/kg/day) did not result in hypothyroxinemia or hypothyroidism. For a daily iodide-sufficient intake of 200 µg/day, the dose of perchlorate required to reduce maternal fT4 levels to a hypothyroxinemic state was estimated at 32.2 µg/kg/day. As iodide intake was lowered to 75 µg/day, the model simulated daily perchlorate dose required to cause hypothyroxinemia was reduced by eightfold. Similarly, the perchlorate intake rates associated with the onset of subclinical hypothyroidism ranged from 54.8 to 21.5 µg/kg/day for daily iodide intake of 250-75 µg/day. This BBDR-HPT axis model for pregnancy provides an example of a novel public health assessment tool that may be expanded to address other endocrine-active chemicals found in food and the environment.

A previously presented physiologically-based pharmacokinetic model for immediate release (IR) methylphenidate (MPH) was extended to characterize the pharmacokinetic behaviors of oral extended release (ER) MPH formulations in adults for the first time. Information on the anatomy and physiology of the gastrointestinal (GI) tract, together with the biopharmaceutical properties of MPH, was integrated into the original model, with model parameters representing hepatic metabolism and intestinal non-specific loss recalibrated against in vitro and in vivo kinetic data sets with IR MPH. A Weibull function was implemented to describe the dissolution of different ER formulations. A variety of mathematical functions can be utilized to account for the engineered release/dissolution technologies to achieve better model performance. The physiological absorption model tracked well the plasma concentration profiles in adults receiving a multilayer-release MPH formulation or Metadate CD, while some degree of discrepancy was observed between predicted and observed plasma concentration profiles for Ritalin LA and Medikinet Retard. A local sensitivity analysis demonstrated that model parameters associated with the GI tract significantly influenced model predicted plasma MPH concentrations, albeit to varying degrees, suggesting the importance of better understanding the GI tract physiology, along with the intestinal non-specific loss of MPH. The model provides a quantitative tool to predict the biphasic plasma time course data for ER MPH, helping elucidate factors responsible for the diverse plasma MPH concentration profiles following oral dosing of different ER formulations.

Bromate (BrO(3)(-)) is a ubiquitous by-product of using ozone to disinfect water containing bromide (Br(-)). The reactivity of BrO(3)(-) with biological reductants suggests that its systemic absorption and distribution to target tissues may display non-linear behavior as doses increase. The intent of this study is to determine the extent to which BrO(3)(-) is systemically bioavailable via oral exposure and broadly identify its pathways of degradation. In vitro experiments of BrO(3)(-) degradation in rat blood indicate a rapid initial degradation immediately upon addition that is >98% complete at concentrations up to 66μM in blood. As initial concentrations are increased, progressively lower fractions are lost prior to the first measurement. Secondary to this initial loss, a slower and predictable first order degradation rate was observed (10%/min). Losses during both phases were accompanied by increases in Br(-) concentrations indicating that the loss of BrO(3)(-) was due to its reduction. In vivo experiments were conducted using doses of BrO(3)(-) ranging from 0.077 to 15.3mg/kg, administered intravenously (IV) or orally (gavage) to female F344 rats. The variable nature and uncertain source of background concentrations of BrO(3)(-) limited derivation of terminal half-lives, but the initial half-life was approximately 10min for all dose groups. The area under the curve (AUC) and peak concentrations (C(t=5')) were linearly related to IV dose up to 0.77mg/kg; however, disproportionate increases in the AUC and C(t=5') and a large decrease in the volume of distribution was observed when IV doses of 1.9 and 3.8mg/kg were administered. The average terminal half-life of BrO(3)(-) from oral administration was 37min, but this was influenced by background levels of BrO(3)(-) at lower doses. With oral doses, the AUC and C(max) increased linearly with dose up to 15.3mgBrO(3)(-)/kg. BrO(3)(-) appeared to be 19-25% bioavailable without an obvious dose-dependency between 0.077 and 1.9mg/kg. The urinary elimination of BrO(3)(-) and Br(-) was measured from female F344 rats for four days following administration of single doses of 8.1mgKBrO(3)/kg and for 15 days after a single dose of 5.0mgKBr/kg. BrO(3)(-) elimination was detected over the first 12h, but Br(-) elimination from BrO(3)(-) over the first 48h was 18% lower than expected based on that eliminated from an equimolar dose of Br(-) (15.5±1.6 vs. 18.8±1.2μmol/kg, respectively). The cumulative excretion of Br(-) from KBr vs. KBrO(3) was equivalent 72h after administration. The recovery of unchanged administered BrO(3)(-) in the urine ranged between 6.0 and 11.3% (creatinine corrected) on the 27th day of treatment with concentrations of KBrO(3) of 15, 60, and 400mg/L of drinking water. The recovery of total urinary bromine as Br(-)+BrO(3)(-) ranged between 61 and 88%. An increase in the fraction of the daily BrO(3)(-) dose recovered in the urine was observed at the high dose to both sexes. The deficit in total bromine recovery raises the possibility that some brominated biochemicals may be produced in vivo and more slowly metabolized and eliminated. This was supported by measurements of dose-dependent increases of total organic bromine (TOBr) that was eliminated in the urine. The role these organic by-products play in BrO(3)(-)-induced cancer remains to be established.

Partition coefficients (PCs) are used in physiologically based pharmacokinetic (PBPK) models to estimate the free concentration of a chemical in specific blood or organs. Biological PC(tissue:blood) (tissue to blood) values were determined for a series of nonvolatile herbicides, insecticides, and fungicides in liver, brain, skin, fat, kidneys, and muscle of male Sprague-Dawley rats using two different analytical methods. The free phase concentration (in phosphate-buffered saline) of a given chemical was measured in the presence and absence of tissue (including blood) and used to calculate the PC, defined as the ratio of the concentration of the chemical in saline to the concentration in the tissue. PCs were determined for 13 compounds with aqueous solubility ranging from 20 to 4100 mg/L, molecular weights from 187.3 to 342.2 g/mol, and log K (ow) values from -0.18 to 3.9. An ultrafiltration high-performance liquid chromatography (HPLC) method was implemented for compounds with log K (ow) near 0.1 or less and a negligible depletion solid-phase microextraction (nd-SPME) method for compounds with higher log K (ow). PC(tissue:saline) coefficients of variation were 0.13 (n = 3 compounds) on average for the HPLC method and 0.29 (n = 10 compounds) for the nd-SPME method. Presented here is one of the most comprehensive data sets of biological partition coefficients for herbicides, insecticides, and fungicides.

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.

Toxicity and exposure evaluations remain the two of the key components of human health assessment. While improvement in exposure assessment relies on a better understanding of human behavior patterns, toxicity assessment still relies to a great extent on animal toxicity testing and human epidemiological studies. Recent advances in computer modeling of the dose-response relationship and distribution of xenobiotics in humans to important target tissues have advanced our abilities to assess toxicity. In particular, physiologically based pharmacokinetic (PBPK) models are among the tools than can enhance toxicity assessment accuracy. Many PBPK models are available to the health assessor, but most are so difficult to use that health assessors rarely use them. To encourage their use these models need to have transparent and user-friendly formats. To this end the Agency for Toxic Substances and Disease Registry (ATSDR) is using translational research to increase PBPK model accessibility, understandability, and use in the site-specific health assessment arena. The agency has initiated development of a human PBPK tool-kit for certain high priority pollutants. The tool kit comprises a series of suitable models. The models are recoded in a single computer simulation language and evaluated for use by health assessors. While not necessarily being state-of-the-art code for each chemical, the models will be sufficiently accurate to use for screening purposes. This article presents a generic, seven-compartment PBPK model for six priority volatile organic compounds (VOCs): benzene (BEN), carbon tetrachloride (CCl(4)), dichloromethane (DCM), perchloroethylene (PCE), trichloroethylene (TCE), and vinyl chloride (VC). Limited comparisons of the generic and original model predictions to published kinetic data were conducted. A goodness of fit was determined by calculating the means of the sum of the squared differences (MSSDs) for simulation vs. experimental kinetic data using the generic and original models. Using simplified solvent exposure assumptions for oral ingestion and inhalation, steady-state blood concentrations of each solvent were simulated for exposures equivalent to the ATSDR Minimal Risk Levels (MRLs). The predicted blood levels were then compared to those reported in the National Health and Nutrition Examination Survey (NHANES). With the notable exception of BEN, simulations of combined oral and inhalation MRLs using our generic VOC model yielded blood concentrations well above those reported for the 95th percentile blood concentrations for the U.S. populations, suggesting no health concerns. When the PBPK tool kit is fully developed, risk assessors will have a readily accessible tool for evaluating human exposure to a variety of environmental pollutants.

A physiologically based pharmacokinetic (PBPK) model was developed for bisphenol A (BPA) in adult rhesus monkeys using intravenous (iv) and oral bolus doses of 100 μg d6-BPA/kg (Doerge et al., 2010). This calibrated PBPK adult monkey model for BPA was then evaluated against published monkey kinetic studies with BPA. Using two versions of the adult monkey model based on monkey BPA kinetic data from Doerge et al. (2010) and Taylor et al. (2011), the aglycone BPA pharmacokinetics were simulated for human oral ingestion of 5 mg d16-BPA per person (Völkel et al., 2002). Völkel et al. were unable to detect the aglycone BPA in plasma, but were able to detect BPA metabolites. These human model predictions of the aglycone BPA in plasma were then compared to previously published PBPK model predictions obtained by simulating the Völkel et al. kinetic study. Our BPA human model, using two parameter sets reflecting two adult monkey studies, both predicted lower aglycone levels in human serum than the previous human BPA PBPK model predictions. BPA was metabolized at all ages of monkey (PND 5 to adult) by the gut wall and liver. However, the hepatic metabolism of BPA and systemic clearance of its phase II metabolites appear to be slower in younger monkeys than adults. The use of the current non-human primate BPA model parameters provides more confidence in predicting the aglycone BPA in serum levels in humans after oral ingestion of BPA.

A physiologically based pharmacokinetic (PBPK) model was developed for the conazole fungicide triadimefon and its primary metabolite, triadimenol. Rat tissue:blood partition coefficients and metabolic constants were measured in vitro for both compounds. Pharmacokinetic data for parent and metabolite were collected from several tissues after intravenous administration of triadimefon to male Sprague-Dawley rats. The model adequately simulated peak blood and tissue concentrations but predicted more rapid clearance of both triadimefon and triadimenol from blood and tissues. Reverse metabolism of triadimenol to triadimefon in the liver was explored as a possible explanation of this slow clearance, with significant improvement in model prediction. The amended model was extrapolated to humans using in vitro metabolic constants measured in human hepatic microsomes. Human equivalent doses (HEDs) were calculated for a rat no observable adverse effect level (NOAEL) dose of 3.4mg/kg/day using area under the concentration curve (AUC) in brain and blood for triadimefon and triadimenol as dosimetrics. All dosimetric-based HEDs were 25-30 fold above the human oral reference dose of 0.03mg triadimefon/kg/day, but did not account for intra-human variability or pharmacodynamic differences. Ultimately, derivations of this model will be able to better predict the exposure profile of these and other conazole fungicides in humans.

Perchlorate can perturb thyroid hormone (TH) homeostasis by competitive inhibition of iodide uptake by the thyroid gland. Until recently, the effects of perchlorate on TH homeostasis were examined by measuring serum concentrations of THs by immunoassay (IA) methods. IA methods are sensitive, but for TH analysis they are compromised by lack of adequate specificity. In this study, we determined the concentrations of six THs: L-thyroxine (T₄), 3,3',5-triiodo-L-thyronine (T₃), 3,3',5'-triiodo-L-thyronine (rT₃), 3,5-diiodo-L-thyronine, 3,3'-diiodo-L-thyronine, and 3-iodo-L-thyronine in the serum of rats administered perchlorate by isotope (¹³C₆-T₄)-dilution liquid chromatography-tandem mass spectrometry. The method recoveries for THs spiked into a serum matrix were between 97.0% and 115%, with a coefficient of variation of 2.1% to 9.4%. Rats were placed on an iodide-deficient or iodide-sufficient diet for 2.5 months, and for the last 2 weeks of that period they were provided drinking water either without or with perchlorate (10 mg/kg body weight/day). No significant differences in serum concentrations of T₃ and T₄ were observed between rats given iodide-deficient and iodide-sufficient diets for 2 or 2.5 months. After 24 h of perchlorate exposure, significantly lower concentrations of T₃ and T₄ were found in the serum of rats administered the iodide-deficient diet but not in rats administered the iodide-sufficient diet. However, after 2 weeks of perchlorate exposure, TH levels in rats fed the iodide-sufficient diet were also significantly lower than those in control rats. Our results suggest that perchlorate affects TH homeostasis and that such effects are more pronounced under iodide-deficient nutrition.

Jet Fuel 8 (JP-8) is a major fuel source used by US and NATO military. JP-8 is a complex mixture of aliphatic and aromatic isomers of hydrocarbons. Tissue/blood partition coefficient (PC) values are chemical-specific parameters used in modeling the kinetic behavior of chemicals. The partition coefficient values for n-alkanes tend to increase with the increasing carbon number, but less is known about the trend for isomers of n-alkanes. PC values were obtained for the n-alkane nonane (C9) and five of its isomers, namely 3-methyloctane, 4-ethylheptane, 2,3-dimethylheptane, 2,2,4-trimethylhexane, 2,2,4,4-tetramethylpentane. The blood:air and tissue:air PC values correlated with the published log octanol/water (O:W) PC values for n-nonane and its isomers. Experimentally determined blood:air and tissue:air PC values for n-nonane with the largest O:W value were greatest and smallest for the isomer 2,2,4,4-tetramethylpentane with the lowest O:W value. As expected the fat tissue had the highest PC values and muscle the lowest for n-nonane and its isomers. For each tissue, a linear relationship was observed between the tissue/blood PC values for the isomers of n-nonane and n-nonane. This suggests that tissue/blood PC values for all isomers of an alkane could be estimated using data collected from only a sub-set of alkanes of equal carbon number. These reported tissue/blood PC values will support the development of a jet fuel physiologically-based pharmacokinetic (PBPK) model.

Bisphenol A (BPA) is a high-production volume industrial chemical used in the manufacture of polycarbonate plastic products and epoxy resin-based food can liners. The presence of BPA in urine of >90% of Americans aged 6-60 is controversial because of the potential for endocrine disruption, particularly during perinatal development, as suggested by in vitro, experimental animal, and epidemiological studies. The current study used LC/MS/MS to measure serum pharmacokinetics of aglycone (active) and conjugated (inactive) BPA in adult and neonatal rhesus monkeys by oral (PND 5, 35, 70) and intravenous injection (PND 77) routes using d6-BPA to avoid sample contamination. The concentration-time profiles observed in adult monkeys following oral administration of 100 μg/kg bw were remarkably similar to those previously reported in human volunteers given a similar dose; moreover, minimal pharmacokinetic differences were observed between neonatal and adult monkeys for the receptor-active aglycone form of BPA. Circulating concentrations of BPA aglycone were quite low following oral administration (< 1% of total), which reflects the redundancy of active UDP-glucuronosyl transferase isoforms in both gut and liver. No age-related changes were seen in internal exposure metrics for aglycone BPA in monkeys, a result clearly different from developing rats where significant inverse age-related changes, based on immaturity of Phase II metabolism and renal excretion, were recently reported. These observations imply that any toxicological effect observed in rats from early postnatal exposures to BPA could over-predict those possible in primates of the same age, based on significantly higher internal exposures and overall immaturity at birth.

Bisphenol A (BPA) is an important industrial chemical used in the manufacture of polycarbonate plastic products and epoxy resin-based food can liners. The presence of BPA in urine of >90% of Americans aged 6-60 suggests ubiquitous and frequent exposure. The current study used LC/MS/MS to measure serum pharmacokinetics of aglycone (active) and conjugated (inactive) BPA in adult and neonatal Sprague-Dawley rats by oral and injection routes. Deuterated BPA was used to avoid issues of background contamination. Linear pharmacokinetics were observed in adult rats treated orally in the range of 0-200 microg/kg bw. Evidence for enterohepatic recirculation of conjugated, but not aglycone, BPA was observed in adult rats. Significant inverse relationships were observed between postnatal age and measures of internal exposures to aglycone BPA and its elimination. In neonatal rats treated orally, internal exposures to aglycone BPA were substantially lower than from subcutaneous injection. The results reinforce the critical role for first-pass Phase II metabolism of BPA in gut and liver after oral exposure that attenuates internal exposure to the aglycone form in rats of all ages. The internal exposures to aglycone BPA observed in adult and neonatal rats following a single oral dose of 100 microg/kg bw are inconsistent with effects mediated by classical estrogen receptors based on binding affinities. However, an impact on alternative estrogen signaling pathways that have higher receptor affinity cannot be excluded in neonatal rats. These findings emphasize the importance of matching aglycone BPA internal dosimetry with receptor affinities in experimental animal studies reporting toxicity.