Physiologically Based Pharmacokinetic Model for Developmental Exposures to TCDD in the Rat

A Physiologically Based Pharmacokinetic (PBPK) Modeling Framework for Mixtures of Dioxin-like Compounds

Humans are exposed to persistent organic pollutants, such as dioxin-like compounds (DLCs), as mixtures. Understanding and predicting the toxicokinetics and thus internal burden of major constituents of a DLC mixture is important for assessing their contributions to health risks. PBPK models, including dioxin models, traditionally focus on one or a small number of compounds; developing new or extending existing models for mixtures often requires tedious, error-prone coding work. This lack of efficiency to scale up for multi-compound exposures is a major technical barrier toward large-scale mixture PBPK simulations. Congeners in the DLC family, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), share similar albeit quantitatively different toxicokinetic and toxicodynamic properties. Taking advantage of these similarities, here we reported the development of a human PBPK modeling framework for DLC mixtures that can flexibly accommodate an arbitrary number of congeners. Adapted from existing TCDD models, our mixture model contains the blood and three diffusion-limited compartments-liver, fat, and rest of the body. Depending on the number of congeners in a mixture, varying-length vectors of ordinary differential equations (ODEs) are automatically generated to track the tissue concentrations of the congeners. Shared ODEs are used to account for common variables, including the aryl hydrocarbon receptor (AHR) and CYP1A2, to which the congeners compete for binding. Binary and multi-congener mixture simulations showed that the AHR-mediated cross-induction of CYP1A2 accelerates the sequestration and metabolism of DLC congeners, resulting in consistently lower tissue burdens than in single exposure, except for the liver. Using dietary intake data to simulate lifetime exposures to DLC mixtures, the model demonstrated that the relative contributions of individual congeners to blood or tissue toxic equivalency (TEQ) values are markedly different than those to intake TEQ. In summary, we developed a mixture PBPK modeling framework for DLCs that may be utilized upon further improvement as a quantitative tool to estimate tissue dosimetry and health risks of DLC mixtures.

Evaluation of Placentation and the Role of the Aryl Hydrocarbon Receptor Pathway in a Rat Model of Dioxin Exposure

BACKGROUND

Our environment is replete with chemicals that can affect embryonic and extraembryonic development. Dioxins, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), are compounds affecting development through the aryl hydrocarbon receptor (AHR).

OBJECTIVES

The purpose of this investigation was to examine the effects of TCDD exposure on pregnancy and placentation and to evaluate roles for AHR and cytochrome P450 1A1 (CYP1A1) in TCDD action.

METHODS

Actions of TCDD were examined in wild-type and genome-edited rat models. Placenta phenotyping was assessed using morphological, biochemical, and molecular analyses.

RESULTS

TCDD exposures were shown to result in placental adaptations and at higher doses, pregnancy termination. Deep intrauterine endovascular trophoblast cell invasion was a prominent placentation site adaptation to TCDD. TCDD-mediated placental adaptations were dependent upon maternal AHR signaling but not upon placental or fetal AHR signaling nor the presence of a prominent AHR target, CYP1A1. At the placentation site, TCDD activated AHR signaling within endothelial cells but not trophoblast cells. Immune and trophoblast cell behaviors at the uterine-placental interface were guided by the actions of TCDD on endothelial cells.

DISCUSSION

We identified an AHR regulatory pathway in rats activated by dioxin affecting uterine and trophoblast cell dynamics and the formation of the hemochorial placenta. https://doi.org/10.1289/EHP9256.

A PBPK model describing the pharmacokinetics of γ-HBCD exposure in mice

The brominated flame retardant, hexabromocyclododecane (HBCD), is added-but not bound-to consumer products and is eventually found in the environment and human tissues. Commercial-grade HBCD mixtures contain three major stereoisomers, alpha (α), beta (β), and gamma (γ), that are typically at a ratio of 12%:6%:82%, respectively. Although HBCD is widely used, the toxicological effects from its exposure in humans are not clearly understood. Using a physiologically based pharmacokinetic (PBPK) model could help improve our understanding of the toxicity of HBCD. The aim of this work was to develop a PBPK model, consisting of five permeability limited compartments (i.e., brain, liver, adipose tissue, blood, and rest of the body), to evaluate the pharmacokinetics of γ-HBCD in C57BL/6 mice. Physiological parameters related to body size, organ weights, and blood flow were taken from the literature. All partition coefficients were calculated based on the log Kow. The elimination in urine and feces was optimized to reflect the percent dose eliminated, as published in the literature. Compared with data from the literature for brain, liver, blood, and adipose tissue, the model simulations accurately described the mouse data set within 1.5-fold of the data points. Also, two examples showing the utility of the PBPK model supplement the information regarding the internal dose that caused the health effects observed during these studies. Although this version of the PBPK model expressly describes γ-HBCD, more efforts are needed to clarify and improve the model to discriminate between the α, β, and γ stereoisomers.

Maternal-Fetal Pharmacology of Drugs: A Review of Current Status of the Application of Physiologically Based Pharmacokinetic Models

