一种中药药代动力学的曲线模拟与预测方法
【技术领域】
[0001] 本发明涉及一种中药药代动力学的曲线模拟与预测方法及其系统，具体而言，涉 及模拟单味药和复方中的一种成分的AUC(时间曲线下面积）-时间曲线，并预测另一种成 分在单味药和复方中的AUC-时间曲线。
【背景技术】
[0002] 中药是复杂的化学体系，存在复杂的药物-药物相互作用，口服单体成分与口服 中药复方制剂，即使其口服剂量完全一样，两者之间会存在很大的差别。然而，ADME相关的 性质预测均是来源于高通量的单体化合物，用它们建模用来预测单体化合物的ADME性质。 但是，因为中药复方存在复杂的药物-药物相互作用，所以这些预测模型对于中药复杂体 系不适用，需要新的研宄思路，建立新的预测模型和方法。
[0003] 公开于该【背景技术】部分的信息仅仅旨在加深对本发明的一般【背景技术】的理解，而 不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
【发明内容】
[0004] 在本发明中，提出一种新的预测策略和预测模型，主要如下：把中药复方化合物， 根据化学结构相似性打分，分为不同的类别，同一类的化合物，被认为结构类似。化合物结 构类似的化合物，可能作用于相同的靶标，表现相似的理化性质、药代动力学性质和药理活 性，为此，提出一个假说，即：口服中药复方制剂，其中，结构类似的化合物，可能相似的代 谢风险，所以，对于中药复方制剂，可以通过获得一个或者少数几个化合物在体内的血药浓 度-时间曲线，从而预测其他化合物在体内的血药浓度-时间曲线，从而获得相关药代动力 学参数，为中药复方药代动力学提供新思路。
[0005] 为此，本发明提供了一种中药药代动力学的曲线模拟与预测方法，可用于单味药 或者中药复方的结构类似化合物的药代动力学的预测，其通过一些步骤，发现和证实化学 结构相似的成分在口服中药方剂的复杂体系后具有相似的药代动力学变化，即相似的r+ 值；然后，在已知成分A在各时刻的r+值，预测其他结构类似成分的auc-时间曲线，从而 建立药代动力学预测新思路和方法。
[0006] 具体地，该方法包括：
[0007] 步骤1.对成分A分别获得其在单味药和复方药中的血药浓度-时间曲线；
[0008] 步骤2.计算成分A在单味药和复方药中各测量时刻的血药浓度-时间曲线下的 面积AUC，进而获得其AUC-时间变化曲线；
[0009] 步骤3.计算成分A在单味药和复方药中各时刻的AUC之和与平均值的比值r+， 其计算公式如下：
[0011] 其中，xjtihyjti)分别为成分A在复方药及单味药中第i个时刻的AUC;
[0012] 步骤4.对成份B做基于化学结构相似性的聚类分析；
[0013] 步骤5.根据聚类分析结果若成分B与成分A同属一类，则可在已知成分A在各时 刻的r+值，成分B在单味药（或复方药）中的AUC值及其在复方药（或单味药）中第3小 时（即t7时刻）的AUC值的情况下，预测成分B在复方药（或单味药）中第4小时（即18 时刻）及以后的AUC-时间曲线，亦即
[0017]其中，xB(ti)、yB(ti)分别为成分B在复方药及单味药中第i个时刻的AUC。
[0018] 在所述步骤1中，通过高效液相色谱质谱/质谱联用来获得单味药和复方中成分 A的血药浓度-时间曲线。
[0019] 在测试血药浓度时，选择甲酸水-甲醇作为流动相，并且采用电喷雾离子化源 ESI，正离子模式检测，以及多反应检测MRM模式对药物离子浓度进行测定。
[0020] 在测试血药浓度时，于样品前处理采用沉淀蛋白方法。
[0021] 所述步骤2中通过梯形法来计算单味药和复方中成分A在各个时间点的血药浓 度-时间曲线下面积，所述梯形法计算公式如下：
[0023] 其中，A(ti)为复方和单味药中成分A在第i时刻的AUC值，C(ti)为复方和单味 药中成分A在第i时刻的血药浓度，并且、=0，々(0) =0，C(0) =0。
[0024] 本发明的方法和装置具有其它特征和优点，这些其它特征和优点将从结合于此的 附图和以下【具体实施方式】中显而易见，或在附图和【具体实施方式】中详细陈述，附图和具体 实施方式共同用于解释本发明的某些原理。
【附图说明】
[0025] 本发明的其它特征和优点将通过下面的参考附图进行的描述而得以显现，在所述 附图中：
[0026] 图1A、图1B和图1C分别显示紫堇碱（I)、a-别隐品碱（II)、延胡索乙素（III)、 四氢小檗碱（IV)、欧前胡素（V)、异欧前胡素（VI)、白当归素（VII)、华法林钠（麗，IS) 和盐酸异丙嗪（IX，IS)的空白血浆样品色谱图、空白血浆标准添加LL0Q色谱图以及给药2 小时后的血浆样品色谱图；
[0027] 图2为SD大鼠口服不同提取物后a-别隐品碱的平均血药浓度-时间曲线图（n =6)；
[0028] 图3为SD大鼠口服不同提取物后延胡索乙素的平均血药浓度-时间曲线图（n= 6);
[0029] 图4为SD大鼠口服不同提取物后四氢小檗碱的平均血药浓度-时间曲线图（n= 6)；
[0030] 图5为SD大鼠口服不同提取物后紫堇碱的平均血药浓度-时间曲线图（n= 6);
[0031] 图6为SD大鼠口服不同提取物后白当归素的平均血药浓度-时间曲线图（n= 6)；
[0032] 图7为SD大鼠口服不同提取物后欧前胡素的平均血药浓度-时间曲线图（n= 6)；
[0033] 图8为SD大鼠口服不同提取物后异欧前胡素的平均血药浓度-时间曲线图（n= 6)；
[0034] 图9显示基于两两之间化学结构相似性打分的元胡止痛方中的7个成分聚类分析 结果；
[0035] 图10显示7个成分在单味药和复方中的AUC-时间变化曲线；
[0036] 图11显示基于AUC-时间曲线变化的元胡止痛方中的7个成分聚类分析的结果； [0037] 图12显示7个成分的r+值的散点图。
【具体实施方式】
[0038] 现在将具体参考本发明的各个实施方案，在附图中和以下的描述中示出这些实施 方案的实例。虽然本发明与示例性实施方案相结合进行描述，但应当了解，本说明书并非旨 在将本发明限制为那些示例性实施方案。相反，本发明旨在不但覆盖这些示例性实施方案， 而且覆盖可以被包括在由所附权利要求所限定的本发明的精神和范围之内的各种替换、修 改、等效形式以及其它实施方案。
[0039] 在本发明的中药药代动力学的曲线模拟与预测方法中，首先建立元胡止痛方中7 个主要活性成分体内分析方法。
[0040] 元胡止痛方主要有延胡索和白芷两味药组成，延胡索中主要成分为生物碱类化合 物，而白芷中主要成分为香豆素类化合物，对两药的药代动力学研宄方面，国内外文献大多 集中在对延胡索乙素的含量测定，但是中药的特点是多成分共同发挥作用，考察其他成分 在体内的含量也是十分有必要的。为了更好的理解元胡止痛方的药理作用机理、体内过程， 有必要建立一种灵敏而又准确的分析方法，以测定生物样品（如血浆、组织匀浆及尿液）中 主要成分的含量。
[0041] 在本发明中将建立同时测定大鼠血浆样品中的a-别隐品碱、延胡索乙素、四氢 小檗碱、紫堇碱、白当归素、欧前胡素、异欧前胡素等7种主要成分的LC-MS-MS方法。
[0042] 本发明所使用的实验仪器：
[0043] WatersAcquityUPLCSystem,WatersQuattroPremierXE，沃特世科技有限公 司）；MIKR0 220R型离心机，德国Hettich公司；
[0044] MX-S型涡旋混合器，美国SCIL0GEX公司；
[0045] XS205型分析天平，瑞士梅特勒-托利多仪器有限公司。
[0046] 本发明所使用的药品与试剂：
[0047] 乙腈（色谱纯，Merck公司）；
[0048] 甲醇（色谱纯，Merck公司）；
[0049]a_别隐品碱（深圳市美荷生物科技有限公司，批号：20100318,纯度：98% ); [0050] 延胡索乙素（中国食品药品检定研宄院，批号：110726-201011，纯度99.9% );
[0051] 紫堇碱（sigma，批号：S449121，纯度：98% );
[0052] 四氢小檗碱（上海融禾医药科技发展有限公司，批号：100516,纯度：98% );
[0053] 白当归素（成都普思生物科技有限公司，批号：RFS-B-100610-04,纯度99. 9% );
[0054] 异欧前胡素（中国食品药品检定研宄院，批号：110827-200407,纯度：98% );
[0055] 欧前胡素（中国食品药品检定研宄院，批号：110826-200712,纯度：99. 9% )。
[0056] 本发明所使用的药材：
[0057] 白]E(AngelicaDahurica(Fisch.exHoffm. )Benth.etHook.f.)产自河北省；
[0058] 延胡索（CorydalisyanhusuoW.T.Wang)产自杭州，
[0059] 均购于河北省安国中药材有限公司，由中国中医科学院中药研宄所何希荣药师鉴 定为药材正品。
[0060] 本发明所使用的实验动物：
[0061] SD大鼠，雄性，体重230±20g，由北京大学医学部实验动物科学部提供，合格证 号：SCXK(京）2009-0017。在温度20-22°C，相对湿度45% -65%，光照/黑暗12h/12h条 件下饲养，自由饮食、饮水，适应性饲养3天后开始实验。
[0062] 发本明所使用的标准品溶液的配制：
[0063] 分别精密称取a_别隐品碱、延胡索乙素、四氢小檗碱、紫堇碱、白当归素、欧前胡 素和异欧前胡素溶于甲醇中，配成储备液浓度分别为lmg/mL，将储备液用甲醇稀释成不同 浓度的混合溶液