Modeling receptor-mediated processes with dioxin: implications for pharmacokinetics and risk assessment

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.

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 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.

Toxicokinetic Modeling as a Tool for Risk Estimation: 2,3,7,8-Tetrachlorodibenzo-P-Dioxin

Concepts of toxicokinetic modeling and the relevance of toxicokinetics for understanding dose-response relationships, species scaling, and risk estimation are broached. A physiological one-compartment model for 2,3,7,8-tetra-chlorodibenzo-p-dioxin (TCDD) is presented in detail. It describes the TCDD burden of the human body, which results from TCDD-contaminated food, in dependence of age. The model was validated using a series of measured values obtained by other authors and this group. They represent lipid-based concentrations of TCDD in liver, blood, adipose tissue, feces, and mother's milk in dependence of age. Special attention was paid to the TCDD burden in infants resulting from feeding with mother's milk or formula. Model simulations demonstrate that TCDD burden can amount to 10 mg/kg of lipids after nursing for 6 months with mother's milk exclusively This is still within the range of the concentrations found in adults. After the nursing period, TCDD burden declines. From the age of 7 years on, there is no longer a difference in the TCDD burden, independently of the food they had received as infants. According to the model, elimination half-life of TCDD from the body is not constant but increases during life starting from a few months in newborns to several years in adults.

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.

My Winding Road: From Microbiology to Toxicology and Environmental Health

I would certainly never have predicted that I would become the director of the National Institute of Environmental Health Sciences (NIEHS) and the National Toxicology Program (NTP) when I was a Jewish girl growing up in Teaneck, New Jersey. My family stressed the importance of education. Yet for a girl there were many not-so-subtle suggestions that the appropriate careers were in teaching or nursing, and the most important thing was to be a wife and mother. Well, I can't disagree with the latter, although I would have to add grandmother to that list of achievements. My parents were both college graduates, but my mom only taught high school English for one year before leaving the field to start our family. My dad returned from World War II and joined his brother in accounting. After my first sister was born, my father joined my mother's family jewelry business and helped to open a second retail store. My mother helped my dad out during the busy times—Christmas and wedding season—but otherwise focused on our growing family of three girls and one boy. This became increasingly challenging when it became clear that my little brother was severely retarded and would require extra care.

A self-regulating cell culture analog device to mimic animal and human toxicological responses

The overall goal of this project is the development of a new methodology for translating advances in molecular level understanding of toxicological responses into a predictive tool for dose response in whole animals and humans exposed to single compounds or mixtures of compounds. The methodology incorporates a mechanistic cellular level model into a PBPK (physiologically based pharmacokinetic) model which simultaneously guides the development of an in vitro cell culture analog (CCA) to the PBPK. Where the PBPK specifies an organ, (e.g., liver) the in vitro or CCA system contains a compartment with the appropriate cell or cell population (e.g., hepatocytes for the liver). The CCA has significant advantages over other in vitro systems and PBPK systems used independently for evaluating metabolic responses to drugs or potentially toxic chemicals where the exchange of metabolites between organs is likely to be important. The CCA system is superior to a PBPK because an a priori description of complete metabolism is not required and secondary, unexpected interactions can be detected. The CCA system, unlike other in vitro systems, gives a dynamic response that realistically simulates in vivo interactions between organs. Furthermore, the CCA allows dosing on the same basis as animal tests (e.g., milligrams per kilogram of body mass equivalent). Because the construction of a CCA is guided by a PBPK, this approach allows extrapolation to low doses and across species, including extrapolation to humans. We have constructed a prototype system and have conducted proof-of-concept experiments using naphthalene as a test chemical. These experiments clearly demonstrate the ability to generate a reactive metabolite in one compartment and detect its effects (on LDH release and glutathione depletion) in a second compartment. However, this prototype device would be expensive to replicate and requires nearly constant supervision from a trained investigator. For this concept to replace animals an inexpensive, self-regulating device is needed. An initial design to accomplish this goal is described as well as the corresponding model using naphthalene as a test compound. (c) 1996 John Wiley & Sons, Inc.

Physiologically based pharmacokinetic/toxicokinetic modeling

Physiologically based pharmacokinetic (PBPK) models differ from conventional compartmental pharmacokinetic models in that they are based to a large extent on the actual physiology of the organism. The application of pharmacokinetics to toxicology or risk assessment requires that the toxic effects in a particular tissue are related in some way to the concentration time course of an active form of the substance in that tissue. The motivation for applying pharmacokinetics is the expectation that the observed effects of a chemical will be more simply and directly related to a measure of target tissue exposure than to a measure of administered dose. The goal of this work is to provide the reader with an understanding of PBPK modeling and its utility as well as the procedures used in the development and implementation of a model to chemical safety assessment using the styrene PBPK model as an example.

Multi-dose-route, multi-species pharmacokinetic models for manganese and their use in risk assessment

Manganese (Mn) is an essential element that may be toxic in conditions of overexposure. Nearly 10 years ago, some of the authors of this article published a proposed methodology to perform a tissue-dose-based risk assessment and a detailed list of data needs necessary to perform the assessment. Since that time, a substantial body of Mn pharmacokinetic (PK) data has been generated in rats and nonhuman primates, allowing for the construction of physiologically based pharmacokinetic (PBPK) models for Mn. This study reviews the development of the Mn PBPK models, reassesses the previously identified data needs, and details potential uses of these models in risk assessment of Mn. Based upon numerous animal experiments, pharmacokinetic (PK) models have effectively simulated tissue kinetics of Mn from both inhaled and oral Mn intake. PK models achieve this by incorporating homeostatic control processes, saturable tissue binding capacities, and preferential fluxes in various tissue regions. While minor data gaps still exist, the models captured the main dose-dependent characteristics of Mn disposition in rodents and monkeys and provide a structure to parameterize an equivalent PK description in humans. These models are organized to contribute to a tissue-dose based risk assessment of Mn that simultaneously considers ingestion and inhalation kinetics of Mn along with homeostatic control of Mn.

Development of a Physiologically Based Pharmacokinetic Model for Decane, a Constituent of Jet Propellent-8

