蔣艷雪，賈宜靜，晁思宏，曹紅斌,*

1. 北京師范大學(xué) 中藥資源保護(hù)與利用北京市重點(diǎn)實(shí)驗(yàn)室，北京 100082 2. 北京師范大學(xué) 地理科學(xué)學(xué)部，北京 100875

基于PBPK模型的大鼠體內(nèi)鎘分布的比較：聯(lián)合暴露及單獨(dú)暴露

蔣艷雪1,2，賈宜靜1,2，晁思宏1,2，曹紅斌1,2,*

1. 北京師范大學(xué) 中藥資源保護(hù)與利用北京市重點(diǎn)實(shí)驗(yàn)室，北京 100082 2. 北京師范大學(xué) 地理科學(xué)學(xué)部，北京 100875

環(huán)境中同時(shí)存在著多種重金屬元素，聯(lián)合暴露與單獨(dú)暴露時(shí)，重金屬在體內(nèi)的蓄積分布情況也可能有所差異。為探究重金屬元素(汞、鉻、砷、鉛)對(duì)鎘(Cd)在體內(nèi)分布的影響，建立了大鼠在Cd暴露下的藥代動(dòng)力學(xué)(PBPK)模型，并進(jìn)行了包括Cd在內(nèi)5種重金屬的聯(lián)合毒性實(shí)驗(yàn)，比較了Cd單獨(dú)給藥與重金屬混合物給藥2種方式下大鼠肝臟、腎臟中的Cd濃度水平。結(jié)果表明，聯(lián)合暴露高(HgCl23.67 mg·kg-1，NaAsO23.67 mg·kg-1，CdCl210.55 mg·kg-1，K2Cr2O76.40 mg·kg-1，Pb(OOCCH3)2·3H2O 133.33 mg·kg-1)、中(HgCl20.367 mg·kg-1，NaAsO20.367 mg·kg-1，CdCl21.055 mg·kg-1，K2Cr2O70.640 mg·kg-1，Pb(OOCCH3)2·3H2O 13.333 mg·kg-1)、低(HgCl20.0367 mg·kg-1，NaAsO20.0367 mg·kg-1，CdCl20.1055 mg·kg-1，K2Cr2O70.0640 mg·kg-1，Pb(OOCCH3)2·3H2O 1.3333 mg·kg-1)劑量組大鼠肝臟中Cd濃度分別為13.37、0.78和0.06 μg·g-1；腎臟中Cd濃度分別為14.41、1.64和0.15 μg·g-1。與對(duì)照組相比，暴露組中Cd濃度有顯著升高，且不同劑量組之間均有顯著性差異。同劑量Cd單獨(dú)暴露的PBPK模擬結(jié)果顯示，肝臟及腎臟中的Cd濃度水平落在聯(lián)合毒性實(shí)驗(yàn)結(jié)果的濃度范圍內(nèi)，初步推斷其他4種重金屬的聯(lián)合暴露并沒(méi)有影響Cd在大鼠腎臟和肝臟中的濃度分布。

鎘；重金屬混合物；PBPK模型；聯(lián)合毒性；體內(nèi)分布

生理藥物代謝動(dòng)力學(xué)模型(Physiologically Based Pharmacokinetic Model, PBPK模型，在毒理學(xué)中一般稱(chēng)之為Physiologically Based Toxicokinetic Model, PBTK模型)是一種通過(guò)模擬藥物或毒物在人類(lèi)或其他動(dòng)物體內(nèi)的動(dòng)力學(xué)過(guò)程、代謝過(guò)程，預(yù)測(cè)其在目標(biāo)靶點(diǎn)有效濃度的工具[1-3]。PBPK模型可以用于描述單一或混合物質(zhì)的藥代/毒代動(dòng)力學(xué)過(guò)程[4-5]，目前被廣泛應(yīng)用于藥物開(kāi)發(fā)、有毒有害物質(zhì)管理等[6-7]?；旌衔镔|(zhì)的PBPK模型的研究對(duì)象多為有機(jī)物[8-11]。

