卡馬西平在癲癇患兒中的群體藥動(dòng)學(xué)研究Δ

焦洋，廖建湘，焦正，黎曙霞（.深圳市兒童醫(yī)院，深圳市5806；.復(fù)旦大學(xué)附屬華山醫(yī)院，上海市 00040；.中山大學(xué)第一附屬醫(yī)院，廣州市 50080）

卡馬西平藥動(dòng)學(xué)復(fù)雜，個(gè)體差異大且治療濃度范圍窄，因此臨床上需要對(duì)該藥進(jìn)行監(jiān)測(cè)和個(gè)體化給藥。非線性混合效應(yīng)模型（Nonlinear Mixed Effect Model，NONMEM）[1]，可對(duì)臨床常規(guī)監(jiān)測(cè)的稀疏數(shù)據(jù)進(jìn)行數(shù)學(xué)分析，得到群體藥動(dòng)學(xué)參數(shù)。該法1999年被美國(guó)食品與藥品管理局（FDA）正式認(rèn)可，作為推薦的藥物研究和開發(fā)手段之一。本研究以單用或合用卡馬西平的中國(guó)兒童癲癇患者為研究對(duì)象，運(yùn)用NONMEM法定量考察卡馬西平清除率的影響因素，估算相對(duì)清除率等群體藥動(dòng)學(xué)參數(shù)，并擬訂給藥方案，協(xié)助臨床用藥。

1 儀器

NONMEM程序:VersionⅥ2.0（美國(guó)加州大學(xué)NONMEM工作組）；計(jì)算機(jī)（Intel Core（TM）263001.86 GHz，內(nèi)存0.99GB，Microsoft中文Windows XP Home SimpChn 672）。

2 方法

2.1 研究對(duì)象

深圳兒童醫(yī)院門診癲癇患者，年齡、性別、民族、身高、體重不限；肝、腎功能和白細(xì)胞、血小板計(jì)數(shù)正常；能按照臨床治療方案服用卡馬西平（商品名:得理多，北京諾華制藥），單用或合并用藥皆可，且連續(xù)用藥4周以上，血藥濃度已達(dá)到穩(wěn)態(tài)者。

2.2 樣本收集

按臨床治療需要，于早晨服藥前空腹采集血樣，并準(zhǔn)確記錄患者生理、病理資料及用藥情況（見表1）。

表1 患兒一般統(tǒng)計(jì)學(xué)資料Tab 1 General information of pediatric patient

2.3 檢測(cè)方法

熒光偏振免疫分析法（TDx，Abbott Co.Irving，Texas，USA），線性范圍為0.5～20 mg·L－1，RSD＜10%。

2.4 模型的建立

2.4.1 藥動(dòng)學(xué)模型:根據(jù)文獻(xiàn)報(bào)道[2]，卡馬西平的藥動(dòng)學(xué)模型采用一級(jí)吸收和消除的一房室開放模型:

公式1、2中，cp為模型預(yù)測(cè)的血藥濃度，c0為初始濃度，e為自然對(duì)數(shù)，Ka為吸收速率常數(shù)，K為消除速率常數(shù)，D為給藥劑量，CL為清除率，Vd為表觀分布體積。Ka、CL、Vd為該模型的基本參數(shù)，且其中的CL和Vd值包括了制劑的生物利用度F值，即CL＝CL/F，Vd＝Vd/F。

由于本研究中所有的數(shù)據(jù)均來自消除相，因此用NONMEM法無法估算Ka，也給估算Vd帶來困難[3～7]。在模型擬合時(shí)，根據(jù)文獻(xiàn)報(bào)道值[2]，將Ka固定在0.5 h－1，然后用NONMEM法進(jìn)行擬合。

