Understanding the linked kinetics of benzo(a)pyrene and 3-hydroxybenzo(a)pyrene biomarker of exposure using physiologically-based pharmacokinetic modelling in rats

3-Hydroxybenzo(a)pyrene (3-OHBaP) in urine has been proposed as a biomarker of occupational exposure to polycyclic aromatic hydrocarbons. However, to reconstruct exposure doses in workers from biomarker measurements, a thorough knowledge of the kinetics of the benzo(a)pyrene (BaP) and 3-OHBaP given different routes of exposure is needed. A rat physiologically-based pharmacokinetic model of BaP and 3-OHBaP was built. Organs (tissues) represented as compartments were based on in vivo experimental data in rats. Tissue: blood partition coefficients, permeability coefficients, metabolism rates, excretion parameters, and absorption fractions and rates for different routes-of-entry were obtained directly from published in vivo time courses of BaP and 3-OHBaP in blood, various tissues and excreta of rats. The latter parameter values were best-fitted by least square procedures and Monte Carlo simulations. Sensitivity analyses were then carried out to ensure the stability of the model and the key parameters driving the overall modeled kinetics. This modeling pointed out critical determinants of the kinetics: (1) hepatic metabolism of BaP and 3-OHBaP elimination rate as the most sensitive parameters; (2) the strong partition of BaP in lungs compared to other tissues, followed by adipose tissues and liver; (3) the strong partition of 3-OHBaP in kidneys; (4) diffusion-limited tissue transfers of BaP in lungs and 3-OHBaP in lungs, adipose tissues and kidneys; (5) significant entero-hepatic recycling of 3-OHBaP. Very good fits to various sets of experimental data in rats from four different routes-of-entry (intravenous, oral, dermal and inhalation) were obtained with the model.