【建模文章解讀 】 采用系統定量毒理學(xué)評估因小分子化合物導致線(xiàn)粒體損傷和膽汁酸毒性，而出現肝毒性呈現種屬差異
原文作者
Grant Generaux1, Vinal V. Lakhani1, Yuching Yang1, Sashi Nadanaciva2, Luping Qiu3, Keith Riccardi4, Li Di4, Brett A. Howell1,, Scott Q. Siler1,, Paul B. Watkins5,6, Hugh A. Barton7, Michael D. Aleo3, Lisl K. M. Shoda1
1 美國DILIsym Services 公司
2 美國輝瑞制藥公司，全球藥物化學(xué)部
3 美國輝瑞制藥公司，藥物安全研究與開(kāi)發(fā)部
4 美國輝瑞制藥公司，醫藥科學(xué)部
5 美國北卡羅來(lái)納大學(xué)教堂山分校，埃舍曼藥學(xué)院
6 美國北卡羅來(lái)納大學(xué)教堂山分校，藥物安全科學(xué)研究所
7 美國輝瑞制藥公司，生物醫學(xué)設計部
導 讀
FDA當前非常鼓勵制藥企業(yè)使用模型的方法評估在研化合物特性，以提高藥物的安全性。本次解讀的案例：采用定量毒理學(xué)軟件DILIsym重現了PF-04895162種屬差異性的肝毒性結果，解釋了化合物在人體產(chǎn)生肝毒性的原因。通過(guò)案例學(xué)習，可了解如何通過(guò)DILIsym定量評估肝損傷。
案例摘要
該案例說(shuō)明如何采用DILIsym軟件(一款評估肝損傷的定量毒理學(xué)計算工具)重現PF-04895162 (ICA-105665)導致的種屬差異性肝毒性。PF-04895162是開(kāi)發(fā)用于治療癲癇癥的藥物，由于在健康志愿者中觀(guān)測到氨基轉移酶的升高而被迫終止(NCT01691274)，然而在臨床前安全性研究中并沒(méi)有發(fā)現肝臟安全性的問(wèn)題。整合表征肝毒性機制以及預測體內肝臟暴露的體外數據到DILIsym軟件中，重現該化合物的臨床肝毒性與未在大鼠身上發(fā)現肝毒性的情況。從多重維度上模擬考察差異性，并發(fā)現人體肝臟暴露要高于大鼠。模擬結果提示人體肝毒性是因為化合物同時(shí)造成了線(xiàn)粒體毒性和膽汁酸轉運體的抑制(假設去除上述一項機理都可以使其肝損傷消除)。模型表明更高的PF-04895162肝臟暴露(比大鼠安全性試驗給藥劑量時(shí)的暴露更高)時(shí)，會(huì )加重線(xiàn)粒體損傷但依然不會(huì )導致大鼠肝毒性，這是因為該種屬不能有效地積累產(chǎn)生細胞毒素的膽汁酸。這些研究提示聯(lián)合體外數據和計算機篩選的方法可以鑒定潛在肝毒性的風(fēng)險，并為類(lèi)似的前瞻性研究奠定的基礎。
推薦理由
很多臨床前看上去很有前景的化合物，但由于臨床的安全性(包括肝臟毒性)問(wèn)題而被終止。本案例聯(lián)合體外數據(表征肝毒性的機制性參數以及預測肝臟暴露的模型數據)和DILIsym軟件(定量評估肝損傷的模型軟件)，重現化合物PF-04895162種屬差異性的肝毒性結果，并解釋了該化合物在人體產(chǎn)生肝毒性的原因。通過(guò)研究工作，說(shuō)明模型可以準確反映觀(guān)測結果，并為其它藥物的肝毒性評估工作奠定了基礎。此外，該工作也可用于詳細學(xué)習DILIsym的模型搭建、模型參數、結果分析等知識。
軟件用途
案例中，利用DILIsym軟件搭建PF-04895162肝損傷評估模型?？疾炝嘶衔镌诟闻K的暴露(通過(guò)生理藥代動(dòng)力學(xué)PBPK模型進(jìn)行預測)、可能產(chǎn)生肝毒性的一種或多種機理(多重維度進(jìn)行分析確認)等相關(guān)因素重現PF-04895162種屬差異性的肝毒性結果。通過(guò)模型的開(kāi)展，發(fā)現化合物同時(shí)造成了線(xiàn)粒體毒性和膽汁酸轉運體的抑制，進(jìn)而產(chǎn)生了人體的肝毒性(去除其中一個(gè)機理都會(huì )使得肝毒性消除)；而大鼠由于不能有效累積膽汁酸，因此即使比安全性更高的肝臟暴露也僅加重線(xiàn)粒體的損失，但不會(huì )直接產(chǎn)生肝毒性。
1 研究背景
FDA提出創(chuàng )新監管科學(xué)計劃(ARSI)作為改善監管科學(xué)質(zhì)量和效率的戰略計劃，此項戰略的一個(gè)方面便是開(kāi)發(fā)計算方法和數學(xué)模型以提高藥品的安全性。DILIsym軟件應用定量系統毒理學(xué)(QST)方法評估藥物性肝損傷，現已成功用于解釋臨床肝信號的數據。此外，DILIsym還能重現大鼠與人肝毒性的種屬差異性。
PF-04895162 (ICA-105665，由Icagen公司發(fā)現)是開(kāi)發(fā)用于治療癲癇病的全新小分子化合物。PF-04895162在多個(gè)臨床前種屬進(jìn)行了安全性評估：大鼠給藥7天的毒性試驗中，丙氨酸氨基轉移酶(ALT)的水平呈現劑量的相關(guān)性，提示有潛在的肝毒性；然而并沒(méi)有檢測到肝損傷的組織學(xué)變化證據，且在7天給藥的重現研究中未再觀(guān)察到ALT的水平變化。大鼠28天和6個(gè)月的毒性試驗也未發(fā)現氨基轉移酶的水平變化和肝臟毒性；獼猴9個(gè)月的毒性研究結果也顯示陰性的肝毒性結果?？傊?，臨床前的研究未發(fā)現肝臟可能是該化合物安全性問(wèn)題的靶器官。
在臨床研究中，健康志愿者單次給藥劑量爬坡至600 mg或多次給藥劑量爬坡至每天兩次給藥200 mg (BID, 共7天)均未發(fā)現有明顯數據(如ALT > 3x ULN)說(shuō)明PF‐04895162會(huì )造成肝損傷。另外一個(gè)7天給藥的研究中，12個(gè)志愿者250 mg BID或300 mg BID給藥后發(fā)現一例ALT的升高。但意外的是，14天給藥的研究，8個(gè)志愿者300 mg BID給藥后有6例出現了ALT的升高，甚至高達5x正常值的上限(ULN)；這項研究也直接導致了該化合物開(kāi)發(fā)計劃的中止。
因此本工作選擇PF‐04895162作為DILIsym軟件評估種屬差異性肝毒性的潛在案例，聯(lián)合肝毒性機制和預測體內暴露的體外數據至DILIsym中以預測肝毒性。特別測定了PF‐04895162與線(xiàn)粒體功能、化合物與膽汁酸轉運體相作用的體外數據；通過(guò)生理藥代動(dòng)力學(xué)(PBPK)模型預測化合物在大鼠及人體內的暴露。模擬結果發(fā)現可以重現該化合物在臨床人群的肝毒性，但非臨床種屬無(wú)肝毒性的結果。此外，模擬發(fā)現人體肝毒性是由于該化合物同時(shí)導致了線(xiàn)粒體損傷和膽汁酸毒性，去除其中一個(gè)機理都不會(huì )使得ALT的升高。最后模型設置的劑量方案以及肝毒性的不確定性也均和臨床與動(dòng)物觀(guān)測結果一致。該案例也說(shuō)明了DILIsym軟件可以用于鑒定潛在的人體肝毒性風(fēng)險。
2 建模數據與處理
2.1 體外表征肝毒性機制的實(shí)驗
為了機制性地模擬PF-04895162介導的肝損傷，收集了體外表征PF-04895162與線(xiàn)粒體毒性、膽汁酸轉運體、氧化應激相互關(guān)系的數據。當前研究所采用的數據來(lái)自先前的報道，且未觀(guān)察到PF-04895162會(huì )產(chǎn)生氧化應激。采用DILIsym的膽汁酸及線(xiàn)粒體子模型進(jìn)行模擬，相關(guān)的毒性數據匯總如下：
膽汁酸轉運體毒性數據：通過(guò)標準的囊泡轉運實(shí)驗測定PF-04895162對人BSEP, MRP3和MRP4以及大鼠Bsep和Mrp3等轉運體的半數抑制濃度(IC50)；假如測定的抑制濃度達不到抑制50%，則外推曲線(xiàn)計算IC50；然后將測定或外推的數據直接用于DILIsym。如果不存在抑制常數(Ki)的數據，DILIsym通常的做法是選擇α=5以表征為混合的膽汁酸抑制特性。實(shí)驗結果發(fā)現PF-04895162對人或大鼠膽汁酸轉運體的作用具有濃度依賴(lài)性：PF-04895162抑制人外排轉運體MRP4、未能直接測到外排轉運體BSEP及 MRP3的IC50 (不能有效抑制)、也未能測到攝取轉運體NTCP的IC50；PF-04895162能夠抑制大鼠Mrp3及Bsep外排轉運體、對攝取轉運體Ntcp的IC50也外排轉運體相近。
線(xiàn)粒體毒性數據：為了表征潛在的大鼠及人線(xiàn)粒體毒性，將不同濃度的PF-04895162與大鼠及人肝細胞孵育24小時(shí)。實(shí)驗結果表明該化合物對基礎耗氧量(OCR)具有微弱的影響，24小時(shí)內OCR基本沒(méi)有下降；當加入羰基氰對三氟甲氧基苯腙(FCCP)后，發(fā)現化合物對肝細胞OCR的作用呈現劑量依賴(lài)性；這些結果提示化合物可能是線(xiàn)粒體電子傳遞鏈(ETC)的輕微抑制劑；為準確描述作用部位處化合物濃度與效應之間的關(guān)系，測定了化合物在肝細胞內的濃度，發(fā)現結果要遠高于理論的平均濃度。MITOsym模塊用于再處理所測定的Seahorse數據(也即體外測定的PF-04895162線(xiàn)粒體毒性的數據)，然后采用一個(gè)轉換因子將MITOym的參數轉為DILIsym的體內數據。
2.2 PBPK建模
采用DILIsym的PBPK子模型重現PF-04895162在大鼠及人觀(guān)測的血漿濃度-時(shí)間曲線(xiàn)。最為簡(jiǎn)化的PBPK模型可以用于表征每個(gè)個(gè)體(人)的PK行為；由于口服給藥后存在腸道飽和吸收，因此需要采用稍微復雜些的PBPK模型表征大鼠的PK特性。大鼠PBPK子模型包含了血液、肝臟、肌肉、腸道和其它組織、人體的模型則僅包含了血液和肝臟兩個(gè)房室(剩余的器官則全部納入系統分布體積中)；兩個(gè)種屬的子模型，均采用血流限速型的組織。相關(guān)設置及優(yōu)化的模型參數如下表所示：
2.3 模擬個(gè)體的基線(xiàn)值
個(gè)體模擬的基線(xiàn)值采用人體測量學(xué)及生物化學(xué)性質(zhì)的平均值，再基于基線(xiàn)值優(yōu)化大鼠及PBPK子模型。
2.4 Human SimCohorts
SimCohorts是比SimPops數據集稍少的虛擬個(gè)體(n = 16)，本工作選擇了Human_ROS_apap_mito_BA_v8A_1_Multi16_A SimCohort，包含了線(xiàn)粒體功能損傷、膽汁酸淤積、氧化應激等個(gè)體敏感的參數，并用于研究肝毒性的機理。模型開(kāi)展多種方式的研究，如果整合PF-04895162所有的毒性機理、僅不存在線(xiàn)粒體功能損傷、僅不存在膽汁酸毒性等情況；當不存在某一項機理時(shí)，查看PF-04895162產(chǎn)生臨床個(gè)體ALT水平升高的情況，并分析對臨床肝毒性的貢獻。
2.5 SimPops
SimPops收集了反映人體測量學(xué)及生物化學(xué)等參數變異的虛擬個(gè)體。本工作所有人體的模擬是基于Human_ROS_apop_mito_BA_v8A_1 (n = 285) SimPops進(jìn)行的，其中健康志愿者的體重、線(xiàn)粒體功能、半胱天冬酶的激活(細胞凋亡)、膽汁酸轉運體的表達、氧化應激等參數的群體分布見(jiàn)支持文件；大鼠的模擬則基于Rat_ROS_apop_mito_BA_v8A_11 (n = 294) SimPops進(jìn)行的，相關(guān)模型參數的群體分布也見(jiàn)支持文件。
2.6 模擬方案
優(yōu)化大鼠PBPK模型的模擬方案為：1 mg/kg IV或10, 30, or 100 mg/kg/day每天口服給藥一次至28天；優(yōu)化人PBPK模型的給藥方案為：300 mg單次口服給藥，并使用300 mg每天口服給藥兩次(BID)測試了建立的模型。對于肝毒性研究，模擬了大鼠100 mg/kg/day口服給藥至28天；人則考察了300 mg BID口服給藥至14天，并緊接著(zhù)一個(gè)14天的清洗期。人體還模擬了200 mg BID口服給藥至7天，以重現早期臨床研究中未發(fā)現肝毒性的結果。
2.7 多重線(xiàn)性回歸分析
為了更好地理解具體哪個(gè)機制參數對血漿ALT升高更有顯著(zhù)意義，采用R進(jìn)行了多重線(xiàn)性回歸分析。由于SimPops中不同參數的數量級處于一個(gè)較大的范圍內，所以將SimPops中的每個(gè)參數值除以其最大值進(jìn)行歸一化處理，使得所有參數的范圍都落在0-1之間。
3 模型結果與分析
3.1 大鼠和人PBPK模型的優(yōu)化
大鼠的PBPK模型是根據靜脈及口服給藥的觀(guān)測結果進(jìn)行優(yōu)化，優(yōu)化過(guò)程中直接使用了實(shí)驗值或者考察一定范圍內變化的影響，最終得到的血漿濃度-時(shí)間曲線(xiàn)和觀(guān)測結果較為一致，如下圖所示：
人體PBPK子模型也根據口服給藥300 mg后的血藥濃度-時(shí)間數據進(jìn)行優(yōu)化，優(yōu)化后的模型結果也可以和單次及14天多次給藥的結果吻合一致，如下圖所示：
此外，還進(jìn)一步得到了大鼠及人肝臟平均的濃度-時(shí)間曲線(xiàn)，如下圖。模擬的大鼠肝臟與全血濃度比與大體體內定量全身放射自顯影技術(shù)(QWBA)的數據進(jìn)行了對比，以確定模擬的肝臟濃度與測定數據一致。當然，沒(méi)有獲得相應人體肝臟的濃度數據，因此預測的肝臟濃度與體外肝細胞實(shí)驗測定細胞內濃度進(jìn)行了對比。
3.2 大鼠肝毒性的模擬
按照制定的模擬方案，整合大鼠PBPK子模型以及兩個(gè)推測的毒性機理(線(xiàn)粒體損傷及膽汁酸轉運體抑制)進(jìn)行了大鼠的群體模擬，模擬結果發(fā)現與臨床前結果一致，未發(fā)現PF‐04895162具有肝毒性，如下圖4 (評估藥物誘導嚴重肝毒性, eDISH，該圖是典型地評估化合物臨床肝臟安全性的方式)。還進(jìn)一步研究了膽汁酸和ATP的動(dòng)力學(xué)變化，發(fā)現肝毒性的膽汁酸鵝去氧膽酸-氨基化合物(CDCA-amide)的水平未發(fā)生較大影響、肝臟ATP的水平僅有輕微的下降，如下圖5。
3.3 人體肝毒性的模擬
人體模擬的eDISH圖顯示PF‐04895162對人具有肝毒性。人體SimPops的預測結果：ALT > 1x ULN的比例為21% (59/285), ALT > 5x ULN的比例為11% (32/285), ALT > 10x ULN的比例為6% (18/285); 膽紅素上升并超出2x ULN的比例為3% (9/285)；而臨床觀(guān)測的結果為：ALT > 1x ULN的比例為75% (6/8), ALT > 5x ULN的比例為12% (1/8), ALT > 10x ULN的比例為0% (0/8)；此外也觀(guān)測到兩例受試者的總膽紅素水平介于1 x 和2x ULN。因此當前大鼠與人體的肝毒性模擬結果與觀(guān)測數值基本吻合；和I期多劑量給藥的少數受試者結果相比，模擬的結果在A(yíng)LT輕微變化階段的比例要少一些，另外嚴重肝損傷的比例要高一些。
模擬發(fā)現有9例個(gè)體(9/285, 3%)將發(fā)生嚴重的肝損傷(ALT > 3x ULN, total bilirubin > 2x ULN)，這個(gè)模擬指標對應臨床嚴重肝損傷的響應為Grade 3 ALT升高(5-10x ULN)，且不伴隨膽紅素的升高。通過(guò)多重線(xiàn)性回歸分析發(fā)現SimPops中的體重是PF-04895162體內暴露的協(xié)變量；此外膽汁酸微管外排能力的降低將顯著(zhù)影響血漿中ALT的水平。
DILIsym模擬的血漿ALT水平峰值時(shí)間(下圖7)和臨床結果也較為一致：臨床試驗中，受試者的ALT水平峰值時(shí)間分布于6.2到21天(大部分介于14到19天)，平均時(shí)間為16天；人SimPops的模擬結果分布于7.7到16.8天，平均時(shí)間為14.2天。
3.4 人體肝毒性的機理探討
進(jìn)一步機制性地研究了人體模型的ATP和膽酸鹽的動(dòng)力學(xué)，發(fā)現與大鼠結果相比，肝臟ATP的水平呈現顯著(zhù)地下降、膽汁酸CDCA-amide也顯著(zhù)上升。結果見(jiàn)下圖8所示：
此外還選擇了小樣本量的SimCohort (n=16)，考察ETC抑制或膽汁酸轉運體抑制關(guān)閉的情況下，對肝毒性預測結果的影響，以具體分析造成人肝毒性的主要機制，結果見(jiàn)表4。比較有趣的是，當去除PF-04895162介導的ETC抑制或轉運體抑制兩者中的一個(gè)機制，模擬結果都無(wú)法展現ALT水平的升高。這些研究表明，PF-04895162在抑制ETC的同時(shí)，也抑制BSEP，從而共同產(chǎn)生肝臟毒性(即使化合物對BSEP的親和力較弱，IC50=311 μmol/L)。
4 模型討論
肝臟暴露對毒性的影響：肝臟的藥物暴露量會(huì )影響肝毒性的預測結果，當前模型預測的人體肝臟暴露高于大鼠。雖然采用了不同的數據對大鼠及人肝臟暴露進(jìn)行了驗證，但似乎也可能人體肝臟實(shí)際暴露會(huì )低于模型預測結果(導致人體模型肝毒性的結果出現假陽(yáng)性)，此種情況存在可能性但又無(wú)法真實(shí)獲知。由于存在大鼠體外肝細胞內的濃度數據，因此搭建了大鼠的PBPK模型(同時(shí)吻合觀(guān)測血漿及肝細胞內的濃度，肝細胞內濃度高于體內QWBA的數值)，并基于此模型考察了大鼠肝臟濃度(高估)以及肝毒性的情況，結果發(fā)現該化合物依然對大鼠沒(méi)有明顯的肝毒性。
人體低劑量給藥后，未發(fā)現肝毒性：模型也考察了200 mg BID，共7天給藥的肝毒性情況(一個(gè)24個(gè)受試者的小臨床試驗采用了該給藥方案，且未觀(guān)察到ALT超過(guò)3x ULN)，人體SimPops的結果提示絕大部分(283/285)都無(wú)肝毒性特征，且兩個(gè)個(gè)體的ALT超過(guò)了3x ULN。此外，170 mg BID，共7天給藥的模擬結果，發(fā)現無(wú)任何個(gè)體的ALT大于3x ULN。這些模擬的結果與臨床觀(guān)測現象基本一致，同時(shí)也說(shuō)明了該化合物劑量依賴(lài)性的肝毒性特征。
模型表明ALT有滯后現象，這反映PF-04895162介導膽汁酸轉運體抑制后，初始不能有效積累膽汁酸以引起毒性。此外，模擬還發(fā)現在不加入自適應免疫響應時(shí)，PF-04895162給藥后將近2周才達到ALT水平的峰值。這些滯后的毒性說(shuō)明PF-04895162可能表現為固有型肝損傷。
模型雖然能重現PF-04895162種屬差異性的肝損傷，以及肝毒性的滯后現象，但ALT變化比例、升高的數值等也有觀(guān)測結果有部分差異。除了當前模型已排除的一些機理，也可能存在部分DILIsym未納入的機理；此外對比的觀(guān)測數據來(lái)自樣本量較少(n=8)小臨床試驗；同時(shí)當前的模型也未添加自適應免疫響應(體外缺少定量表征該機理的數據)；最后當前的模擬也未考慮臨床監測的情況：不管是否出現了肝損傷的信號，模型都將繼續運行直到肝細胞的損失量超過(guò)80%以上。
5 總結
該工作采用PK數據和考察肝毒性機理的體外數據，重現了PF-04895162種屬差異性的肝毒性結果，包括表征ALT水平升高的滯后現象。機制性的模型說(shuō)明該化合物抑制膽汁酸轉運體引起細胞內膽汁酸的積累，進(jìn)而誘導線(xiàn)粒體的毒性；此外還連同了ETC抑制破壞線(xiàn)粒體的功能(兩者共同作用，其中一個(gè)都無(wú)法造成ALT水平的變化)。該工作也為后續類(lèi)似的案例模擬提供了研究基礎，通過(guò)DILIsym軟件這種QST的機制性模型可以很好地評估藥物肝毒性的特征。
6 應用軟件與模塊
該案例應用的軟件是DILIsym (v8A)，并基于此版本設定大鼠和人體人群(SimPops)的參數。
關(guān)于DILIsym軟件
DILIsym是基于定量系統毒理學(xué)的計算機模型軟件，運用微分方程來(lái)概括藥物可能導致肝損傷的關(guān)鍵過(guò)程，通過(guò)改變參數來(lái)創(chuàng )建模擬的患者群體或動(dòng)物，以反映遺傳和非遺傳性變異，現已廣泛用于評估藥物性肝損傷的特性。軟件整合多個(gè)子模型(如肝細胞生命周期、線(xiàn)粒體功能損傷和毒性、膽汁酸的處置、生物標記物的釋放等)到模擬過(guò)程中。此外，DILIsym還引入了生理藥代動(dòng)力學(xué)(PBPK)模型以表征藥物的代謝與處置過(guò)程，反映體內暴露情況。DILIsym軟件是由DILI-sim Initiative開(kāi)發(fā)出來(lái)的，該群體是一個(gè)開(kāi)放性、民辦合作組織，匯集了學(xué)術(shù)、制藥企業(yè)及US FDA等業(yè)內的科學(xué)家。DILIsym已經(jīng)在多個(gè)上市前及上市藥物中進(jìn)行了測試，其在80%以上的預測病例中都是準確的。
參考文獻
Grant Generaux, Vinal V. Lakhani, Yuching Yang, et al., Quantitative systems toxicology (QST) reproduces species differences in PF‐04895162 liver safety due to combined mitochondrial and bile acid toxicity. Pharmacol Res Perspect. 2019 Oct 9; 7(6): e00523. IF=2.052
延伸閱讀
1. 視頻 | 采用DILIsym軟件評估藥物性肝損傷DILI