Decane, a 10-carbon n-alkane and one of the highest vapor phase constituents of jet propellent-8 (JP-8), was selected to represent the semivolatile fraction for the initial development of a physiologically based pharmacokinetic (PBPK) model for JP-8. Rats were exposed to decane vapors at time-weighted average concentrations of 1200, 781, or 273 ppm in a 32-L Leach chamber for 4 h. Time-course samples for 1200 ppm and end-of-exposure samples for 781 and 273 ppm decane exposures were collected from blood, brain, liver, fat, bone marrow, lung, skin, and spleen. The pharmacokinetics of decane could not be described by flow-limited assumptions and measured in vitro tissue/air partition coefficients. A refined PBPK model for decane was then developed using flow-limited (liver and lung) and diffusion-limited (brain, bone marrow, fat, skin, and spleen) equations to describe the uptake and clearance of decane in the blood and tissues. Partition coefficient values for blood/air and tissue/blood were estimated by fitting end-of-exposure pharmacokinetic data and assumed to reflect the available decane for rapid exchange with blood. A portion of decane is speculated to be sequestered in "deep" pools in the body, unavailable for rapid exchange with blood. PBPK model predictions were adequate in describing the tissues and blood kinetics. For model validation, the refined PBPK model for decane had mixed successes at predicting tissue and blood concentrations for lower concentrations of decane vapor, suggesting that further improvements in the model may be necessary to extrapolate to lower concentrations.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

Physiologically based pharmacokinetic (PBPK) models are particularly useful for simulating exposures to environmental toxicants for which, unlike pharmaceuticals, there is often little or no human data available to estimate the internal dose of a putative toxic moiety in a target tissue or an appropriate surrogate. This article reviews the current state of knowledge and approaches for application of PBPK models in the process of deriving reference dose, reference concentration, and cancer risk estimates. Examples drawn from previous U.S. Environmental Protection Agency (EPA) risk assessments and human health risk assessments in peer-reviewed literature illustrate the ways and means of using PBPK models to quantify the pharmacokinetic component of the interspecies and intraspecies uncertainty factors as well as to conduct route to route, high dose to low dose and duration extrapolations. The choice of the appropriate dose metric is key to the use of the PBPK models for the various applications in risk assessment. Issues related to whether uncertainty factors are most appropriately applied before or after derivation of human equivalent dose (or concentration) continue to be explored. Scientific progress in the understanding of life stage and genetic differences in dosimetry and their impacts on variability in susceptibility, as well as ongoing development of analytical methods to characterize uncertainty in PBPK models, will make their use in risk assessment increasingly likely. As such, it is anticipated that when PBPK models are used to express adverse tissue responses in terms of the internal target tissue dose of the toxic moiety rather than the external concentration, the scientific basis of, and confidence in, risk assessments will be enhanced.

A possible role of multidrug resistance-associated protein 2 (Mrp2) in hepatic excretion of PCB126, an environmental contaminant: PBPK/PD modeling

3,3',4,4',5'-Pentachlorobiphenyl (PCB126) is a carcinogenic environmental pollutant and its toxicity is mediated through binding with aryl hydrocarbon receptor (AhR). Earlier, we found that PCB126 treated F344 rats had 110-400 times higher PCB126 concentration in the liver than in the fat. Protein binding was suspected to be a major factor for the high liver concentration of PCB126 despite its high lipophilicity. In this research, we conducted a combined pharmacokinetic/pharmacodynamic study in male F344 rats. In addition to blood and tissue pharmacokinetics, we use the development of hepatic preneoplastic foci (glutathione-S-transferase placental form [GSTP]) as a pharmacodynamic endpoint. Experimental data were utilized for building a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model. PBPK/PD modeling was consistent with the experimental PK and PD data. Salient features of this model include: (1) bindings between PCB126 and hepatic proteins, particularly the multidrug resistance-associated protein (Mrp2), a protein transporter; (2) Mrp2-mediated excretion; and (3) a relationship between area under the curve of PCB126 in the livers and % volume of GSTP foci. Mrp2 involvement in PCB126 pharmacokinetics is supported by computational chemistry calculation using a three-dimensional quantitative structure-activity relationship model of Mrp2 developed by S. Hirono et al. (2005, Pharm. Res. 22, 260-269). This work, for the first time, provided a plausible role of a versatile hepatic transporter for drugs, Mrp2, in the disposition of an important environmental pollutant, PCB126.

Pharmacokinetic modeling of manganese in the rat IV: Assessing factors that contribute to brain accumulation during inhalation exposure

A recently published physiologically based pharmacokinetic (PBPK) model successfully accounted for steady-state tissue manganese (Mn) concentration seen with normal dietary intakes and for biphasic, whole-body time-course profiles observed with tracer (54Mn) dosing. In this present study, PBPK modeling was used to evaluate Mn kinetics and brain concentrations in rats exposed to Mn both in their diet and by inhalation. Three published studies were used: (1) rats fed on diets ranging from 2 to 100 ppm, (2) rats on 125 ppm in diet and exposed via inhalation at 0.0 to 3.00 mg Mn/m3 each day for 14 d, and (3) rats to 0.1 or 0.5 mg Mn/m3 for 6 h/d, 5 d/wk over a 90-d period. The original model structure with well-mixed and "deep" compartments for each tissue could not describe rapid increases in tissue concentrations and rapid declines seen in high concentration inhalation studies. A second structure was developed that included (1) saturable, high-affinity binding of Mn in all tissues and (2) asymmetric diffusion from blood into brain (i.e., transport into and out of specific brain regions such as the striatum was described with different diffusion constants). This second model was consistent with liver and striatum experimental data. Preferential increases in some brain regions were predicted for exposures above 0.2 mg/m3 and had a rapid (i.e., 1 or 2 wk) return to steady-state levels. Multi-dose-route PBPK models for Mn based on this alternative model structure can be readily scaled to evaluate tissue Mn kinetics in other species and for human populations. Once validated across test animals, these PBPK models will be useful in tissue-dose based risk assessment with manganese.

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.

Development of a physiologically based pharmacokinetic model for deltamethrin in the adult male Sprague-Dawley rat

Deltamethrin (DLT) is a type II pyrethroid insecticide widely used in agriculture and public health. DLT is a potent neurotoxin that is primarily cleared from the body by metabolism. To better understand the dosimetry of DLT in the central nervous system, a physiologically based pharmacokinetic (PBPK) model for DLT was constructed for the adult, male Sprague-Dawley rat that employed both flow-limited (brain, gastrointestinal [GI] tract, liver, and rapidly perfused tissues) and diffusion-limited (fat, blood/plasma, and slowly perfused tissues) rate equations. The blood was divided into plasma and erythrocytes. Cytochrome P450-mediated metabolism was accounted for in the liver and carboxylesterase (CaE)-mediated metabolism in plasma and liver. Serial blood, brain, and fat samples were taken for DLT analysis for up to 48 h after adult rats received 2 or 10 mg DLT/kg po. Hepatic biotransformation accounted for approximately 78% of these administered doses. Plasma CaEs accounted for biotransformation of approximately 8% of each dosage. Refined PBPK model forecasts compared favorably to the 2- and 10-mg/kg po blood, plasma, brain, and fat DLT profiles, as well as profiles subsequently obtained from adult rats given 1 mg/kg iv. DLT kinetic profiles extracted from published reports of oral and iv experiments were also used for verification of the model's simulations. There was generally good agreement in most instances between predicted and the limited amount of empirical data. It became clear from our modeling efforts that there is considerably more to be learned about processes that govern GI absorption and exsorption, transport, binding, brain uptake and egress, fat deposition, and systemic elimination of DLT and other pyrethroids. The current model can serve as a foundation for construction of models for other pyrethroids and can be improved as more definitive information on DLT kinetic processes becomes available.

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.

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.

