数理模型在中药药性及配伍规律研究中的应用概况

应用现代科学语言，如量化思维、数学模型等，表征和阐释其科学内涵是中医药理论传承、创新行之有效的方法之一 ^{［ 5 ］} 。数理模型现在已被广泛应用于中药药性、配伍机制和剂量等研究领域 ^{［ 6-7 ］} 。本文通过归纳、总结近年来基于数理模型的中药药性理论和复方配伍研究，以期为用现代科学语言解释深奥的中药药性理论和复方配伍，为深入阐释其科学内涵，并为中医药国际化和现代化提供研究思路。

1 数理模型在中药药性理论研究中的应用

1.1 四性

在中国传统医学理论中， “ 四性 ” 是指药物的寒、热、温、凉，是相对于与所治疾病的 “ 寒 ”“ 热 ” 性质而言 ^{［ 8-9 ］} 。在中医临床治疗中，常以病症的 “ 寒 ”“ 热 ” 作为用药依据，正如《神农本草经》指出 “ 药有寒热温凉四气（性），疗寒以热药，疗热以寒药 ” ^{［ 9 ］} 。《脾胃论》中也明确记载 “ 一物之内，气味兼有，一药之中，理性具焉，主对治疗，由是而出 ” 。因此，当机体发生病变时，中医的辨证需要辨别出是寒证还是热证才能做出正确的诊断 ^{［ 10 ］} 。

贺福元等 ^{［ 9 ］} 依据热力学第一定律、 Hess 定律建立了 “ 状态函数 ” 数学模型，并对其进行了实验验证。给予实验大鼠黄连、吴茱萸以及甘草，记录给药前机体的热焓以及药物消除后的热焓，并计算差值表明药物的寒热药性，若差值大于 0 ，则表明该药物是寒性药，反之则为热性药。研究结果表明，黄连和吴茱萸分别属于寒性和热性药。该研究结果表明，中药的 “ 四性 ” 可体现在给药后中药代谢热焓与机体热焓改变值之差。

杨岩涛等 ^{［ 11 ］} 依据热力学第一定律建立了病理状态下的中药 “ 四性 ” 数学模型。应用 HR3000 数显热量计测量小柴胡汤加减方及寒证小鼠模型服药前后体内代谢产物的燃烧焓。并依据燃烧焓判断药物的寒热属性。若小柴胡汤配伍组与平性药物之差大于 0 ，则表明该药物为寒性药，反之则为热性药。研究结果表明，芩 - 柴组、芩 - 柴 - 甘草组为寒性药物；小柴胡汤组、参 - 姜 - 夏 - 枣组、参 - 姜 - 枣 - 夏 - 甘草组、芩柴 - 参 - 姜 - 枣 - 夏组为热性药物；甘草组为平性药物。该结果符合中药 “ 四性学说 ” 的临床应用规律 ^{［ 11 ］} 。研究结果也表明 “ 燃烧焓 ” 模型可以用于表征中药 “ 四性 ” 属性。

中药 “ 四性 ” 与其微量元素和无机元素的含量密切相关 ^{［ 8 ， 12 ］} 。吴文莉等 ^{［ 8 ］} 对 105 种中药中 42 种元素建立了中药 “ 四性 ” 的 4 类 Fisher 判别方程，并应用其方程探寻中药 “ 四性 ” 与其元素含量的关系。该研究不仅将中药的寒热温凉 “ 四性 ” 进行定量化分析，并对 “ 四性 ” 与化学元素含量进行了关联。研究结果充分显示了微量元素是中药寒热温凉性质的一种关键因素，其含量能在一定程度上反映中药的性能。研究结果为中药物质基础的深入研究提供了指导。

