The importance of biological realism in dioxin risk assessment models

A mechanism-based mathematical model of aryl hydrocarbon receptor-mediated CYP1A induction in rats using beta-naphthoflavone as a tool compound

β-Naphthoflavone (BNF) is a synthetic flavone that selectively and potently induces CYP1A enzymes via aryl hydrocarbon receptor activation. Mechanism-based mathematical models of CYP1A enzyme induction were developed to predict the time course of enzyme induction and quantitatively evaluate the interrelationship between BNF plasma concentrations, hepatic CYP1A1 and CYP1A2 mRNA levels, and CYP1A enzyme activity in rats in vivo. Male Sprague-Dawley rats received a continuous intravenous infusion of vehicle or 1.5 or 6 mg · kg(-1) · h(-1) BNF for 6 h, with blood and liver sampling. Plasma BNF concentrations were determined by liquid chromatography-tandem mass spectrometry. Hepatic mRNA levels of CYP1A1 and CYP1A2 were determined by TaqMan. Ethoxyresorufin O-deethylation was used to measure the increase in CYP1A enzyme activity as a result of induction. The induction of hepatic CYP1A1/CYP1A2 mRNA and CYP1A activity occurred within 2 h after BNF administration. This caused a rapid increase in metabolic clearance of BNF, resulting in plasma concentrations declining during the infusion. Overall, the enzyme induction models developed in this study adequately captured the time course of BNF pharmacokinetics, CYP1A1/CYP1A2 mRNA levels, and increases in CYP1A enzyme activity data for both dose groups simultaneously. The model-predicted degradation half-life of CYP1A enzyme activity is comparable with previously reported values. The present results also confirm a previous in vitro finding that CYP1A1 is the predominant contributor to CYP1A induction. These physiologically based models provide a basis for predicting drug-induced toxicity in humans from in vitro and preclinical data and can be a valuable tool in drug development.

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.

No-threshold dose-response curves for nongenotoxic chemicals: findings and applications for risk assessment

We tested the hypothesis that no threshold exists when estradiol acts through the same mechanism as an active endogenous estrogen. A Michaelis-Menten (MM) equation accounting for response saturation, background effects, and endogenous estrogen level fit a turtle sex-reversal data set with no threshold and estimated the endogenous dose. Additionally, 31 diverse literature dose-response data sets were analyzed by adding a term for nonhormonal background; good fits were obtained but endogenous dose estimations were not significant due to low resolving power. No thresholds were observed. Data sets were plotted using a normalized MM equation; all 178 data points were accommodated on a single graph. Response rates from approximately 1% to >95% were well fit. The findings contradict the threshold assumption and low-dose safety. Calculating risk and assuming additivity of effects from multiple chemicals acting through the same mechanism rather than assuming a safe dose for nonthresholded curves is appropriate.

Use of sensitivity analysis to assess reliability of metabolic and physiological models

Because ethical considerations often preclude directly determining the human health effects of treatments or interventions by experimentation, such effects are estimated by extrapolating reactions predicted from animal experiments. Under such conditions, it must be demonstrated that the reliability of the extrapolated predictions is not excessively affected by inherent data limitations and other components of model specification. This is especially true of high-level models composed of ad hoc algebraic equations whose parameters do not correspond to specific physical properties or processes. Models based on independent experimental data restricting the numerical space of parameters that do represent actual physical properties can be represented at a more detailed level. Sensitivities of the computed trajectories to parameter variations permit more detailed attribution of uncertainties in the predictions to these low-level properties. S-systems, in which parameters are estimated empirically, and physiological models, whose parameters can be estimated accurately from independent data, are used to illustrate the applicability of trajectory sensitivity analysis to lower-level models.

Physiological modeling of butadiene disposition in mice and rats

The earliest physiological models of 1,3-butadiene disposition reproduced uptake of the gas from closed chambers but over-predicted steady-state circulating concentrations of the mutagenic intermediates 1,2-epoxybut-3-ene and 1,2:3,4-diepoxybutane. A preliminary model based on the observation of a transient complex between cytochrome P450 and microsomal epoxide hydrolase on the endoplasmic reticulum membrane reproduced the blood epoxide concentrations as well as the chamber uptake data. This model was enhanced by the addition of equations for the production and detoxication of 3,4-epoxybutane-1,2-diol in the liver, lungs, and kidneys. The model includes flow-restricted delivery of butadiene and its metabolites to compartments for lungs, liver, fat, kidneys, gastrointestinal tract, other rapidly perfused tissues, and other slowly perfused tissues. Blood was distributed among compartments for arterial, venous, and tissue capillary spaces. Channeling of the three bound epoxides to epoxide hydrolase and their release from the endoplasmic reticulum are competing processes in this model. Parameters were estimated to fit data for chamber uptake of butadiene and epoxybutene, steady-state blood concentrations of epoxybutene and diepoxybutane, and the fractions of the inhaled dose of butadiene that appears as various excreted metabolites. The optimal values of the apparent K(m)s of membrane-bound epoxides for epoxide hydrolase were only 5% of the values for the cytosolic substrate, consistent with the observation of a transient complex between epoxide hydrolase and the cytochrome P450 that produces the epoxide. This proximity effect corresponds to the notion that epoxides produced in situ have privileged access to epoxide hydrolase. The model also predicts considerable accumulation of epoxybutanediol, in agreement with the observation that most of the DNA adducts in animals exposed to butadiene arise from this metabolite.

Environmental contaminants and biochemical response in eel exposed to Po river water

Groups of eels (Anguilla anguilla) were exposed to Po river water under control conditions at "la Casella" Fluvial Hydrobiology Station for 30-day periods between November and May. At the end of the exposure period, fish were sacrificed and cytochrome P4501A was evaluated in liver microsomes using both catalytic (EROD) and semiquantitative immunodetection (ELISA) assays. At the same time, water samples were taken for chemical analyses of PCBs, PAHs and pesticides. Eel was chosen as bio-indicator on the bases of experiments on cytochrome P450 expression and activity under basal or induced conditions and in consideration of the ecological and economic relevance and ease of handling. Low levels of cytochrome P450 during Winter were followed by a step increase during Spring, associated with a substantial concentration increase of agrochemicals. A good correlation was found between measurements of cytochrome P4501A by EROD and ELISA in dose-effect experiments, however, ELISA showed a higher sensitivity. Immunochemical techniques may be used in addition to enzyme activity measurements both to detect lower levels of cytochrome P450 induction (where they proved more sensitive) and as quality control. These results suggest that measurement of cytochrome P4501A under controlled conditions, using eel as bio-indicator, can be a useful tool in monitoring Po river ecosystems.