有毒物質(zhì)與疾病登記處聯(lián)合機(jī)構(gòu)(ATSDR)基于二元的證據(jù)權(quán)重方法，研究了鉛(Pb)、鎘(Cd)、鉻(Cr)和砷(As)之間的兩兩交互作用[12]。但是，3種以上重金屬?gòu)?fù)合暴露的毒性實(shí)驗(yàn)數(shù)據(jù)和流行病學(xué)調(diào)查數(shù)據(jù)則十分有限，且不同重金屬在生物體內(nèi)的半衰期相差太大，模型的基本形式和數(shù)學(xué)結(jié)構(gòu)也不同，因此重金屬類(lèi)混合物質(zhì)的PBPK模型尚未建立。

重金屬混合物除了在毒性作用模式(mode of action)方面會(huì)相互影響外，在吸收、代謝、分布及排泄過(guò)程中也會(huì)相互影響。如，Yáez等[13]在Cd和As同時(shí)暴露的急性毒性試驗(yàn)中發(fā)現(xiàn)，As的存在會(huì)降低肝臟中Cd的濃度。Garcia等[14]發(fā)現(xiàn)，與Pb(300 mg·L-1)單獨(dú)存在時(shí)相比，Cd(10 mg·L-1)的存在會(huì)減少血液中Pb的含量。在一項(xiàng)研究Pb、Cd對(duì)腎臟中金屬硫蛋白表達(dá)影響的實(shí)驗(yàn)中發(fā)現(xiàn)，Pb和Cd在腎臟中的濃度不受彼此的影響，在血清中的濃度也不受彼此影響[15]。

本文針對(duì)重金屬單獨(dú)暴露與復(fù)合暴露產(chǎn)生的內(nèi)暴露劑量存在差異的狀況，建立重金屬Cd在大鼠體內(nèi)吸收、分布、代謝、排泄(absorption, distribution, metabolism and excretion, ADME)的PBPK模型(Cd-PBPK model)，模擬Cd在大鼠不同器官中的分布。同時(shí)，進(jìn)行含Cd重金屬混合物的聯(lián)合毒性實(shí)驗(yàn)，將實(shí)驗(yàn)大鼠肝臟、腎臟中Cd濃度與模型模擬結(jié)果比較，分析Cd及其他重金屬在ADME過(guò)程中的相互影響。

1 材料與方法(Materials and methods)

1.1 PBPK模型建立方法

利用Matlab/Simulink建立生理藥代動(dòng)力學(xué)模型，并進(jìn)行調(diào)試，用于模擬Cd在大鼠體內(nèi)的代謝、分布過(guò)程。在利用一套實(shí)驗(yàn)數(shù)據(jù)進(jìn)行參數(shù)優(yōu)化之后，將另一組毒性實(shí)驗(yàn)數(shù)據(jù)輸入已建立的模型進(jìn)行驗(yàn)證。

1.1.1 模型結(jié)構(gòu)

假設(shè)Cd及其代謝產(chǎn)物在大鼠體內(nèi)的分布轉(zhuǎn)運(yùn)是血流限速過(guò)程，建立血流限速型PBPK模型。該模型滿(mǎn)足以下假設(shè)：

1)Cd的透膜速度極快，其在各個(gè)器官的分布由血液進(jìn)入該器官的血流速率決定，即鎘進(jìn)入到組織細(xì)胞中的速率遠(yuǎn)大于其在血管中的運(yùn)輸速率；

2)從各個(gè)器官中流出的血液中的Cd及其代謝產(chǎn)物的濃度已經(jīng)與該器官組織細(xì)胞內(nèi)液中的濃度達(dá)到了平衡。

整體模型和血流循環(huán)過(guò)程如圖1所示。模型納入了主要的靶器官——肝臟和腎臟，一共包括6個(gè)器官(組織)房室和2個(gè)血液室。

按照物質(zhì)平衡原理，圖1中肺和其他組織是非消除器官，鎘在非消除器官中含量的變化速率可用微分方程式(1)表示：

(1)

圖1 生理藥物代謝動(dòng)力學(xué)模型(PBPK模型)模擬的鎘的血流循環(huán)Fig. 1 Blood circulation of Cd simulated by Physiologically Based Pharmacokinetic Model (PBPK model)

