中藥復方藥物動力學總量統計矩法的實驗驗證研究

賀福元,鄧凱文,劉文龍,石繼連,伍勇,劉偉,賀慶平,李博

[摘要] 目的：對已建立的藥物動力學總量統計矩法用補陽還五湯中黃芪甲苷、芍藥苷、川芎嗪3個成分進行驗證，建立中藥復方多成分體系的藥物動力學實驗方法。方法：采用反相高效液相色譜法測定，其條件為C18色譜柱（4.6 mm×250 mm，5 μm）；黃芪甲苷ELSD測定；流動相乙腈1%乙酸溶液（35∶65），流速1 mL·min-1；柱壓 （15.0±0.408） MPa；溫度30 ℃；芍藥苷、川芎嗪的檢測波長254 nm，流動相甲醇水（33∶67），流速1 mL·min-1；柱壓（20.71±0.42） MPa。單個成分房室模型藥物動力學參數采用DAS軟件處理，總量統計矩按公式計算。結果：3個單成分在大鼠體內藥物動力學都遵循二室模型（P<0.01），與疊加總濃度比，單個成分各參數相差最大上萬倍，而總量統計矩參數的RSD為3.5%，3個成分藥物動力學的總量統計矩參數AUCt （119.8±27.20） g·min·L^-1；MRTt （210.0±54.49） min；VRTt （5.608±2.723）×104 min2；CLt為（0.319 6±0.068 8） mL·min-1·kg-1； Vt為（64.12±8.243） mL·kg-1； t0.95t為（588.9±149.4） min。這就說明3個成分總量的半衰期為（145.5±37.76） min，3藥在0～674.2 min內代謝95%的濃度。結論：中藥復方多成分藥物動力學可采用總量統計矩法進行研究，該法能表征多成分體系總量的量時變化的動力學規律。

[關鍵詞] 總量統計矩；藥物動力學；譜動學；補陽還五湯；黃芪甲苷；芍藥苷；川芎嗪

中藥及復方藥物動力學研究是闡明中藥復方產生作用物質基礎的關鍵問題，是闡明機體對中藥及復方量時效關系的基礎，也是中藥復方制成控緩釋制劑的理論基礎[1]。由于中藥復方為多成分體系，如按單成分藥物動力學研究將會存在以下問題[25]：①怎樣建立適宜多成分藥物動力學研究的數學模型，并對參數進行測定；②在對照品缺少的情況下，中藥復方多成分怎樣進行藥物動力學研究；③怎樣建立宏觀綜合效應與總成分、單成分濃度對應關系。針對這些問題，作者前期創立了總量統計矩法，并有論文發表和成果獲獎[6]，但由于公式推導較復雜，參數的理論與實驗測定問題沒有深入闡明清楚，可能會影響同行在進行中藥復方藥物動力學研究時采用難度，為此，本論文先從第一個問題著手，詳細闡明建模的思路及參數的意義及測定方法，再逐步報道其相關問題[7]。

補陽還五湯系清代王清任《醫林改錯》首載，該方由黃芪、當歸、芍藥、川芎、桃仁、紅花、地龍組成，具有補氣、通絡的功能。在抗血栓、抗腦栓血、抗衰老、抗血脂及免疫功能方面都有廣泛的作用[8]。本室對該方進行了大量的抗腦血栓作用機制研究，發現該方具有重要的抗血栓活性作用。其主要有效成分部位為總苷、總苷元、總生物堿和總多糖，其代表性成分為黃芪甲苷、芍藥苷、川芎嗪等等，本文選這3個成分作為代表性成分來闡明中藥復方多成分如何進行以總量統計矩數學模型為基礎的藥物動力學研究。

1 多成分藥物動力學的建模數學思路