Pregnancy and the postpartum period are associated with several physiological changes that can alter the pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. For certain drugs, dosing changes may be required during pregnancy and postpartum to achieve drug exposures comparable to what is observed in non-pregnant subjects. There is very limited data on fetal exposure of drugs during pregnancy, and neonatal exposure through transfer of drugs via human milk during breastfeeding. Very few systematic clinical pharmacology studies have been conducted in pregnant and postpartum women due to ethical issues, concern for the fetus safety as well as potential legal ramifications. Over the past several years, there has been an increase in the application of modeling and simulation approaches such as population PK (PopPK) and physiologically based PK (PBPK) modeling to provide guidance on drug dosing in those special patient populations. Population PK models rely on measured PK data, whereas physiologically based PK models incorporate physiological, preclinical, and clinical data into the model to predict drug exposure during pregnancy. These modeling strategies offer a promising approach to identify the drugs with PK changes during pregnancy to guide dose optimization in pregnancy, when there is lack of clinical data. PBPK modeling is also utilized to predict the fetal exposure of drugs and drug transfer via human milk following maternal exposure. This review focuses on the current status of the application of PBPK modeling to predict maternal and fetal exposure of drugs and thereby guide drug therapy during pregnancy.

Physiologically based toxicokinetic models and in silico predicted partition coefficients to estimate tetrachlorodibenzo-p-dioxin transfer from feed into growing pigs

Tetrachlorodibenzo-p-dioxin (TCDD) is a ubiquitous, toxic, persistent and bioaccumulative organic pollutant. TCDD can potentially enter the food chain through contaminated food of animal origin as a consequence of feed contamination. Prediction of the TCDD transfer from feed into animal products is thus important for human health risk assessment. Here, we develop several physiologically based toxicokinetic (PBTK) models of TCDD transfer from contaminated feed into growing pigs (Sus scrofa) exposed to doses ranging from 24.52 to 3269.25 ng of TCDD. We test the consequences of explicit dose-dependent absorption (DDA) versus the indirect effects of a self-induced liver metabolism (SIM). The DDA and SIM models showed similar fit to experimental data, although currently it is not possible to unequivocally make statement on a mechanistic preference. The performance of both toxicokinetic models was successfully evaluated using the 1999 Belgian case of contaminated fats for feeding. In combination with toxicokinetic models of other dioxin congeners, they can be used to formulate maximum allowance levels of dioxins in feedstuffs for pigs. Additionally, the implementation of in silico-predicted partition coefficients was explored as a useful alternative to predict TCDD tissue distribution in low-dose scenarios without recurring to animal experiments.

Application of qualitative and quantitative uncertainty assessment tools in developing ranges of plausible toxicity values for 2,3,7,8-tetrachlorodibenzo-p-dioxin

Increasing interest in characterizing risk assessment uncertainty is highlighted by recent recommendations from the National Academy of Sciences. In this paper we demonstrate the utility of applying qualitative and quantitative methods for assessing uncertainty to enhance risk-based decision-making for 2,3,7,8-tetrachlorodibenzo-p-dioxin. The approach involved deconstructing the reference dose (RfD) via evaluation of the different assumptions, options, models and methods associated with derivation of the value, culminating in the development of a plausible range of potential values based on such areas of uncertainty. The results demonstrate that overall RfD uncertainty was high based on limitations in the process for selection (e.g., compliance with inclusion criteria related to internal validity of the co-critical studies, consistency with other studies), external validity (e.g., generalizing findings of acute, high-dose exposure scenarios to the general population), and selection and classification of the point of departure using data from the individual studies (e.g., lack of statistical and clinical significance). Building on sensitivity analyses conducted by the US Environmental Protection Agency in 2012, the resulting estimates of RfD values that account for the uncertainties ranged from ~1.5 to 179 pg/kg/day. It is anticipated that the range of RfDs presented herein, along with the characterization of uncertainties, will improve risk assessments of dioxins and provide important information to risk managers, because reliance on a single toxicity value limits the information needed for making decisions and gives a false sense of precision and accuracy.

Statistical Approach to Identify Threshold and Point of Departure in Dose-Response Data

The study presents an integrated, rigorous statistical approach to define the likelihood of a threshold and point of departure (POD) based on dose-response data using nested family of bent-hyperbola models. The family includes four models: the full bent-hyperbola model, which allows for transition between two linear regiments with various levels of smoothness; a bent-hyperbola model reduced to a spline model, where the transition is fixed to a knot; a bent-hyperbola model with a restricted negative asymptote slope of zero, named hockey-stick with arc (HS-Arc); and spline model reduced further to a hockey-stick type model (HS), where the first linear segment has a slope of zero. A likelihood-ratio test is used to discriminate between the models and determine if the more flexible versions of the model provide better or significantly better fit than a hockey-stick type model. The full bent-hyperbola model can accommodate both threshold and nonthreshold behavior, can take on concave up and concave down shapes with various levels of curvature, can approximate the biochemically relevant Michaelis-Menten model, and even be reduced to a straight line. Therefore, with the use of this model, the presence or absence of a threshold may even become irrelevant and the best fit of the full bent-hyperbola model be used to characterize the dose-response behavior and risk levels, with no need for mode of action (MOA) information. Point of departure (POD), characterized by exposure level at which some predetermined response is reached, can be defined using the full model or one of the better fitting reduced models.

