劉遠(yuǎn)，王中鈺，陳景文，喬顯亮

大連理工大學(xué)環(huán)境學(xué)院，工業(yè)生態(tài)與環(huán)境工程教育部重點(diǎn)實(shí)驗(yàn)室，大連116024

三氯乙烯(TCE)作為脫脂劑和清洗劑被廣泛應(yīng)用于五金、電鍍和電子等行業(yè)[1]。2012年TCE在中國的總產(chǎn)量約57萬t[2]。由于被大量長期使用，TCE成為空氣、地下水和食品中常檢出的污染物[3]。由于具有很強(qiáng)的揮發(fā)性，在使用脫脂劑或清洗劑的電子、五金和電鍍等工作場所的空氣中，TCE的檢出濃度可以達(dá)到18～683 mg·m-3[4-5]。TCE的職業(yè)暴露會產(chǎn)生過敏綜合癥[6]和致癌風(fēng)險[7]。2012年，TCE被美國環(huán)保局(US EPA)和國際癌癥研究機(jī)構(gòu)列為1類致癌物[7-8]。研究表明，TCE在體內(nèi)的代謝產(chǎn)物與肺部、肝臟和腎臟腫瘤均存在相關(guān)性[1,9]。因此，有必要評估TCE及代謝產(chǎn)物在職業(yè)暴露中對人體健康的致癌風(fēng)險。

傳統(tǒng)的人體暴露風(fēng)險評價一般基于環(huán)境中污染物的外暴露水平[10-11]。污染物進(jìn)入人體后會在不同組織和器官中進(jìn)行分配，在肝臟等器官中還會存在代謝過程[12]。基于生理學(xué)的藥代動力學(xué)(physiologically based pharmacokinetic, PBPK)模型可以用于模擬預(yù)測外源物質(zhì)在體內(nèi)的動態(tài)分布和代謝過程，已被用于芘[13]、氯仿[14]和尼古丁[15]等物質(zhì)在人體內(nèi)的分布、代謝和健康風(fēng)險研究。本研究采用吸附管采樣-熱脫附/氣相色譜-質(zhì)譜法分析了大連市某企業(yè)車間生產(chǎn)工況下空氣中TCE濃度，基于PBPK模型預(yù)測了呼吸暴露途徑下人體組織中TCE的動態(tài)變化、代謝產(chǎn)物生成情況和致癌風(fēng)險，并對模型結(jié)果進(jìn)行了敏感性和不確定性分析。本研究的發(fā)現(xiàn)可以為TCE的職業(yè)暴露致癌風(fēng)險提供參考。

1 材料與方法(Materials and methods)

1.1 車間空氣中TCE樣品采集和分析

基于本課題組建立的工廠車間空氣中揮發(fā)性有機(jī)物的測定方法[16]，于2017年8月份使用吸附管(Tenax, 6 mm×4 mm×200 mm)采樣-熱脫附/氣相色譜-質(zhì)譜法(Markes Series 2 Unity TD/Agilent 7890A-5975C GCMS)分析了大連市某企業(yè)車間生產(chǎn)工況下空氣中TCE濃度，將該企業(yè)辦公區(qū)域設(shè)置為對照，車間和辦公區(qū)分別設(shè)置3次重復(fù)。

1.2 TCE在人體內(nèi)的代謝產(chǎn)物

根據(jù)前人研究，TCE在人體內(nèi)主要被細(xì)胞色素P450酶氧化或者被谷胱甘肽(GSH)轉(zhuǎn)移酶催化代謝[7,12,17]。P450酶氧化主要在肝臟中進(jìn)行。進(jìn)入肝臟的TCE被其氧化為三氯乙醛(CHL)；CHL被醇脫氫酶轉(zhuǎn)化為三氯乙醇(TCOH)。TCOH既可被糖酯化生成相應(yīng)的葡糖苷酸(TCOG)，也可被醛脫氫酶氧化為三氯乙酸(TCA)，還可被還原為二氯乙酸(DCA)。TCOG可通過肝腸循環(huán)轉(zhuǎn)化為TCOH，TCA可脫氯轉(zhuǎn)化為DCA。GSH酶催化途徑大致為肝臟中大量谷胱甘肽S-轉(zhuǎn)移酶與TCE發(fā)生親核取代生成S-二氯乙烯基-L-谷胱甘肽(DCVG)。DCVG在腎臟中轉(zhuǎn)化為S-二氯乙烯基-L-半胱氨酸(DCVC)，最后以N-乙酰-S-(1,2-二氯乙烯)-L-氨酸(NacDCVC)排出體外。GSH酶催化途徑產(chǎn)生的DCVC對人類腎癌具有重要貢獻(xiàn)[7]。CHL、TCA和DCA也具有毒理學(xué)效應(yīng)[1]。