其中，下標(biāo)i指不同器官；V為器官體積(L)；C為器官或組織中的Cd濃度(μmolL-1)；t為模擬時(shí)間；Q為器官的血流量(L)；CA為動(dòng)脈血中鎘濃度(μmolL-1)；P為鎘在器官i中的組織-血液分配系數(shù)。

圖1中肝和腎是消除器官，鎘在消除器官中含量的變化速率可用微分方程式(2)表示：

(2)

其中，Kout為排泄系數(shù)。

在血液室中，根據(jù)鎘在血漿和血紅細(xì)胞中的轉(zhuǎn)移，列為2個(gè)微分方程，如式(3)和式(4)所示。

血漿：

(3)

血紅細(xì)胞：

(4)

其中，下標(biāo)plas和rbc分別指血漿和血紅細(xì)胞；Qc為肺血流量(L)；Ppr為鎘從血漿進(jìn)入血紅細(xì)胞的滲透系數(shù)；Prp為鎘從血紅細(xì)胞進(jìn)入血漿的滲透系數(shù)。

1.1.2 模型參數(shù)

模型中各器官的生理參數(shù)在表1中列出。大鼠的基本身體參數(shù)，如器官體積、心輸出量、各器官血流分?jǐn)?shù)等都來(lái)自文獻(xiàn)[16]。Cd在各個(gè)組織器官中的分配系數(shù)來(lái)自參數(shù)優(yōu)化，基于Cha等[17]的實(shí)驗(yàn)數(shù)據(jù)，使模型模擬結(jié)果最大限度地逼近實(shí)驗(yàn)結(jié)果，以確定參數(shù)最優(yōu)取值。

由于Cd在體內(nèi)的半衰期較長(zhǎng)，達(dá)數(shù)年甚至數(shù)十年。因此，在模型模擬時(shí)，時(shí)間步長(zhǎng)按照天來(lái)計(jì)算。

1.1.3 整體仿真模型建立

在Simulink仿真環(huán)境下，根據(jù)各個(gè)器官房室的微分方程，設(shè)計(jì)各個(gè)器官的子模型，且將其封裝成為整體模型中的子模塊。將各個(gè)基本模塊根據(jù)整體模型和血流循環(huán)規(guī)律進(jìn)行連接，并且配以產(chǎn)生并且存儲(chǔ)系統(tǒng)時(shí)間的模塊，使之成為完整的實(shí)時(shí)仿真模型。Simulink環(huán)境下建立的整體模型結(jié)構(gòu)如圖2所示。

1.2 實(shí)驗(yàn)方法與材料

1.2.1 實(shí)驗(yàn)方法

8周齡SD大鼠80只，雌、雄各半，購(gòu)回后馴養(yǎng)一周使其適應(yīng)實(shí)驗(yàn)室環(huán)境，平均體重為(127.8±2.2) g。根據(jù)實(shí)驗(yàn)設(shè)計(jì)隨機(jī)分為4組，分別為：對(duì)照組(A)，高劑量組(B)，中劑量組(C)，低劑量組(D)。每個(gè)劑量組共20只大鼠，雌、雄各10只，每5只飼養(yǎng)于一只籠中。將5種重金屬藥品混合，溶于經(jīng)蒸汽滅菌的水中，配制后供不同實(shí)驗(yàn)組自由飲用，對(duì)照組飲用未添加重金屬的滅菌水。由于實(shí)驗(yàn)規(guī)模所限，在進(jìn)行實(shí)驗(yàn)動(dòng)物的分組給藥時(shí)，沒(méi)有分出單獨(dú)給藥的組，即沒(méi)有單一重金屬的對(duì)照組。各組大鼠均自由攝取食物，室溫為(25±1) ℃，室內(nèi)濕度50%±10% RH。

圖2 大鼠體內(nèi)Cd的PBPK模型整體結(jié)構(gòu)Fig. 2 Whole structure of Cd-PBPK model in rats

表1 大鼠體內(nèi)鎘PBPK模型的參數(shù)取值Table 1 Parameters values in Cd-PBPK model in rats

每天記錄飲水量，每周稱(chēng)量體重一次，根據(jù)體重變化調(diào)整給藥量，連續(xù)暴露3個(gè)月。

