《生物药剂学与药物动力学》课程教学资源（文献资料）中药复方药物动力学总量统计矩法的实验验证研究_贺福元

第38卷第2期Vol 38,Issue 2中国中药杂志2013年1月January, 2013·药代动力学·中药复方药物动力学总量统计矩法的实验验证研究贺福元12,3*，邓凯文*，刘文龙"，石继连"，伍勇，刘伟"，贺庆平"，李博（1.湖南中医药大学药学院，湖南长沙410208；2.中药药性与药效国家中医药管理局重点实验室，湖南长沙4102083.湖南中医药大学现代中药制技术与评价实验室，湖南长沙4102084.湖南中医药大学第一附属医院，湖南长沙410007）[摘要】目的：对已建立的药物动力学总量统计矩法用补阳还五汤中黄芪甲苷、芍药苷、川芎嗪3个成分进行验证，建立中药复方多成分体系的药物动力学实验方法。方法：采用反相高效液相色谱法测定，其条件为Cls色谱柱（4.6mm×250mm，5μm）：黄芪甲苷ELSD测定：流动相乙→1%乙酸溶液（35:65），流速1mL：min*：柱压（15.0±0.408）MPa：温度30℃：荐药苷、川芎嗪的检测波长254nm，流动相甲醇-水（33:67），流速1mL：min：柱压（20.71±0.42）MPa。单个成分房室模型药物动力学参数采用DAS软件处理，总量统计矩按公式计算。结果：3个单成分在大鼠体内药物动力学都遵循二室模型（P<0.01），与叠加总浓度比，单个成分各参数相差最大上万倍，而总量统计矩参数的RSD为3.5%，3个成分药物动力学的总量统计矩参数AUC（119.8±27.20）gmin·L-;MRT,（210.0±54.49）min；VRT,（5.608±2.723）×10*min2：CL为(0.3196±0.0688）mLmin-·kg：V,为（64.12±8.243）mL·kg：to.95为（588.9±149.4）min。这就说明3个成分总量的半衰期为（145.5±37.76）min，3药在0~674.2min内代谢95%的浓度。结论：中药复方多成分药物动力学可采用总量统计矩法进行研究，该法能表征多成分体系总量的量时变化的动力学规律。[关键词】总量统计矩：药物动力学：谱动学：补阳还五汤：黄芪甲苷：芍药苷：川芎嗪中药及复方药物动力学研究是阐明中药复方产采用难度，为此，本论文先从第一个问题着手，详细生作用物质基础的关键问题，是阐明机体对中药及阐明建模的思路及参数的意义及测定方法，再逐步报道其相关问题。复方量时效关系的基础，也是中药复方制成控缓释制剂的理论基础。由于中药复方为多成分体补阳还五汤系清代王清任《医林改错》首载，该方由黄芪、当归、芍药、川芎、桃仁、红花、地龙组成，系，如按单成分药物动力学研究将会存在以下问题2-]：①怎样建立适宜多成分药物动力学研究的数具有补气、通络的功能。在抗血栓、抗脑栓血、抗衰老、抗血脂及免疫功能方面都有广泛的作用。本学模型，并对参数进行测定：②在对照品缺少的情况下，中药复方多成分怎样进行药物动力学研究：③怎室对该方进行了大量的抗脑血栓作用机制研究，发样建立宏观综合效应与总成分、单成分浓度对应关现该方具有重要的抗血栓活性作用。其主要有效成系。针对这些问题，作者前期创立了总量统计矩法，分部位为总、总苷元、总生物碱和总多糖，其代表性并有论文发表和成果获奖，但由于公式推导较复成分为黄芪甲苷、芍药苷、川芎嗪等等，本文选这3个杂，参数的理论与实验测定问题没有深入阐明清楚，成分作为代表性成分来阐明中药复方多成分如何进行可能会影响同行在进行中药复方药物动力学研究时以总量统计矩数学模型为基础的药物动力学研究。1多成分药物动力学的建模数学思路[稿件编号］20120524002多成分的药物动力学研究是回答各成分浓度随[基金项目】国家自然科学基金项目（81073142，30901971）：湖南时间量变规律，总体上有2种思路，一种是根据药物省自然科学基金重点项目（11JJ2055）：湖南省教育厅“十二五”药学浓度具有叠加性，可将各单成分浓度直接相加构成重点学科项目[通信作者]贺福元，教投，博士生导师，主要从事中药药理学、谱总量药物浓度，然后再按单成分药物动力学方法进动学、谱效学、谱效动力学、谱量学、中药药剂学及中医药信息数学研行参数求算，这种方法能获得总量药物动力学参数，究工作，Tel:(0731)85381372，13787213681，E-mail:pharmsharking@但会丢失各成分的药物动力学参数，没有建立起总tom. com·253

·药代动力学· 中药复方药物动力学总量统计矩法的实验验证研究 贺福元1，2，3* ，邓凯文4 ，刘文龙1 ，石继连1 ，伍勇1 ，刘伟1 ，贺庆平1 ，李博1 ( 1. 湖南中医药大学 药学院，湖南 长沙 410208; 2. 中药药性与药效国家中医药管理局重点实验室，湖南 长沙 410208; 3. 湖南中医药大学 现代中药制技术与评价实验室，湖南 长沙 410208; 4. 湖南中医药大学 第一附属医院，湖南 长沙 410007) ［摘要］ 目的: 对已建立的药物动力学总量统计矩法用补阳还五汤中黄芪甲苷、芍药苷、川芎嗪 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% 的浓度。结论: 中药复方多成分药物动力学可采用总量 统计矩法进行研究，该法能表征多成分体系总量的量-时变化的动力学规律。 ［关键词］ 总量统计矩; 药物动力学; 谱动学; 补阳还五汤; 黄芪甲苷; 芍药苷; 川芎嗪 ［稿件编号］ 20120524002 ［基金项目］ 国家自然科学基金项目( 81073142，30901971) ; 湖南 省自然科学基金重点项目( 11JJ2055) ; 湖南省教育厅“十二五”药学 重点学科项目 ［通信作者］ * 贺福元，教授，博士生导师，主要从事中药药理学、谱 动学、谱效学、谱效动力学、谱量学、中药药剂学及中医药信息数学研 究工作，Tel: ( 0731) 85381372，13787213681，E-mail: pharmsharking@ tom. com 中药及复方药物动力学研究是阐明中药复方产 生作用物质基础的关键问题，是阐明机体对中药及 复方量-时-效关系的基础，也是中药复方制成控缓 释制剂的理论基础［1］。由于中药复方为多成分体 系，如按单成分药物动力学研究将会存在以下问 题［2-5］: ①怎样建立适宜多成分药物动力学研究的数 学模型，并对参数进行测定; ②在对照品缺少的情况 下，中药复方多成分怎样进行药物动力学研究; ③怎 样建立宏观综合效应与总成分、单成分浓度对应关 系。针对这些问题，作者前期创立了总量统计矩法， 并有论文发表和成果获奖［6］，但由于公式推导较复 杂，参数的理论与实验测定问题没有深入阐明清楚， 可能会影响同行在进行中药复方药物动力学研究时 采用难度，为此，本论文先从第一个问题着手，详细 阐明建模的思路及参数的意义及测定方法，再逐步 报道其相关问题［7］。 补阳还五汤系清代王清任《医林改错》首载，该 方由黄芪、当归、芍药、川芎、桃仁、红花、地龙组成， 具有补气、通络的功能。在抗血栓、抗脑栓血、抗衰 老、抗血脂及免疫功能方面都有广泛的作用［8］。本 室对该方进行了大量的抗脑血栓作用机制研究，发 现该方具有重要的抗血栓活性作用。其主要有效成 分部位为总苷、总苷元、总生物碱和总多糖，其代表性 成分为黄芪甲苷、芍药苷、川芎嗪等等，本文选这 3 个 成分作为代表性成分来阐明中药复方多成分如何进行 以总量统计矩数学模型为基础的药物动力学研究。 1 多成分药物动力学的建模数学思路 多成分的药物动力学研究是回答各成分浓度随 时间量变规律，总体上有 2 种思路，一种是根据药物 浓度具有叠加性，可将各单成分浓度直接相加构成 总量药物浓度，然后再按单成分药物动力学方法进 行参数求算，这种方法能获得总量药物动力学参数， 但会丢失各成分的药物动力学参数，没有建立起总 ·253· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue2中国中药杂志2013年1月January,2013量药物动力学参数与单个成分药物动力学参数的关.系：第二方法是进行函数曲线转换后再叠加处理，先-MRT,(3)2.将单成分药物动力学曲线转变成一种公用函数曲线，再进行曲线叠加，最常用的方法就是泰勒由上诸式可知：设M表示r室模型第i项e的（Tatlor）转换叠加：傅利叶（Fourier）转换叠加等指数项前面的系数，α，表示第i项e的幂指系数。但对转换后的叠加函数不能直接给有价的药物动力一、二、三室成分的各总量统计矩参数意义及计算参学参数，因此较宜采用统计矩转换叠加构成总量统见文献[6]。非线性模型的统计矩参数：按统计矩计矩参数，与单成分统计矩参数一样，能表征多成分原理可得到非线性单成分零、一、二阶统计矩参数通宏观的药物动力学行为。由于单成药物动力学的统式，为式(5)，(6)，(7)。计矩模型及参数已被熟悉，其总量统计矩参数也易AUC.t+(5)为大家所接受，为此可以先将单成分的药物动力学V+2V.统计矩化，亦正态分布曲线化，然后再进行叠加构成2cg+9kgco+12kz(6)MRT, =总量统计矩参数，亦总量正态分布曲线化，由此可得V_(12k+6c)总量零、一、二阶矩参数。因此总量统计矩参数既能2cg +85kgco +20k.cg +180k.co +144kmVRT, =3(7)36V(2k+c)2反映总量药物浓度的动力学行为，又能与单个成分的统计矩参数关联，可用于房室、非房室及不规则成各总量统计矩参数意义及计算参见文献6]。分的药物动力学研究，是整合多成分药物动力学较不规则模型的统计矩参数：按统计矩参数定义，宜的数学处理方法。用梯形法或抛物线法计算零阶矩（AUC），按文献的1.1单成分的统计矩将单成分的药物动力学曲先计算AUMC（时间与曲线乘积下的面积），再计算MRT，再计算AUVC（时间平方与曲线乘积下的面线（参数）转变成统计矩曲线（参数）。其思路见图积），再计算VRT。1。单成分药时曲线的统计矩化过程实际上就是将1.2多成分体系的药动学总量统计矩参数的设定药物动力学曲线转变成正态分布曲线，其中重要参数有零阶矩（AUC）、一阶矩（MRT）和二阶矩及与单个成分指标的关系根据统计矩原理进行加(VRT)。合处理，见图2，图中分别为2个成分的统计矩曲线，经统计矩原理加合计算后可得图中b的总量统计矩曲线，由此可得总量统计矩参数：总量零阶矩AUC，、总量一阶矩MRT、总量二阶矩VRT。假设中药复方为一个由n个单体成分群构成体系，其中m个成分服从线性药物动力学。第i个成分服从ri室模型。S.个成分服从非线性药物动力学。S,个为不a药物动力学曲线：b.统计矩曲线。图1药物动力学曲线（参数）转变成统计矩曲线（参数）规则模型，该中药复方成分总量的零阶矩、一阶矩、Fig.1Statistical moment cuvers (parameters)translated from二阶矩及由此衍生的一系列动力学参数概念如下。pharmacokinetic cuvers (parameters)线性房室模型的统计矩参数：按统计矩原理可MRT,VRT得到单成分零、一、二阶统计矩参数通式，为式（1），(2),(3)。tM.2(1)AUC,=a整合前2药的总量统计矩正态曲线b.整合后2药的总量统计矩正态曲线。(2)MRT, =图2统计矩的加合性Fig.2Adding property for statistical moment·254·

量药物动力学参数与单个成分药物动力学参数的关 系; 第二方法是进行函数曲线转换后再叠加处理，先 将单成分药物动力学曲线转变成一种公用函数曲 线，再进行曲线叠加，最常用的方法就是 泰 勒 ( Tatlor) 转换叠加; 傅利叶( Fourier) 转换叠加等［9］， 但对转换后的叠加函数不能直接给有价的药物动力 学参数，因此较宜采用统计矩转换叠加构成总量统 计矩参数，与单成分统计矩参数一样，能表征多成分 宏观的药物动力学行为。由于单成药物动力学的统 计矩模型及参数已被熟悉，其总量统计矩参数也易 为大家所接受，为此可以先将单成分的药物动力学 统计矩化，亦正态分布曲线化，然后再进行叠加构成 总量统计矩参数，亦总量正态分布曲线化，由此可得 总量零、一、二阶矩参数。因此总量统计矩参数既能 反映总量药物浓度的动力学行为，又能与单个成分 的统计矩参数关联，可用于房室、非房室及不规则成 分的药物动力学研究，是整合多成分药物动力学较 宜的数学处理方法［6］。 1. 1 单成分的统计矩 将单成分的药物动力学曲 线( 参数) 转变成统计矩曲线( 参数) 。其思路见图 1。单成分药时曲线的统计矩化过程实际上就是将 药物动力学曲线转变成正态分布曲线，其中重要参 数有 零 阶 矩 ( AUC ) 、一 阶 矩 ( MRT ) 和 二 阶 矩 ( VRT) 。 a. 药物动力学曲线; b. 统计矩曲线。 图 1 药物动力学曲线( 参数) 转变成统计矩曲线( 参数) Fig. 1 Statistical moment cuvers( parameters) translated from pharmacokinetic cuvers( parameters) 线性房室模型的统计矩参数: 按统计矩原理可 得到单成分零、一、二阶统计矩参数通式，为式( 1) ， ( 2) ，( 3) 。 AUCri = ∑ ri i = 1 Mi αi ( 1) MRTri = ∑ ri i = 1 Mi α2 i ∑ ri i = 1 Mi αi ( 2) MRTri = 2 ∑ ri i = 1 Mi α3 i ∑ ri i = 1 Mi αi － MRT2 ri ( 3) 由上诸式可知: 设 Mi表示 ri室模型第 i 项 e 的 指数项前面的系数，αi 表示第 i 项 e 的幂指系数。 一、二、三室成分的各总量统计矩参数意义及计算参 见文献［6］。非线性模型的统计矩参数: 按统计矩 原理可得到非线性单成分零、一、二阶统计矩参数通 式，为式( 5) ，( 6) ，( 7) 。 AUCn = km c0 Vm + c 2 0 2Vm ( 5) MRTn = 2c 2 0 + 9km c0 + 12k 2 m Vm ( 12km + 6c0 ) ( 6) VRTn = 2c 4 0 + 85k 2 m c 2 0 + 20km c 3 0 + 180k 3 m c0 + 144k 4 m 36V2 m ( 2km + c0 ) 2 ( 7) 各总量统计矩参数意义及计算参见文献［6］。 不规则模型的统计矩参数: 按统计矩参数定义， 用梯形法或抛物线法计算零阶矩( AUC) ，按文献的 先计算 AUMC( 时间与曲线乘积下的面积) ，再计算 MRT，再计算 AUVC( 时间平方与曲线乘积下的面 积) ，再计算 VRT ［2］。 1. 2 多成分体系的药动学总量统计矩参数的设定 及与单个成分指标的关系 根据统计矩原理进行加 合处理，见图 2，图中 a 分别为 2 个成分的统计矩曲 线，经统计矩原理加合计算后可得图中 b 的总量统 计矩曲线，由此可得总量统计矩参数: 总量零阶矩 AUCt、总量一阶矩 MRTt、总量二阶矩 VRTt。假设中 药复方为一个由 n 个单体成分群构成体系，其中 m 个成分服从线性药物动力学。第 i 个成分服从 ri 室 模型。S1个成分服从非线性药物动力学。S2个为不 规则模型，该中药复方成分总量的零阶矩、一阶矩、 二阶矩及由此衍生的一系列动力学参数概念如下。 a. 整合前 2 药的总量统计矩正态曲线; b. 整合后 2 药的总量统计矩 正态曲线。 图 2 统计矩的加合性 Fig. 2 Adding property for statistical moment ·254· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue2中国中药杂志January,20132013年1月(13)t12. = 0.693MRT..1.3中药复方的总量统计矩参数包括总量统计零、总体表观消除平衡常数倒数为MRT倒数，为、二阶矩以及由此衍生的其他总量统计矩参数。式(14)。总量零阶矩（AUC）：根据统计矩的加合性可计算中药复方多成分的总量统计矩，为下式（8）。10.693(14)k.=MRT..11/2.91M+2 2km.fo.i + cal +AUC, =AUC, =MA总体表观清除率为式（15）。台2Vm.lX.oCL, = -(15)T(8)Car.gdtAIC台J式中m表示线性成分的个数，r为第i成分代总体表观分布容积为式（16）。谢的室数，s,为非线性成分的个数，s,为产生多峰不X...)×AUMC规则成分个数。(16)V, = CI, × MRT, =(SAUC)总量一阶矩（MRT）：即为总量平均驻留时间，亦为各成分的总体均值，可定义如下式（9）。95%总量与累积95%总量浓度时间区间：由于≥MRT, ×AUC,中药复方为多成分体系，用零阶矩只表征其总成分MRT,=(9)的生物利用度，用一阶矩只表征诸成分平均驻留时BAUC,间中心，这还不能完整表征多成分的药物动力学行总量二阶矩（VRT，）：亦可采用相同方法处理，为，因为平均驻留时间中心相同的一组成分其药物的可得式（10）。种类与离散程度可以不同，也就是95%总量与累积5[(MRT + VRT,)× AUC]95%时间区间不同，因此将一、二阶矩结合可表征95%台- MRT?(10)VRT,=总量浓度时间区间，亦为式（17)。2.AuCc.95%总量浓度时间区间=[MRT,-1.96VRT,，MRT,+由上诸式可知，一个复方的总体药物动力学趋1.96VRT,J(17)势由单个成分的动力学参数决定。其总量零阶矩等累积总量95%浓度时间区间=［-，a](18)于各单成分的零价矩之和：总量一阶矩等于各单个其中to.gs按累积正态分布函数求出。成分一阶矩对零阶矩的算术平均值：总量二阶矩等(r-MRT,)20.95195%=2VRT,dt(19)2元VRT/于各单个成分的一阶矩的平方与二阶矩之和对零阶矩的算术平均值，再与总量一阶矩平方之差。这些这样，式（8）～~（19）构成了中药复方多成分药性质是由总量统计矩的加合性决定的，亦有式物动力学总量统计矩数学模型及参数体系6.10]。(11),(12)。2材料与方法MRT AUG MRT, 12.1仪器与试药Agilent-200液相色谱仪：All-Z AUC, + MRT,eAUC-o+tech3300ELSD检测与UV测定。瑞士Precisa公司(11)的20A型电子天平：德国Heraeus公司的高速冷冻AUC, +VRT, + MRT ZAUC, = VRT + MRT离心机及SZ-1型快速混匀器。黄芪甲苷（批号台0110781-200512）、芍药苷（批号110736-200629）、川AUCVRT + MRT(12)营嗪（批号110817-200305）（购于国家药品生物制业品检定所），甲醇、乙睛为色谱纯，蒸馏水为二次重式中MRT，为1~i成分的总量一阶矩，VRT，为对应的总量二阶矩为：MRT，为i~n成分的总量蒸馏水，其他试剂均为分析纯。其复方药材来源为一阶矩，VRT，为对应的总量二阶矩。由总量零、一、黄芪为豆科植物蒙古黄芪Astragalusmembranaceus二阶矩可衍生出其它总量药物动力学参数。Bge.var.monghokicus（Bgo）Hsiao干燥根：当归为总量表观半衰期、清除率和表观分布容积：根据伞形科植物当归Angelicasinensis（Oliv）Diels的干统计矩与半衰期有关系，中药复方多成分的总量表燥根：川芎为伞形科植物川芎Ligusticumchuanxiong观半衰期为式（13）。Hort的干燥根：赤为为毛科植物芍药Paeonia.255