多成分的藥物動力學研究是回答各成分濃度隨時間量變規律，總體上有2種思路，一種是根據藥物濃度具有疊加性，可將各單成分濃度直接相加構成總量藥物濃度，然后再按單成分藥物動力學方法進行參數求算，這種方法能獲得總量藥物動力學參數，但會丟失各成分的藥物動力學參數，沒有建立起總量藥物動力學參數與單個成分藥物動力學參數的關系；第二方法是進行函數曲線轉換后再疊加處理，先將單成分藥物動力學曲線轉變成一種公用函數曲線，再進行曲線疊加，最常用的方法就是泰勒（Tatlor）轉換疊加；傅利葉（Fourier）轉換疊加等[9]，但對轉換后的疊加函數不能直接給有價的藥物動力學參數，因此較宜采用統計矩轉換疊加構成總量統計矩參數，與單成分統計矩參數一樣，能表征多成分宏觀的藥物動力學行為。由于單成藥物動力學的統計矩模型及參數已被熟悉，其總量統計矩參數也易為大家所接受，為此可以先將單成分的藥物動力學統計矩化，亦正態分布曲線化，然后再進行疊加構成總量統計矩參數，亦總量正態分布曲線化，由此可得總量零、一、二階矩參數。因此總量統計矩參數既能反映總量藥物濃度的動力學行為，又能與單個成分的統計矩參數關聯，可用于房室、非房室及不規則成分的藥物動力學研究，是整合多成分藥物動力學較宜的數學處理方法[6]。

1.1 單成分的統計矩 將單成分的藥物動力學曲線（參數）轉變成統計矩曲線（參數）。其思路見圖1。單成分藥時曲線的統計矩化過程實際上就是將藥物動力學曲線轉變成正態分布曲線，其中重要參數有零階矩（AUC）、一階矩（MRT）和二階矩（VRT）。

線性房室模型的統計矩參數：按統計矩原理可得到單成分零、一、二階統計矩參數通式，為式（1），（2），（3）。

由上諸式可知：設M_i表示r_i室模型第i項e的指數項前面的系數，αi表示第i項e的冪指系數。一、二、三室成分的各總量統計矩參數意義及計算參見文獻[6]。非線性模型的統計矩參數：按統計矩原理可得到非線性單成分零、一、二階統計矩參數通式，為式（5），（6），（7）。

各總量統計矩參數意義及計算參見文獻[6]。

不規則模型的統計矩參數：按統計矩參數定義，用梯形法或拋物線法計算零階矩（AUC），按文獻的先計算AUMC（時間與曲線乘積下的面積），再計算MRT，再計算AUVC（時間平方與曲線乘積下的面積），再計算VRT[2]。

1.2 多成分體系的藥動學總量統計矩參數的設定及與單個成分指標的關系 根據統計矩原理進行加合處理，見圖2，圖中a分別為2個成分的統計矩曲線，經統計矩原理加合計算后可得圖中b的總量統計矩曲線，由此可得總量統計矩參數：總量零階矩AUCt、總量一階矩MRTt、總量二階矩VRTt。假設中藥復方為一個由n個單體成分群構成體系，其中m個成分服從線性藥物動力學。第i個成分服從ri室模型。S1個成分服從非線性藥物動力學。S2個為不規則模型，該中藥復方成分總量的零階矩、一階矩、二階矩及由此衍生的一系列動力學參數概念如下。

1.3 中藥復方的總量統計矩參數 包括總量統計零、一、二階矩以及由此衍生的其他總量統計矩參數。

總量零階矩（AUCt）：根據統計矩的加合性可計算中藥復方多成分的總量統計矩，為下式（8）。

式中m表示線性成分的個數，r_i 為第i成分代謝的室數，s_i為非線性成分的個數，s_i為產生多峰不規則成分個數。

總量一階矩（MRTt）：即為總量平均駐留時間，亦為各成分的總體均值，可定義如下式（9）。

總量二階矩（VRTt）：亦可采用相同方法處理，可得式（10）。

由上諸式可知，一個復方的總體藥物動力學趨勢由單個成分的動力學參數決定。其總量零階矩等于各單成分的零價矩之和；總量一階矩等于各單個成分一階矩對零階矩的算術平均值；總量二階矩等于各單個成分的一階矩的平方與二階矩之和對零階矩的算術平均值，再與總量一階矩平方之差。這些性質是由總量統計矩的加合性決定的，亦有式（11），（12）。

