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

Abstract

De nombreux travailleurs sont exposés aux hydrocarbures aromatiques polycycliques (HAP). Le benzo(a)pyrène (BaP) fait partie de ce groupe de polluants. Cette substance a été classée cancérogène reconnu chez l’humain. Pour évaluer l'exposition aux HAP cancérogènes, plusieurs chercheurs ont proposé d’utiliser la mesure du 3-hydroxybenzo(a)pyrène (3-OHBaP) dans l’urine des travailleurs exposés. Dans le cadre du présent projet, deux approches de modélisation ont été développées et appliquées pour permettre une meilleure compréhension de la toxicocinétique du BaP et son biomarqueur d’intérêt actuel, le 3-OHBaP, et pour aider à interpréter les résultats de surveillance biologique. Un modèle toxicocinétique à plusieurs compartiments a été développé sur la base des données préalablement obtenues sur le rat par notre groupe. Selon le modèle, le BaP injecté par voie intraveineuse est rapidement distribué du sang vers les tissus (t½ ≈ 4 h), avec une affinité particulière pour les poumons et les composantes lipidiques des tissus. Le BaP est ensuite distribué vers la peau et le foie. Au foie, le BaP est promptement métabolisé et le 3-OHBaP est formé avec une demi-vie de ≈ 3 h. Le métabolisme pulmonaire du BaP a également été pris en compte, mais sa contribution à la cinétique globale du BaP a été jugée négligeable. Une fois formé, le 3-OHBaP est distribué vers les différents organes presque aussi rapidement que la molécule mère (t½ ≈ 2 h). Le profil temporel du 3-OHBaP dans le rein montre une accumulation transitoire en raison de la différence observée entre le taux d’entrée (t½ = 28 min) et le taux de sortie (t½ = 4,5 h). La clairance totale de 3-OHBaP du corps est principalement gouvernée par le taux de transfert de la bile vers le tractus gastro-intestinal (t½ ≈ 4 h). Le modèle toxicocinétique à plusieurs compartiments a réussi à simuler un ensemble indépendant de profils urinaires publiés sur le 3-OHBaP. Ce modèle toxicocinétique à compartiments s'est avéré utile pour la determination des facteurs biologiques déterminants de la cinétique du BaP et du 3-OHBaP. Par la suite, un modèle pharmacocinétique à base physiologique (PCBP) reproduisant le devenir du BaP et du 3-OHBaP chez le rat a été construit. Les organes (ou tissus) représentés comme des compartiments ont été choisis en fonction de données expérimentales obtenues in vivo chez le rat. Les coefficients de partition, les coefficients de perméabilité, les taux de métabolisation, les paramètres d'excrétion, les fractions absorbées et les taux d'absorption pour différentes voies d’exposition ont été obtenus directement à partir des profils sanguins, tissulaires, urinaires et fécaux du BaP et du 3-OHBaP. Les valeurs de ces derniers paramètres ont été calculées par des procédures Monte-Carlo. Des analyses de sensibilité ont ensuite été réalisées pour s’assurer de la stabilité du modèle et pour établir les paramètres les plus sensibles de la cinétique globale. Cette modélisation a permis d’identifier les facteurs déterminants de la cinétique: 1) la sensibilité élevée des paramètres de la métabolisation hépatique du BaP et du 3-OHBaP ainsi que du taux d'élimination; 2) la forte distribution du BaP dans les poumons par rapport à d'autres tissus; 3) la distribution considérable du BaP dans les tissus adipeux et le foie; 4) la forte distribution du 3-OHBaP dans les reins; 5) le transfert limité du BaP par la diffusion tissulaire dans les poumons; 6) le transfert limité du 3-OHBaP par la diffusion tissulaire dans les poumons, les tissus adipeux et les reins; 7) la recirculation entéro-hépatique significative du 3-OHBaP. Suite à des analyses de qualité des ajustements des équations du modèle aux données observées, les probabilités que les simulations reproduisent les données expérimentales par pur hasard se sont avérées toujours inférieures à 10% pour les quatre voies d’exposition : intraveineuse, orale, cutanée et respiratoire. Nous avons extrapolé les modèles cinétiques du rat à l’humain afin de se doter d’un outil permettant de reconstituer les doses absorbées chez des travailleurs exposés dans diverses industries à partir de mesures de l'évolution temporelle du 3-OHBaP dans leur urine. Les résultats de ces modélisations ont ensuite été comparés à ceux de simulations obtenues avec un modèle toxicocinétique à compartiment unique pour vérifier l’utilité comparative d’un modèle simple et complexe. Les deux types de modèle ont ainsi été construits à partir de profils sanguins, tissulaires, urinaires et fécaux du BaP et du 3-OHBaP sur des rats exposés. Ces données ont été obtenues in vivo par voie intraveineuse, cutanée, respiratoire et orale. Ensuite, les modèles ont été extrapolés à l’humain en tenant compte des déterminants biologiques essentiels des différences cinétiques entre le rat et l’humain. Les résultats ont montré que l'inhalation n'était pas la principale voie d'exposition pour plusieurs travailleurs étudiés. Les valeurs de concentrations de BaP dans l’air utilisées afin de simuler les profils d’excrétion urinaire chez les travailleurs étaient différentes des valeurs de concentrations de BaP mesurées dans l’air. Une exposition au BaP par voie cutanée semblait mieux prédire les profils temporels observés. Finalement, les deux types de modélisation se sont avérés utiles pour reproduire et pour interpréter les données disponibles chez des 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.