1. 3 中药复方的总量统计矩参数 包括总量统计零、 一、二阶矩以及由此衍生的其他总量统计矩参数。 总量零阶矩( AUCt ) : 根据统计矩的加合性可计 算中药复方多成分的总量统计矩，为下式( 8) 。 AUCt = ∑ n i = 1 AUCi = ∑ m i = 1 ∑ ri j = 1 Mi，j αi，j + ∑ s1 l = 1 2km，lc0，l + c 2 0，l 2Vm，l + ∑ s2 q = 1 ∫ ∞ 0 cnr，qdt ( 8) 式中 m 表示线性成分的个数，ri为第 i 成分代 谢的室数，si为非线性成分的个数，si为产生多峰不 规则成分个数。 总量一阶矩( MRTt ) : 即为总量平均驻留时间， 亦为各成分的总体均值，可定义如下式( 9) 。 MRTt = ∑ n i = 1 MRTi × AUCi ∑ n i = 1 AUCi ( 9) 总量二阶矩( VRTt ) : 亦可采用相同方法处理， 可得式( 10) 。 VRTt = ∑ n i = 1 ［( MRT2 i + VRTi ) × AUCi ］ ∑ n i = 1 AUCi － MRT2 t ( 10) 由上诸式可知，一个复方的总体药物动力学趋 势由单个成分的动力学参数决定。其总量零阶矩等 于各单成分的零价矩之和; 总量一阶矩等于各单个 成分一阶矩对零阶矩的算术平均值; 总量二阶矩等 于各单个成分的一阶矩的平方与二阶矩之和对零阶 矩的算术平均值，再与总量一阶矩平方之差。这些 性质是由总量统计矩的加合性决定的，亦 有 式 ( 11) ，( 12) 。 MRTt∑ n i = 1 AUCi = MRTp1∑ j l = 1 AUCl + MRTp2∑ n s = j+1 AUCs ( 11) VRTt + MRT2 t∑ n i = 1 AUCi = VRTp1 + MRT2 p1 ∑ j l = 1 AUCl + VRTp2 + MRT2 p2∑ n s = j+1 AUCs ( 12) 式中 MRTp1 为 1 ～ i 成分的总量一阶矩，VRTp1 为对应的总量二阶矩为; MRTp2为 i ～ n 成分的总量 一阶矩，VRTp2为对应的总量二阶矩。由总量零、一、 二阶矩可衍生出其它总量药物动力学参数。 总量表观半衰期、清除率和表观分布容积: 根据 统计矩与半衰期有关系，中药复方多成分的总量表 观半衰期为式( 13) 。 t1 /2，t = 0. 693MRTt，iv ( 13) 总体表观消除平衡常数倒数为 MRTt，iv倒数，为 式( 14) 。 kt，iv = 1 MRTt，iv = 0. 693 t1 /2，t ( 14) 总体表观清除率为式( 15) 。 Clt = Xo，t( i) ∑ n i = 1 AUCt ( 15) 总体表观分布容积为式( 16) 。 Vt = Clt × MRTt = Xo，t，( v) × AUMCt ( ∑ n i = 1 AUCt ) 2 ( 16) 95% 总量与累积 95% 总量浓度时间区间: 由于 中药复方为多成分体系，用零阶矩只表征其总成分 的生物利用度，用一阶矩只表征诸成分平均驻留时 间中心，这还不能完整表征多成分的药物动力学行 为，因为平均驻留时间中心相同的一组成分其药物的 种类与离散程度可以不同，也就是 95% 总量与累积 95%时间区间不同，因此将一、二阶矩结合可表征 95% 总量浓度时间区间，亦为式( 17) 。 95% 总量浓度时间区间 =［MRTt － 1. 96VRTt 1 2 ，MRTt + 1. 96VRTt 1 2 ］ ( 17) 累积总量 95% 浓度时间区间 =［－ ∞ ，a］ ( 18) 其中 t0. 95按累积正态分布函数求出。 95% = ∫ t0. 95 －∞ 1 槡2πVRT1 /2 t e － ( t－MRTt ) 2 2VRTt dt ( 19) 这样，式( 8) ～ ( 19) 构成了中药复方多成分药 物动力学总量统计矩数学模型及参数体系［6，10］。 2 材料与方法 2. 1 仪器与试药 Agilent-1200 液相色谱仪，All￾tech 3300 ELSD 检测与 UV 测定。瑞士 Precisa 公司 的 20 A 型电子天平; 德国 Heraeus 公司的高速冷冻 离心机 及 SZ-1 型 快 速 混 匀 器。黄 芪 甲 苷 ( 批 号 110781-200512) 、芍药苷( 批号 110736-200629) 、川 芎嗪( 批号 110817-200305) ( 购于国家药品生物制 品检定所) ，甲醇、乙腈为色谱纯，蒸馏水为二次重 蒸馏水，其他试剂均为分析纯。其复方药材来源为: 黄芪为豆科植物蒙古黄芪 Astragalus membranaceus Bge. var. monghokicus ( Bgo) Hsiao 干燥根; 当归为 伞形科植物当归 Angelica sinensis ( Oliv) Diels 的干 燥根; 川芎为伞形科植物川芎 Ligusticum chuanxiong Hort 的干 燥 根; 赤芍为毛莨科植物芍药 Paeonia ·255· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue 2中国中药杂志2013年1月January, 2013lactloraPall的干燥根茎；桃仁为蔷薇科植物物桃统计矩参数。Prunuspersica（Linn）Batsch的干燥成熟种子：红花3结果为菊科植物红花CarthamustinctoriusL的干燥花；3.1保留时间黄芪甲苷：血样经预处理后，在黄地龙为环节动物门钜蚓科动物参环毛蚓Pheretima芪甲苷色谱条件下进样测定，得血浆色谱图，见图aspergillum（Perrier）的干燥主体。处方黄芪60g，3。大鼠血浆中黄芪甲苷峰在保留时间为15min左当归9g，川芎6g赤芍9g，桃仁9g，红花9g，地龙右可完全分离。药着写川芎嗪：在符药耆、川写嗪9g。Wistar大白鼠，雌雄各半，体重220~300g，由对照品液的色谱条件下测得为药苷、川芎嗪的血浆湖南中医药大学实验动物中心提供，许可证号色谱图。大鼠血浆补阳还五汤中药苷峰在保留时SCXK（湘）2001-2004。间为11.5min左右，川芎嗪峰在保留时间为21.5min左右可完全分离。2.2色谱条件色谱柱UltimateXB-Cis柱（4.6mm×250mm，5μm），黄芪甲苷采用Agi-3.2标准曲线的制备及线性关系的考察分别取试管6只，加入802.1mg·L-黄芪甲苷甲醇溶液lent1200-Alltech3300ELSD测定，流动相乙睛-%对照品液10，25，50，100，200，400μL，水浴蒸干，乙酸溶液（35:65），流速1mLmin-，柱压（15.0±0.408）MPa；温度30℃。芍药苷与川芎嗪采用Agi-再加入空白血浆400μL，经样品预处理后，残渣用50μL甲醇溶解，使血浆中黄芪甲苷的质量浓度分别为lentUV测定，由于两者同时以指纹图谱的形式测定，可在原测定波长的基础上调整为254nm12-3，20.05~802.1mg·L-；另取试管6只，分别加入320.4mg·L为药苷甲醇溶液对照品液1,15,30,60,75，100流动相甲醇-水（33：67），流速1mL·min-，柱压μL，水浴蒸干，再加空白血浆400μL，经样品预处理后，(20.71±0.42）MPa；温度35℃。残渣用50μL甲醇溶解，使血浆中芍药苷的质量浓度2.3给药与取血将补阳还五汤水提醇沉制得的分别为0.801~80.10mg·L-；又取试管6只，分别加复方注射液（分别含黄芪甲苷2.160g·L-，苟药苷入5.140mgL-川芎嗪甲醇溶液对照品液2，5，10，2.578g·L-、川芎嗪25.12mg.L-）。实验大鼠1520，40μL，水浴蒸干，另加入空白血浆400μL，经样尾静脉给药2mL，分别于给药后0.25，0.5，1.25，1.75,2.5,4.5,5.5，8，17.5h，眼眶取血浆400μL，品预处理后，残渣用50证甲醇溶解，使血浆中川弯嗪的质量浓度分别为25.70~514.0μg·L-。在前述色并在1h内进行样品预处理。谱条件下进行测定，以黄芪甲者浓度为纵坐标，黄芪甲2.4样品预处理将补阳还五汤待测血浆400μL苷峰面积为横坐标进行线性回归处理，得回归方程为置手尖底塑料刻度离心管中，加入三氯醋酸液1C=0.18324-0.01,结果表明黄芪甲苷在20~800mgmL，涡旋混合5min，5000r·min-离心10min，取·L-呈良好线性（r=0.9997），黄芪甲苷的最低检测三氯醋酸层，常温下氮气流吹干，残渣用甲醇50μL浓度为0.6mg·L-（S/N>3）：以芍药苷面积为纵坐溶解，取20μL进样。标，药苷浓度为横坐标进行线性回归处理，得回归方2.5空白血浆处理于实验前，取大鼠血样按样品程为A=178.49612C-3.361结果表明荷药昔在处理项下同法制得空白血浆样。0.80~80mg·L-呈良好线性（r=0.9999），芍药苷的2.6对照品液制备精密称取黄芪甲苷对照品最低检测浓度为0.12mg·L-（S/N>3）：以川芎嗪面8.0mg，用甲醇定容至10mL，制得每1mL含黄芪积为纵坐标，川弯嗪浓度为横坐标进行线性回归处理，甲苷0.8mg的对照品液：精密称取芍药苷对照品得回归方程为A=64.70791C+86.32575，结果表明川3.2mg，用甲醇定容至10mL，制得每1mL含芍药芎嗪在25.70~514.0ugL-呈良好线性（r=苷0.32mg的对照品液：精密称取川芎嗪对照品0.9997），川嗪的最低检测质量浓度为8μg·L-15.0mg，用甲醇定容至100mL，再稀释10倍制得每(S/N>3)。1mL含川芎嗪0.005mg的对照品液。3.3精密度和回收率精密度按文献分低、中、高2.7数据处理与分析结果用DAS2.0软处理，分3个浓度进行日内和日间精密度、加样回收率试验别得单个成分、3成分叠加及平均的经典药物动力研究。学参数，按自建的统计矩法分析得药物动力学总量·256·

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( 湘) 2001-2004。 2. 2 色 谱 条 件 色 谱 柱 Ultimate XB-C18 柱 ( 4. 6 mm × 250 mm，5 μm) ，黄芪甲苷［11］ 采用 Agi￾lent 1200-Alltech 3300 ELSD 测定，流动相乙腈-1% 乙酸溶液( 35∶ 65) ，流速 1 mL·min － 1 ，柱压( 15. 0 ± 0. 408) MPa; 温度 30 ℃。芍药苷与川芎嗪采用 Agi￾lent UV 测定，由于两者同时以指纹图谱的形式测 定，可在原测定波长的基础上调整为 254 nm［12-13］， 流动相甲醇-水( 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 个浓度进行日内和日间精密度、加样回收率试验 研究。 ·256· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue2中国中药杂志2013年1月January,20130.004.919.814.7414.864.939.860.004.959.9175101251517.52022.5 25275305751012315175202252527530 25 57.5 10 123 15 17.5: 20 225 25 275 3025575101251517.520225252753a黄芪甲苷色谱图（1.黄芪甲苷对照品：2.血浆样品：3.含药血浆样品）：h.芍药苷与川芎嗪色谱图（1.药药苷对照品：2.川嗪对照品：3.含药血浆样品：4.空白血浆样品）。图3补阳还五汤血浆样品中黄芪甲苷、药苷与川芎嗪HPLC图Fig 3HPLC Chromatography of astragaloside IV, paeonflorin and tetramethylpyrazine in Buyang Huanwu injection精密度试验：取大鼠空白血浆400μL，加入黄个浓度的样品，经预处理后分别测定日内精密度和芪甲苷、芍药苷、川芎嗪对照品液，配制低、中、高3日间精密度，结果见表1。表1血浆中黄芪甲苷、药苷、川嗪的精密度测定（n=3)Table1Accuracy of determination of astragaloside V,paeonflorinin, tetramethylpyrazine in plasma (n =3)日内精密度日间精密度加入质量浓度成分/mg -L-t检测量/mg·L-1RSD/%检测量/mg·L-!RSD/%黄芪甲苷40. 1138.1340. 612. 12. 22.72.7401.1397.3403.61.9795.71.6802.1793.72.21.5苟药苔1.6021.5621.6522.940.0540.3839.202.62. 4 80.1078.4079. 392. 90. 025 72. 40. 036 72. 8川芎嗪0. 024 22.50.25700.24870.24552. 32.62.90.514 00.507 80.514 3回收率试验：取空白血浆，加入一定量的黄芪甲实验前夜禁食12h，尾静脉给补阳还五汤注射液2苷、芍药苷、川嗪对照品液，分别配制高、中、低3mL，定时取血与处理，可得药-时曲线，见图4。结果个不同浓度的样品，经预处理后测定，按回归方程计用DAS2.0软处理，分别得单个成分、3个成分叠加算测得量，得到的平均方法回收率列为表2。由表及平均的药物动力学参数，见表3，按自建的统计矩1~2可知，本法能用于3个成分的药物动力学研究法分析得药物动力学总量统计矩参数，见表4。中的浓度测定。3.5补阳还五汤中黄芪甲苷、苟药苷、川芎嗪药物3.4药物动力学测定取大鼠100只，随机分为常规药物动力学及总量统计矩参数计算先用10组，每组10只，每一组为一实验点，一组为空白。DAS软件进行曲线拟合，得表3中常见的药物动力.257

a. 黄芪甲苷色谱图( 1. 黄芪甲苷对照品; 2. 血浆样品; 3. 含药血浆样品) ; b. 芍药苷与川芎嗪色谱图( 1. 芍药苷对照品; 2. 川芎嗪对照品; 3. 含药血浆样品; 4. 空白血浆样品) 。 图 3 补阳还五汤血浆样品中黄芪甲苷、芍药苷与川芎嗪 HPLC 图 Fig. 3 HPLC Chromatography of astragalosideⅣ，paeonflorin and tetramethylpyrazine in Buyang Huanwu injection 精密度试验: 取大鼠空白血浆 400 μL，加入黄 芪甲苷、芍药苷、川芎嗪对照品液，配制低、中、高 3 个浓度的样品，经预处理后分别测定日内精密度和 日间精密度，结果见表 1。 表 1 血浆中黄芪甲苷、芍药苷、川芎嗪的精密度测定( n = 3) Table 1 Accuracy of determination of astragaloside Ⅳ，paeonflorinin，tetramethylpyrazine in plasma ( n = 3) 成分 加入质量浓度 /mg·L － 1 日内精密度 日间精密度 检测量/mg·L － 1 RSD/% 检测量/mg·L － 1 RSD/% 黄芪甲苷 40. 11 40. 61 2. 1 38. 13 2. 2 401. 1 397. 3 2. 7 403. 6 2. 7 802. 1 793. 7 1. 9 795. 7 1. 6 芍药苷 1. 602 1. 562 2. 2 1. 652 1. 5 40. 05 40. 38 2. 9 39. 20 2. 6 80. 10 78. 40 2. 4 79. 39 2. 9 川芎嗪 0. 025 7 0. 024 2 2. 4 0. 036 7 2. 8 0. 257 0 0. 248 7 2. 5 0. 245 5 2. 3 0. 514 0 0. 507 8 2. 6 0. 514 3 2. 9 回收率试验: 取空白血浆，加入一定量的黄芪甲 苷、芍药苷、川芎嗪对照品液，分别配制高、中、低 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 中常见的药物动力 ·257· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue2中国中药杂志2013年1月January,2013表2血浆中黄芪甲苷、药苷、川芎嗪的回收率（±s，n=3)学参数，复方中黄芪甲苷、药苷、川芎嗪3个成分Table2Recovery of astragaloside IV,paeonflorinin,tetram-在大鼠体内药物动力学都遵循二室模型（P<ethylpyrazine in plasma ( ± s,n =3)0.01），再依式（1）~式（4）的计算公式对黄芪甲苷、加入质量浓度检测量RSD平均回收率芍药苷、川芎嗪单个成分的统计矩参数。如没有进成分/mg · L-1I/mg · L-11%1%行曲线拟合，也可直接用梯形法或抛物线法按有关39.2597. 86黄芪甲蓉40. 112.1统计矩公式进行计算，然后运用式（5）～式（11）将401.1398. 199. 262. 33药的决量统计矩进行加合计算得总量药物动力学802.1795.999.221.4芍药苷1.6021. 56797.792.4统计矩参数，在计算过程一定要注意只有AUC（零40. 0541. 43103. 43.3阶矩）、AUMC（时间与曲线乘积下的面积）、AUVC80.1080.77100.81.7（时间平方与曲线乘积下的面积）能进行相加，而0. 025 798. 562.7川芎嗪0. 025 33MRT，VRT不能相加，只能按公式MRT=AUMC/0.25700.248296.584.70.514 00.503 898. 011.4AUC，VRT=AUVC/AUC－MRT计算，其他参数在70800:01700t600.66000.55004.1-800.4童i0.30200.2K1091100O5001 500100050020040060001.0001 50090080800800+5药行700+考药谷700一川喝子60600600—川嗪4-黄花甲#50黄装甲貨宣40043成分含加403减分加1003002001000001 5005001 0001500500200400,600-80010001200/mia.黄芪甲苷：h.药苷：c川嗪：d.补阳还五汤中黄芪甲苷、药药苷与川芎嗪3个成分叠加：e.黄芪甲苷、芍药苷与川萼嗪3个成分叠加总量（全图）黄芪甲苷、药药苷与川草嗪3个成分叠加总量（放大）。图4补阳还五汤中黄芪甲苷、药药苷、川草嗪及3个成分叠加的药时曲线Fig.4Pharmacokinetic profiles of the astragalosideIV,paeonflorinin and tetramethylpyrazine in the BYHWD injection and their con-centration overlapedMRT参数基础上进行计算。用上法算得补阳还五3.6补阳还五汤中黄芪甲苷、药苷、川嗪3个汤中3个成分的统计矩参数列表4，其中总量统计成分95%总量浓度与累积95%总量浓度代谢时间矩参数AUC为（119.8±27.20）g：min·L-1;MRT区间按式（14）计算得3药的95%总量浓度在[210.0-1.96(5.608×10*)1/2,210.0 +1.96(5.608x为（210.0±54.49）min；VRT为（5.608±2.723）×10*）1/2]min，亦在-254.2~674.2min时间内代谢10*min2;tun为（145.5±37.76）min;CL,为（0.319695%的浓度。按式（16）算得tog为588.9min，故累积±0.06882）mLmin-1:kg-";V,为（64.12±95%总量浓度代谢时间区间为［-0，588.9]min。8.243)mLkgl;to.95.为(588.9±149.4)min。·258·