1.3 PBPK模型結(jié)構(gòu)和計(jì)算方法

TCE及代謝產(chǎn)物的人體PBPK模型由包含血?dú)饨粨Q系統(tǒng)、支氣管組織(TBT)、脂肪組織(Fat)、充分灌注室(RPT)、非充分灌注室(PPT)、腸組織(Gut)、肝臟組織(Liver)、動脈血(AB)和靜脈血(VB)的生理學(xué)腔室組成，其中腎臟包含于RPT(圖1)。模型初始參數(shù)中TCE血液/空氣分配系數(shù)(Pb)為9.2，其變異系數(shù)為50%[9]。體重(BW)、心輸出血流量(QC)和肺換氣流量(QP)分別為70 kg、315 L·h-1和436 L·h-1，其變異系數(shù)分別為79%、40%和30%[9]。與生理學(xué)腔室有關(guān)的其他參數(shù)及變異情況見表1～3。對于TCE進(jìn)入人體的吸收、分布、代謝和排泄過程簡單表述如下，詳細(xì)公式見Clewell等[17]的報道。

吸收過程：當(dāng)空氣中TCE濃度為CInh(mg·L-1)時，經(jīng)血?dú)饨粨Q系統(tǒng)后TCE在動脈血中的濃度CAB(mg·L-1)如式(1)所示，

(1)

式中，Pb表示TCE的血液/空氣分配系數(shù)。

分布過程：TCE在靜脈血和各組織靜脈血中的濃度CVB(mg·L-1)和CVTissue(mg·L-1)分別如式(2)和(3)所示，其中ATissue(mg)為各組織內(nèi)目標(biāo)化合物質(zhì)量，

(2)

(3)

代謝過程：以TCE在代謝組織中的轉(zhuǎn)化速率RAMTissue(mg·h-1)和腎臟中DCVC的轉(zhuǎn)化為例，其分別如式(4)和(5)所示，

(4)

(5)

式中，MW表示化合物分子量。

排泄過程：以NacDCVC的排泄為例，其如式(6)所示，

(6)

圖1 三氯乙烯(TCE)及代謝產(chǎn)物的人體基于生理學(xué)的藥代動力學(xué)(PBPK)模型Fig. 1 Human physiologically based pharmacokinetic (PBPK) model for trichloroethylene (TCE) and its metabolites

根據(jù)質(zhì)量守恒原理，將描述TCE在人體各組織中藥代動力學(xué)過程的微分方程聯(lián)立方程組，車間空氣中TCE濃度設(shè)為CInh(mg·L-1)，各組織中TCE初始濃度均為零，以1×10-5h為步長，職業(yè)工人工作時間為暴露時長(8 h)，模擬24 h。模型參數(shù)初始值如表1～3所示，使用Python(3.6)科學(xué)計(jì)算庫SciPy中integrate模塊下的odeint函數(shù)對方程組進(jìn)行計(jì)算求解，結(jié)果得到人體各組織中TCE及代謝產(chǎn)物濃度。其中代謝產(chǎn)物TCA、DCA、TCOH和TCOG對應(yīng)的體內(nèi)表觀分布體積所占體重的比例分別為0.1、0.1、0.65和0.65，其變異系數(shù)均為79%[9]。

1.4 外暴露和內(nèi)暴露方法評價職業(yè)暴露風(fēng)險

TCE職業(yè)工人外暴露致癌風(fēng)險基于US EPA推薦的吸入風(fēng)險評估模型和職業(yè)暴露參數(shù)[10,18]。致癌風(fēng)險(Risk)如式(7)所示，

Risk=IUR×EC

(7)