1.2.2 實(shí)驗(yàn)試劑及劑量

依據(jù)SD大鼠急性毒性實(shí)驗(yàn)結(jié)果，取HgCl2(Amresco公司，分析純)、Cr2K2O7(Sigma公司，分析純)、CdCl2(Alfa Aesar公司，分析純)、NaAsO2(阿拉丁公司，分析純)以及Pb(OOCCH3)2·3H2O(Alfa Aesar公司，分析純)5種藥品對(duì)大鼠半數(shù)致死劑量(LD50)的1/15，分別為HgCl23.67 mg·kg-1，NaAsO23.67 mg·kg-1，CdCl210.55 mg·kg-1，K2Cr2O76.40 mg·kg-1，以及Pb(OOCCH3)2·3H2O 133.33 mg·kg-1作為高劑量，充分混合后稀釋10倍和100倍分別作為中劑量和低劑量，即各組間距為10倍。

1.2.3 器官重金屬含量測(cè)定

大鼠處死后，將組織器官稱(chēng)重，立刻保存在-20 ℃冰箱，檢測(cè)肝臟、腎臟中重金屬含量(由天津海世達(dá)檢測(cè)技術(shù)有限公司測(cè)定)。采用電感耦合等離子體質(zhì)譜儀(ICP-MS)(安捷倫7500a，美國(guó))檢測(cè)。所用標(biāo)準(zhǔn)物質(zhì)為凍干人尿鉛標(biāo)準(zhǔn)物質(zhì)(GBW09105e，中國(guó)疾病預(yù)防控制中心職業(yè)衛(wèi)生與中毒控制所)，凍干牛血鉛、鎘成份分析標(biāo)準(zhǔn)物質(zhì)(GBW09140g，中國(guó)疾病預(yù)防控制中心職業(yè)衛(wèi)生與中毒控制所)和質(zhì)量控制物質(zhì)凍干牛血汞(ZK021-2，中國(guó)疾病預(yù)防控制中心職業(yè)衛(wèi)生與中毒控制所)，國(guó)家標(biāo)準(zhǔn)物質(zhì)豬肉(GBW08552，中國(guó)科學(xué)院上海原子核研究所)，豬肝(GBW10015，中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所)，牛肝(GBW(E)080193，解放軍總醫(yī)院臨床基礎(chǔ)學(xué)研究所)。

1.3 數(shù)據(jù)統(tǒng)計(jì)與分析

采用SPSS 22.0軟件(IBM SPSS Inc.)進(jìn)行方差齊性檢驗(yàn)和單因素方差分析(One Way ANOVA)。若方差齊性，采用SNK檢驗(yàn)(Student-Newman-Keuls法)；若方差不齊，采用Dunnett’s T3檢驗(yàn)。P < 0.05、P < 0.01表示差異顯著。

2 結(jié)果(Results)

2.1 模型驗(yàn)證

根據(jù)建立的生理藥代動(dòng)力學(xué)模型，以Asagba等[18]的毒性實(shí)驗(yàn)為基礎(chǔ)數(shù)據(jù)，進(jìn)行模型的驗(yàn)證，結(jié)果如圖3所示。

圖3表明，在暴露1個(gè)月和暴露3個(gè)月時(shí)，肝臟和腎臟中Cd的模擬濃度與實(shí)驗(yàn)測(cè)得的濃度一致，模擬精度較好。上述建立的模型可用于模擬鎘單獨(dú)暴露時(shí)在大鼠體內(nèi)的吸收、分布、代謝、排泄過(guò)程。

圖3 單一Cd暴露的大鼠肝、腎Cd濃度PBPK模型模擬值與文獻(xiàn)實(shí)驗(yàn)值的比較Fig. 3 Comparison between PBPK simulated and reported measured concentrations of Cd in liver and kidney of rat exposed to single Cd

2.2 5種重金屬聯(lián)合暴露時(shí)大鼠體內(nèi)Cd的分布