式中MRT _p1為1～i成分的總量一階矩，VRT _p1為為對應的總量二階矩為；MRT _p2為i～n成分的總量一階矩，VRT _p2為為對應的總量二階矩。由總量零、一、二階矩可衍生出其它總量藥物動力學參數。

總量表觀半衰期、清除率和表觀分布容積：根據統計矩與半衰期有關系，中藥復方多成分的總量表觀半衰期為式（13）。

95%總量與累積95%總量濃度時間區間：由于中藥復方為多成分體系，用零階矩只表征其總成分的生物利用度，用一階矩只表征諸成分平均駐留時間中心，這還不能完整表征多成分的藥物動力學行為，因為平均駐留時間中心相同的一組成分其藥物的種類與離散程度可以不同，也就是95%總量與累積95%時間區間不同，因此將一、二階矩結合可表征95%總量濃度時間區間，亦為式（17）。

這樣，式（8）～（19）構成了中藥復方多成分藥物動力學總量統計矩數學模型及參數體系[6，10]。

2 材料與方法

2.1 儀器與試藥 Agilent1200液相色譜儀，Alltech 3300 ELSD檢測與UV測定。瑞士Precisa公司的 20 A型電子天平；德國Heraeus公司的高速冷凍離心機及SZ1型快速混勻器。黃芪甲苷（批號110781200512）、芍藥苷（批號110736200629）、川芎嗪（批號110817200305）（購于國家藥品生物制品檢定所），甲醇、乙腈為色譜純，蒸餾水為二次重蒸餾水，其他試劑均為分析純。其復方藥材來源為：黃芪為豆科植物蒙古黃芪Astragalus membranaceus Bge.var. monghokicus （Bgo） Hsiao 干燥根；當歸為傘形科植物當歸Angelica sinensis （Oliv） Diels的干燥根；川芎為傘形科植物川芎Ligusticum chuanxiong Hort的干燥根；赤芍為毛莨科植物芍藥Paeonia lactlora Pall的干燥根莖；桃仁為薔薇科植物物桃Prunus persica （Linn） Batsch的干燥成熟種子；紅花為菊科植物紅花Carthamus tinctorius L.的干燥花；地龍為環節動物門鉅蚓科動物參環毛蚓Pheretima aspergillum （Perrier）的干燥主體。處方黃芪60 g，當歸9 g，川芎6 g 赤芍9 g，桃仁9 g， 紅花9 g，地龍9 g。Wistar大白鼠，雌雄各半，體重220～300 g，由湖南中醫藥大學實驗動物中心提供，許可證號SCXK（湘）20012004。

2.2 色譜條件 色譜柱Ultimate XBC18柱（4.6 mm×250 mm，5 μm），黃芪甲苷[11]采用Agilent 1200 Alltech 3300 ELSD測定，流動相乙腈1%乙酸溶液（35∶65），流速1 mL·min-1，柱壓（15.0±0.408） MPa；溫度30 ℃。芍藥苷與川芎嗪采用Agilent UV測定，由于兩者同時以指紋圖譜的形式測定，可在原測定波長的基礎上調整為254 nm[1213]，流動相甲醇水（33∶67），流速1 mL·min-1，柱壓（20.71±0.42） MPa；溫度35 ℃。

2.3 給藥與取血 將補陽還五湯水提醇沉制得的復方注射液（分別含黃芪甲苷2.160 g·L^-1，芍藥苷2.578 g·L^-1、川芎嗪25.12 mg·L^-1）。實驗大鼠尾靜脈給藥2 mL，分別于給藥后0.25，0.5，1.25，1.75，2.5，4.5，5.5，8，17.5 h，眼眶取血漿400 μL，并在1 h內進行樣品預處理。