刘进等 ^{［ 12 ］} 以 K 、 Ca 、 Mg 、 Mn 、 Fe 、 Cu 、 Zn 7 种无机元素为指标，采用支持向量机和留一法识别 193 种中药的药性。作者首先研究这 7 种元素对平性药物和非平性药物的影响；然后，探究这 7 种元素对温热药以及寒凉药的影响。研究发现， K 、 Ca 、 Mg 、 Mn 、 Fe 等元素含量对平性药和非平性药的判别影响较大； Ca 、 Fe 的含量对温热药的识别影响较大。研究表明，无机元素的含量与中药 “ 四性 ” 属性有一定的联系。且基于此联系，可应用特定无机元素的含量判别药物的寒热温凉。

1.2 五味

人们在防治疾病中发现许多药物 “ 味 ” 与其 “ 效 ” 密切相关，因此，用 “ 味 ” 来描述药物的某一或某些功能特性，并应用该属性指导临床合理用药 ^{［ 13 ］} 。中药 “ 五味 ” 是指辛、酸、苦、甘、咸，大多数中药的 “ 五味 ” 是由口尝而得。常用的评价方法是由有经验的药工通过口尝鼻闻判断的。但该方法主观性太强。因此，需要客观、可量化的手段对中药 “ 五味 ” 进行科学的定量表征 ^{［ 14 ］} 。

张培等 ^{［ 15 ］} 基于 197 种中药饮片的现代药理和临床数据，建立了中医 “ 五味 ” 的贝叶斯网络模型，并利用网络拓扑图和条件概率表预测苦味、辛味和甘味的药理活性，明确了苦味、辛味和甘味与相关药理指标的联系，实现了根据药理药效指标判定中药 “ 五味 ” 属性。研究结果表明，贝叶斯网络模型不仅能有效地预测中药的 “ 五味 ” 属性，还可以直观地表达 “ 五味 ” 与药效作用的复杂关系。

刘杨茜 ^{［ 14 ］} 采集了 40 种中药饮片和 33 种对照药材的全息化学指纹图谱，计算中药饮片与对照药材的欧氏距离相似度，依据该相似度对其药味进行量化分析，若中药饮片与对照药材的欧氏距离越近，表明药味相似。作者还应用多元回归分析方法建立了中药主成分 “ 五味 ” 量化指标，绘制了中药 “ 五味 ” 五行图，建立了可同时定性、定量表征中药 “ 五味 ” 的预测模型。研究结果为进一步阐明中药 “ 五味 ” 的化学物质基础和完善 “ 五味 ” 理论研究提供新思路。

1.3 归经

中药归经学说是中药药性理论的核心内容之一，是中医药理论体系的重要组成部分 ^{［ 16 ］} 。 “ 归经 ” 指药物对脏腑、经络病变部位的选择性作用，是历代中医药学家在长期与疾病作斗争的医疗实践中，总结出来的一种认识药物作用部位的客观规律。随着中医药现代化、国际化的进程加快，传统的中药归经体系尚需进一步完善和拓展 ^{［ 17-19 ］} ，其科学内涵也亟需应用现代科学语言进行诠释和表征。

李芳等 ^{［ 20 ］} 将方剂归经量化并结合方剂中各药物所占的比例，建立方剂归经系数模型。共研究了 21 个经典方剂的归经系数。以麻黄汤为例，计算得到了麻黄汤主要归肺经，与其 “ 发汗解表、宣肺平喘 ” 的功效一致。所建立的模型为中药归经的现代化研究提供了方法学参考、奠定了基础。

李芳等 ^{［ 21 ］} 运用层次分析法进一步研究方剂的归经强度，以药物归经量化和相对药量为依据建立了判断矩阵，构建了方剂归经权重模型。以四君子汤为例进行研究，发现四君子汤中归经权重值最大的是人参，且四君子汤主要归脾经。研究结果与四君子汤主治脾胃气虚相吻合。所建模型即可以量化方剂的归经强度，还能表征各药物在方剂结构中的地位和意义。研究结果对探索方用、药用的规律性等具有重要意义。

1.4 毒性研究