式中，IUR為TCE的吸入單位風(fēng)險4.10×10-6(μg·m-3)-1[19]，EC為職業(yè)工人在工作車間內(nèi)的暴露濃度(μg·m-3)，如式(8)所示：

(8)

式中，CA為空氣中TCE濃度(μg·m-3)，ET為暴露時間(8 h·d-1)，EF為暴露頻率(250 d·a-1)，ED為持續(xù)暴露階段(20 a)，AT為平均壽命(70 a)。

TCE職業(yè)工人的內(nèi)暴露致癌風(fēng)險基于PBPK模型預(yù)測的最大內(nèi)暴露濃度與相應(yīng)的內(nèi)暴露單位風(fēng)險的乘積[20,21]。TCE及代謝產(chǎn)物在組織中的內(nèi)暴露單位風(fēng)險計(jì)算如下，以IURTCE對應(yīng)的空氣濃度(2.44×10-1μg·m-3)[22]做為CInh，模型計(jì)算過程如1.3節(jié)所述，得到安全閾值下TCE及代謝產(chǎn)物對應(yīng)的最大內(nèi)暴露濃度，基于風(fēng)險閾值1.0×10-6將最大內(nèi)暴露濃度轉(zhuǎn)換為相應(yīng)的內(nèi)暴露單位風(fēng)險[23]。以肝臟中TCE和組織中TCE代謝產(chǎn)物DCVC為例，最大內(nèi)暴露濃度分別為1.97 μg·m-3和8.27×10-4μg·m-3；其內(nèi)暴露單位風(fēng)險分別為5.08×10-7(μg·m-3)-1和1.21×10-3(μg·m-3)-1。

由于個體間PBPK模型參數(shù)存在變異性，模型生理參數(shù)值和代謝參數(shù)值分別服從正態(tài)分布和對數(shù)正態(tài)分布[9,17]，其范圍用變異系數(shù)(CV)表示，CV = 100%×標(biāo)準(zhǔn)偏差/均值，如表1～3所示?；赑BPK模型結(jié)合Monte Carlo方法可以體現(xiàn)模型參數(shù)變化對預(yù)測結(jié)果的影響，從而可以提高模型預(yù)測結(jié)果的準(zhǔn)確性[24]。預(yù)測TCE及代謝產(chǎn)物的內(nèi)暴露濃度時，利用Monte Carlo方法從PBPK模型參數(shù)的概率分布范圍(均值±標(biāo)準(zhǔn)偏差)內(nèi)隨機(jī)抽取10萬次并代入模型，計(jì)算過程如1.3節(jié)所述，得到TCE及代謝產(chǎn)物對應(yīng)的最大內(nèi)暴露濃度。由于PBPK模型參數(shù)存在較大的變異性[9,17]，為了降低參數(shù)極值對預(yù)測結(jié)果的影響，選擇模型預(yù)測的最大內(nèi)暴露濃度5%～95%分位數(shù)的平均值來計(jì)算致癌風(fēng)險。

1.5 模型參數(shù)的敏感性分析及不確定性分析

敏感性系數(shù)(SC)常用于判斷PBPK模型預(yù)測結(jié)果對模型參數(shù)的敏感性[25-26]。SC計(jì)算為模型參數(shù)初始值θ變化±5%時模型預(yù)測結(jié)果f(θ)的相對變化程度，如式(9)所示，

(9)

其可分為低敏感(0 < |SC| ≤ 0.5)、中敏感(0.5 < |SC| ≤ 1)和高敏感(|SC| > 1)[26]。Monte Carlo方法常用來評價PBPK模型預(yù)測結(jié)果的不確定性[24]。以PBPK模型預(yù)測肝臟中TCE和組織中TCE代謝產(chǎn)物DCVC的內(nèi)暴露濃度為例，結(jié)合Monte Carlo方法在模型參數(shù)的概率分布范圍(均值±標(biāo)準(zhǔn)偏差)內(nèi)隨機(jī)抽取后代入模型，計(jì)算過程如1.3節(jié)所述，得到不同參數(shù)值下TCE及代謝產(chǎn)物最大內(nèi)暴露濃度的不確定性范圍[17]。

2 結(jié)果與討論(Results and discussion)

2.1 職業(yè)工人體內(nèi)TCE及代謝產(chǎn)物的動力學(xué)模擬