2.4 樣品預處理 將補陽還五湯待測血漿400 μL置于尖底塑料刻度離心管中，加入三氯醋酸液1 mL，渦旋混合5 min，5 000 r·min-1離心10 min，取三氯醋酸層，常溫下氮氣流吹干，殘渣用甲醇50 μL溶解，取20 μL進樣。

2.5 空白血漿處理 于實驗前，取大鼠血樣按樣品處理項下同法制得空白血漿樣。

2.6 對照品液制備 精密稱取黃芪甲苷對照品8.0 mg，用甲醇定容至10 mL，制得每1 mL含黃芪甲苷0.8 mg的對照品液；精密稱取芍藥苷對照品3.2 mg，用甲醇定容至10 mL，制得每1 mL含芍藥苷0.32 mg的對照品液；精密稱取川芎嗪對照品5.0 mg，用甲醇定容至100 mL，再稀釋10倍制得每1 mL含川芎嗪0.005 mg的對照品液。

2.7 數據處理與分析 結果用DAS 2.0軟處理，分別得單個成分、3成分疊加及平均的經典藥物動力學參數，按自建的統計矩法分析得藥物動力學總量統計矩參數。

3 結果

3.1 保留時間 黃芪甲苷：血樣經預處理后，在黃芪甲苷色譜條件下進樣測定，得血漿色譜圖，見圖3。大鼠血漿中黃芪甲苷峰在保留時間為15 min左右可完全分離。芍藥苷與川芎嗪：在芍藥苷、川芎嗪對照品液的色譜條件下測得芍藥苷、川芎嗪的血漿色譜圖。大鼠血漿補陽還五湯中芍藥苷峰在保留時間為11.5 min左右，川芎嗪峰在保留時間為21.5 min左右可完全分離[2]。

3.2 標準曲線的制備及線性關系的考察 分別取試管6只，加入802.1 mg·L^-1黃芪甲苷甲醇溶液對照品液10， 25， 50， 100， 200， 400 μL，水浴蒸干，再加入空白血漿400 μL，經樣品預處理后，殘渣用50 μL甲醇溶解，使血漿中黃芪甲苷的質量濃度分別為20.05～802.1 mg·L^-1；另取試管6只，分別加入320.4 mg·L^-1芍藥苷甲醇溶液對照品液1，15，30，60，75，100 μL，水浴蒸干，再加空白血漿400 μL，經樣品預處理后，殘渣用50 μL甲醇溶解，使血漿中芍藥苷的質量濃度分別為0.801～80.10 mg·L^-1；又取試管6只，分別加入5.140 mg·L^-1川芎嗪甲醇溶液對照品液2，5，10，15，20，40 μL，水浴蒸干，另加入空白血漿400 μL，經樣品預處理后，殘渣用50 μL甲醇溶解，使血漿中川芎嗪的質量濃度分別為25.70～514.0 μg·L^-1。在前述色譜條件下進行測定，以黃芪甲苷濃度為縱坐標，黃芪甲苷峰面積為橫坐標進行線性回歸處理，得回歸方程為C=0.183 2A-0.01，結果表明黃芪甲苷在20～800 mg·L^-1呈良好線性（r=0.999 7），黃芪甲苷的最低檢測濃度為0.6 mg·L^-1 （S/N>3）；以芍藥苷面積為縱坐標，芍藥苷濃度為橫坐標進行線性回歸處理，得回歸方程為A=178.496 12C-3.361，結果表明芍藥苷在0.80～80 mg·L^-1呈良好線性（r=0.999 9），芍藥苷的最低檢測濃度為0.12 mg·L^-1（S/N>3）；以川芎嗪面積為縱坐標，川芎嗪濃度為橫坐標進行線性回歸處理，得回歸方程為A=64.707 91C+86.325 75，結果表明川芎嗪在25.70～514.0 μg·L^-1呈良好線性（r=0.999 7），川芎嗪的最低檢測質量濃度為8 μg·L^-1 （S/N>3）。