2.4.2 固定效應(yīng)模型:首先建立不含任何固定效應(yīng)參數(shù)且各參數(shù)間相對(duì)獨(dú)立的最簡(jiǎn)模型（即視所有個(gè)體的清除率都相等的模型）。根據(jù)NONMEM對(duì)最簡(jiǎn)模型的擬合結(jié)果，觀察各影響因素與CL關(guān)系的散點(diǎn)圖，用加法、乘法或指數(shù)模型等逐步添加解釋性變量（Covariate）。由于目標(biāo)函數(shù)的差值（Difference in Objective Functions，DOBJ）近似χ2分布，當(dāng)χ2＞3.84時(shí)（df＝1，P＜0.05），統(tǒng)計(jì)學(xué)上具顯著性差異，逐一添加解釋性變量可得全量回歸模型（Full regression model，F(xiàn)RM）。用似然比檢驗(yàn)法（Likelihood ratio test）對(duì)各解釋性變量進(jìn)行檢驗(yàn)，方法如下:定檢驗(yàn)水準(zhǔn)α＝0.01，固定FRM中某一固定效應(yīng)參數(shù)，得限制性回歸模型（Restricted regression model，RRM）。RRM和FRM之差DOBJ＞6.64時(shí)，則該固定效應(yīng)參數(shù)對(duì)CL的群體典型值有顯著性影響，反之為干擾因素，予以剔除。對(duì)FRM重新整理、簡(jiǎn)化后可得最終回歸模型。

2.4.3 統(tǒng)計(jì)學(xué)模型:統(tǒng)計(jì)學(xué)模型中的個(gè)體間變異模型采用比例模型，公式為:

個(gè)體自身變異模型（亦稱殘差模型）采用加法模型，公式為:

公式3中CLj為第j個(gè)患者的清除率；CL為群體典型值，即能代表患者群體特性的清除率；ηCL，j是均數(shù)為0、方差為ω2CL的個(gè)體間變異。公式4中c0ij為血藥濃度實(shí)測(cè)值；cij為血藥濃度預(yù)測(cè)值；εij是均數(shù)為0、方差為σ2的殘差變異。

2.5 模型的驗(yàn)證

本研究采用自舉驗(yàn)證法（Bootstrap）[4]，該法特別適合那些難以用常規(guī)方法導(dǎo)出對(duì)參數(shù)的區(qū)間估計(jì)、假設(shè)檢驗(yàn)等問題。當(dāng)缺少驗(yàn)證數(shù)據(jù)時(shí)或樣本量較少（如兒科患者或因倫理學(xué)原因）時(shí)，Bootstrap具有獨(dú)特的優(yōu)點(diǎn)，目前已被美國(guó)FDA采用。

3 結(jié)果

3.1 最終回歸模型

通過模型假設(shè)檢驗(yàn)（表2、表3），確定體重、服藥劑量為卡馬西平的清除率的影響因子；性別、合用丙戊酸鈉不影響卡馬西平清除率。最終回歸模型如下:

各參數(shù)的運(yùn)算結(jié)果見表4（95%置信區(qū)間＝參數(shù)值±1.96·殘差），最終模型為公式:

上式中各參數(shù)的含義如下:

CL:清除率（單位:L·h－1）；WT:體重（單位:kg）；Dose:日劑量（單位:g·kg－1）。

表2 混合效應(yīng)和目標(biāo)函數(shù)值Tab 2 Fixed effect and objective value

表3 RRM和FRM的目標(biāo)函數(shù)值之差Tab 3 Difference of the objective function between RRM and FRM

表4 最終模型的各群體藥動(dòng)學(xué)參數(shù)值Tab 4 Population pharmacokinetic parameters of carbamazepine by finally model

群體預(yù)測(cè)值與觀測(cè)值的相關(guān)性見圖1；個(gè)體預(yù)測(cè)值與觀測(cè)值的相關(guān)性見圖2。

圖1 群體預(yù)測(cè)值與觀測(cè)值的相關(guān)性Fig 1 Relationship between PRED and DV

圖2 個(gè)體預(yù)測(cè)值與觀測(cè)值的相關(guān)性Fig 2 Relationship between IPRE and DV

3.2 驗(yàn)證

利用NONMEM程序組中已有的Bootstrap軟件寫以下程序命令:NMBS CBZ4DA 2000。其意義為Bootstrap檢驗(yàn)再抽樣共2000次。

再抽樣2000次計(jì)算各參數(shù)的中位數(shù)與FRM的參數(shù)相比相對(duì)誤差如果＜5%[6]，則表明，F(xiàn)RM計(jì)算的參數(shù)是穩(wěn)定可靠的。檢驗(yàn)結(jié)果見表5、表6（其中，相對(duì)誤差（Relative error）＝│FRM-Bootstrap│/FRM×100%）。

表5 自舉驗(yàn)證法結(jié)果Tab 5 Result of bootstrap validation

表6 最終模型參數(shù)值與檢驗(yàn)值比較Tab 6 Comparison of FRM and bootstrap

結(jié)果表明，F(xiàn)RM計(jì)算的參數(shù)是穩(wěn)定可靠的。

4 討論