重金屬聯(lián)合暴露毒性實(shí)驗(yàn)中，暴露3個(gè)月后，各個(gè)劑量組腎臟、肝臟中Cd的濃度結(jié)果如表2所示?？梢钥闯觯?、雄大鼠肝臟及腎臟中Cd的濃度均隨給藥濃度的增加而增加。在肝臟和腎臟中Cd濃度均為雌性大鼠>雄性大鼠。對(duì)于雌、雄大鼠肝臟及腎臟中的Cd濃度，各劑量組與對(duì)照組間，以及各劑量組之間均有顯著差異(P<0.01)。

2.3 PBPK模型模擬結(jié)果與毒性實(shí)驗(yàn)結(jié)果的比較

將重金屬聯(lián)合暴露實(shí)驗(yàn)結(jié)果(肝臟、腎臟中Cd的濃度)與利用Cd的PBPK模型模擬得到的Cd在不同器官中的濃度進(jìn)行比較，分析在有其他重金屬存在的情況下，Cd在不同器官中的分布是否會(huì)受到影響。結(jié)果表明，高、中、低劑量組的肝臟和腎臟中Cd濃度的模擬結(jié)果基本落在重金屬聯(lián)合暴露毒性實(shí)驗(yàn)樣本點(diǎn)的濃度范圍之內(nèi)(圖4)。

3 討論(Discussion)

目前，大部分對(duì)重金屬聯(lián)合暴露的體內(nèi)分布實(shí)驗(yàn)研究多為2～3種重金屬，4種及以上的研究較少。Cobbina等[19]分析了4種重金屬(NaAsO20.003 mg·kg-1，CdCl20.68 mg·kg-1，Pb(OOCCH3)20.34 mg·kg-1，HgCl21.01 mg·kg-1)連續(xù)暴露120 d的ICR小鼠的大腦、肝臟和腎臟中的重金屬含量。結(jié)果表明，與Cd單獨(dú)暴露相比，多種重金屬聯(lián)合暴露(Pb+Cd、Pb+Hg+Cd、Pb+Cd+As+Hg)并沒(méi)有使大腦、肝臟和腎臟中的Cd濃度升高。Rai等[20]研究了As(NaAsO2，3.80 mg·kg-1)、Cd(CdCl2，0.98 mg·kg-1)、Pb(Pb(OOCCH3)2，2.22 mg·kg-1)對(duì)大鼠從孕期到產(chǎn)后大腦發(fā)育的聯(lián)合毒性。結(jié)果表明，在As、Cd、Pb聯(lián)合暴露與Cd單獨(dú)暴露情境下，大腦皮層中Cd含量無(wú)顯著差異。其研究所用的As、Cd、Pb劑量分別與本實(shí)驗(yàn)的高劑量組(3.67 mg·kg-1)、中劑量組(1.06 mg·kg-1)、低劑量組(1.33 mg·kg-1)接近。

Pandya等[21]在Cd和Pb聯(lián)合毒性實(shí)驗(yàn)中發(fā)現(xiàn)，Pb、Cd聯(lián)合暴露時(shí)大鼠體內(nèi)肝的Cd含量較Cd單獨(dú)暴露時(shí)顯著下降，其實(shí)驗(yàn)劑量遠(yuǎn)低于本實(shí)驗(yàn)的低劑量組，暴露方式也不同，為腹腔注射方法；但Massó等[22]研究了飲水途徑Cd和Pb聯(lián)合暴露對(duì)圍產(chǎn)期大鼠肝臟氧化損傷的影響，所用Cd、Pb劑量分別約為本實(shí)驗(yàn)中劑量組的2倍和6倍，結(jié)果表明肝臟中Cd含量在聯(lián)合暴露時(shí)并無(wú)顯著變化；Shukla和Chandra[23]研究了Cd、Pb復(fù)合暴露情境下大鼠體內(nèi)各器官中重金屬分布情況，Cd通過(guò)腹腔注射CdCl2

圖4 雄性大鼠肝、腎Cd濃度的單一Cd暴露PBPK模擬值與重金屬聯(lián)合暴露實(shí)驗(yàn)值的比較(A肝臟，B腎臟)Fig. 4 Comparison of Cd concentration in liver (A) or kidney (B) of male rats between PBPK simulation of single Cd exposure and toxicity experiment of metal mixtures