3.3 精密度和回收率精密度 按文獻分低、中、高3個濃度進行日內和日間精密度、加樣回收率試驗研究。

精密度試驗：取大鼠空白血漿400 μL，加入黃芪甲苷、芍藥苷、川芎嗪對照品液，配制低、中、高3個濃度的樣品，經預處理后分別測定日內精密度和日間精密度，結果見表1。

回收率試驗：取空白血漿，加入一定量的黃芪甲苷、芍藥苷、川芎嗪對照品液，分別配制高、中、低3個不同濃度的樣品，經預處理后測定，按回歸方程計算測得量，得到的平均方法回收率列為表2。由表1～2可知，本法能用于3個成分的藥物動力學研究中的濃度測定。

3.4 藥物動力學測定 取大鼠100只，隨機分為10組，每組10只，每一組為一實驗點，一組為空白。實驗前夜禁食12 h，尾靜脈給補陽還五湯注射液2 mL，定時取血與處理，可得藥時曲線，見圖4。結果用DAS 2.0軟處理，分別得單個成分、3個成分疊加及平均的藥物動力學參數，見表3，按自建的統計矩法分析得藥物動力學總量統計矩參數，見表4。

3.5 補陽還五湯中黃芪甲苷、芍藥苷、川芎嗪藥物常規藥物動力學及總量統計矩參數計算 先用DAS軟件進行曲線擬合，得表3中常見的藥物動力學參數，復方中黃芪甲苷、芍藥苷、川芎嗪3個成分在大鼠體內藥物動力學都遵循二室模型（P<0.01），再依式（1）～式（4）的計算公式對黃芪甲苷、芍藥苷、川芎嗪單個成分的統計矩參數。如沒有進行曲線擬合，也可直接用梯形法或拋物線法按有關統計矩公式進行計算，然后運用式（5）～式（11）將3藥的決量統計矩進行加合計算得總量藥物動力學統計矩參數，在計算過程一定要注意只有AUC（零階矩）、AUMC（時間與曲線乘積下的面積）、AUVC（時間平方與曲線乘積下的面積）能進行相加，而MRT，VRT不能相加，只能按公式MRT=AUMC/AUC，VRT=AUVC/AUC-MRT2計算，其他參數在MRT參數基礎上進行計算。用上法算得補陽還五湯中3個成分的統計矩參數列表4，其中總量統計矩參數AUCt為（119.8±27.20） g·min·L^-1；MRTt為（210.0±54.49） min；VRTt為（5.608±2.723）×104 min2；t1/2為（145.5±37.76） min；CLt為（0.319 6±0.068 82） mL·min-1·kg-1；Vt為（64.12±8.243） mL·kg-1；t0.95t為（588.9± 149.4） min。

3.6 補陽還五湯中黃芪甲苷、芍藥苷、川芎嗪3個成分95%總量濃度與累積95%總量濃度代謝時間區間 按式（14）計算得3藥的95%總量濃度在[210.0-1.96（5.608×104）1/2，210.0+1.96（5.608×104）1/2] min，亦在-254.2～674.2 min時間內代謝95%的濃度。按式（16）算得t0.95t為588.9 min，故累積95%總量濃度代謝時間區間為[-∞，588.9] min。

3.7 補陽還五湯中黃芪甲苷、芍藥苷、川芎嗪的藥動學參數的比較 以3個成分疊加擬合而得的藥物動力學參數（P<0.01）為標準對照組，表3中其他組與之的百分比值列于表5。

表3諸成分參數比較可采用等級資料的χ²檢驗，其χ²=279 8χ²_0.01（36）=58.619，說明諸成分參數差異顯者（P<0.01），由表5可知，以川芎嗪的參數B的比值最小，為0.050%；以芍藥苷的K₁₀的比值最大，為2 883%，兩者相差6萬倍，黃芪甲苷、芍藥苷、川芎嗪組的RSD分別為18.32%，172.4%，169.7%，這說明中藥復方中各成分的藥物動力學參數相差懸殊，彼此之間參數不能相互代替。平均組與疊加組的χ²=165.7χ²_0.01（9）=21.666，其RSD也為117.9%，遠大于一般藥物動力學實驗誤差5%，因此不能用各成分藥物動力學參數的平均值來代表總量的藥物動力學參數，也就是說總量的藥物動力學參數不是各成分藥物動力學參數的簡單相加。