大連市某企業(yè)生產(chǎn)工況下車間空氣中TCE濃度為(39.2±24.4) μg·m-3，辦公區(qū)對照點(diǎn)的濃度為(0.08±0.06) μg·m-3?？諝庵蠺CE濃度在車間區(qū)高于辦公區(qū)與其作為清洗劑的使用有關(guān)。當(dāng)職業(yè)工人在車間內(nèi)工作8 h后離開，基于PBPK模型模擬24 h內(nèi)TCE經(jīng)呼吸暴露在人體內(nèi)的分布和代謝過程。模擬結(jié)果表明TCE及代謝產(chǎn)物在人體組織中的濃度隨暴露時間增長逐漸升高。預(yù)測的組織中TCE最大濃度呈現(xiàn)出脂肪>腸>充分灌注室>支氣管>非充分灌注室>肝臟>靜脈血>動脈血的趨勢(圖2a)。相對于其他組織，TCE易于在脂肪組織中富集，可能與TCE的脂肪/血液分配系數(shù)PFat= 73高有關(guān)[17]。前人發(fā)現(xiàn)芘(PFat= 11.8)[13]和甲基叔丁基醚(PFat= 6.5)[27]等也具有易在脂肪內(nèi)富集的特征。

模擬結(jié)果表明，TCE代謝產(chǎn)物在人體組織中的最大濃度呈現(xiàn)出CTCA> CTCOH> CTCOG> CDCA> CCHL> CDCVC的趨勢(圖2b)。代謝產(chǎn)物濃度分布與其相應(yīng)的生成途徑有一定關(guān)聯(lián)，表現(xiàn)為P450酶氧化途徑產(chǎn)物濃度大于GSH酶催化途徑產(chǎn)物濃度。實(shí)驗(yàn)表明TCE在大鼠體內(nèi)的代謝也表現(xiàn)出類似的趨勢[28]。TCE在嚙齒動物和人體內(nèi)具有相似的代謝過程，這說明該模型可以模擬TCE在人體內(nèi)的分布和代謝過程。

2.2 外暴露和內(nèi)暴露方法評估TCE職業(yè)暴露風(fēng)險

當(dāng)該車間職業(yè)工人在(39.2±24.4) μg·m-3TCE濃度水平下，暴露8 h·d-1，連續(xù)暴露20年時，采用外暴露方法評估致癌風(fēng)險均值為1.31×10-5，高于可接受的風(fēng)險閾值1.0×10-6。污染物進(jìn)入人體后會在不同組織和器官中進(jìn)行分配，在肝臟等器官中還會存在代謝過程[12]。通過PBPK模型計(jì)算進(jìn)入靶組織的具有致癌活性物質(zhì)的劑量，可以基于內(nèi)暴露劑量來評估致癌風(fēng)險[11]。在本研究監(jiān)測的濃度水平下，通過PBPK模型預(yù)測TCE內(nèi)暴露濃度計(jì)算的致癌風(fēng)險與基于外暴露濃度計(jì)算的風(fēng)險基本一致(表4)。前人研究表明，TCE的主要代謝產(chǎn)物如肝臟中的TCA和腎臟中的DCVC也具有致癌風(fēng)險[1,7]。本研究采用PBPK模型預(yù)測的TCA和DCVC內(nèi)暴露濃度分別計(jì)算獲得的致癌風(fēng)險值是基于TCE外暴露濃度風(fēng)險值的1.73和1.17倍(表4)。考慮到本次調(diào)查車間空氣TCE濃度較低(μg·m-3)，本研究選取文獻(xiàn)報道的TCE高暴露水平(119 mg·m-3)來分析暴露水平對風(fēng)險評價結(jié)果的影響[5]。如表4所示，高濃度水平下，基于TCE、TCA和DCVC內(nèi)暴露的致癌風(fēng)險是外暴露風(fēng)險的1.05～1.95倍。以上分析說明TCE的暴露水平對內(nèi)暴露方法與外暴露方法評價結(jié)果存在一定影響，暴露水平越高，內(nèi)、外暴露的風(fēng)險評價結(jié)果差異越大，這與前人的研究發(fā)現(xiàn)一致[20,23]。綜上，評價TCE的人體健康風(fēng)險非常需要考慮其在人體內(nèi)的分布和代謝過程。