Carcinogenic risks of dioxin: Mechanistic considerations

Dioxins and dioxin-like chemicals demonstrate high affinity binding to the aryl hydrocarbon receptor (AhR), a ligand activated transcription factor, which mediates most, if not all, of the toxic responses of these agents. Since dioxins are not directly genotoxic their carcinogenic effect is likely the result of their tumor promoting activity produced by activation of the AhR. For the purpose of risk assessment extrapolation from effects in the observable high dose range to background dietary exposure is necessary. In the present review, we discuss various aspects of low-dose-response of receptor-mediated processes in general, including threshold phenomena with regard to tumor promotion during multi-stage carcinogenesis. In this connection the reversibility of tumor promotion plays an important role but this may not be valid for dioxins due to their long half life. The relevance of cytochrome P 4501 A-induction as biomarker for prediction of carcinogenic effects of dioxins at low doses is considered. Dioxins may act in concert with endogenous ligands of the AhR, an effect which becomes particularly relevant at low toxicant concentrations. At present, however, the nature and role of these postulated ligands are unknown. Furthermore, it is unclear whether dioxins produce synergistic tumor promotional effects with non-dioxin-like chemicals to which humans are also exposed. Dioxins and, e.g., non-dioxin-like PCBs act through different receptors and there is, albeit yet limited, experimental evidence from experimental studies to suggest that they may act on different target cell populations within the same target organ. From the available data the existence of a (physiological) threshold of effects cannot be proven and may not even exist. For regulatory purposes the application of a so called "practical threshold" for the carcinogenic effect of dioxins is proposed. Further mechanistic studies should be conducted to get insight into the dose-response characteristics of relevant events of dioxin-like and non-dioxin-like agents and into the consequences of potential interactions between both group of compounds during carcinogenesis.

Dioxin health risk to infants using simulated tissue concentrations

Dioxin concentrations in infant and child were simulated using physiologically based pharmacokinetic (PBPK) models developed for these groups. The infant model was validated by comparing the simulated concentration with the measured concentration from the literature, and they showed good agreement. Simulations with our PBPK model showed temporal patterns in concentrations in various tissues. For risk assessment, estimated concentrations of 29 dioxins in the liver were summed up in a toxic equivalency (TEQ) basis to be compared with actual 2,3,7,8-TCDD concentrations in rat liver associated with toxicity. Maximum liver concentrations in breast-fed and formula-fed infants were 16.8pg TEQ/g and 3.5pg TEQ/g, respectively. The level in breast-fed infant liver was approximately 1/300 of the level associated with hepatocellular carcinoma and 1/5 of the level found in maternal rat liver associated with alterations in reproductive organs in the next generation. Based on our analysis, the present contamination level is not safe enough, but further dose-response data is required for a quantitative risk assessment.

Mechanistic approaches for mixture risk assessments-present capabilities with simple mixtures and future directions

Mechanistic studies with simple mixtures have provided insights into the nature of interactions among chemicals that lead to non-additive effects and have elucidated the exposure conditions under which interactions are likely to occur. This paper discusses studies on four mixtures: (1) 1,1-dichloroethylene and trichloroethylene, (2) carbon tetrachloride and Kepone, (3) hexane and methyl-n-butylketone, and (4) coplanar and non-coplanar polychlorinated biphenyls. These mechanistic studies show that interactions should be described at the level of target tissue dose and are best categorized as either pharmacokinetic (PK) or pharmacodynamic (PD) interactions. In PK interactions the presence of a second chemical alters the kinetics such that a unit of administered dose no longer produces a unit of dose at the target tissue. In PD interactions, the presence of other compounds alters the PDs such that a unit tissue dose no longer produces a unit of response. Physiologically based pharmacokinetic (PBPK) models for mixtures have become important tools for predicting conditions under which interactions are likely to alter the assumption of additivity and have permitted calculation of interaction thresholds with more confidence. New cumulative risk assessment approaches have provided opportunities to classify compounds on the basis of similar chemistry-based modes of action (cholinesterase inhibitors) or similar physiological modes of action (diverse chemicals that alter a common biological outcome, such as defeminization of the developing nervous system). The latter examples present challenges for expanding our risk assessment paradigm to focus on the biology of responses more than on the kinetics of the xenobiotics. Some of the future advances in mixture research will depend on progress in systems biology, a discipline that integrates information across multiple level of biological organization producing PD models of normal function and assessing conditions under which exposures to chemicals lead to the perturbations sufficiently great to produce toxicity and disease. We describe briefly the elements of a systems biology approach for assessing the interactions between various PCB congeners.

Simulation of dioxin accumulation in human tissues and analysis of reproductive risk

Polychlorinated dibenzo-p-dioxin, polychlorinated dibenzofuran and dioxin-like polychlorinated biphenyl concentrations in human liver, kidney, fat, blood, muscle, richly perfused tissue (brain, lung etc.) and skin were simulated to assess the health risk for Japanese fetuses. A 40-year time course of dioxin accumulation via food ingestion was simulated using a physiologically based pharmacokinetic (PBPK) model. In richly perfused tissue, the concentration estimated by the PBPK model showed better agreement with measured concentrations than that calculated by the one-compartment model. Fetal dioxin concentration was simulated based on the assumption that the fetal concentration was almost equal to the concentration in the mother's richly perfused tissue. To assess the reproductive risk, the estimated concentration in human fetus was compared with that in rat fetus in which reproductive function showed signs of alteration by 2,3,7,8-TCDD in previous reports [Toxicol. Appl. Pharmacol. 114 (1992) 118; 146 (1997) 11; Toxicol. Sci. 53 (2000) 411; 57 (2000) 275]. The present daily intake of 2,3,7,8-TCDD is approximately 1/50 of the amount that leads to possible reproductive toxicity in the next generation. However, when 29 kinds of dioxin congeners are considered, the present level is 1/5 of the hazardous levels. For species extrapolation of dioxin risk, further study on tissue concentration versus toxicity is required.

Toxicokinetic modeling and its applications in chemical risk assessment

