Modélisation toxicocinétique du benzo(a)pyrène et 3-hydroxybenzo(a)pyrène pour l’interprétation des données de surveillance biologique de l’exposition chez les travailleurs

Many workers are exposed to polycyclic aromatic hydrocarbons (PAHs). Benzo(a) pyrene (BaP) is part of this group of pollutants. This substance has been classified as a known carcinogen in humans. To assess exposure to carcinogenic PAHs, several researchers have proposed using the measurement of 3-hydroxybenzo(a)pyrene (3-OHBaP) in the urine of exposed workers. In this project, two modeling approaches were developed and applied to enable a better understanding of the toxicokinetics of BaP and its biomarker of current interest, 3-OHBaP, to help interpret the results of biological monitoring. A multi-compartment toxicokinetic model was developed based on the data previously obtained in rats by our group of research. According to the model, BaP injected intravenously is rapidly distributed from blood to tissues (t½ ≈ 4 h), with a particular affinity for lungs and lipid components of tissues. Subsequently, BaP is distributed to the liver and the skin. Once in the liver, BaP is promptly metabolized and 3-OHBaP is formed with a half-life of about 3 h. Pulmonary biotransformation of BaP was also taken into account, but its contribution to the overall kinetics of BaP was considered negligible. Once formed, 3-OHBaP is distributed to various organs almost as fast as the parent compound (t½ ≈ 2 h). An accumulation of 3-OHBaP profile is present in the kidneys because of the difference between the uptake rate (t½ = 28 min) and the ouput rate (t½ = 4.5 h). Total clearance of 3-OHBaP from the blood stream is primarily governed by the rate of transfer of the bile to the gastrointestinal tract (t ½ ≈ 4 h). The multi-compartment toxicokinetic model was able to simulate an independent set of published 3-OHBaP urinary profiles. This toxicokinetic compartmental model has proved useful for the determination of the main biological features of the kinetics of BaP and 3-OHBaP. Thereafter, a physiological pharmacokinetic model (PBPK) reproducing the fate of BaP and 3-OHBaP rats was built. Organs (or tissues) represented as compartments were selected based on experimental data obtained in vivo in rats. Partition coefficients, coefficients of permeability, biotransformation rates, excretion parameters, and absorption fraction for different exposure routes were obtained directly from the profiles of BaP and 3-OHBaP in blood, various tissues and excreta. The values of these parameters were calculated by Monte Carlo procedures. Sensitivity analyses were then performed to ensure the stability of the model and to determine the most sensitive parameters. This modeling has identified the following features: 1) a high sensitivity of hepatic metabolism and elimination rates of BaP and 3-OHBaP; 2) a large distribution of BaP in the lungs compared to other tissues; 3) a considerable distribution of BaP in adipose tissues and liver; 4) a significant distribution of 3-OHBaP in the kidneys; 5) a diffusion-limited transfer of BaP in the lungs, 6) a diffusion-limited transfer of 3-OHBaP in lungs, adipose tissues and kidneys; and 7) a significant entero-hepatic recycling of 3-OHBaP. Following a series of analysis of goodness of fit, the probabilites that the model simulations reproduced the experimental data due to pure chance were always below 10%, for the four routes of exposure: intravenous, oral, dermal and respiratory. Subsequently, we have extrapolated the kinetic models from rats to humans in order reproduce the temporal evolution of 3-OHBaP biomarker of exposure in the urine of workers occupationally expose. Results of these models were then compared to simulations obtained with a single compartment toxicokinetic model to verify the comparative usefulness of simple and complex model. Both types of models have been constructed from blood, tissue, urinary and faecal profiles of BaP and 3-OHBaP in rats. These data were obtained in vivo by intravenous, subcutaneous, oral and respiratory exposure. The models were extrapolated to humans taking into account the essential biological determinants of kinetic differences between rats and humans. Results showed that inhalation was not the primary route of exposure for many workers studied. The values of air concentrations of BaP used to simulate the urinary excretion profiles were different from those measured in the air. Dermal exposure to BaP seemed to better predict the temporal patterns observed. Finally, the two types of modeling have been proved useful to reproduce and to interpret experimental data obtained in workers.

Keyphrases: Kinetic modeling, Biological monitoring, Reverse dosimetry