The Seveso accident: A look at 40 years of health research and beyond

A 1976 chemical factory explosion near Seveso, Italy exposed residents to high levels of 2,3,7,8-tetracholorodibenzo-p-dioxin (TCDD or dioxin). Dioxin is a known human carcinogen and potent endocrine disruptor. It is highly lipophilic and has a long half-life in humans. Much of what we know and can learn about the risks of dioxin exposure on human health arose from the tragic circumstances of Seveso. This review aims to describe the Seveso accident, summarize the results of 40 years of research on the health of the Seveso population since the accident, and discuss next-stage research on the health of Seveso residents, their children, and grandchildren.

Release and toxicity of adipose tissue-stored TCDD: Direct evidence from a xenografted fat model

BACKGROUND

Persistent organic pollutants (POPs) are known to accumulate in adipose tissues (AT). This storage may be beneficial by diverting POPs from other sensitive tissues or detrimental because of chronic release of pollutants as indirectly suggested during weight loss. The aim is to study the biological and/or toxic effects that chronic POP release from previously contaminated grafted AT could exert in a naïve mouse.

METHODS

C57BL/6J male mice were exposed intraperitoneally to 2,3,7,8-tetrachlorodibenzo-p-doxin (TCDD); their epididymal fat pads were collected and grafted on the back skin of uncontaminated recipient mice whose brain, liver, and epididymal ATs were analyzed (TCDD concentration, relevant gene expression). Kinetics of release and redistribution were modeled using Physiologically Based PharmacoKinetics (PBPK).

RESULTS

The grafts released TCDD over a period of 10 weeks with different kinetics of distribution in the three organs studied. A PBPK model was used to simulate the AT releasing process and the incorporation of TCDD into the major organs. At three weeks post-graft, we observed significant changes in gene expression in the liver and the host AT with signatures reminiscent of inflammation, gluconeogenesis and fibrosis as compared to the control.

CONCLUSIONS

This study confirms that AT-stored TCDD can be released and distributed to the organs of the recipient hence leading to distinct changes in gene expression. This original model provides direct evidence of the potential toxic-relevant effects when endogenous sources of contamination are present.

Risk for animal and human health related to the presence of dioxins and dioxin-like PCBs in feed and food

The European Commission asked EFSA for a scientific opinion on the risks for animal and human health related to the presence of dioxins (PCDD/Fs) and DL-PCBs in feed and food. The data from experimental animal and epidemiological studies were reviewed and it was decided to base the human risk assessment on effects observed in humans and to use animal data as supportive evidence. The critical effect was on semen quality, following pre- and postnatal exposure. The critical study showed a NOAEL of 7.0 pg WHO2005-TEQ/g fat in blood sampled at age 9 years based on PCDD/F-TEQs. No association was observed when including DL-PCB-TEQs. Using toxicokinetic modelling and taking into account the exposure from breastfeeding and a twofold higher intake during childhood, it was estimated that daily exposure in adolescents and adults should be below 0.25 pg TEQ/kg bw/day. The CONTAM Panel established a TWI of 2 pg TEQ/kg bw/week. With occurrence and consumption data from European countries, the mean and P95 intake of total TEQ by Adolescents, Adults, Elderly and Very Elderly varied between, respectively, 2.1 to 10.5, and 5.3 to 30.4 pg TEQ/kg bw/week, implying a considerable exceedance of the TWI. Toddlers and Other Children showed a higher exposure than older age groups, but this was accounted for when deriving the TWI. Exposure to PCDD/F-TEQ only was on average 2.4- and 2.7-fold lower for mean and P95 exposure than for total TEQ. PCDD/Fs and DL-PCBs are transferred to milk and eggs, and accumulate in fatty tissues and liver. Transfer rates and bioconcentration factors were identified for various species. The CONTAM Panel was not able to identify reference values in most farm and companion animals with the exception of NOAELs for mink, chicken and some fish species. The estimated exposure from feed for these species does not imply a risk.

Physiologically Based Pharmacokinetic Modeling in Pregnancy: A Systematic Review of Published Models

During recent years there has been a surge in developing and applying physiologically based pharmacokinetic (PBPK) models in pregnant women to better understand and predict changes in drug pharmacokinetics throughout pregnancy. As a consequence, the number of publications focusing on pregnancy PBPK models has increased substantially. However, to date these models, especially across various platforms, have not been systematically evaluated. Hence, this review aims to assess published PBPK models in pregnancy used for therapeutic purposes.

Prenatal contribution of 2, 2', 4, 4'-tetrabromodiphenyl ether (BDE-47) to total body burden in young children