In recent years physiologically based pharmacokinetic (PBPK) modeling has found frequent application in risk assessments where PBPK models serve as important adjuncts to studies on modes of action of xenobiotics. In this regard, studies on mode of action provide insight into both the sites/mechanisms of action and the form of the xenobiotic associated with toxic responses. Validated PBPK models permit calculation of tissue doses of xenobiotics and metabolites for a variety of conditions, i.e. at low-doses, in different animal species, and in different members of a human population. In this manner, these PBPK models support the low-dose and interspecies extrapolations that are important components of current risk assessment methodologies. PBPK models are sometimes referred to as physiological toxicokinetic (PT) models to emphasize their application with compounds causing toxic responses. Pharmacokinetic (PK) modeling in general has a rich history. Data-based PK compartmental models were developed in the 1930's when only primitive tools were available for solving sets of differential equations. These models were expanded in the 1960's and 1970's to accommodate new observations on dose-dependent elimination and flow-limited metabolism. The application of clearance concepts brought many new insights about the disposition of drugs in the body. In the 1970's PBPK/PT models were developed to evaluate metabolism of volatile compounds of occupational importance, and, for the first time, dose-dependent processes in toxicology were included in PBPK models in order to assess the conditions under which saturation of metabolic and elimination processes lead to non-linear dose response relationships. In the 1980's insights from chemical engineers and occupational toxicology were combined to develop PBPK/PT models to support risk assessment with methylene chloride and other solvents. The 1990's witnessed explosive growth in risk assessment applications of PBPK/PT models and in applying sensitivity and variability methods to evaluate model performance. Some of the compounds examined in detail include butadiene, styrene, glycol ethers, dioxins and organic esters/aids. This paper outlines the history of PBPK/PT modeling, emphasizes more recent applications of PBPK/TK models in health risk assessment, and discusses the risk assessment perspective provided by modern uses of these modeling approaches.

Toxicokinetics and physiologically based toxicokinetics in toxicology and risk assessment

Toxicokinetics is the study of kinetics of absorption, distribution, metabolism, and excretion of a xenobiotic under the conditions of toxicity evaluation. Conventional toxicokinetics uses the hypothetical compartments, and the model is composed of rate equations that describe the time course of drug and chemical disposition. The utility of toxicokinetics in toxicity evaluation and interpretation of animal toxicology data is emerging as an important tool in product discovery and development. With implementation of the International Conference on Harmonization (ICH) guidelines on systemic exposure and dose selection, toxicokinetics have been integrated in routine toxicity evaluations. Although traditional compartmental/noncompartmental models are generally adequate for assessing systemic exposure, they are unable to the predict time course of drug disposition in target tissues and often fail to relate systemic drug levels to a biological response. Physiologically based toxicokinetic (PB-TK) models address this deficiency of traditional compartmental models. PB-TK models are the kinetic models of the uptake and disposition of chemicals based on rates of biochemical reactions, physiological and anatomical characteristics. These models, when developed appropriately, can predict target organ drug distribution in different species under variety of conditions. This minireview discusses the basic principles, and applications of traditional compartmental toxicokinetic and physiologically based toxicokinetics (PB-TK) models in drug development and risk assessment. Special emphasis will be placed on discussion related to interpretation of the ICH guidelines related to toxicokinetics and the utility of toxicokinetics data in dose selection for toxicity and carcinogenicity studies. The utility of PB-TK models in risk assessment of methylene chloride, vinyl chloride, retinoic acid, dioxin, and inhaled organic esters is discussed.

Population toxicokinetic analysis of 2,3,7,8-tetrachlorodibenzo-p-dioxin using Bayesian techniques

Understanding the kinetics of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) concentrations in humans is an important step for TCDD cancer risk assessment. In this paper longitudinal series of serum TCDD concentration measurements on U.S. veterans of the Vietnam war, who were exposed to dioxin during herbicide-spraying operations, are studied. The overall aim is to use these data to infer the dynamics of TCDD concentrations in humans. This is done by identifying a kinetic model describing the dioxin time course at the individual level. The individual toxicokinetic model is then expanded into a population model within a Bayesian hierarchical framework which allows residual variations across subjects that cannot be explained by observed covariates. Other complications in the data, such as unknown exposure histories, are also resolved implicitly through the hierarchical model. Moreover, the choice of a Bayesian approach enables the accumulation of external source of information in the form of prior distributions. The model is subjected to various diagnostic checks and analyses of sensitivity to distributional assumptions showing a good fit in terms of both the population and the kinetic features.

Determination of tissue-blood partition coefficients for a physiological model for humans, and estimation of dioxin concentration in tissues

The tissue-blood partition coefficients for a physiologically based pharmacokinetic (PBPK) model were determined, and the concentrations of 17 congeners of polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/Fs) in tissues in Japanese people were estimated using the model. According to the PBPK model established by Lawrence and Gobas [Chemosphere 35 (1997) 427-452], we assumed a steady-state fugacity model for Japanese people in general, and set the route of PCDD/Fs exposure only from food intake. The required partition coefficients for liver, kidney, adipose, muscle, skin, bile, gut and viscera (richly perfused tissue) were calculated using available autopsy data from eight Japanese men and women who were not accidentally exposed to PCDD/Fs. For validation of the partition coefficients, estimated PCDD/F concentrations in liver, kidney, fat, blood and muscle using the model were compared to other two sets of measured concentration data in Japanese tissues. Good agreement was obtained between estimated data and measured data, and most of the measured data were within the simulated concentration range in liver, kidney, blood and muscle. From these results, our model and calculated partition coefficients seem applicable for the estimation of congener-specific concentrations in human tissues.

The use of quantitative histological and molecular data for risk assessment and biologically based model development

In organs with diverse cell populations, it is not uncommon for one type of cell to respond while others are spared. Even in an organ with common cell types, such as hepatocytes within the liver, the population of cells may respond with different sensitivities for injury or for biochemical responses to toxicants. In the liver, many tumor promoters induce cytochrome P450 enzymes and other proteins in centrilobular cells at much lower doses than required to cause induction in periportal cells. In addition, these induction responses appear to occur at the level of individual cells--a 50% response of the liver for induction does not represent 50% induction in all cells. Instead, half of the cells are fully induced and half are unaffected. Cells "switch" from one phenotypic state to another. Over the past 10 years, several attempts have been made to model these cellular switches and to understand their relevance for hepatic tumor promotion and risk assessment. The data used for analyzing these switches include responses of the entire liver (total induction), responses of individual cells in the liver (regional induction), and cellular responses such as proliferation and apoptosis. This brief overview describes the development of biologically based, dose-response (BBDR) models for protein induction and tumor promotion in liver by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) with emphasis on the role of specific types of histological and molecular data in providing insights about mechanisms for cellular switches and their implications for tumor promotion. As the biological basis of these switches become unraveled and incorporated into the models, these BBDR models should eventually serve to improve risk assessments with a variety of liver tumor promoters with receptor-based modes of action.

Mode of action and tissue dosimetry in current and future risk assessments