同理，以疊加成的3個成分統計矩參數為對照組，表4中其他的總量統計矩參數與之百分比值列于表6。

表4諸成分的統計矩參數的χ²=1 962χ²_0.01（25）=44.314，說明諸成分與方法的總量統計矩參數差異顯著（P<0.01）。由表6可知，以川芎嗪的參數AUCt的比值最小，為0.054 6%；以芍藥苷的CLt的比值最大，為1 080%，兩者相差1.7萬倍，黃芪甲苷、芍藥苷、川芎嗪組的RSD分別為31.63%，174.2%，113.6%，這也說明中藥復方中各成分的藥物動力學總量統計矩參數相差懸殊，彼此參數之間不能相互代替。平均組與疊加組的χ²=212.3χ²_0.01（5）=15.086，其RSD也為108.3%，也大于一般藥物動力學實驗誤差5%，因此也不能用各成分統計矩參數的平均值來代表總量統計矩參數，但按式（6）～（13）總量統計矩公式計算的各參數的RSD為3.510%，小于一般藥物動力學實驗誤差5%，同時兩者的χ²=1.006 1χ²_0.01（5）=15.086，因此只有按總量統計矩計算方法對各成分的統計矩參數進行加合計算，所得結果才能與成分疊加的總量的統計矩參數一致，也就是說本文所建立的式（1）～（13）的總量統計矩的數學模型能客觀真實地反映各單個成分總量的藥物動力學行為。

4 分析與討論

4.1 中藥及復方多成分藥物動力學研究 中藥及復方為多成分復雜體系，其產生的整體效應為諸有效成（輔助成分）在體內作用的綜合表現。對于中藥多成分的藥物動力學研究目前仍采用單一成分的房室模型研究方法，這必然產生前述3個問題，其中最關鍵的問題是怎樣表達多成分藥物動力學的數學模型和多成分與效應的關聯性問題[14]，前文解決中藥及復方藥物動力學數學模型問題，本文采用實驗進行驗證。根據統計矩原理可知[4，15]，統計矩為非房室模型，能表達線性房室、非線性MichaelisMenten和不規則曲線模型，能用已知曲線通過統計矩公式直接計算，也可通過各取血點用梯形法或拋物線法進行計算，故能用來表達中藥及復方單成分與多成分藥物動力學整體行為：由總量零、一、二階矩及衍生的總體表觀清除率、總體表觀分布容積構成主要總量統計矩參數體系。由表5，6的對比數據可知，只有按總量統計矩方法對各成分的藥物動力學參數進行加合計算，所得總量統計矩參數結果才能與諸成分疊加的總量的統計矩參數一致，才能反映多成分的整體濃度量變的動力學規律。

4.2 總量統計矩的特點 由統計矩原理和式（5）～（13）可知：總量統計矩具有2大特點[15]：①具加合性，單個成分經統計矩處能構成總量統計矩參數；反之，總量統計矩又能分解成單個成分統計矩；總量統計矩減去一部分成分總量統計矩會得另一部分總量統計矩，利用這一原理可進行血樣空白干擾的消除；②具偶聯性，能與多維變量曲線偶聯，建立起多維曲線下的總量統計矩參數。如與HPLC指紋圖譜與藥動學偶聯得保留時間、藥物動力學代謝時間二維變量對藥物濃度曲線總量統計矩參數；HPLC指紋圖譜與效應偶聯得保留時間、藥物濃度（響應值）二維變量對藥物效應曲線總量統計矩參數，前者稱譜動學、后者稱譜效學[1617]，如HPLC指紋圖譜與藥動學、效應偶聯得保留時間、藥物動力學代謝時間、藥物濃度（峰響應值）三維變量對效應曲線總量統計矩參數，稱譜效動力學。因此，總量統計計矩可用于任何多成分變量的中心和離散性分析，是研究中藥多成分體系動態量變、效應變化的有力的數學工具。