表2 重金屬聯(lián)合暴露后大鼠肝臟及腎臟中Cd的濃度(μg·g-1)Table 2 Concentrations of Cd in liver and kidney of rats exposed to metal mixtures (μg·g-1)

注：平均值±標(biāo)準(zhǔn)偏差(mean±SD)；采用Dunnett’ s T3方法進(jìn)行多重比較，a, b, c, d, e上標(biāo)字母相同時(shí)表示在1%水平上差異無(wú)統(tǒng)計(jì)學(xué)意義。高劑量為HgCl23.67 mg·kg-1，NaAsO23.67 mg·kg-1，CdCl210.55 mg·kg-1，K2Cr2O76.40 mg·kg-1，Pb(OOCCH3)2·3H2O 133.33 mg·kg-1；中劑量為HgCl20.367 mg·kg-1，NaAsO20.367 mg·kg-1，CdCl21.055 mg·kg-1，K2Cr2O70.640 mg·kg-1，Pb(OOCCH3)2·3H2O 13.333 mg·kg-1；低劑量為HgCl20.0367 mg·kg-1，NaAsO20.0367 mg·kg-1，CdCl20.1055 mg·kg-1，K2Cr2O70.0640 mg·kg-1，Pb(OOCCH3)2·3H2O 1.3333 mg·kg-1。

Note: values were shown as mean±standard deviation. The Dunnett’s T3 method was used for multiple comparisons. Within the same columns, values with different superscripts letters (a, b, c, d, e) were significantly different (P < 0.01). High dose, HgCl23.67 mg·kg-1, NaAsO23.67 mg·kg-1, CdCl210.55 mg·kg-1, K2Cr2O76.40 mg·kg-1, Pb(OOCCH3)2·3H2O 133.33 mg·kg-1; medium dose, HgCl20.367 mg·kg-1, NaAsO20.367 mg·kg-1, CdCl21.055 mg·kg-1, K2Cr2O70.640 mg·kg-1, Pb(OOCCH3)2·3H2O 13.333 mg·kg-1; low dose, HgCl20.0367 mg·kg-1, NaAsO20.0367 mg·kg-1, CdCl20.1055 mg·kg-1, K2Cr2O70.0640 mg·kg-1, Pb(OOCCH3)2·3H2O 1.3333 mg·kg-1.

給藥，兩組劑量分別與本實(shí)驗(yàn)中低劑量組接近，Pb通過(guò)飲水添加Pb(OOCCH3)2給藥，劑量與本實(shí)驗(yàn)低劑量相當(dāng)，結(jié)果也表明Cd在肝臟中的含量未受到Pb的影響。

Shukla和Chandra[23]還發(fā)現(xiàn)與Cd單一暴露組相比，Pb+Cd顯著增加了腎臟中的Cd含量。而Wang等[15]使用與本實(shí)驗(yàn)高劑量組接近的劑量，使大鼠通過(guò)自由飲水進(jìn)行復(fù)合暴露，發(fā)現(xiàn)在血清和腎皮質(zhì)中Cd濃度在Pb+Cd組與Cd組間無(wú)顯著差異，Cd主要富集在腎臟中，而Pb主要分布在血液中。

根據(jù)本研究結(jié)果，結(jié)合已有的文獻(xiàn)報(bào)道，僅做出初步推斷，即在肝臟和腎臟中，其他4種重金屬的共同作用未對(duì)Cd在大鼠肝臟和腎臟中的濃度分布產(chǎn)生影響。

[1] Barton H A, Chiu W A, Setzer R W, et al. Characterizing uncertainty and variability in physiologically based pharmacokinetic models: State of the science and needs for research and implementation [J]. Toxicological Sciences, 2007, 99(2): 395-402

[2] Bouvier D M, Prieto P, Blaauboer B J, et al. Physiologically-based Kinetic Modelling (PBK Modelling): Meeting the 3Rs agenda. The report and recommendations of ECVAM Workshop 63. [J]. Alternatives to Laboratory Animals: ATLA, 2008, 35(6): 661-671