Two fundamental concepts have emerged to organize contemporary approaches to chemical risk assessment - mode of action and tissue dosimetry. Mode of action specifies the nature of the interactions between the chemical and the body that lead to toxic responses and should, under optimal circumstances, also specify the form of the tissue dose that leads to these effects. This paper highlights recent development of biologically based dose response (BBDR) models for specific toxic endpoints that use knowledge on mode of action to specify measures of dose. These dose measures then are used to support low dose and interspecies extrapolations. We first focus on a series of dose response models developed for several compounds that produce nasal toxicity. These examples demonstrate a range of model structures from simple dosimetry models (methylmethacrylate) to linkage of dosimetry with specific biological processes involved in carcinogenesis (formaldehyde). Two BBDR models with dioxin illustrate the organization of biological and dosimetry information into specific testable hypotheses that could distinguish these different models and lead to a more uniform approach to risk assessment for this compound. A final section discusses the impact of molecular biology and the genomic revolution in relation to development of BBDR models for specific toxic endpoints.

Physiological modeling of a proposed mechanism of enzyme induction by TCDD

A physiological model was previously constructed to facilitate extrapolation of surrogates for the effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in rat liver to doses comparable to human environmental exposures. The model included induction of P450 isozymes and suggested the presence of multiple binding sites with different affinities for the TCDD-liganded Ah receptor at CYP1A1 dioxin responsive elements. The model also indicated that protein synthesis on the mRNA template exhibited saturation kinetics with respect to message levels. In the present work the earlier model was revised to include the increased proteolysis of the Ah receptor on binding TCDD, more realistic representations of gene transcription and mRNA translation, and different stability for each mRNA. The revised model includes multiple TCDD-liganded Ah receptor binding sites for CYP1A1 and CYP1B1 genes, a lag of 0.2 day for production of mRNA and induced proteins, and stabilization of mRNA by a poly(A) tail. The model reproduced the transient depletion of the Ah receptor subsequent to binding ligand and the dose-response of the receptor in rats treated with biweekly oral doses of TCDD in corn oil. The model reproduced tissue TCDD concentrations observed for several dosing scenarios. Such robustness indicates the utility of the model in estimating internal dose. The model also reproduced the observed dose-response patterns for mRNA and protein for CYP1A1, CYP1A2, and CYP1B1 after repeated dosing. Neither of the two dissociation constants for the Ah receptor bound to the CYP1B1 gene is negligible, supporting the assumption of multiple response elements for this gene. The poorer induction of CYP1B1 was predicted to be due to lower affinity of the dioxin responsive elements for binding the liganded Ah receptor, suggesting the involvement of other regulatory factors, and a shorter poly(A) tail on CYP1B1 mRNA, leading to a shorter lifetime. Saturation in the kinetics of protein synthesis was linked to the limited number of ribosomes that could bind to each message molecule, resulting in fewer ribosomes bound per message at higher doses. Predicted induction at low doses was found to vary widely with the assumptions used in the construction of a model. More detailed descriptions of biological processes might provide more reliable predictions of enzyme induction.

Contribution of planar (0-1 ortho) and nonplanar (2-4 ortho) fractions of Aroclor 1260 to the induction of altered hepatic foci in female Sprague-Dawley rats

The hepatic tumor promoting activity of the planar 0-1 ortho ( approximately 9.7% w/w) and the nonplanar 2-4 ortho ( approximately 90.3% w/w) fraction of the commercial PCB mixture Aroclor 1260 was studied using a medium-term two-stage initiation/promotion bioassay in female Sprague-Dawley rats. Fractionation was carried out on an activated charcoal column. The composition of the effluent from the column was tested by GC-ECD. The absence of planar compounds in the 2-4 ortho fraction was confirmed by GC-MS analysis. The dioxin-like toxic potency of the fractions was determined with the DR-CALUX assay. The animal experiment was started with the initiation procedure (diethylnitrosamine injection, 30 mg/kg body wt ip, 24 h after (2)/(3) hepatectomy), followed 6 weeks later by the promotion treatment, which consisted of a weekly subcutaneous injection during 20 weeks. Exposure groups (n = 10) received the following treatments (dose/kg body wt/week): Aroclor 1260 (10 mg), 0-1 ortho fraction (0.97 mg), 2-4 ortho fraction (1, 3, or 9 mg), a reconstituted 0-4 ortho fraction (9.97 mg), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB 153; 1 or 9 mg), 2,3,7,8-TCDD (1 microg; positive control) or corn oil (1 ml; vehicle control). One group did not receive a promotion treatment. All exposure groups exhibited a significantly increased volume fraction of the liver occupied by hepatic foci positive for the placental form of glutathione-S-transferase-p compared to the corn oil control, except for the groups treated with 0-1 ortho fraction and 1 mg PCB 153/kg body wt/week. Approximately 80% of the total tumor promoting capacity of the reconstituted 0-4 ortho fraction could be explained by the 2-4 ortho PCB fraction while the 0-1 ortho fraction had only a negligible contribution. These results suggest that the majority of the tumor promotion potential of PCB mixtures resides in the non-dioxin-like fraction, which is not taken into account in the toxic equivalency factor (TEF) approach for risk assessment of PCBs. This may result in an underestimation of the tumor promotion potential of environmental PCB mixtures.

Relative contributions of affinity and intrinsic efficacy to aryl hydrocarbon receptor ligand potency

Models of receptor action are valuable for describing properties of ligand-receptor interactions and thereby contribute to mechanism-based risk assessment of receptor-mediated toxic effects. In order to build such a model for the aryl hydrocarbon receptor (AHR), binding affinities and CYP1A induction potencies were measured in PLHC-1 cells and were used to determine intrinsic efficacies for 10 halogenated aromatic hydrocarbons (HAH): 2,3,7, 8-tetrachlorodibenzo-p-dioxin (TCDD), 2,3,7, 8-tetrachlorodibenzofuran (TCDF), and eight polychlorinated biphenyls (PCB). TCDD, TCDF, and non-ortho-substituted PCBs 77, 81, 126, and 169 behaved as full agonists and displayed high-intrinsic efficacy. In contrast, the mono- and di-ortho-substituted PCBs bound to the AHR but displayed lower or no intrinsic efficacy. PCB 156 was a full agonist, but with an intrinsic efficacy 10- to 50-fold lower than non-ortho-substituted PCBs. PCB 118 was a very weak partial agonist. PCBs 105 and 128 were shown to be competitive antagonists in this system. The model was then used to predict CYP1A induction by binary mixtures. These predictions were tested with binary mixtures of PCB 126, 128, or 156 with TCDD. Both PCB 156 (a low-intrinsic efficacy agonist) and PCB 128 (a competitive antagonist) inhibited the response to TCDD, while the response to TCDD and PCB126 was additive. These data support the following conclusions: 1) only 1-2% of the receptors in the cell need be occupied to achieve 50% of maximal CYP1A induction by one of the high-intrinsic efficacy agonists, demonstrating the existence of "spare" receptors in this system; 2) the insensitivity of fish to ortho-substituted PCBs is due to both reduced affinity and reduced intrinsic efficacy compared to non-ortho-substituted PCBs; 3) PCB congeners exhibit distinct structure-affinity and structure-efficacy relationships. Separation of AHR ligand action into the properties of affinity and intrinsic efficacy allows for improved prediction of the behavior of complex mixtures of ligands, as well as mechanistic comparisons across species and toxic endpoints.