4.3 總量統計矩藥物動力學模型使用時應注意的問題 由式（1）～式（11）的推導過程可知，總量統計矩藥物動力學模型是建立在藥物濃度化學動力學量變的模型的基礎上，其適用性在于能用本模型表達體內各成分的濃度整體量變。不能簡單地用本模型計算出來的總量來推斷藥物生物效應（藥效與毒效）的大小，原因在于各成分的效應與毒效應系數不同，因此多成分的總量效應并不簡單地與總濃度呈相關性，不能用單個成分的濃度相加來表達。總成分濃度量變統計矩參數與總效應統計矩參數并不重迭，只有將各成分的濃度量變乘上效應系數后，將濃度轉化為效應后，才能用總量統計矩計算得出總效應的總效應統計矩參數，才能用來評價藥物總效應，這時總量統計矩對單個成分效應的整合與對單個成分藥物動力學整合一樣，其總效應統計矩參數能指導臨床用，這方而的詳細闡明請讀者參考相關文獻[7]。中藥多成的藥物動力學與效應動力學不對應是有別于西藥單一成分的藥物動力學與藥效動力學對應最為明顯的地方。

本文詳細闡明了中藥復方多成分藥物動力學總量統計矩數學模型的建立原理、實驗與計算方法，適用于已知能測準的中藥代表性成分，但中藥絕大多數成分的沒有對照品或為未知成分，因此怎樣利用指紋圖譜測定中藥復方多成分技術，并與藥物動力學關聯構成的譜動學來測定中藥復方表觀藥物動力學參數，亦譜動學參數將另文詳述。

[參考文獻]

[1] Loew D，Kaszhin .Approaching the problem of bioequivalence of herbal medicinal products[J]. Phytother Res，2002，16（8）：705.

[2] 賀福元.中藥復方藥物動力學數學模型的建立及對補陽還五湯的研究[D]. 成都：成都中醫藥大學，2006：6.

[3] Ludbrook J. Statistical techniques for comparing measurers and method of measurement：a critical review[J].Clin Exp Pharmacol Physiol，2002，29（7）：527.

[4] 魏樹禮.生物藥劑學與藥物動力學[M]. 北京：北京醫科大學&中國協和醫科大學聯合出版社，1997：1204.

[5] Vengpedersin .Noncoupartmentalybased pharmacokinetics modeling[J]. Adv Deliv Rev， 2001，48（2/3）：265.

[6] 賀福元，羅杰英，鄧凱文. 中藥復方動力學數學模型總量統計矩法的研究[J]. 世界科學技術——中醫藥現代化，2006， 9（6）：13.

[7] 賀福元，鄧凱文，鄒歡，等. 中藥復方譜動學與譜效動力學差異性的研究[J]. 中國中藥雜志，2011，36（2）：136.

[8] 陳冬，楊潔紅. 補陽還五湯抗腦缺血作用的研究進展[J].中華中醫藥學刊，2010，28（1）： 72.

[9] 數學手冊編寫組.數學手冊[M]. 8版.北京：高等教育出版社， 2000：219，531.

[10] 段曉鵬，賀福元，周晉，等.補陽還五湯指紋圖譜加合性的研究[J]. 中國中藥雜志，2011，36（23）： 3247.

[11] 李曉亮，邵旭，牛海軍，等. 高效液相色譜——蒸發光散射檢測法測定補陽還五苷酮膠囊中黃芪甲苷的含量[J].安徽醫藥，2010，14（11）：1286.

[12] 鄧凱文，賀福元，劉文龍. RPHPLC法測定大鼠血漿中補陽還五湯及總苷中芍藥苷的含量及藥代動力學研究[J].中藥材，2010，33（11）：1768.