[3] Loizou G, Spendiff M, Barton H A, et al. Development of good modelling practice for physiologically based pharmacokinetic models for use in risk assessment: The first steps [J]. Regulatory Toxicology and Pharmacology, 2008, 50(3): 400-411

[4] Haddad S, Béliveau M, Tardif R, et al. A PBPK modeling-based approach to account for interactions in the health risk assessment of chemical mixtures [J]. Toxicological Sciences, 2001, 63(1): 125-131

[5] Yang R S, El-Masri H A, Thomas R S, et al. Chemical mixture toxicology: From descriptive to mechanistic, and going on to in silico toxicology [J]. Environmental Toxicology and Pharmacology, 2004, 18(2): 65-81

[6] Zhao P, Zhang L, Grillo J A, et al. Applications of physiologically based pharmacokinetic (PBPK) modeling and simulation during regulatory review [J]. Clinical Pharmacology and Therapeutics, 2011, 89(2): 259-267

[7] Edginton A N, Theil F, Schmitt W, et al. Whole body physiologically-based pharmacokinetic models: Their use in clinical drug development [J]. Expert Opinion on Drug Metabolism and Toxicology, 2008, 4(9): 1143-1152

[8] Jasper M N, Martin S A, Oshiro W M, et al. Application of biologically based lumping to investigate the toxicokinetic interactions of a complex gasoline mixture [J]. Environmental Science and Technology, 2016, 50(6): 3231-3238

[9] Olie J D N, Bessems J G, Clewell H J, et al. Evaluation of semi-generic PBTK modeling for emergency risk assessment after acute inhalation exposure to volatile hazardous chemicals [J]. Chemosphere, 2015, 132: 47-55

[10] LeFew W, El-Masri H. Computational estimation of errors generated by lumping of physiologically-based pharmacokinetic (PBPK) interaction models of inhaled complex chemical mixtures [J]. Inhalation Toxicology, 2012, 24(1): 36-46

[11] Knaak J B, Dary C C, Okino M S, et al. Parameters for Carbamate Pesticide QSAR and PBPK/PD Models for Human Risk Assessment in Reviews of Environmental Contamination and Toxicology [M]. Springer, 2008: 53-212

[12] ATSDR. Interaction profile for arsenic, cadmium, chromium and lead [R]. Atlanta, GA: U.S. Department of Health and Human Services, Public Health Service, 2004

[14] Garcia T A, Corredor L. Biochemical changes in the kidneys after perinatal intoxication with lead and/or cadmium and their antagonistic effects when coadministered [J]. Ecotoxicology and Environmental Safety, 2004, 57(2): 184-189

[15] Wang L, Chen D, Wang H, et al. Effects of lead and/or cadmium on the expression of metallothionein in the kidney of rats [J]. Biological Trace Element Research, 2009, 129(1): 190-199

[16] Brown R P, Delp M D, Lindstedt S L, et al. Physiological parameter values for physiologically based pharmacokinetic models [J]. Toxicology and Industrial Health, 1997, 13(4): 407-484

[17] Cha C W. A study on the effect of garlic to the heavy metal poisoning of rat [J]. Journal of Korean Medical Science, 1987, 2(4): 213-224

[18] Asagba S. Comparative effect of water and food-chain mediated cadmium exposure in rats [J]. BioMetals, 2010, 23(6): 961-971

[19] Cobbina S J, Chen Y, Zhou Z, et al. Low concentration toxic metal mixture interactions: Effects on essential and non-essential metals in brain, liver, and kidneys of mice on sub-chronic exposure [J]. Chemosphere, 2015, 132: 79-86

[20] Rai A, Maurya S K, Khare P, et al. Characterization of developmental neurotoxicity of As, Cd, and Pb mixture: Synergistic action of metal mixture in glial and neuronal functions [J]. Toxicological Sciences, 2010, 118(2): 586-601

[21] Pandya C D, Pillai P P, Gupta S S. Lead and cadmium co-exposure mediated toxic insults on hepatic steroid metabolism and antioxidant system of adult male rats [J]. Biological Trace Element Research, 2010, 134(3): 307-317