Many scientists made estimates of the body burden of PBDEs from breastmilk and house dust. Interestingly, they have not included the prenatal contribution to the body burden in young children after birth. In order to address how the prenatal contribution is important in the risk assessment of PBDEs in infants up to five years old, we used the median measurements of BDE-47 as a model chemical in 108 neonates in Korea, and made simulations of its disposition out of body from birth to five years. During the simulation periods, the environmental exposure was considered for house dust, babyfood, breastmilk consumption, etc., with assumption of typical exposure scenario applicable to general infants in Korea. About 22% of the total amounts of BDE-47 in newborn remained up to 5years after birth. The relative amounts of BDE-47 from the prenatal source were 20%, 14%, 10%, 8%, 6%, and 4% of the total body burden for 1-, 2-, 3-, 4- and 5-year after birth, respectively. The contribution from breastfeeding was 95.2% and 92.2% of the total postnatal exposure amounts at 6-month and 1-year after birth, respectively. After cease of breastfeeding at 1-yr, house dust and food were the important sources of exposure up to 5-yr; however, their contributions to the bodyburden were negligible with consideration of the remaining amounts of the analytes from the breastmilk and prenatal exposure. Suggestively, the innate amounts and pharmacokinetics should be counted in estimating bodyburden of BDE-47.

Physiologically based pharmacokinetic toolkit to evaluate environmental exposures: Applications of the dioxin model to study real life exposures

Chlorinated dibenzo-p-dioxins (CDDs) are a series of mono- to octa-chlorinated homologous chemicals commonly referred to as polychlorinated dioxins. One of the most potent, well-known, and persistent member of this family is 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). As part of translational research to make computerized models accessible to health risk assessors, we present a Berkeley Madonna recoded version of the human physiologically based pharmacokinetic (PBPK) model used by the U.S. Environmental Protection Agency (EPA) in the recent dioxin assessment. This model incorporates CYP1A2 induction, which is an important metabolic vector that drives dioxin distribution in the human body, and it uses a variable elimination half-life that is body burden dependent. To evaluate the model accuracy, the recoded model predictions were compared with those of the original published model. The simulations performed with the recoded model matched well with those of the original model. The recoded model was then applied to available data sets of real life exposure studies. The recoded model can describe acute and chronic exposures and can be useful for interpreting human biomonitoring data as part of an overall dioxin and/or dioxin-like compounds risk assessment.

An assessment of dioxin exposure across gestation and lactation using a PBPK model and new data from Seveso

On July 10, 1976, an explosion at a chemical plant in Seveso, Italy, released up to 30kg of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-the most potent dioxin congener. Twenty years later, the Seveso Women's Health Study (SWHS) initiated a follow-up assessment of a cohort of female Seveso residents. Researchers collected serial blood, measured for TCDD levels, and recorded information about the women's medical history after the explosion. The study's aims were to: 1) modify the human PBPK model for TCDD (Emond et al. 2004; Emond et al. 2005; NCEA-USEPA, 2010) to include repetitive gestation and lactation; 2) simulate TCDD blood concentrations during different life stages including pregnancy and lactation, under different exposure scenarios; and 3) use this PBPK model to compare the influence of gestation and lactation on elimination of TCDD. After optimization of the model, it was assessed using data from the SWHS cohort. The 23 women in Subcohort A, were 4-39years old and in Subcohort B, the 18 women were 3-17years old when the explosion occurred. The model accurately predicted the blood concentrations during the 20years post-exposure, including periods of pregnancy and lactation. The model was also used to analyze the contribution of gestation and lactation to the mother's elimination of TCDD. The results suggest that gestation and lactation do not significantly impact TCDD blood elimination. Future efforts will focus on using additional data to evaluate the PBPK model and improving the mathematical descriptions of lactation and multiple gestations.

Uncertainty analysis in 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cancer dose-response for three occupational cohorts

While the U.S. EPA has issued a draft report with a 1% TCDD effective dose (ED01) of 87.9pg/kg/day based on continuous integration of key scientific evidence, a detailed and comprehensive uncertainty analysis has not been well documented. In this study, a new estimate for ED01 was derived based on uncertainty analysis by quantitatively assessing the potential bias arising from the selection of kinetic models, dose-response models and cohorts. The cumulative serum lipid concentration (CSLC) and cumulative body burden (CBB) were reconstructed as dose metrics using a concentration- and age-dependent pharmacokinetic model (CADM), physiologically based pharmacokinetic model (PBPK), and age-dependent half-life model (FV), and the reconstructed dose metrics based on CADM and PBPK were generally higher than those based on the FV model. Three dose-response curves (linear, multiplicative and power) were used to link dose metrics and cancer risk to estimate ED01, and the linear model resulted in the lowest ED01, followed by the power model and multiplicative model, for the same cohort. Meanwhile, ED01 based on the CADM model was the highest, followed by those based on the PBPK model and first-order model. Finally, the ED01 was estimated to be 17.03±7.83pg/kg/day by statistically analyzing the distribution of ED01 values based on various kinetic models, cohorts and dose-response models. The study presented here strengthens the scientific basis for understanding the potential health implications of TCDD exposure.

ITE and TCDD differentially regulate the vascular remodeling of rat placenta via the activation of AhR

