A physiologically based pharmacokinetic (PBPK) model for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) was combined with a five-compartment geometric model of hepatic zonation to predict both total and regional induction of CYP450 proteins within the liver. Three literature studies on TCDD pharmacokinetics and protein induction in female rats were analyzed. In simulating low-dose behavior for mRNA in whole liver and, particularly, in representing immunohistochemical observations, the five-compartment model was more successful than conventional homogeneous one-compartment liver models, The five-compartment liver model was used with the affinity of TCDD for the Ah receptor (AhR) held constant across all the liver (K-b = 0.2 nM). The presumed affinities of the AhR-TCDD complex for TCDD responsive elements in the CYP1A1 (K-d1) and CYP1A2 (K-d2) genes varied between adjacent compartments by a factor of 3, This parameterization leads to predicted 81-fold differences in affinities between the centrilobular and the periportal regions. The affinities used for AhR-TCDD complex binding to TCDD response elements for CYP1A2 in compartment 3 (the midzonal area) ranged from 0.08 to 1.0 nM in the three studies modeled. For CYP1A1 the corresponding dissociation constant in compartment 3 varied from 0.6 to 2.0 nM. In each compartment, the Hill coefficient for induction had to be 4 or greater to match the immunohistochemical results. This multicompartment liver model is consistent with data on protein and mRNA induction throughout the liver and on the regional distribution of these proteins, No previous model has incorporated regional variations in induction, The PBPK analysis based on the multicompartment liver model suggests that the low-dose behavior for hepatic CYP1A1/CYP1A2 induction by TCDD is highly nonlinear. (C) 1997 Academic Press.