4.1 劑量對(duì)卡馬西平清除率的影響

本研究與多數(shù)研究的結(jié)果一致:劑量是影響清除率最主要的因素之一（P＜0.001），且隨著卡馬西平劑量的增加，清除率以非線性形式增加。其原因可能與卡馬西平的水溶性差，大劑量時(shí)卡馬西平的吸收程度下降有關(guān)[6]。

4.2 體重對(duì)卡馬西平清除率的影響

一般認(rèn)為，體重是影響表觀分布體積的因素，但體重是否是清除率的獨(dú)立影響因素，不同藥物間各不相同。本研究的結(jié)果顯示，卡馬西平的清除率和體重有關(guān)，且隨著體重的增加以非線性方式增加[6]。

4.3 合用丙戊酸對(duì)卡馬西平清除率的影響

在癲癇的臨床治療中卡馬西平和丙戊酸的聯(lián)用很常見。一般認(rèn)為，卡馬西平能顯著降低丙戊酸的清除率，而丙戊酸對(duì)卡馬西平代謝的影響，各國(guó)學(xué)者看法不一。用NONMEM法對(duì)卡馬西平的群體藥動(dòng)學(xué)研究結(jié)果也不完全相同，有日本學(xué)者認(rèn)為丙戊酸會(huì)稍增加卡馬西平的清除率[2]，國(guó)內(nèi)重慶[7]、上海[6]也認(rèn)為丙戊酸會(huì)稍增加卡馬西平的清除率，而西班牙等[3]的研究顯示丙戊酸不改變卡馬西平的清除率。

本研究的結(jié)果表明，丙戊酸對(duì)卡馬西平的清除率無明顯影響。本研究的119例患兒中合用丙戊酸的有27例，占22.7%。在最簡(jiǎn)模型與丙戊酸的散點(diǎn)圖中，二者無明顯趨勢(shì)性變化。在添加丙戊酸因素計(jì)算目標(biāo)函數(shù)值時(shí)，僅減少2.154。故認(rèn)為丙戊酸對(duì)卡馬西平清除率的無明顯影響。

本文建立了Bootstrap驗(yàn)證法:采用完整數(shù)據(jù)建立模型，然后用再抽樣技術(shù)（Resampling），生成大量驗(yàn)證數(shù)據(jù)對(duì)每一組數(shù)據(jù)計(jì)算其參數(shù)，最后根據(jù)這些參數(shù)觀察模型的可靠性。當(dāng)缺少驗(yàn)證數(shù)據(jù)時(shí)或樣本量較少時(shí)，如兒科患者或因倫理學(xué)原因，自舉法具有獨(dú)特的優(yōu)點(diǎn)，目前已在群體藥動(dòng)學(xué)中作為驗(yàn)證方法，并被美國(guó)FDA采用。

綜上所述，本研究揭示了體重、劑量為卡馬西平清除率的影響因素，性別、合用丙戊酸不干擾卡馬西平的代謝，并且最終模型的各參數(shù)值與文獻(xiàn)報(bào)道接近。

[1]Sheiner LB，Rosenberg B，Marathe VV.Estimation of population characteristicsofpharmacokineticparameters from routine elinieal date[J].Pharmaeokinet Biopharm，1997，25（2）:445.

[2]Yukawa E，Aoyama T.Detection of Carbamazepine Drug Interaction by multiple peak approach screening using routine clinical pharmacokinetic data[J].Clin Pharmacol，1996，36（9）:752.

[3]Hartley R，F(xiàn)orsythe WI，Mclain B，et al.Daily variation in steady state plasma concentrations of carbamazepine and its metabolites in epileptic children[J].Clin Pharmacokinet，1991，20（3）:237.

[4]Efron B，Tibshirani RJ.An Introduction to the Bootstrap[M].London:Chapman＆ Hall.1993.

[5]Efron B，Gong G.Aleisurely look at the bootstrap.The jackknife，and cross-validation[J].Am Stat，1983，37（1）:36.

[6]焦 正，鐘明康.NONMEM法估算中國(guó)癲癇患者卡馬西平的清除率[J].中國(guó)藥理學(xué)通報(bào)，2003，19（3）:314.

[7]王 剛，劉 彬，黃曉英，等.非線性混合效應(yīng)模型法研究卡馬西平在癲癇兒童中的群體藥動(dòng)學(xué)[J].中國(guó)醫(yī)院藥學(xué)雜志，2007，27（7）:1001.