Vascular remodeling in the placenta is essential for normal fetal development. The previous studies have demonstrated that in utero exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD, an environmental toxicant) induces the intrauterine fetal death in many species via the activation of aryl hydrocarbon receptor (AhR). In the current study, we compared the effects of 2-(1'H-indole-3'-carbonyl)-thiazole-4-carboxylic acid methyl ester (ITE) and TCDD on the vascular remodeling of rat placentas. Pregnant rats on gestational day (GD) 15 were randomly assigned into 5 groups, and were exposed to a single dose of 1.6 and 8.0 mg/kg body weight (bw) ITE, 1.6 and 8.0 µg/kg bw TCDD, or an equivalent volume of the vehicle, respectively. The dams were sacrificed on GD20 and the placental tissues were gathered. The intrauterine fetal death was observed only in 8.0 µg/kg bw TCDD-exposed group and no significant difference was seen in either the placental weight or the fetal weight among all these groups. The immunohistochemical and histological analyses revealed that as compared with the vehicle-control, TCDD, but not ITE, suppressed the placental vascular remodeling, including reduced the ratio of the placental labyrinth zone to the basal zone thickness (at least 0.71 fold of control), inhibited the maternal sinusoids dilation and thickened the trophoblastic septa. However, no marked difference was observed in the density of fetal capillaries in the labyrinth zone among these groups, although significant differences were detected in the expression of angiogenic growth factors between ITE and TCDD-exposed groups, especially Angiopoietin-2 (Ang-2), Endoglin, Interferon-γ (IFN-γ) and placenta growth factor (PIGF). These results suggest ITE and TCDD differentially regulate the vascular remodeling of rat placentas, as well as the expression of angiogenic factors and their receptors, which in turn may alter the blood flow in the late gestation and partially resulted in intrauterine fetal death.

Bayesian Analysis of a Lipid-Based Physiologically Based Toxicokinetic Model for a Mixture of PCBs in Rats

A lipid-based physiologically based toxicokinetic (PBTK) model has been developed for a mixture of six polychlorinated biphenyls (PCBs) in rats. The aim of this study was to apply population Bayesian analysis to a lipid PBTK model, while incorporating an internal exposure-response model linking enzyme induction and metabolic rate. Lipid-based physiologically based toxicokinetic models are a subset of PBTK models that can simulate concentrations of highly lipophilic compounds in tissue lipids, without the need for partition coefficients. A hierarchical treatment of population metabolic parameters and a CYP450 induction model were incorporated into the lipid-based PBTK framework, and Markov-Chain Monte Carlo was applied to in vivo data. A mass balance of CYP1A and CYP2B in the liver was necessary to model PCB metabolism at high doses. The linked PBTK/induction model remained on a lipid basis and was capable of modeling PCB concentrations in multiple tissues for all dose levels and dose profiles.

Distribution of carbon-14 labeled C60 ([14C]C60) in the pregnant and in the lactating dam and the effect of C60 exposure on the biochemical profile of urine

This study was conducted to determine the distribution of [(14)C]C60 in the pregnant rat and fetuses, and in the lactating rat and offspring. Pregnant rats were dosed on gestation day (gd) 15 and lactating rats were dosed on postnatal day (pnd) 8 via tail vein injection with a suspension of approximately 0.3 mg [(14)C]C60 kg(-1) body weight prepared in polyvinylpyrrolidone (PVP), or with PVP alone. Tissues were collected at 24 and 48 h after dosing. The largest portion of the administered dose was detected in the liver (approximately 43%, pregnant dam; approximately 35%, lactating dam) and lung (approximately 25%, lactating dam). Radioactivity (approximately 6%) was distributed to the reproductive tract, placenta and fetuses of the pregnant dam. Lactating rats had radioactivity distributed to the milk (3140 dpm g(-1) tissue, 24 h; 1620 dpm g(-1) tissue, 48 h), and to the pups' GI tract (2.8%, 24 h; 4.4% 48 h) and liver (<1%). Blood radioactivity was significant at 24 h (14-19%) and at 48 h (7%) after dosing; largely accounted for in the plasma fraction. Less that 4% of the dose was recovered in the maternal spleen, heart, brain, urine or feces. Metabolomics analysis of urine indicated that dams exposed to [(14)C]C60 had decreased metabolites derived from the Krebs cycle and increased metabolites derived from the urea cycle or glycolysis, as well as alterations in the levels of some sulfur-containing amino acids and purine/pyrimidine metabolites. This study demonstrated that [(14)C]C60 crosses the placenta and is transmitted to offspring via the dam's milk and subsequently systemically absorbed.

Interpretation of studies on the developmental reproductive toxicology of 2,3,7,8-tetrachlorodibenzo-p-dioxin in male offspring