表2 血浆中黄芪甲苷、芍药苷、川芎嗪的回收率( x珋± s，n = 3) Table 2 Recovery of astragaloside Ⅳ，paeonflorinin，tetram￾ethylpyrazine in plasma( x珋± s，n = 3) 成分 加入质量浓度 /mg·L － 1 检测量 /mg·L － 1 平均回收率 /% RSD /% 黄芪甲苷 40. 11 39. 25 97. 86 2. 1 401. 1 398. 1 99. 26 2. 3 802. 1 795. 9 99. 22 1. 4 芍药苷 1. 602 1. 567 97. 79 2. 4 40. 05 41. 43 103. 4 3. 3 80. 10 80. 77 100. 8 1. 7 川芎嗪 0. 025 7 0. 025 33 98. 56 2. 7 0. 257 0 0. 248 2 96. 58 4. 7 0. 514 0 0. 503 8 98. 01 1. 4 学参数，复方中黄芪甲苷、芍药苷、川芎嗪 3 个成分 在大鼠 体 内 药 物 动 力 学 都 遵 循 二 室 模 型 ( P ＜ 0. 01) ，再依式( 1) ～ 式( 4) 的计算公式对黄芪甲苷、 芍药苷、川芎嗪单个成分的统计矩参数。如没有进 行曲线拟合，也可直接用梯形法或抛物线法按有关 统计矩公式进行计算，然后运用式( 5) ～ 式( 11) 将 3 药的决量统计矩进行加合计算得总量药物动力学 统计矩参数，在计算过程一定要注意只有 AUC( 零 阶矩) 、AUMC( 时间与曲线乘积下的面积) 、AUVC ( 时间平方与曲线乘积下的面积) 能进行相加，而 MRT，VRT 不 能 相 加，只 能 按 公 式 MRT = AUMC / AUC，VRT = AUVC /AUC － MRT2 计算，其他参数在 a. 黄芪甲苷; b. 芍药苷; c. 川芎嗪; d. 补阳还五汤中黄芪甲苷、芍药苷与川芎嗪 3 个成分叠加; e. 黄芪甲苷、芍药苷与川芎嗪 3 个成分叠加总 量( 全图) ; f. 黄芪甲苷、芍药苷与川芎嗪 3 个成分叠加总量( 放大) 。 图 4 补阳还五汤中黄芪甲苷、芍药苷、川芎嗪及 3 个成分叠加的药时曲线 Fig. 4 Pharmacokinetic profiles of the astragalosideⅣ，paeonflorinin and tetramethylpyrazine in the BYHWD injection and their con￾centration overlaped 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。 ·258· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol. 38, Issue2中国中药杂志2013年1月January,2013chinaJoumalot ChicMetic补阳还五汤复方静脉给药大鼠血液中黄芪甲苷（2.160g·L-）芍药苷（2.578g·L-）、川芎嗪（25.12mg·L-）的药表3物动力学参数（±s,n=10）Pharmacokinetic parameters of astragalosideV (2.160 g·L-), paeonflorinin(2.578 g·L), tetramethylpyrazine(25.12Table3mg L-") in white rat after iw Buyanghuanwu injection ( ±s,n =10)A/g·L-fupa /min"t成分B/g L-1α/min-!β/min-黄芪甲苷366. 3 ±45. 38452.5 ±57. 330. 042 78 ±0. 008 3510. 004 509 ±0. 001 12416. 77 ±3. 243芍药苷1 309 ± 159. 323. 65 ±3.3850. 390 9 ±0. 081 520. 011 07 ±0. 001 7871. 837 ±0.349 8川芎嗪0. 352 1 ±0. 082 000. 218 0 ±0. 056 290. 027 43 ±0.005 1620. 004 387 ±0. 000 760 526. 11 ±5. 092平均558.4 ± 58.71158. 8 ± 19. 290. 153 7 ±0. 027 550. 006 654 ±0. 000 913 314.91 ±2. 166成分叠加拟合435. 8 ±61. 48468. 2 ± 39. 080. 047 83 ±0. 009 2460. 004 638 ±0. 001 12515. 05 ± 3. 223fu/2g / min~lKg /min-IKan /min -!Kro /min -Ve /mL. - kg-!P成分x0.01(9)=21.666，其RSD也值最大，为2883%，两者相差6万倍，黄芪甲苷、芍为117.9%，远大于一般药物动力学实验误差5%，·259

表 3 补阳还五汤复方静脉给药大鼠血液中黄芪甲苷( 2. 160 g·L － 1 ) 、芍药苷( 2. 578 g·L － 1 ) 、川芎嗪( 25. 12 mg·L － 1 ) 的药 物动力学参数( x珋± s，n = 10) Table 3 Pharmacokinetic parameters of astragalosideⅣ( 2. 160 g·L － 1 ) ，paeonflorinin( 2. 578 g·L － 1 ) ，tetramethylpyrazine( 25. 12 mg·L － 1 ) in white rat after iv Buyanghuanwu injection ( x珋± s，n = 10) 成分 A/g·L －1 B/g·L －1 α/min －1 β/min －1 t1/2α /min －1 黄芪甲苷 366. 3 ±45. 38 452. 5 ±57. 33 0. 042 78 ±0. 008 351 0. 004 509 ±0. 001 124 16. 77 ±3. 243 芍药苷 1 309 ±159. 3 23. 65 ±3. 385 0. 390 9 ±0. 081 52 0. 011 07 ±0. 001 787 1. 837 ±0. 349 8 川芎嗪 0. 352 1 ±0. 082 00 0. 218 0 ±0. 056 29 0. 027 43 ±0. 005 162 0. 004 387 ±0. 000 760 5 26. 11 ±5. 092 平均 558. 4 ±58. 71 158. 8 ±19. 29 0. 153 7 ±0. 027 55 0. 006 654 ±0. 000 913 3 14. 91 ±2. 166 成分叠加拟合 435. 8 ±61. 48 468. 2 ±39. 08 0. 047 83 ±0. 009 246 0. 004 638 ±0. 001 125 15. 05 ±3. 223 成分 t1/2β /min －1 K12 /min －1 K21 /min －1 K10 /min －1 VC /mL·kg －1 P 黄芪甲苷 162. 5 ±39. 85 0. 014 142 ±0. 003 646 0. 025 67 ±0. 005 254 0. 007 471 ±0. 001 633 5. 331 ±0. 578 7 ＜0. 01 芍药苷 63. 94 ±9. 265 0. 142 1 ±0. 038 88 0. 017 82 ±0. 003 097 0. 242 1 ±0. 043 76 3. 921 ±0. 470 8 ＜0. 01 川芎嗪 162. 3 ±28. 14 0. 009 549 ±0. 002 912 0. 013 17 ±0. 002 583 0. 009 095 ±0. 001 430 91. 68 ±18. 86 ＜0. 01 平均 129. 6 ±19. 64 0. 055 27 ±0. 012 68 0. 018 89 ±0. 001 932 0. 086 21 ±0. 0153 4 33. 64 ±6. 486 成分叠加拟合 157. 8 ±37. 60 0. 017 38 ±0. 004 962 0. 026 97 ±0. 004 945 0. 008 114 ±0. 001 546 10. 62 ±0. 927 9 ＜0. 01 表 4 补阳还五汤中单成分及 3 个成分总量统计矩参数( x珋± s，n = 10) Table 4 Total quantum statistical moment parameters of astragaloside Ⅳ，paeonflorinin，tetramethylpyrazine in Buyanghuanwu injec￾tion and their their concentration overlaped ( x珋± s，n = 10) 成分 总量统计矩参数 AUCt /g·min·L － 1 MRTt /min VRTt /min2 CLt /mL·min － 1· kg － 1 Vt /mL·kg － 1 t0. 95t /min 黄芪甲苷 114. 2 ± 26. 99 218. 3 ± 55. 93 5. 688 × 104 ± 2. 743 × 104 0. 152 7 ± 0. 034 25 31. 83 ± 4. 222 600. 1 ± 150. 9 芍药苷 5. 567 ± 0. 566 6 37. 40 ± 5. 909 5 299 ± 1 510 3. 594 ± 0. 344 6 134. 4 ± 25. 27 156. 0 ± 23. 46 川芎嗪 0. 063 08 ± 0. 012 47 193. 8 ± 35. 92 5. 154 × 104 ± 1. 909 × 104 3. 157 ± 0. 542 1 606. 5 ± 130. 4 561. 7 ± 102. 6 平均 39. 93 ± 9. 065 149. 8 ± 25. 91 3. 791 × 104 ± 1. 315 × 104 2. 301 ± 0. 221 7 257. 6 ± 39. 51 465. 9 ± 79. 98 总量统计矩计算 119. 8 ± 27. 20 210. 0 ± 54. 49 5. 608 × 104 ± 2. 723 × 104 0. 319 6 ± 0. 068 82 64. 12 ± 8. 243 588. 9 ± 149. 4 成分叠加拟合 115. 6 ± 27. 52 210. 7 ± 52. 43 5. 327 × 104 ± 2. 514 × 104 0. 332 7 ± 0. 075 18 66. 83 ± 6. 831 580. 7 ± 142. 2 3. 7 补阳还五汤中黄芪甲苷、芍药苷、川芎嗪的药 动学参数的比较 以 3 个成分叠加拟合而得的药物 动力学参数( P ＜ 0. 01) 为标准对照组，表 3 中其他 组与之的百分比值列于表 5。 表 5 补阳还五汤静脉给药大鼠血液中黄芪甲苷、芍药苷、川芎嗪药物动力学参数对比 Table 5 Comparison of pharmacokinetic parameters of astragaloside Ⅳ，paeonflorinin，tetramethylpyrazine in white rat after iv Buyan￾ghuanwu injection 成分 A B α β t1/2α t1/2 β K12 K21 K10 VC RSD/% 黄芪甲苷 84. 05 96. 65 89. 44 97. 22 111. 4 103. 0 81. 37 95. 18 92. 08 50. 20 18. 32 芍药苷 300. 4 5. 050 817. 3 238. 7 12. 21 40. 52 817. 6 66. 07 2 983 36. 92 172. 4 川芎嗪 0. 080 0 0. 050 0 57. 35 94. 59 173. 5 102. 9 54. 94 48. 83 112. 1 863. 3 169. 7 平均 128. 1 33. 92 321. 4 143. 5 99. 07 82. 13 318. 0 70. 04 1062 316. 8 117. 9 注: 成分叠加拟合均为 100. 0( 表 6 同) 。 表 3 诸成分参数比较可采用等级资料的 χ 2 检 验，其 χ 2 = 279 8χ 2 0. 01( 36) = 58. 619，说明诸成分参 数差异显者( P ＜ 0. 01) ，由表 5 可知，以川芎嗪的参 数 B 的比值最小，为 0. 050% ; 以芍药苷的 K10的比 值最大，为2 883% ，两者相差 6 万倍，黄芪甲苷、芍 药苷、川芎嗪组的 RSD 分别为 18. 32% ，172. 4% ， 169. 7% ，这说明中药复方中各成分的药物动力学参 数相差悬殊，彼此之间参数不能相互代替。平均组 与叠加组的 χ 2 = 165. 7χ 2 0. 01( 9) = 21. 666，其 RSD 也 为 117. 9% ，远大于一般药物动力学实验误差 5% ， ·259· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

Vol 38, Issue2第38卷第2期中国中药杂志January,20132013年1月因此不能用各成分药物动力学参数的平均值来代表同理，以叠加成的3个成分统计矩参数为对照总量的药物动力学参数，也就是说总量的药物动力学组，表4中其他的总量统计矩参数与之百分比值列参数不是各成分药物动力学参数的简单相加。于表6。表6补阳还五汤有效单成分及总量统计矩参数分析Table 6Analyses of total quantum statistical moment parameters of the single ingredients in Buyanghuanwu injection and their concen-tration overlapedV.成分AUC,MRT,VRT,CLRSD/%fo.95r103. 6黄芪甲苷98. 79106. 845. 9047, 63103.331, 631080药药苷4.81617.759.947201.126.86174. 2川芎嗪0. 054691.9896. 75948.9907.596.73113. 6平均34, 5471. 1071.17691.6385.580.23108.3103.699.67105.396. 0695. 94101.43.510总量统计矩计算进行验证。根据统计矩原理可知4,15，统计矩为非表4诸成分的统计矩参数的×2=1962》xo.01(25）=44.314,说明诸成分与方法的总量统计矩房室模型，能表达线性房室、非线性Michaelis-Ment-参数差异显著（P<0.01）。由表6可知，以川芎嗪en和不规则曲线模型，能用已知曲线通过统计矩公的参数AUC，的比值最小，为0.0546%；以芍药苷的式直接计算，也可通过各取血点用梯形法或抛物线CL，的比值最大，为1080%，两者相差1.7万倍，黄法进行计算，故能用来表达中药及复方单成分与多芪甲者、为药昔、川芎嗪组的RSD分别为31.63%，成分药物动力学整体行为：由总量零、一、二阶矩及174.2%，113.6%，这也说明中药复方中各成分的药衍生的总体表观清除率、总体表观分布容积构成主物动力学总量统计矩参数相差悬殊，彼此参数之间要总量统计矩参数体系。由表5，6的对比数据可不能相互代替。平均组与叠加组的×=212.3知，只有按总量统计矩方法对各成分的药物动力学Xa.o1(5)=15.086，其RSD也为108.3%，也大于一般参数进行加合计算，所得总量统计矩参数结果才能药物动力学实验误差5%，因此也不能用各成分统与诸成分叠加的总量的统计矩参数一致，才能反映计矩参数的平均值来代表总量统计矩参数，但按式多成分的整体浓度量变的动力学规律。4.2总量统计矩的特点由统计矩原理和式（6）~（13）总量统计矩公式计算的各参数的RSD(5）（13）可知：总量统计矩具有2大特点：①为3.510%，小于一般药物动力学实验误差5%，同时两者的×=1.0061≤X0.01(s)=15.086，因此只有具加合性，单个成分经统计矩处能构成总量统计矩按总量统计矩计算方法对各成分的统计矩参数进行参数：反之，总量统计矩又能分解成单个成分统计加合计算，所得结果才能与成分叠加的总量的统计矩：总量统计矩减去一部分成分总量统计矩会得另矩参数一致，也就是说本文所建立的式（1）～（13）一部分总量统计矩，利用这一原理可进行血样空白的总量统计矩的数学模型能客观真实地反映各单个于扰的消除：②具偶联性，能与多维变量曲线偶联，成分总量的药物动力学行为。建立起多维曲线下的总量统计矩参数。如与HPLC4分析与讨论指纹图谱与药动学偶联得保留时间、药物动力学代4.1中药及复方多成分药物动力学研究中药及谢时间二维变量对药物浓度曲线总量统计矩参数：复方为多成分复杂体系，其产生的整体效应为诸有HPLC指纹图谱与效应偶联得保留时间、药物浓度效成（辅助成分）在体内作用的综合表现。对于中（响应值）二维变量对药物效应曲线总量统计矩参数，前者称谱动学、后者称谱效学16-7，如HPLC指药多成分的药物动力学研究目前仍采用单一成分的房室模型研究方法，这必然产生前述3个问题，其中纹图谱与药动学、效应偶联得保留时间、药物动力学最关键的问题是怎样表达多成分药物动力学的数学代谢时间、药物浓度（峰响应值）三维变量对效应曲模型和多成分与效应的关联性问题4，前文解决中线总量统计矩参数，称谱效动力学。因此，总量统计药及复方药物动力学数学模型问题，本文采用实验计矩可用于任何多成分变量的中心和离散性分析，·260

因此不能用各成分药物动力学参数的平均值来代表 总量的药物动力学参数，也就是说总量的药物动力学 参数不是各成分药物动力学参数的简单相加。 同理，以叠加成的 3 个成分统计矩参数为对照 组，表 4 中其他的总量统计矩参数与之百分比值列 于表 6。 表 6 补阳还五汤有效单成分及总量统计矩参数分析 Table 6 Analyses of total quantum statistical moment parameters of the single ingredients in Buyanghuanwu injection and their concen￾tration overlaped 成分 AUCt MRTt VRTt CLt Vt t0． 95t RSD/% 黄芪甲苷 98. 79 103. 6 106. 8 45. 90 47. 63 103. 3 31. 63 芍药苷 4. 816 17. 75 9. 947 1080 201. 1 26. 86 174. 2 川芎嗪 0. 0546 91. 98 96. 75 948. 9 907. 5 96. 73 113. 6 平均 34. 54 71. 10 71. 17 691. 6 385. 5 80. 23 108. 3 总量统计矩计算 103. 6 99. 67 105. 3 96. 06 95. 94 101. 4 3. 510 表 4 诸成分的统计矩参数的 χ 2 = 1 962  χ 2 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% ，这也说明中药复方中各成分的药 物动力学总量统计矩参数相差悬殊，彼此参数之间 不能相互代替。平均组与叠加组的 χ 2 = 212. 3  χ 2 0. 01( 5) = 15. 086，其 RSD 也为 108. 3% ，也大于一般 药物动力学实验误差 5% ，因此也不能用各成分统 计矩参数的平均值来代表总量统计矩参数，但按式 ( 6) ～ ( 13) 总量统计矩公式计算的各参数的 RSD 为 3. 510% ，小于一般药物动力学实验误差 5% ，同 时两者的 χ 2 = 1. 006 1χ 2 0． 01( 5) = 15. 086，因此只有 按总量统计矩计算方法对各成分的统计矩参数进行 加合计算，所得结果才能与成分叠加的总量的统计 矩参数一致，也就是说本文所建立的式( 1) ～ ( 13) 的总量统计矩的数学模型能客观真实地反映各单个 成分总量的药物动力学行为。 4 分析与讨论 4. 1 中药及复方多成分药物动力学研究 中药及 复方为多成分复杂体系，其产生的整体效应为诸有 效成( 辅助成分) 在体内作用的综合表现。对于中 药多成分的药物动力学研究目前仍采用单一成分的 房室模型研究方法，这必然产生前述 3 个问题，其中 最关键的问题是怎样表达多成分药物动力学的数学 模型和多成分与效应的关联性问题［14］，前文解决中 药及复方药物动力学数学模型问题，本文采用实验 进行验证。根据统计矩原理可知［4，15］，统计矩为非 房室模型，能表达线性房室、非线性 Michaelis-Ment￾en 和不规则曲线模型，能用已知曲线通过统计矩公 式直接计算，也可通过各取血点用梯形法或抛物线 法进行计算，故能用来表达中药及复方单成分与多 成分药物动力学整体行为: 由总量零、一、二阶矩及 衍生的总体表观清除率、总体表观分布容积构成主 要总量统计矩参数体系。由表 5，6 的对比数据可 知，只有按总量统计矩方法对各成分的药物动力学 参数进行加合计算，所得总量统计矩参数结果才能 与诸成分叠加的总量的统计矩参数一致，才能反映 多成分的整体浓度量变的动力学规律。 4. 2 总量统计矩的特点 由统计矩原理和式 ( 5) ～ ( 13) 可知: 总量统计矩具有 2 大特点［15］: ① 具加合性，单个成分经统计矩处能构成总量统计矩 参数; 反之，总量统计矩又能分解成单个成分统计 矩; 总量统计矩减去一部分成分总量统计矩会得另 一部分总量统计矩，利用这一原理可进行血样空白 干扰的消除; ②具偶联性，能与多维变量曲线偶联， 建立起多维曲线下的总量统计矩参数。如与 HPLC 指纹图谱与药动学偶联得保留时间、药物动力学代 谢时间二维变量对药物浓度曲线总量统计矩参数; HPLC 指纹图谱与效应偶联得保留时间、药物浓度 ( 响应值) 二维变量对药物效应曲线总量统计矩参 数，前者称谱动学、后者称谱效学［16-17］，如 HPLC 指 纹图谱与药动学、效应偶联得保留时间、药物动力学 代谢时间、药物浓度( 峰响应值) 三维变量对效应曲 线总量统计矩参数，称谱效动力学。因此，总量统计 计矩可用于任何多成分变量的中心和离散性分析， ·260· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue 2中国中药杂志January,20132013年1月五汤的研究[D]．成都：成都中医药大学，2006:6是研究中药多成分体系动态量变、效应变化的有力[3]Ludbrook J. Statistical techniques for comparing measurers and的数学工具。method of measurement:a critical review D. Clin Exp Pharmacol4.3总量统计矩药物动力学模型使用时应注意的Physiol,2002,29(7):527问题由式（1）～式（11）的推导过程可知，总量统[4]魏树礼。生物药剂学与药物动力学[M]．北京：北京医科大计矩药物动力学模型是建立在药物浓度化学动力学学&中国协和医科大学联合出版社，1997：1204.[5]Veng-pedersin , Noncoupartmentaly-based pharmacokinetics mod-量变的模型的基础上，其适用性在于能用本模型表eling DJ Adv Deliv Rev, 2001,48 (2 /3) :265.达体内各成分的浓度整体量变。不能简单地用本模[6] 贺福元，罗杰英，邓凯文。中药复方动力学数学模型-总量统型计算出来的总量来推断药物生物效应（药效与毒计矩法的研究，世界科学技术——中医药现代化，2006，效)的大小，原因在于各成分的效应与毒效应系数9 (6) :13.不同，因此多成分的总量效应并不简单地与总浓度[7] 贺福元，邓凯文，邹欢，等，中药复方谱动学与谱效动力学差呈相关性，不能用单个成分的浓度相加来表达。总异性的研究中国中药杂志，2011,36(2)：136.[8]陈冬，杨洁红．补阳还五汤抗脑缺血作用的研究进展.中成分浓度量变统计矩参数与总效应统计矩参数并不华中医药学刊2010,28（1)：72.重迭，只有将各成分的浓度量变乘上效应系数后，将[9]数学手册编写组。数学手册[M].8版，北京：高等教育出版浓度转化为效应后，才能用总量统计矩计算得出总社，2000:219，531效应的总效应统计矩参数，才能用来评价药物总效[10]段晓鹏，贺福元，周晋，等，补阳还五汤指纹图谱加合性的研应，这时总量统计矩对单个成分效应的整合与对单究，中国中药杂志，2011,36（23）：3247.[i]李晓亮，邵旭，牛海军，等：高效液相色谱潜—蒸发光散射检个成分药物动力学整合一样，其总效应统计矩参数测法测定补阳还五苷酮胶囊中黄芪甲苷的含量．安徽医能指导临床用，这方而的详细阐明请读者参考相关药，2010,14(11)：1286文献]。中药多成的药物动力学与效应动力学不[12]邓凯文，贺福元，刘文龙、RP-HPLC法测定大鼠血浆中补阳对应是有别于西药单一成分的药物动力学与药效动还五汤及总苷中芍药苷的含量及药代动力学研究.中药力学对应最为明显的地方。材，2010,33(11)：1768[13]贺福元，刘文龙，贺庆平，等，反相高效液相色谱法测定大鼠本文详细阐明了中药复方多成分药物动力学总血浆中补阳还五汤总生物碱中川芎嗪的含量及药代动力学量统计矩数学模型的建立原理、实验与计算方法，适研究。中国药物与临床，2005,5（12）：920.用于已知能测准的中药代表性成分，但中药绝大多[14]贺福元，罗杰英。中药复方药物动力学研究方法现状及研究数成分的没有对照品或为未知成分，因此怎样利用方法探讨中草药，2005，35（10）：1582指纹图谱测定中药复方多成分技术，并与药物动力[15]贺福元，周宏灏，邓凯文，等，指纹图谱的一种定性定量研究新方法：总量统计矩分析法D).药学学报，2008，43（2）：学关联构成的谱动学来测定中药复方表观药物动力195.学参数，亦谱动学参数将另文详述。[6]】贺福元，罗杰英，邓凯文，等。中药谱效学研究方法初探D][参考文献]世界科学技术—中医药现代化，2004,6（6）：44[1] Loew D, Kaszhin. Approaching the problem of bioequivalence of[17]贺福元，邓凯文，刘文龙，等：中药与人体的作用本质“网通herbal medicinal products UJ : Phytother Res, 2002, 16 (8) :虹势”的多重遗传谱效动力学.中国实验方剂学杂志，705.2011,17(2) :240.[2］贺福元。中药复方药物动力学数学模型的建立及对补阳还·261 :

是研究中药多成分体系动态量变、效应变化的有力 的数学工具。 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］ Veng-pedersin . Noncoupartmentaly-based pharmacokinetics mod￾eling［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］ 邓凯文，贺福元，刘文龙 . RP-HPLC 法测定大鼠血浆中补阳 还五汤及总苷中芍药苷的含量及药代动力学研究［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. ·261· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013

第38卷第2期Vol 38,Issue 2中国中药杂志January, 20132013年1月Experimental studies on pharmacokinetics of three components inBuyanghuanwu injection on base of total quantum statistical momentHE Fu-yuan'-2, 3, DENG Kai-wen', LIU Wen-ong', SHI JiHian',WU Yong, LIU Wei', HE Qing-ping', LI Bo'(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 MedicineHunan University of Tradition Chinese Medicine, Changsha 410208,China;4. The First Afinity 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, tet-ramethylpyrazine in Buyanghuanwu injection, in order to establish a pharmacokinetic experimental method with multi-component tradi-tional Chinese medicine(TCM) compound system. Method: The RP-HPLC was adopted, with the chromatographic column of Crs,4.6 mm ×250 mm, 5 μm.As for astragaloside IV, the ELSD detector was adopted with acetonitrile-water(35: 65) as the mobile phaseat1mL·min'; the pressure of column was (15.0±0.408)MPa, the column temperature was 30 C.Regarding paeoniflorin andtetramethylpyrazine, the detection of wavelengths was 254 nm, with acetonitrile-water (35: 65) as the mobile phase at 1 mL·minthe column pressure of (15. 17 ±0.41) MPa The pharmacokinetic parameters for single component were dealt with DAS and the totalquantum statistical moment(TQSM) parameters were calculated using formulations. Result: All of the three components followed thetwo compartmental pharkacokinetic model (P<0. O1) 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 TQSMpharmacokinetic parameters of the three components in Buyanghuanwu injection showed that AUC,, MRT,, VRT,, CL,,V, were(119.8±27.20)gminL-l(210.0±54.49)min，(5.608±2.723)×10min2,(0.3196±0.0688)mLmin-.kg*and(64.12 ±8.243) mL .kg*, respectively, suggesting that the half-life time for the three components were (145.5 ±37.76) min and95% of them were metabolized within 0-674.2 min.Conclusion: The TQSM can be used to study pharmacokinetic parameters ofmulti-component TCM compound, because the method can characterize the pharmacokinetic regularity of quantum-ime change in amulti-component system[Key words]total quantum statistical moment; pharmacokinetics; pharmacokinetics with chromatographic fingerprint; Buyanghuanwu injection; astragaloside IV: paeoniflorin; tetramethylpyrazinedoi:10.4268/cjcmm20130222[贵任编辑陈玲]·262·

Experimental studies on pharmacokinetics of three components in Buyanghuanwu injection on base of total quantum statistical moment HE Fu-yuan1，2，3* ，DENG Kai-wen4 ，LIU Wen-long 1 ，SHI Ji-lian1 ， WU Yong 1 ，LIU Wei 1 ，HE Qing-ping 1 ，LI Bo 1 ( 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，tet￾ramethylpyrazine in Buyanghuanwu injection，in order to establish a pharmacokinetic experimental method with multi-component tradi￾tional Chinese medicine( TCM) compound system. Method: The RP-HPLC 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 acetonitrile-water( 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 acetonitrile-water ( 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 compo￾nent 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 half-life time for the three components were ( 145. 5 ± 37. 76) min and 95% of them were metabolized within 0-674. 2 min. Conclusion: The TQSM can be used to study pharmacokinetic parameters of multi-component TCM compound，because the method can characterize the pharmacokinetic regularity of quantum-time change in a multi-component system. ［Key words］ total quantum statistical moment; pharmacokinetics; pharmacokinetics with chromatographic fingerprint; Buyang￾huanwu injection; astragaloside IV; paeoniflorin; tetramethylpyrazine doi: 10． 4268 /cjcmm20130222 ［责任编辑 陈玲］ ·262· 第 38 卷第 2 期 2013 年 1 月 Vol. 38，Issue 2 January，2013