Dose-response relationships for induction of CYP1A1 and CYP1A2 enzyme activity in liver, lung, and skin in female mice following subchronic exposure to polychlorinated biphenyls

The Toxic Equivalency Factor (TEF) method is used to estimate potential health risks associated with exposure to dioxin-like chemicals. The TEF method is a relative potency (REP) scheme that assumes dose additivity, that the chemicals produce the same response through the same mechanism, and that the REP of a chemical is equivalent for all responses. The present study estimates the REP for five PCBs with dioxin-like activity. Mice were exposed to either 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 2,3,3',4, 4'-pentachlorobiphenyl (105), 2,3',4,4',5-pentachlorobiphenyl (118), 3,3',4,4',5-pentachlorobiphenyl (126), or 2,3,3',4,4'-, 5-hexachlorobiphenyl (156), five days/week for 13 weeks by oral gavage in a corn oil vehicle. Three days after the last dose, animals were euthanized and the ethoxyresorufin-O-deethylase activity was determined in liver, lung, and skin. Acetanilide-4-hydroxylase activity was determined in liver. In addition, liver and skin disposition of the chemicals were determined. REPs were estimated using a statistical method previously described (DeVito et al., Toxicol. Appl. Pharmacol.147, 267-280, 1997). For any given compound, the REP generally varied by less than a factor of four across endpoints when calculated based on an administered dose. However, typically there was one response for every chemical in which the REP was different by an order of magnitude or more from the other responses. There was some evidence that the REPs may be dose-dependent. While, in general, these data support the use of a single point estimate of the REP, the issue of dose-dependency requires targeted investigation.

Liver tumor-promoting effect of beta-naphthoflavone, a strong CYP 1A1/2 inducer, and the relationship between CYP 1A1/2 induction and Cx32 decrease in its hepatocarcinogenesis in the rat

Interrelationships among induction of cytochrome P-450 (CYP) 1A1/2, decrease in connexin 32 (Cx32), and liver tumor-promoting activity by beta-naphthoflavone (BNF) in the promotion stage were examined in a 2-stage liver carcinogenesis model. A total of 20 male Fischer 344 rats were initiated with a single intraperitoneal injection of 150 mg/kg of diethylnitrosamine (DEN) or were given the saline vehicle alone. Starting 2 weeks later, they were fed a diet containing 2%, 1%, or 0% BNF for 6 weeks. All animals were subjected to a two-thirds partial hepatectomy at week 3 and were sacrificed at week 8. Absolute and relative liver weights were significantly increased in the DEN+BNF groups as compared to the DEN-alone group. Diffuse hepatocellular hypertrophy with cytoplasmic eosinophilia, sometimes accompanied by development of adenoma-like hepatic foci, was observed in the BNF-treated rats. Remarkable induction of cytochrome CYP 1A1/2 and significant increase in CYP 2E1 were noted in the DEN+BNF groups, and positive immunohistochemical staining for both was observed diffusely. The areas of Cx32-positive spots per hepatocyte in the centrilobular areas of livers of the BNF-treated rats were significantly decreased, but no changes were observed in periportal areas. The numbers and areas of foci positive for glutathione S-transferase placental form were increased in the BNF-treated groups. These results suggest that BNF is a liver tumor promoter that, unlike phenobarbital, does not induce CYP 2B1/2 isozymes, and there seems to be no direct relationship between CYP 1A1/2 induction and Cx32 reduction in BNF hepatocarcinogenesis.

Mechanistic aspects--the dioxin (aryl hydrocarbon) receptor

The Ah receptor mediates the toxicological responses of 2,3,7,8-TCDD and related compounds. Receptor-deficient animals were shown to be resistant to the toxic effects of dioxin, although there is also evidence for the existence of a receptor-independent pathway for dioxin-induced toxicity. In the cytosol the receptor is present in a non-activated ligand binding conformation. Association with Arnt in the nucleus turns the receptor complex into a ligand activated form. The physiological role of the receptor is not yet understood; however, the conservation of the receptor in a wide range of animal species (including humans) suggests a fundamental role in cellular physiology.

Toxicokinetics

The toxicokinetic determinants of dioxin and related chemicals depend on three major properties: lipophilicity, metabolism, and binding to CYP1A2 in the liver. The induction of CYP1A2 is partially under the control of the aryl hydrocarbon receptor (AhR). Lipophilicity increases with more chlorination and controls absorption and tissue partitioning. Metabolism is the rate-limiting step for elimination. Induction of CYP1A2 leads to hepatic sequestration of TCDD. Binding to this inducible hepatic protein results in non-linear dose-dependent tissue distribution: with increasing doses, the relative concentration in extra-hepatic tissues decreases while that in liver increases. The induction of this protein occurs in both animals and humans and results in an increase in the liver to fat ratio of these compounds. Humans have similar sensitivities to rodents for dioxin-like compounds when using tissue concentration (from in vitro studies), body burden, average lifetime serum lipid concentration, or lifetime area-under-the-curve concentration based on both low dose (biochemical) and high dose (cancer) driven endpoints. To reach the same tissue concentration in humans as rodents however, humans need a lower daily intake than rodents based on differences in pharmacokinetic behaviour. This clearly indicates that physiologically based pharmacokinetic models should be explored for the estimation of the daily intake of dioxin-like compounds in humans based on tissue dose levels or derivatives of those.

Risk ranges for various endpoints following exposure to 2,3,7,8-TCDD

Conducting a dose-response analysis for an environmental contaminant requires a careful evaluation of most of the available data focusing on both the magnitude of the effect and the possible ranges of dose-response shapes which fit the data. This paper calculates potency values (1% response exposures) for human and animal data on cancer, non-cancer and biological effect endpoints for TCDD and finds that a reasonable estimate for 1% excess cancer would be between 1 and 50 pg/kg/day. The paper also evaluates the adequacy of linear and non-linear models for fitting these data and concludes that the assumption of a threshold dose-response is not fully supported by these data.

Relative lipid content as the sole mechanistic determinant of the adipose tissue:blood partition coefficients of highly lipophilic organic chemicals