传统医学概念中，广义上的 “ 毒 ” 指的是药物的偏性；狭义的 “ 毒 ” 是药物对人体产生的不良影响及损害，即药物毒性和副反应 ^{［ 22 ］} 。随着中药的广泛应用，国内外关于中药引起不良反应的报道逐渐增多，中药安全性越来越受到人们的重视 ^{［ 23 ］} 。

刘红杰等 ^{［ 23 ］} 通过文献检索，将肾毒性和非肾毒性中药纳入数据库，利用神经网络模型预测分析中药肾毒性，并将其与 “ 四气 ”“ 五味 ”“ 归经 ” 等属性因素进行相关性检验，筛选相关性变量因素，构建神经网络模型，并用受试者工作特征曲线进行验证。研究结果表明，热性相对于其他因素对中药肾毒性的影响更加明显，远高于其他因素。所建模型即可为中药肾毒性研究提供一定的参考，也可以用于中药肾毒性的预测。

刘红杰等 ^{［ 24 ］} 进一步用 Logistic 回归预测模型对 111 种肾毒性和 398 种非肾毒性中药的 “ 四性 ”“ 五味 ”“ 归经 ” 进行模型构建并用受试者工作特征曲线进行验证。研究发现，中药毒性与中药 “ 四性 ”“ 五味 ” 等属性相关，但与 “ 归经 ” 没有明显关联。同时，还发现 “ 热性 ” 和 “ 苦味 ” 是中药肾毒性的危险因素，而 “ 平性 ” 和 “ 甘味 ” 是中药肾毒性的保护因素。研究结果表明，所建立的 Logistic 回归预测模型可有效预测中药肾毒性，为深入研究中药肾毒性作用机制和科学内涵奠定了基础。

为了考察辛味中药挥发油成分的皮肤细胞毒性，杨文国等 ^{［ 25 ］} 运用非线性典型相关分析建立了一种综合中药药性特征与挥发油皮肤细胞毒性相关的数学模型。通过研究 33 种辛味药材的 “ 四性 ”“ 五味 ”“ 归经 ” 与细胞毒性的关系，研究发现，中药的药性特征尤其是 “ 归经 ” 属性会影响辛味中药挥发油成分的皮肤细胞毒性。该研究及其研究结果为从中药 “ 属性 ” 角度深入揭示中药细胞毒性奠定了基础。

2 数理模型在中药复方配伍研究中的应用

2.1 单一数理模型

2.1.1 粗糙集 在中药复方研究中，粗糙集把中药复方看作一种知识，运用复杂的数学原理对其进行挖掘和学习，处理其所包含的大量非线性数据 ^{［ 26 ］} 。粗糙集在研究复方中各药物重要性时侧重于各味药的贡献，并且可以区分出该药物所对应疾病的证候类型。因此，粗糙集为中药配伍规律的研究提供了新视角。

刘娟等 ^{［ 27 ］} 运用粗糙集对治疗乙肝的复方进行数据挖掘，结果表明补气药在治疗乙肝的复方中发挥重要作用。其研究还发现当归是治疗乙肝的重要药物且集中治疗气滞血瘀型乙肝。应用粗糙集研究中药复方配伍不仅可以筛选出复方中的关键药物，而且还能区分出该复方主治疾病的证候类型。但粗糙集的使用需要结合先验知识对药物的属性进行约简。

2.1.2 网络药理学 基于网络药理学的组分中药研究策略很好地体现组分中药配伍的合理性、有效组分间的协同、相加、拮抗等药物相互作用特点。

于金高 ^{［ 28 ］} 通过动物实验和细胞模型发现水通道靶点可能是甘草 - 芜花之间的禁忌作用的主要靶点。同时，作者通过网络药理学和分子对接模型验证了水通道靶点在甘草 - 芜花配伍禁忌中的重要作用。网络药理学通过构建 “ 药材 - 成分 - 靶点 ” 的相互作用网络，推断复方化学成分与疾病靶点之间的作用关系，从而预测复方中的关键化学成分及作用靶点。同时，网络药理学还可以通过分子对接模型与药理实验的结果相互佐证。因此，网络药理学研究模式与中医药理论的有效结合为中药复方配伍研究提供了新的研究模式和视角。但网络药理学的研究结论需要临床和药理实验的进一步验证。