There have been several studies on the maternal administration of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and effects in the reproductive tract of male offspring, subsequent to risk assessments undertaken in 2001. This review compares the methodology and results to examine key methodological features, and consistency in reported outcomes. Maternal dosing at >0.8 microg TCDD/kg causes lethality and weight loss, and it is difficult to distinguish between direct and indirect effects of TCDD at these dose levels. Statistically significant effects of maternal doses of <1 microg TCDD/kg (i.e. the dose levels relevant for risk assessment) on prostate weight or epididymal sperm counts in offspring were reported in the minority of studies. The pharmacokinetics of TCDD differs considerably between acute and chronic dosing, and with dose level of TCDD. On the basis of body burden, TCDD had different potency at inducing adverse effects in the only comparison study between acute and chronic dosing. Understanding of the pharmacokinetics of TCDD and relationship to adverse effects in offspring is required. These analyses identify key features of TCDD developmental toxicity in male offspring, and identify data needs for future risk assessment.

Physiologically based pharmacokinetic modeling of deltamethrin: development of a rat and human diffusion-limited model

Mirfazaelian et al. developed a physiologically based pharmacokinetic (PBPK) model for the pyrethroid pesticide deltamethrin in the rat. This model describes gastrointestinal (GI) tract absorption as a saturable process mediated by phase III efflux transporters which pump deltamethrin out of the intestinal enterocytes into the GI tract lumen, resulting in minimal net absorption at low concentrations and increasing absorption at higher concentrations. In the present study, the dose dependency in absorption of deltamethrin was examined in male Long Evans rats using po exposures predicted by the Mirfazaelian model to yield different po bioavailability values. No difference in the bioavailability from single po doses of 0.3 and 3.0 mg/kg deltamethrin was observed. Based on this finding, the Mirfazaelian PBPK model was modified to exclude a saturable absorption process. Other changes to the Mirfazaelian model included describing all tissue compartments with diffusion-limited kinetics and a single blood compartment. These changes improved model predictions of deltamethrin tissue concentration data from the present study and the literature. The rat model was then scaled to humans. The model predicted a twofold greater peak deltamethrin brain concentration and threefold greater area under the curve (AUC(0-48 h)) for humans following an po exposure of 1 mg/kg. Based on this model, humans would have greater distribution of deltamethrin to the brain for the same administered po dose compared to rats. The relative sensitivity to deltamethrin between rats and humans depends on both pharmacokinetic and pharmacodynamic differences. Species differences in the pharmacodynamic responses to deltamethrin between rats and humans remain uncharacterized.

A stochastic whole-body physiologically based pharmacokinetic model to assess the impact of inter-individual variability on tissue dosimetry over the human lifespan

Physiologically based pharmacokinetic (PBPK) models have proven to be successful in integrating and evaluating the influence of age- or gender-dependent changes with respect to the pharmacokinetics of xenobiotics throughout entire lifetimes. Nevertheless, for an effective application of toxicokinetic modelling to chemical risk assessment, a PBPK model has to be detailed enough to include all the multiple tissues that could be targeted by the various xenobiotics present in the environment. For this reason, we developed a PBPK model based on a detailed compartmentalization of the human body and parameterized with new relationships describing the time evolution of physiological and anatomical parameters. To take into account the impact of human variability on the predicted toxicokinetics, we defined probability distributions for key parameters related to the xenobiotics absorption, distribution, metabolism and excretion. The model predictability was evaluated by a direct comparison between computational predictions and experimental data for the internal concentrations of two chemicals (1,3-butadiene and 2,3,7,8-tetrachlorodibenzo-p-dioxin). A good agreement between predictions and observed data was achieved for different scenarios of exposure (e.g., acute or chronic exposure and different populations). Our results support that the general stochastic PBPK model can be a valuable computational support in the area of chemical risk analysis.

On the translation of uncertainty from toxicokinetic to toxicodynamic models--the TCDD example

When estimating human health risks from exposure to TCDD using toxicokinetic and toxicodynamic models, it is important to understand how model choice and assumptions necessary for modeling add to the uncertainty of risk estimates. Several toxicokinetic models have been proposed for the risk assessment of dioxins, in particular the elimination kinetics in humans has been a matter of constant debate. For a long time, a simple linear elimination kinetics has been common choice. Thus, it was used for the statistical analysis of the largest occupationally exposed cohort, the German Boehringer cohort. We challenge this assumption by considering, amongst others, a nonlinear modified Michaelis-Menten-type elimination kinetics, the so-called Carrier kinetics. Using the area under the lipid TCDD concentration time curve as dose metrics, we model the time to cancer-related death using the Cox proportional hazards model as toxicodynamic model. This risk assessment set-up was simulated in order to quantify uncertainty of both the dose (TCDD body burden) and the risk estimates, depending on the use of the kinetic model, variations of carcinogenic effect of TCDD and variations of latency period (lag time). If past exposure is estimated assuming a linear elimination kinetics although a Carrier kinetics actually holds, then high exposures in reality will be underestimated through statistical analysis and low exposures will be overestimated, respectively. This bias will carry over on the estimated individual concentration-time curves and the therefrom derived TCDD dose metric values. Using biased dose values when estimating a dose-response relationship will finally lead to biased risk estimates. The extent of bias and the decrease of precision are quantified in selected scenarios through this simulation approach. Our findings are in concordance with recent results in the field of dioxin risk assessment. They also reinforce the general demand for the scheduled uncertainty assessments in risk analyses.