The adipose tissue:blood partition coefficient (PCat:b) refers to the ratio of chemical concentration or solubility in adipose tissue and blood. The solubility of a chemical in adipose tissue or whole blood is equal to the sum total of its solubility in lipid and water fractions of these matrices. For highly lipophilic organic chemicals (HLOCs, i.e., chemicals with log n-octanol:water partition coefficients (PCo:w) greater than four), their solubility in the water fractions of both tissue and blood is negligible, and therefore their solubility in lipid fractions of tissue and blood alone determines PCat:b. Since the numerical value representing chemical solubility in lipids is likely to be the same for both blood lipids and adipose tissue lipids, the PCat:b values should be hypothetically, equal to the ratio of lipid content of adipose tissue and blood. The objective of the present study was therefore to verify whether the PCat:bs of HLOCs (volatile organics, dioxins, PCBs, PBBs, DDT) are equal to the ratio of adipose tissue and blood lipid levels. The data on lipid content of rat and human blood and adipose tissues were obtained from the literature. The calculated tissue:blood lipid ratios were comparable to the human and rat PCat:b of volatile organic chemicals, dioxins, PCBs, PBBs and/or DDT obtained from the literature. These results then suggest that, regardless of the identity and PCo:w of HLOCs, their PCat:b is equal to the ratio of lipid in adipose tissues and blood.

Uncertainty in estimating exposure using a toxicokinetic model. The example of 2,3,7,8-tetrachlorodibenzo-p-dioxin

This paper deals with sources of uncertainty in the use of a minimal physiological toxicokinetic model to obtain dose estimates for a dose-response analysis of cancer in an occupational cohort. Toxicokinetic models make it possible to construct exposure parameters that are more closely related to the individual dose than traditional measures of exposures to toxic agents. However, the process introduces a wide array of sources of uncertainty. Selecting a model structure to describe the kinetics of a toxic agent implies necessarily making simplifications and assumptions that influence the range of applicability of the model. Once a model has been selected, the value of certain model parameters (constants) must be assigned, for example, from anthropometric data. The question then arises of how sensitive the model predictions are to variations in the values of these constants. Other model parameters, typically those describing the kinetics of the agent, are next estimated from actual data. There may be limitations in the data concerning, for example, sparseness (too few observations per subject) or missing values. The methods used for parameter estimation carry their own set of assumptions that need to be appropriate to the situation at hand. In summary, the dioxin example is used to characterize the sources of uncertainty at different levels, such as model structure, methods and data used for parameter estimation, estimation of occupational exposure, and imputation of missing values in exposure indices derived from the kinetic model.

Uncertainty in biomonitoring and kinetic modeling

Uncertainty in exposure assessment and uncertainty in kinetic models of early effects after exposure to a toxin are addressed in this paper. Sources of uncertainty in the determination of exposure of workers in chemical industry exposed to dioxins are exhibited and a simple kinetic model for biomonitor measurements of the concentrations from occupational exposure is derived. Model uncertainty, and uncertainty in the model parameters of physiologically-based pharmacokinetic models (PBPK models) are addressed when these models are used to estimate the effective dose in risk assessment. Uncertainty in the model parameters originating from the use of different statistical analysis methods is exhibited for Hill type nonlinear kinetics of enzyme induction mediated by a toxin.

Effects of TCDD on thyroid hormone homeostasis in the rat

A physiological dosimetric model was constructed to describe the effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) on circulating thyroid hormones in the rat and to test the hypothesis that these hormonal changes cause chronically elevated serum thyrotropin (thyroid stimulating hormone, TSH), which mediates growth promotion and may lead to thyroid tumors in TCDD-treated rats. The model included diffusion restricted distribution of TCDD among compartments for liver, kidney, white fat, slowly and rapidly perfused tissues, and the thyroxine-sensitive tissues brown fat, pituitary, and thyroid. Blood was distributed among major vessels and the capillary beds of the tissues. Metabolism of TCDD was limited to the liver. Secretion of 3,5,3'-triiodothyronine (T3) and thyroxine (3,5,3',5'-tetraiodothyronine, T4) from the thyroid was modeled as stimulated by circulating TSH, whose release from the pituitary was regulated by the hypothalamic peptides thyrotropin releasing hormone (activating) and somatostatin (inhibiting). Release of these peptides was represented as inhibited and activated, respectively, by circulating T4. Binding proteins for T3 and T4 and metabolism of the hormones by deiodination were included in thyroxine-sensitive tissues. Induction of hepatic UDP-glucuronosyltransferase-1*6 (UGT), the enzyme which glucuronidates T4, was modeled as induced by the complex formed between TCDD and the aryl hydrocarbon receptor. The computed extent of deiodination, primacy of the thyroid in generating T3 from T4, dependence of liver and kidney on locally produced T3, and export of T3 formed in the pituitary agreed with experimental observations. The model reproduced the observed decrease in circulating T4 and elevated serum TSH following chronic administration of TCDD. The altered levels were attributed to the increased clearance of T4 by the induced UGT and the consequent modification of feedback control of hormone releases. These results are consistent with the hypothesis of growth stimulation by elevated TSH, but measured values of this hormone in blood of rats vary over a large range, and the change induced by TCDD is often small. Measured UGT levels are less variable and the increase in this protein is much greater, suggesting that this response may be a more reliable biomarker for effects of TCDD on the thyroid.

Induction of altered hepatic foci by a mixture of dioxin-like compounds with and without 2,2',4,4',5,5'-hexachlorobiphenyl in female Sprague-Dawley rats

The hepatic tumor-promoting activity of a mixture of polyhalogenated aromatic hydrocarbons (PHAHs) was studied in a medium term two-stage initiation/promotion bioassay in female Sprague-Dawley rats. The PHAH mixture contained 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 1, 2,3,7,8-pentachlorodibenzo-p-dioxin (PeCDD), 2,3,4,7, 8-pentachlorodibenzofuran (PeCDF), 3,3',4,4',5-pentachlorobiphenyl (PCB 126), 2,3',4,4',5-pentachlorobiphenyl (PCB 118), 2,3,3',4,4', 5-hexachlorobiphenyl (PCB 156), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB 153) and covered >90% of the total toxic equivalents (TEQ) present in Baltic herring. To determine possible interactive effects of di-ortho-substituted PCBs, the PHAH mixture was tested with (PHAH+) and without (PHAH-) PCB 153. Rats were initiated by a diethylnitrosamine injection (30 mg/kg body wt i.p.) 24 h after a partial 23 hepatectomy. Six weeks after initiation, the PHAH mixtures were administered once a week by subcutaneous injections for 20 weeks. Treatment with the PHAH mixtures caused liver enlargement and an increased activity of the hepatic cytochrome P4501A1/2 and P4502B1/2. All PHAH exposure groups exhibited an increased occurrence of hepatic foci positive for the placental form of glutathione-S-transferase. In the PHAH-group dosed 1 microgram TEQ/kg body wt/week, the volume fraction of the liver occupied by foci was significantly lower compared to the TEQ equivalent dosed TCDD group (3.8 vs 8.7%). The volume fraction was significantly increased in the groups treated with 0.5, 1, or 2 micrograms TEQ/kg body wt/week of the PHAH+ mixture (4.5, 5.2, and 6.6%, respectively) compared to the corn oil group (2.0%), but to a lower extent than expected on basis of the TEQ doses. Overall, the TEQ-based administered dose overestimated the observed tumor-promoting effects of this PHAH mixture. The applicability of the toxic equivalency factor concept, the role of differences in toxicokinetic properties and interactive effects of PCB 153 on hepatic deposition of the dioxin-like congeners are discussed.

