基于贝叶斯的苯并(α)芘PBPK模型优化与健康风险评估应用  

Optimization of benzo (α) pyrene PBPK model based on bayes and its application in health risk assessment

作　　者：孟祥晖 王宏洋[2] 孙宇巍 张明慧[2] 朱光灿[1,3] 沈亚琴[2] 柳慧芳 MENG Xiang-hui;WANG Hong-yang;SUN Yu-wei;ZHANG Ming-hui;ZHU Guang-can;SHEN Ya-qin;LIU Hui-fang(Southeast University School of Energy and Environment,Nanjing 210096,China;State Key Laboratory of Environment Criteria and Risk Assessment,Chinese Research Academy of Environmental Sciences,Beijing 100012,China;Key Laboratory of Environmental Medical Engineering,Ministry of Education,Southeast University,Nanjing 210096,China)

机构地区：[1]东南大学能源与环境学院,江苏南京210096 [2]中国环境科学研究院,环境基准与风险评估国家重点实验室,北京100012 [3]东南大学,环境医学工程教育部重点实验室,江苏南京210096

出　　处：《中国环境科学》2022年第5期2370-2378,共9页China Environmental Science

基　　金：国家重点研发计划项目(2021YFE0106600)。

摘　　要：应用基于生理的药代动力学(PBPK)模型预测苯并(α)芘(BaP)暴露的人体内部剂量,基于贝叶斯的马尔科夫链蒙特卡洛模拟(MCMC)方法对模型参数进行校准和优化,最后运用已优化的模型对BaP内暴露基准值进行推导.研究发现,基于贝叶斯的MCMC方法对模型后验参数校准后,模型精度明显提高,两个数据集验证结果显示残差平方和分别降低了72%和94%.PBPK模型以BaP和子代谢物3-羟基苯并(α)芘(3-OHBaP)的体内动力学过程为结构基础,模拟BaP体内浓度分布大小为脂肪>肾脏>皮肤>缓慢灌注组织>快速灌注组织>静脉血>肝脏;3-OHBaP体内浓度分布大小为肾脏>快速灌注组织>脂肪>肺>静脉血>缓慢灌注组织>肝脏>皮肤.敏感性分析显示,快速灌注组织-血分配系数对模型输出影响最大,灵敏度系数超过了200%;排泄系数影响最小,只有肾小球过滤率K_(BR)的灵敏度系数超过了1%.以美国国家环境保护局推荐的参考浓度2.0×10^(-6)mg/m^(3)为外暴露安全基准值,基于PBPK模型推导了职业暴露的BaP生物监测当量(BE),结果显示BE值为0.405pmol/mol肌酐(尿液3-OHBaP平均浓度),为基于人体内暴露剂量水平进行定量健康风险评估奠定了基础.The physiologically based pharmacokinetic(PBPK)model was used to predict the human internal dose of benzo(α)pyrene(BaP)exposure.The parameters of the model were optimized based on bayes statistics and Markov Chain Monte Carlo simulation(MCMC),and the optimized model was adopted to derive the reference value regarding the internal dose of BaP.It was found that the accuracy of the model was significantly improved after calibrating the posterior parameters by Monte Carlo simulation,and the validation results of two datasets showed that the sums of squared residuals were reduced by 72%and 94%.The PBPK model was based on the pharmacokinetic of BaP and its metabolite 3-hydroxybenzo(α)pyrene(3-OHBaP).The internal concentration distribution of BaP followed the order of fat>kidney>skin>slowly perfused tissue>richly perfused tissue>venous blood>liver,while that of 3-OHBaP was in the order of kidney>richly perfused tissue>fat>lung>venous blood>slowly perfused tissue>liver>skin.Besides,the sensitivity analysis indicated that the rich perfused tissue-blood distribution coefficients showed the strongest influence on the model output,which sensitivity coefficient exceeded 200%.While the coefficients related to excretion showed the weakest influence,and only sensitivity coefficient of glomerular filtration rate K_(BR) exceeded 1%.According to the reference concentration 2.0×10^(-6)mg/m^(3) recommended by U.S.Environmental Protection Agency,the biomonitoring equivalent of BaP was derived based on the optimized PBPK model.The results showed that the reference value for the occupational populations was 0.405pmol/mol creatinine(i.e.,the mean concentration of 3-OHBaP in urine),which lays a foundation for the quantitative health risk assessment based on the human internal dose.

关 键 词：苯并(α)芘 贝叶斯统计 马尔科夫链蒙特卡洛方法(MCMC) 基于生理的药代动力学(PBPK)模型 3-OHBaP 生物监测当量 

分 类 号：X592[环境科学与工程—环境工程]