Use of a physiologically based pharmacokinetic model for rats to study the influence of body fat mass and induction of CYP1A2 on the pharmacokinetics of TCDD

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a highly lipophilic chemical that distributes into adipose tissue, especially at low doses. However, at high doses TCDD sequesters in liver because it induces cytochrome P450 1A2 (CYP1A2) that binds TCDD. A physiologically based pharmacokinetic (PBPK) model was developed that included an inducible elimination rate of TCDD in the Sprague-Dawley rat. Objectives of this work were to characterize the influence of induction of CYP1A2 and adipose tissue mass fraction on the terminal elimination half-life (t1/2) of TCDD using this PBPK model. When the model assumes a fixed elimination of TCDD, t1/2 increases with dose, due to hepatic sequestration. Because experimental data indicate that the t1/2 of TCDD decreases with dose, the model was modified to include an inducible elimination rate. The PBPK model was then used to compare the t1/2 after an increase of adipose tissue mass fraction from 6.9 to 70%. The model suggests that at low exposures, increasing adipose tissue mass increases the terminal t1/2. However, at higher exposures, as CYP1A2 is induced, the relationship between adipose tissue mass and t1/2 reaches a plateau. This demonstrates that an inducible elimination rate is needed in a PBPK model in order to describe the pharmacokinetics of TCDD. At low exposures these models are more sensitive to parameters related to partitioning into adipose tissue.

Estimated cancer risk of dioxins to humans using a bioassay and physiologically based pharmacokinetic model

The health risk of dioxins and dioxin-like compounds to humans was analyzed quantitatively using experimental data and mathematical models. To quantify the toxicity of a mixture of three dioxin congeners, we calculated the new relative potencies (REPs) for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 1,2,3,7,8-pentachlorodibenzo-p-dioxin (PeCDD), and 2,3,4,7,8- pentachlorodibenzofuran (PeCDF), focusing on their tumor promotion activity. We applied a liver foci formation assay to female SD rats after repeated oral administration of dioxins. The REP of dioxin for a rat was determined using dioxin concentration and the number of the foci in rat liver. A physiologically based pharmacokinetic model (PBPK model) was used for interspecies extrapolation targeting on dioxin concentration in liver. Toxic dose for human was determined by back-estimation with a human PBPK model, assuming that the same concentration in the target tissue may cause the same level of effect in rats and humans, and the REP for human was determined by the toxic dose obtained. The calculated REPs for TCDD, PeCDD, and PeCDF were 1.0, 0.34, and 0.05 for rats, respectively, and the REPs for humans were almost the same as those for rats. These values were different from the toxic equivalency factors (TEFs) presented previously (Van den Berg, M., Birnbaum, L., Bosveld, A.T.C., Brunstrom, B., Cook, P., Feeley, M., Giesy, J.P., Hanberg, A., Hasegawa, R., Kennedy, S.W., Kubiak, T., Larsen, J.C., Rolaf van Leeuwen, F.X., Liem, A.K.D., Nolt, C., Peterson, R.E., Poellinger. L., Safe, S., Schrenk, D., Tillitt, D, Tysklind, M., Younes, M., Waern, F., Zacharewski, T., 1998. Toxic equivalency factors (TEFs) for PCBs, PCDDs, PCDFs for humans and wildlife. Environ. Health Perspect. 106, 775-792). The relative risk of excess liver cancer for Japanese people in general was 1.7-6.5 x 10(-7) by TCDD only, and 2.9-11 x 10(-7) by the three dioxins at the present level of contamination.

Adverse health effects in humans exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)

The environmental contaminant 2,3,7,8-tetrachlordibenzo-p-dioxin (TCDD) belongs to the category of highly toxic, persistent organic pollutants that accumulate in animal fat and plant tissues. Today, background TCDD levels in human fat are showing a decreasing trend. The food chain is the main source of exposure in the human population. TCDD regulates the expression of a wide range of drug-metabolizing enzymes and has an impact on a large number of biological systems. The most pronounced effects have occurred in occupational settings following the uncontrolled formation of TCDD after industrial accidents, as well as in rare intentional intoxications. Although the acute effects of TCDD exposure are well described in the literature, the long-term consequences have been underevaluated. The most well-known symptoms of severe acute intoxication are chloracne, porphyria, transient hepatotoxicity, and peripheral and central neurotoxicity. Because of the long-term persistence of TCDD in the human body, atherosclerosis, hypertension, diabetes, vascular ocular changes, and signs of neural system damage, including neuropsychological impairment, can be present several decades after massive exposure. Such chronic effects are nonspecific, multifactorial, and may be causally linked to TCDD only in heavily intoxicated subjects. This opinion is supported by the dose-dependent effect of TCDD found in exposed workers and by experimental animal studies.