Toxicokinetics of an environmentally relevant mixture of dioxin-like PHAHs with or without a non-dioxin-like PCB in a semi-chronic exposure study in female Sprague Dawley rats

Female Sprague Dawley rats were treated subcutaneously for 20 weeks with an environmentally relevant mixture of dioxin-like PHAHs with (PHAH+) or without (PHAH-) 2,2',4,4',5,5'-hexachlorobiphenyl. The hepatic retention (% of given dose) of the various PHAH congeners differed considerably and in the following order: 2,3,4,7,8-pentachlorodibenzofuran (30.5-43.1%), 3,3',4,4',5-pentachlorobiphenyl (12.8-17.6%), 1,2,3,7,8-pentachlorodibenzo-p-dioxin (6.9-10.8%), 2,3,7,8-tetrachlorodibenzo-p-dioxin (3.2-4.5%), 2,3,3',4,4',5-hexachlorobiphenyl (1.0-1.7%), 2,2',4,4',5,5'-hexachlorobiphenyl (0.5-0.8%) and 2,3',4,4',5-pentachlorobiphenyl (0.2-0.4%). A decrease of the hepatic retention of 2,3,7,8-TCDD, 1,2,3,7,8-PeCDD and 2,3,4,7,8-PeCDF was found at increasing doses of the PHAH+ mixture. 2,2',4,4',5,5'-Hexachlorobiphenyl increased the hepatic retention (1.3-2 times) of all congeners in the PHAH+ group, compared to the TEQ equivalent dosed PHAH- group. No interactions were observed on the ethoxyresorufin-O-deethylase activity.

A pharmacodynamic analysis of TCDD-induced cytochrome P450 gene expression in multiple tissues: dose- and time-dependent effects

The ability of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD, dioxin) to alter gene expression and the demonstration that the induction of CYP1A2 is responsible for hepatic TCDD sequestration suggest that both pharmacokinetic and pharmacodynamic events must be incorporated for a quantitative description of TCDD disposition. In this paper, a biologically based pharmacodynamic (BBPD) model for TCDD-induced biochemical responses in multiple tissues was developed. The parameters responsible for tissue response were estimated simultaneously with a refined physiologically based pharmacokinetic (PBPK) model developed by Wang et al. (1997a), by using the time-dependent effects of TCDD on induced CYP1A1/CYP1A2 gene expression in multiple target tissues (liver, lungs, kidneys, and skin) of female Sprague-Dawley rats treated with 10 microgram TCDD/kg for 30 min, 1, 3, 8, or 24 h, or 7, 14, or 35 days. This refined BBPD model developed based on the time-course of TCDD-induced CYP1A1/CYP1A2 protein expression, and associated enzymatic activities well described the dose-dependent effects of TCDD on cytochrome P450 protein expression and associated enzyme activities in the multiple tissues of female Sprague-Dawley rats at 3 days following a single exposure to TCDD (0.01-30.0 micromgram TCDD/kg). This is the first BBPD model to quantitatively describe the time- and dose-dependent effects of TCDD on induced CYP1A1/CYP1A2 protein expression and associated enzyme activities in multiple target tissues for TCDD-induced biochemical responses.

Physiologically-based pharmacokinetic modeling of pyrene in the rat

The objective of the present study was to develop a physiologically-based model to simulate the oral and i.v. pharmacokinetics of pyrene in the rat. The physiologically-based pharmacokinetic (PBPK) model for pyrene consisted of the following tissue compartments: liver, lungs, adipose tissue, slowly perfused tissues, and richly perfused tissues interconnected with arterial and venous blood pools. The tissue:blood partition coefficients required for the pyrene PBPK model were estimated by equilibrium dialysis. Using perfusion-limited descriptions for tissue uptake and previously determined in vitro-derived hepatic metabolism rate constants (V(max) and K(m)), the PBPK model predicted a faster clearance of pyrene than that suggested by the experimental data. The biological basis of PBPK model then provided an opportunity to refine the estimate of V(max), and to explore and uncover additional mechanistic determinants of pyrene disposition in vivo. Accordingly, the in vitro V(max) had to be lowered by about a factor of 10 to adequately simulate experimental data on pyrene pharmacokinetics. Further, the model simulations could be matched with the experimental data on tissue concentrations of pyrene only with the considerations of (i) diffusion-limited uptake in slowly perfused tissues and adipose tissue, and (ii) binding to proteins in metabolizing tissues (lungs and liver). The present study successfully integrated the available data on oral and i.v. pharmacokinetics of pyrene using a physiological model framework, and identified several mechanistic data gaps that should be addressed by future research efforts.

Determination of parameters responsible for pharmacokinetic behavior of TCDD in female Sprague-Dawley rats

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is the most toxic member of a class of planar and halogenated chemicals. Improvements in exposure assessment of TCDD require scientific information on the distribution of TCDD in target tissues and cellular responses induced by TCDD. Since 1980, several physiologically based pharmacokinetic (PBPK) models for TCDD and related compounds have been reported. Some of these models incorporated the induction of a hepatic binding protein in response to interactions of TCDD, the Ah receptor, and DNA binding sites and described the TCDD disposition in a biological system for certain data sets. Due to the limitations of the available experimental data, different values for the same physical parameters of these models were obtained from the different studies. The inconsistencies of the parameter values limit the application of PBPK models to risk assessment. Therefore, further refinement of previous models is necessary. This paper develops an improved PBPK model to describe TCDD disposition in eight target tissues. The interaction of TCDD with the Ah receptor and with hepatic inducible CYP1A2 were also incorporated into the model. This model accurately described the time course distribution of TCDD following a single oral dose of 10 microg/kg, as well as the TCDD concentration on Day 3 after six different doses, 0.01, 0.1, 0.3, 1, 10, and 30 microg TCDD/kg, in target tissues. This study extends previous TCDD models by illustrating the validity and the limitation of the model and providing further confirmation of the potential PBPK model for us in optimal experimental design and extrapolation across doses and routes of exposure. In addition, this study demonstrated some critical issues in PBPK modeling.