[22] Massó E L, Corredor L, Antonio M T. Oxidative damage in liver after perinatal intoxication with lead and/or cadmium [J]. Journal of Trace Elements in Medicine and Biology, 2007, 21(3): 210-216

[23] Shukla G S, Chandra S V. Concurrent exposure to lead, manganese, and cadmium and their distribution to various brain regions, liver, kidney, and testis of growing rats [J]. Archives of Environmental Contamination and Toxicology, 1987, 16(3): 303-310

◆

Comparison of Cadmium Distribution in Rats between Exposed to Individual Cadmium and Heavy Metals Mixture Based on PBPK Model

Jiang Yanxue1,2, Jia Yijing1,2, Chao Sihong1,2, Cao Hongbin1,2,*

1. Beijing Area Major Laboratory of Protection and Utilization of Traditional Chinese Medicine, Beijing Normal University, Beijing 100082, China 2. Falculty of Geographical Science, Beijing Normal University, Beijing 100875, China

Received 4 May 2016 accepted 29 June 2016

The distribution in vivo of metals mixture might be different from that of individual metals. In order to explore the influence of four metals (mercury, chromium, arsenic and lead) on the distribution of cadmium (Cd) in vivo, a physiologically based pharmacokinetic model (PBPK) was constructed to simulate the absorption, distribution, metabolism and excretion (ADME) of Cd in rats exposed to individual Cd. The levels of Cd in liver and kidney exposed to metals mixture were obtained from a joint toxicity experiment. The experiment results showed that the concentrations of Cd in liver in high (HgCl23.67 mg·kg-1, NaAsO23.67 mg·kg-1, CdCl210.55 mg·kg-1, K2Cr2O76.40 mg·kg-1, Pb(OOCCH3)2·3H2O 133.33 mg·kg-1), medium (HgCl20.367 mg·kg-1, NaAsO20.367 mg·kg-1, CdCl21.055 mg·kg-1, K2Cr2O70.640 mg·kg-1, Pb(OOCCH3)2·3H2O 13.333 mg·kg-1) and low (HgCl20.0367 mg·kg-1, NaAsO20.0367 mg·kg-1, CdCl20.1055 mg·kg-1, K2Cr2O70.0640 mg·kg-1, Pb(OOCCH3)2·3H2O 1.3333 mg·kg-1) dose groups were 13.37, 0.78 and 0.06 μg·g-1, respectively; and those in kidney were 14.41, 1.64 and 0.15 μg·g-1, respectively. There were significant differences between the control group and exposure groups. The simulated concentrations of Cd in liver and kidney were in the range of the corresponding experimental concentrations. It is preliminarily inferred that the co-exposure with other four heavy metals had no obvious influence on Cd distribution in liver and kidney of rats.

cadmium; metals mixture; PBPK model; joint toxicity; distribution in vivo

國(guó)家科技支撐計(jì)劃“村鎮(zhèn)環(huán)境評(píng)價(jià)與預(yù)警平臺(tái)建設(shè)”基金資助項(xiàng)目(2012BAJ24B04)

蔣艷雪(1992-)，女，博士研究生，研究方向?yàn)榄h(huán)境健康，E-mail: 201431190020@mail.bnu.edu.cn；

*通訊作者(Corresponding author)， E-mail: caohongbin@bnu.edu.cn

10.7524/AJE.1673-5897.20160504001

2016-05-04 錄用日期：2016-06-29

1673-5897(2016)6-121-07

X171.5

A

曹紅斌(1966—)，女，工學(xué)博士，教授，主要研究方向?yàn)榄h(huán)境風(fēng)險(xiǎn)評(píng)價(jià)與管理、天然藥物安全性評(píng)價(jià)。

蔣艷雪, 賈宜靜, 晁思宏, 等. 基于PBPK模型的大鼠體內(nèi)鎘分布的比較：聯(lián)合暴露及單獨(dú)暴露[J]. 生態(tài)毒理學(xué)報(bào)，2016, 11(6): 121-127

Jiang Y X, Jia Y J, Chao S H, et al. Comparison of cadmium distribution in rats between exposed to individual cadmium and heavy metals mixture based on PBPK model [J]. Asian Journal of Ecotoxicology, 2016, 11(6): 121-127 (in Chinese)