2.1.3 神经网络模型 神经网络用于中药复方配伍的药性特征预测、药对配伍和组效关系 3 方面研究具有无可比拟的优势：其参数比统计方法的自变量多；无论变量是何种类型，是否满足正态性、独立性等条件均可应用 ^{［ 29 ］} 。

芦霞 ^{［ 30 ］} 用 50 种治疗感冒的中药复方建立了感冒中药药性预测的神经网络模型，并用另外 10 种感冒复方对所建立的神经网络模型进行验证。结果显示，所建立的神经网络模型可以根据中药的基本药性（如性、味、归经等）对感冒复方的药物性能做出准确的预测，如神经网络模型预测出板蓝根在 “ 清热解毒，润肺止咳 ” 方面发挥疗效，与板蓝根的疗效相符。该研究将中药药性和神经网络数理模型结合，对于中药配伍的研究具有很好的指导意义。神经网络是处理非线性数据的最佳工具，而且其自变量相对于其他模型有着无法比拟的优势。神经网络的特点、能力与中药复方的 “ 多成分、多靶点、多途径 ” 的特点相吻合。但由于神经网络的算法不能解释影响变量的因素，无法检验变量的权重系数 ^{［ 23 ］} ，仍需与其他数学方法和模型结合，进行更全面、深入的分析。

2.1.4 贝叶斯网络 贝叶斯网络用概率分布表示依赖关系的强弱，是不确定性知识的表达、推理和建模最有效的方法之一 ^{［ 29 ］} 。

王晓燕 ^{［ 31 ］} 构建了贝叶斯网络，并应用此网络分析了中药寒、热药性与药性特征标记之间的配伍模式。选用中医公认的典型寒、热药物，对其化学成分以及采集正常大鼠、虚热证、虚寒证大鼠给药后的药效指标，探讨寒热药性、药效与物质成分的关系。结果表明，寒热中药不仅能使机体病理状态下的寒热得到纠正，作用于正常机体后也能引起正常机体寒热效应，并且相同药效中药的物质组分在配伍模式上存在共性。贝叶斯网络依据已知的中药信息预测中药配伍中隐含的复杂关系。应用这一数理模型可以挖掘经典中药复方和中药大品种的临床价值，指导临床实践，提高临床疗效，并为创制中药新药奠定基础。但其仍需药理实验、临床试验结果的验证。

2.1.5 “ 阴阳球 - 八纲三级结构系统 ” 数理模型 郭述强 ^{［ 32 ］} 基于对中医阴阳太极图的比较和认识，构建了 “ 阴阳球 ” 模型，并将其与中医理论中的八纲（寒、热、虚、实、表、里、阴、阳）衔接在一起，构建了 “ 阴阳球 - 八纲三级结构系统数理模型 ” 。同时，应用所建立的模型探究《脾胃论》所记载的中药复方蕴含的数学规律：首先，根据复方的 “ 药力 ” 构建 “ 药力 ” 数据库， “ 药力 ” 即该中药对于它所对应的 “ 证 ” 的作用力大小；然后，将《脾胃论》中筛选出的药方导入 “ 阴阳球 - 八纲三级结构系统数理模型 ” 的输入板块，以空间定位的方式研究复方的药、方、证在三维空间的变化轨迹。此研究应用数学思维模式对传统中医学的辨证论治思想进行清晰、明确地阐述 ^{［ 32 ］} 。

2.2 复合数理模型

单一数理模型的分析结果虽然有助于明确中药复方中的重要药材、配伍禁忌的关键靶点以及寒热药性的配伍，但无法实现量化表征。因此，不同数理模型或方法的结合是实现量化分析行之有效的手段之一，有助于深层次地、全面地阐释中药复方配伍的科学内涵。

2.2.1 正交设计结合模型 曾勇等 ^{［ 33 ］} 将正交设计与路径分析结合，建立了