[13] 賀福元，劉文龍，賀慶平，等.反相高效液相色譜法測定大鼠血漿中補陽還五湯總生物堿中川芎嗪的含量及藥代動力學研究[J].中國藥物與臨床，2005，5（12）：920.

[14] 賀福元，羅杰英.中藥復方藥物動力學研究方法現狀及研究方法探討[J]. 中草藥，2005，35（10）：1582.

[15] 賀福元，周宏灝，鄧凱文，等.指紋圖譜的一種定性定量研究新方法：總量統計矩分析法[J].藥學學報，2008，43（2）：195.

[16] 賀福元，羅杰英，鄧凱文，等.中藥譜效學研究方法初探[J] 世界科學技術——中醫藥現代化，2004，6（6）：44.

[17] 賀福元，鄧凱文，劉文龍，等. 中藥與人體的作用本質“網通虹勢”的多重遺傳譜效動力學[J].中國實驗方劑學雜志，2011，17（2）：240.

Experimental studies on pharmacokinetics of three components in

Buyanghuanwu injection on base of total quantum statistical moment

HE Fuyuan1， 2， 3*， DENG Kaiwen4， LIU Wenlong1， SHI Jilian1，

WU Yong1， LIU Wei1， HE Qingping1， LI Bo1

（1. Department of Pharmaceutics， Hunan University of Traditional Chinese Medicine， Changsha 410208， China；

2. Key laboratory of Property and Pharmacodynamics， State Administration of Traditional Chinese Medicine，

Changsha 410208， China； 3. Pharmaceutical Preparation Technology and Evaluation Laboratory of Traditional Chinese Medicine，

Hunan University of Tradition Chinese Medicine， Changsha 410208， China；

4. The First Affinity Hospital， Hunan University of Tradition Chinese Medicine， Changsha 410007， China）

[Abstract] Objective： To verify established the total quantum statistic moments model with astragaloside IV， paeoniflorin， tetramethylpyrazine in Buyanghuanwu injection， in order to establish a pharmacokinetic experimental method with multicomponent traditional Chinese medicine（TCM） compound system. Method： The RPHPLC was adopted， with the chromatographic column of C18， 4.6 mm×250 mm， 5 μm. As for astragaloside IV， the ELSD detector was adopted with acetonitrilewater（35∶65） as the mobile phase at 1 mL·min-1； the pressure of column was （15.0±0.408） MPa， the column temperature was 30 ℃. Regarding paeoniflorin and tetramethylpyrazine， the detection of wavelengths was 254 nm， with acetonitrilewater （35∶65） as the mobile phase at 1 mL·min-1， the column pressure of （15.17±0.41） MPa. The pharmacokinetic parameters for single component were dealt with DAS and the total quantum statistical moment（TQSM） parameters were calculated using formulations. Result： All of the three components followed the two compartmental pharkacokinetic model （P<0.01） in rats. Compared with the superimposed total concentration， each single component showed difference in parameters up to 10 000 times at most， whereas the RSD of TQSM parameters was 3.510%. The TQSM pharmacokinetic parameters of the three components in Buyanghuanwu injection showed that AUCt， MRTt， VRTt， CLt，Vt were （119.8±27.20） g·min·L^-1， （210.0±54.49） min， （5.608±2.723）×104 min2， （0.319 6±0.068 8） mL·min-1·kg-1 and （64.12±8.243） mL·kg-1， respectively， suggesting that the halflife time for the three components were （145.5±37.76） min and 95% of them were metabolized within 0674.2 min. Conclusion： The TQSM can be used to study pharmacokinetic parameters of multicomponent TCM compound， because the method can characterize the pharmacokinetic regularity of quantumtime change in a multicomponent system.

[Key words] total quantum statistical moment； pharmacokinetics； pharmacokinetics with chromatographic fingerprint； Buyanghuanwu injection； astragaloside IV； paeoniflorin； tetramethylpyrazine

doi：10.4268/cjcmm20130222

[責任編輯 陳玲]