表1 與組織類型有關(guān)的人體PBPK模型生理參數(shù)值aTable 1 Physiological parameter values of human PBPK model related to types of tissue a

注：a參數(shù)值來自Clewell等[17]的報道，體積分?jǐn)?shù)值和血流量分?jǐn)?shù)值需分別乘以體重0.75和心輸出血量轉(zhuǎn)換。b變異系數(shù)值來自Hack等[9]的報道。c包含于充分關(guān)注組織參數(shù)值中。

Note:aParameter values from Clewell et al[17], fractional volumes and fractional blood flows scaled by 0.75 power of the body weight and cardiac output, respectively.bThe CV values from Hack et al[9].cValue belongs to QRPT.

表2 肝臟和腸組織中的人體PBPK模型代謝參數(shù)值aTable 2 Metabolic parameter values of human PBPK model in liver and gut tissue a

注：a參數(shù)值來自Clewell等[17]的報道，代謝速率需乘以體重0.75轉(zhuǎn)換，動力學(xué)常數(shù)和排泄常數(shù)需乘以體重-0.25轉(zhuǎn)換。變異系數(shù)值除b為79%和c為490%，其余均為200%，變異系數(shù)值都來自Hack等[9]的報道。d無參數(shù)。TCA表示三氯乙酸，DCA表示二氯乙酸，TCOH表示三氯乙醇，TCOG表示三氯乙醇葡糖苷酸。PO表示比例常數(shù)，ku表示排泄常數(shù)，kehr表示膽汁排泄常數(shù)。

Note:aParameter values from Clewell et al[17]; values of Vmaxscaled by 0.75 power of the body weight; values of kinetics and excretion scaled by -0.25 power of the body weight. The CV values were 200% except forbof 79% andcof 490%, and the CV values from Hack et al[9].dNo parameters. TCA, trichloroacetate; DCA, dichloroacetate; TCOH, trichloroethanol; TCOG, trichloroethanol glucuronide; PO, percent of TCE; ku, urinary excretion; kehr, biliary excretion.

表3 腎臟和支氣管組織中的人體PBPK模型代謝參數(shù)值aTable 3 Metabolic parameter values of human PBPK model in kidney and trachea-bronchial tissue a

注：a參數(shù)值來自Clewell等[17]的報道，代謝速率需乘以體重0.75轉(zhuǎn)換，動力學(xué)常數(shù)需乘以體重-0.25轉(zhuǎn)換，變異系數(shù)值均為200%來自Hack等[9]的報道。b無參數(shù)。DCVC表示S-二氯乙烯基-L-半胱氨酸，CHL表示三氯乙醛。

Note:aParameter values from Clewell et al[17]; values of Vmaxand kinetics scaled by 0.75 power of and -0.25 power of the body weight, respectively. All the CV values were 200% and from Hack et al[9].bNo parameters. DCVC, S-(1,2-dichlorovinyl)-L-cysteine; CHL, chloral.

圖2 人體PBPK模型預(yù)測組織中TCE(a)及代謝產(chǎn)物(b)的濃度分布注：VB表示靜脈血，AB表示動脈血。Fig. 2 Levels of TCE (a) and its metabolites levels (b) in tissues predicted by human PBPK modelNote: VB, venous blood; AB, arterial blood.

2.3 模型參數(shù)對預(yù)測TCE及代謝產(chǎn)物濃度的影響