The early embryo loss caused by 2,3,7,8-tetrachlorodibenzo-p-dioxin may be related to the accumulation of this compound in the uterus

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a reproductive and developmental toxicant that can alter endocrine status, leading to decreased fertility and altered embryonic development; however, there are limited reports on TCDD toxicity during early pregnancy. In the present study, pregnant and pseudopregnant NIH mice were exposed to TCDD orally (2, 50 and 100 ng/kg body weight) during early gestation (days 1-8), pre-implantation stages (days 1-3), and peri-implantation to early post-implantation stages (days 4-8). TCDD concentration in uterus, liver, kidney, brain and fat on day 9 of pregnancy was monitored by an aryl hydrocarbon receptor (Ah receptor, AhR)-mediated LacZ reporter system in yeast. Results showed that the number of implanted embryos was significantly reduced on day 9 of gestation by 50 and 100 ng/kg TCDD exposure. The number of implantation sites was lower for animals exposed to TCDD on days 1-3 versus those exposed during days 4-8. Decidualization in pseudopregnant mice was also inhibited by TCDD exposure. TCDD concentrations as low as 2 ng/kg significantly decreased serum progesterone levels but had no effect on serum estradiol. TCDD level in the uterus was equal to levels in the liver, but lower than the fat tissue. These results suggest that TCDD sensitivity might be attributed its local accumulation in the uterus.

Comparison of the use of a physiologically based pharmacokinetic model and a classical pharmacokinetic model for dioxin exposure assessments

In epidemiologic studies, exposure assessments of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) assume a fixed elimination rate. Recent data suggest a dose-dependent elimination rate for TCDD. A physiologically based pharmacokinetic (PBPK) model, which uses a body-burden-dependent elimination rate, was developed previously in rodents to describe the pharmacokinetics of TCDD and has been extrapolated to human exposure for this study. Optimizations were performed using data from a random selection of veterans from the Ranch Hand cohort and data from a human volunteer who was exposed to TCDD. Assessment of this PBPK model used additional data from the Ranch Hand cohort and a clinical report of two women exposed to TCDD. This PBPK model suggests that previous exposure assessments may have significantly underestimated peak blood concentrations, resulting in potential exposure misclassifications. Application of a PBPK model that incorporates an inducible elimination of TCDD may improve the exposure assessments in epidemiologic studies of TCDD.

Comparison of the concentrations of polychlorinated dibenzo-p-dioxins, dibenzofurans, and dioxin-like polychlorinated biphenyls in maternal and fetal blood, amniotic and allantoic fluids in cattle

To characterize the maternal-fetal transport of lipophilic endocrine disrupting chemicals, concentrations of polychlorinated (2,3,7,8-substituted) dibenzo-p-dioxins (PCDDs), dibenzofurans (PCDFs), and dioxin-like polychlorinated biphenyls (PCBs) were measured in maternal and fetal blood, and amniotic and allantoic fluids in cattle. Total toxicity equivalent quantity (TEQ) was highest in amniotic fluid on a fat-weight basis, whereas it was highest in maternal blood on a total weight basis. TEQ was lowest in allantoic fluid on either basis; 26 of 29 congeners analyzed in this experiment were detected in one or more samples. The largest number of congeners was detected in amniotic fluid. O8CDD, 2,3,4,7,8-P5CDF and 2,3',4,4',5-P5CB were the major congeners in PCDDs, PCDFs and PCBs, respectively. The O8CDD concentration was higher in fetal blood than in maternal blood on a fat-weight basis, whereas concentrations of other congeners were lower in fetal blood than in maternal blood. Furthermore, on a fat-weight basis, the O8CDD concentration was considerably higher in allantoic fluid compared with other samples. Concentrations of major PCB congeners were higher in amniotic fluid than in maternal and fetal blood on a fat-weight basis. In conclusion, it is suggested that lipophilic endocrine-disrupting chemicals contained in maternal blood are all transferred to the fetal circulation via the placenta in cattle. Furthermore, the results of this experiment imply that O8CDD has different transportation systems from other dioxins in the circulation, and that a considerable amount of PCBs is excreted and accumulated in amniotic fluid during the fetal stage in cattle.

Computational pharmacokinetics during developmental windows of susceptibility

Computational modeling has an increasing role in analyses of biological effects, including how the body handles chemicals (i.e., pharmacokinetics or toxicokinetics) and how the body responds to chemicals (i.e., pharmacodynamics or toxicodynamics). Pharmacokinetic models increasingly describe not just adult humans and animals, but also changes with age and life stage (e.g., pregnancy and fetal exposures, lactational exposures, and childhood growth). Physiologically based pharmacokinetic models provide an important route to estimate the potential changes in internal dose that may occur throughout the life cycle. These models require inputs describing changes in physiology, metabolism, and exposure with age and life stage. A particular challenge exists when the "equivalent" developmental period in the rodents and humans differs (e.g., early postnatal in rats and in utero in humans) such that the "equivalent" window of susceptibility to toxic effects of the chemical may involve substantially different exposures (e.g., lactational versus placental transfer). Pharmacodynamic modeling could similarly address changes with age, but few such models currently exist. The growth of systems biology is anticipated to change this over the coming decade.