為了考察不同模型參數(shù)對PBPK模型輸出結(jié)果的影響[25]，本研究進(jìn)行了模型敏感性分析。由于DCVC與腎癌存在關(guān)系，且DCVC的生成主要在肝臟和腎臟中進(jìn)行[7]，以模型預(yù)測的體內(nèi)DCVC最大濃度、肝臟和充分灌注室(包括腎臟)中TCE最大濃度為輸出結(jié)果考察敏感性。根據(jù)前人的研究[26]，基于敏感性系數(shù)(SC)可將三氯乙烯的人體PBPK模型參數(shù)分為低敏感(0 < |SC|≤ 0.5)、中敏感(0.5 < |SC| ≤ 1)和高敏感(|SC| > 1)。模擬條件為職業(yè)工人在該企業(yè)車間空氣TCE監(jiān)測水平下工作8 h，模擬24 h。本研究得到具有中敏感性和高敏感性的參數(shù)如表5所示。模型預(yù)測結(jié)果對心輸出血流量(QC)和充分灌注室血流量(QRPT)具有高敏感性，這與血流量對于TCE及代謝產(chǎn)物在不同組織中的輸送和分配具有重要影響有關(guān)[29]。模型預(yù)測結(jié)果對DCVC生成常數(shù)(KF)、腎臟體積(VKidney)、肝臟體積(VLiver)以及肝臟/血液分配系數(shù)(PLiver)等具有中敏感性，這可能與TCE由血液進(jìn)入主要代謝器官肝臟后代謝為DCVG，DCVG再進(jìn)入腎臟代謝為DCVC有關(guān)[12]。

表4 基于PBPK模型預(yù)測的TCE及代謝產(chǎn)物組織內(nèi)暴露致癌風(fēng)險值Table 4 Cancer risks of internal tissue levels of TCE and its metabolites based on PBPK model

注：a基于組織最大濃度5%～95%分位數(shù)的均值。b空氣濃度為39.2 μg·m-3。c空氣濃度為119 mg·m-3。

Note:aBased on mean value of 5thto 95thpercentile maximal tissue concentration.bAir level of 39.2 μg·m-3.cAir level of 119 mg·m-3.

表5 TCE及代謝產(chǎn)物人體PBPK模型的中敏感性和高敏感性系數(shù)Table 5 Medium and high sensitivity coefficients in human PBPK model for TCE and its metabolites

注：QC表示心輸出血流量，QTissue表示組織血流量，QP表示肺換氣流量，BW表示體重。

Note: QC, cardiac output; QTissue, blood flow to tissue; QP, pulmonary ventilation; BW, body weight.

由于人體生理參數(shù)和代謝參數(shù)在個體層面存在較大的變異[9,25]，這些參數(shù)的變化會影響PBPK模型輸出結(jié)果。本研究基于表1～3中各參數(shù)的變異性，以體內(nèi)TCE及代謝產(chǎn)物的最大濃度為輸出結(jié)果，采用Monte Carlo方法對PBPK模型的不確定性進(jìn)行了分析。在本研究監(jiān)測的TCE濃度水平下，以肝臟中TCE濃度和組織DCVC濃度為例，結(jié)果表明，10萬次Monte Carlo模擬時其最大濃度值趨于收斂，說明模型輸出結(jié)果趨于穩(wěn)定。CTCE-Liver對應(yīng)的5%、50%和95%分位數(shù)分別為8.31×10-5、2.85×10-4和8.07×10-4mg·L-1；CDCVC對應(yīng)的5%、50%和95%分位數(shù)分別為1.85×10-8、1.10×10-7和6.25×10-7mg·L-1。不確定性分析表明，模型參數(shù)變化會顯著地影響PBPK模型的輸出結(jié)果，但輸出結(jié)果差異一般小于1個數(shù)量級，處于模型預(yù)測的可接受水平。

采用吸附管采樣-熱脫附/氣相色譜-質(zhì)譜法測定了某企業(yè)車間生產(chǎn)工況下TCE濃度，為(39.2±24.4) μg·m-3。PBPK模型預(yù)測表明，工人暴露8 h后TCE在不同組織中的最大濃度呈現(xiàn)出脂肪>腸>充分灌注室>支氣管>非充分灌注室>肝臟>靜脈血>動脈血的趨勢。體內(nèi)會生成多種與致癌有關(guān)的代謝產(chǎn)物，最大濃度表現(xiàn)為TCA > DCA > CHL > DCVC?；赥CE、TCA和DCVC內(nèi)暴露濃度與TCE外暴露濃度計(jì)算的致癌風(fēng)險評價結(jié)果存在顯著差異。敏感性分析表明血流量參數(shù)對模型輸出結(jié)果有重要影響，不確定性分析表明模型參數(shù)變化會顯著地影響PBPK模型輸出結(jié)果，但變異在可接受水平。TCE的人體健康風(fēng)險評價需要考慮其在體內(nèi)的分布和代謝過程。