浙中医：基于UPLC-Q-TOF-MS代谢组学技术研究淫羊藿-川芎对骨关节炎大鼠的调控作用

骨关节炎（osteoarthritis，OA）是一种以关节软骨的退行性病变和继发性骨质增生为主要病变的慢性关节疾病，以关节疼痛、滑膜炎症和活动受限为主要表现[1]。目前化学药治疗以缓解疼痛、抗炎或改善功能状态为主，但不良反应明显，不宜长期应用，因此尚缺乏切实有效的治疗手段和药物[2]。中医临床多以“补肾活血”为治疗准则[3]，常以出自《太平圣惠方》的补肾活血药对淫羊藿-川芎配伍使用。淫羊藿补肾壮阳而治肾虚之根本，入肝肾祛风湿强筋骨；川芎活血祛瘀、祛风止痛治疗血瘀气滞[4-5]。

OA疾病发展过程及其治疗与代谢密切相关，代谢组学的研究可用于发现疾病的发病机制及探究药对干预下的调控机制[6-8]。研究表明，在OA发展过程中，软骨细胞代谢发生改变，基质金属蛋白酶-13（matrix metalloproteinase-13，MMP-13）参与软骨的降解，滑膜细胞和软骨细胞局部合成的促炎细胞因子白细胞介素-1β（interleukin-1β，IL-1β）和肿瘤坏死因子-α（tumor necrosis factor-α，TNF-α）造成关节软骨持续损伤[9-11]。药理研究表明，淫羊藿可调节软骨细胞代谢，促进软骨基质合成和细胞增殖[12-13]。川芎的有效成分可抑制炎症因子IL-1β、TNF-α、MMP-13的表达[14-15]。

淫羊藿 - 川芎对 OA 的作用机制尚不清楚，因此本研究基于超高效液相色谱 - 四极杆 - 飞行时间质谱联用仪（ ultra performance liquid chromatography-quadrupole time-of-flight mass spectrometry ， UPLC-Q-TOF-MS ）的代谢组学方法考察淫羊藿 - 川芎配伍及其单药治疗 OA 的内源性代谢产物变化及代谢轮廓，明确其发挥药效的代谢调控通路。

1 材料

1.1 动物

清洁级 SD 大鼠 66 只，体质量（ 190 ± 10 ） g ，雌雄各半，购自上海斯莱克实验动物有限公司，生产许可证号 SCXK （沪） 2013-0018 ，使用许可证号 SCXK （浙） 2013-0184 。温度和相对湿度分别维持在（ 22 ± 1 ）℃和（ 60 ± 10 ） % 。动物实验经浙江中医药大学伦理委员会批准（批准号 ZSLL-2018-009 ）。

1.2 药材

淫羊藿（批号 141117 ）和川芎（批号 140527 ）均购自华东医药股份有限公司，经浙江中医药大学药学院张汝松教授分别鉴定为小檗科植物淫羊藿 Epimedium brevicornu Maxim. 的干燥叶、伞形科植物川芎 Ligusticum chuanxiong Hort. 的干燥根茎。

1.3 药品与试剂

塞来昔布胶囊（批号 R55865 ）购自美国辉瑞制药有限公司；盐酸氨基葡萄糖胶囊（批号 1504174 ）购自海南澳美华制药有限公司； MMP-13 ELISA 试剂盒（批号 20151201A ）购自杭州乐峰生物科技有限公司；色谱纯甲醇和乙腈购自德国 Merck 公司；分析级乙酸（批号 200903006 ）购自杭州化学试剂有限公司；水为超纯水；青霉素（批号 F4056202 ）购自华北制药股份有限公司。

1.4 仪器

Waters Acquity UPLC 、 Waters Q-TOF-Premier （美国 Waters 公司）； DMI3000B 型显微镜（德国 Leica 公司）； Power Wave 340 型酶标仪（美国 Bio-Tek 公司）； WH-866 型涡旋震荡器（华利达公司）； SENCOR 系列旋转蒸发仪（上海申生科技有限公司）； 5804R 型高速冷冻离心机、移液枪（德国 Eppendorf 公司）。

2 方法

2.1 给药 溶液的制备

2.1.1 淫羊藿溶液的制备与分析 根据课题组前期研究确定最佳提取工艺进行药液制备 [16] ，取淫羊藿饮片 5000 g ，置于多功能提取罐中，加 12 倍量 50% 乙醇，室温浸泡 1 h ，加热回流提取 1 h ，提取 2 次，合并药液，滤过，浓缩成 2 g 生药 /mL 的药液，置 4 ℃冰箱保存备用。浓缩的药液稀释 4 倍，进样检测分析 [17] 。

2.1.2 川芎溶液的制备与分析 取川芎饮片 5000 g ， 按“ 2.1.1 ”项下方法制备与分析。

2.1.3 淫羊藿 - 川芎溶液的制备与分析 取淫羊藿饮片和川芎饮片各 2500 g ，按“ 2.1.1 ”项下方法制备与分析。

2.2 动物分组、造模、给药及取材

66 只大鼠适应性饲养 1 周后，按体质量随机分为正常组（ n ＝ 8 ）、假手术组（ n ＝ 8 ）、模型组（ n ＝ 8 ）、病理经时评价组（ n ＝ 10 ）、淫羊藿组（ n ＝ 8 ）、川芎组（ n ＝ 8 ）、淫羊藿 - 川芎组（ n ＝ 8 ）、阳性药物组（ n ＝ 8 ）。除正常组和假手术组外，其余各组均采用交叉韧带 - 半月板损伤术造模。假手术组切开关节囊，打开关节腔，不切断前交叉韧带和损伤半月板。术后所有手术大鼠 im 青霉素 [ 40 000 U/( 只 · d) ] 。

造模 6 周后，各组 ig 相应药物， 1 次 /d ，连续 8 周。淫羊藿组、川芎组、淫羊藿 - 川芎组分别 ig 相应给药溶液（ 2 g 生药 /mL ），阳性药物组 ig 盐酸氨基葡萄糖（ 0.006 g/mL ）和塞来昔布（ 0.05 g/mL ），正常组和模型组 ig 等体积蒸馏水（ 20 mL/kg ）。

病理经时评价组于第 0 、 2 、 4 、 6 、 8 周随机取 2 只大鼠，颈椎脱臼处死，其余各组大鼠于第 14 周处死。取关节软骨，福尔马林溶液中固定， EDTA 脱钙、脱水、石蜡固定包埋，切片，用于苏木素 - 伊红（ hematoxylin-eosin staining ， HE ）和番红 O- 固绿（ safranin O/fast green staining ， SO ）染色。

2.3 样品的采集与制备

于第 0 、 2 、 4 、 6 、 8 周，采集假手术组血液样本，于第 6 、 8 、 10 、 12 、 14 周采集各给药组血液样本，于第 0 、 2 、 4 、 6 、 8 、 10 、 12 、 14 周正常组、模型组血液样本置于抗凝管中，摇匀， 4 ℃、 3000 r/min 离心 10 min ，分离得到血浆，置 − 80 ℃冰箱保存备用，后续用于 UPLC-Q-TOF-MS 分析和 ELISA 检测。

取血浆样品，解冻至室温，移取 200 µL 血浆样品，加入 600 µL 乙腈。涡旋 1 min ， 4 ℃、 12 000 r/min 离心 10 min ，取上清液上机检测。吸取各样品上清液 50 μL 混合成质控样品。

2.4 色谱和质谱条件

2.4.1 色谱条件 Waters AcquityTM C18 色谱柱（ 100 mm × 2.1 mm ， 1.7 µm ），流动相为含 0.1% 醋酸的水溶液（ A ） - 含 0.1% 醋酸的甲醇溶液（ B ），梯度洗脱： 0 ～ 0.5 min ， 3% B ； 0.5 ～ 8 min ， 3% ～ 80% B ； 8 ～ 16 min ， 80% ～ 98% B ； 16 ～ 21 min ， 98% B ； 21 ～ 21.5 min ， 98% ～ 100% B ； 21.5 ～ 24 min ， 100% B ； 24 ～ 25 min ， 100% ～ 3% B ； 25 ～ 30 min ， 3% B 。体积流量为 300 µL/min ；柱温为 35 ℃；进样体积为 5 µL ；自动进样器温度为 4 ℃。

2.4.2 质谱条件 采用电喷雾离子源（ ESI ），正离子扫描模式，毛细管电压为 3.0 kV ，锥孔电压为 40 V ，碰撞电压为 10 ～ 60 V ，离子源温度为 120 ℃，脱溶剂气温为 350 ℃，脱溶剂气体积流量为 450 L/h ，碰撞气体为氩气。扫描范围 m /z 50 ～ 1000 ，扫描时间 0.3 s ，扫描间隔 0.02 s 。

TOF 离子飞行方式采用“ V ”模式。使用 0.5 ng/µL 的亮氨酸 - 脑啡肽溶液进行相对分子质量测定校正。

2.5 数据处理与统计学分析

2.5.1 代谢组学数据处理 将原始数据导入 MassLynx v4.1 and MarkerLynx 4.1 软件，经过峰识别、峰对齐和归一化等处理后，得到由保留时间、质荷比和峰面积组成的数据矩阵，导入 SIMCA-P 12.0 进行无监督的主成分分析（ principal component analysis ， PCA ）和有监督的正交偏最小二乘 - 判别分析（ orthogonal partial least squares-discriminant analysis ， OPLS-DA ）。根据 S-plot 图中的变量投影重要性（ variable importance in projection ， VIP ）＞ 1 和独立样本 t 检验的 P ＜ 0.05 ，筛选差异代谢物作为潜在生物标志物。根据选取物质精确的质荷比，检索 HMDB （ http://www.hmdb.ca/ ）和 METLIN （ http://metlin.scripps.edu ）对代谢物进行匹配鉴定。将鉴定所得潜在生物标志物导入京都基因与基因组百科全书（ Kyoto encyclopedia of genes and genomes ， KEGG ）数据库（ http://www.kegg.ip ），明确相关代谢途径，并分析其调控作用。

2.5.2 统计学分析 数据均用 表示，采用 SPSS 25.0 软件进行分析，两组间比较采用 t 检验。

3 结果

3.1 淫羊藿、川芎及药对的主要成分含量测定

淫羊藿有效成分绿原酸、朝藿定 A 、朝藿定 B 、朝藿定 C 、淫羊藿苷、宝藿苷 I ，其质量分数分别为 0.618 0 、 1.480 、 1.968 、 2.283 、 4.494 、 0.840 0 mg/g 。川芎有效成分绿原酸、藁本内酯，其质量分数分别为 0.502 2 、 11.39 mg/g 。淫羊藿 - 川芎中绿原酸、阿魏酸、朝藿定 A 、朝藿定 B 、朝藿定 C 、淫羊藿苷、宝藿苷 I 、藁本内酯质量分数分别为 1.318 、 0.502 6 、 1.562 、 2.104 、 1.988 、 5.869 、 0.858 9 、 9.216 mg/g 。

3.2 病理结果

3.2.1 OA 模型大鼠关节软骨组织病理学经时评价 如图 1 所示，正常组及假手术组关节面光滑平整，颜色红润鲜艳， SO 着色均匀，软骨细胞排列整齐，四层结构及潮线清晰可见。造模第 2 周模型组大鼠关节软骨较正常组无明显变化，仅在软骨层表面有较小、浅的裂隙；第 4 周关节软骨层变薄，表面粗糙，局部软骨出现片状脱落，部分软骨基质不着色；第 6 周关节软骨层严重破坏，各层软骨细胞簇集生长，排列紊乱；第 8 周病变加剧，可见软骨细胞坏死，潮线更加不规则或完全消失，甚至可见软骨下骨暴露。提示第 6 周 OA 模型建立成功。

3.2.2 淫羊藿 - 川芎对 OA 大鼠膝关节组织病理学的影响 如图 2 所示，正常组关节面光滑平整，颜色红润，软骨细胞排列整齐，四层结构及潮线清晰可见， SO 着色均匀，颜色鲜艳。模型组关节软骨层明显变薄，软骨部分脱落，软骨细胞簇集，排列紊乱，潮线多不完整或完全消失， SO 着色重度减弱，可见坏死崩解细胞，甚至可见软骨下骨质外露。各给药组软骨细胞排列有序，软骨细胞聚集生长现象逐渐消失，软骨层逐渐恢复。淫羊藿 - 川芎组和阳性药物组关节软骨层光滑，无明显软骨细胞聚集， SO 着色均匀，但阳性药物组软骨细胞内可见大量脂肪滴。淫羊藿组软骨层表面有轻微裂缝，未见明显软骨脱落，软骨细胞排列较为整齐，但存在细胞聚集；川芎组软骨层薄且脱落明显，淫羊藿组和川芎组均存在失染现象，但淫羊藿组较川芎组软骨层厚，软骨细胞簇集程度较弱。说明淫羊藿、川芎及其配伍、阳性药物均对 OA 有治疗作用。淫羊藿 - 川芎配伍治疗效果与阳性药物相似，较单用更佳。

3.3MMP-13结果

3.3.1 OA 模型大鼠血浆中 MMP-13 水平经时评价 如表 1 所示，与正常组比较，假手术组 MMP-13 水平无明显差异，模型组造模 2 周后 MMP-13 水平显著升高（ P ＜ 0.05 、 0.01 ）。

3.3.2 单药及配伍对 OA 大鼠血浆中 MMP-13 水平的影响 各给药组于第 12 、 14 周的 MMP-13 水平较同时期模型组有显著性差异（ P ＜ 0.01 ）。淫羊藿 - 川芎组和阳性药物组血浆中 MMP-13 水平显著低于淫羊藿组和川芎组（ P ＜ 0.05 、 0.01 ）。说明阳性药物组、淫羊藿、川芎与淫羊藿 - 川芎均有治疗 OA 的作用，且配伍使用效果较好。

3.4 代谢组学分析

3.4.1 OA 大鼠模型代谢组学分析 各组大鼠 0 ～ 8 周血浆代谢动态反映了 OA 的发展，如图 3 所示，第 0 、 2 周， 3 组血浆样本代谢无明显变化，第 4 、 6 、 8 周， 3 组样本均呈现分离趋势，假手术组的趋势更接近模型组。正常组大鼠未进行任何手术， 6 周的代谢变化只反映了大鼠正常衰老的情况；而假手术组模仿 OA 模型组大鼠实施了手术，体现了创伤对机体内源性成分的影响，其代谢变化由大鼠正常衰老和手术创伤共同作用产生。因此，在进行 OA 代谢组学研究时，以假手术组作为空白对照组更为合理。

如图 4 所示，正常组呈现一定的代谢变化趋势，可能是大鼠正常衰老的生理变化。假手术组较正常大鼠 0 ～ 2 周分离趋势较明显， 4 ～ 8 周变化趋于缓和，可能与手术创伤及机体自身修复有关。模型组造模前后样本分离明显，随着时间的推移，第 4 ～ 8 周呈现有规律的分离， OA 大鼠内源性代谢产物较正常组有明显变化。第 4 ～ 6 周变化最明显，第 6 ～ 8 周趋于稳定，说明第 6 周 OA 模型建立成功。

为了使组间差异最大化，消除组内变异对分类的负面影响，建立模型组和假手术组的 OPLS-DA 模型（图 5 ），模型组和假手术组能明显分开，两者在代谢产物方 面具有明显差异。在 200 次迭代中应用排列检验，说明没有过拟合。

3.4.2 给药组大鼠代谢组学分析 如图 6-A ～ H 所示，各给药组于第 10 ～ 14 周与模型组达到分离，并逐渐向正常组靠拢，在第 14 周趋势最为明显。说明各给药组均可调控其向正常状态转变，其中淫羊藿 - 川芎组和阳性药物组调控作用明显优于单药。各给药组 OPLS-DA 图（图 6-I 、 J ）显示，淫羊藿 - 川芎和阳性药物组疗效优于淫羊藿和川芎单药，进行了模型验证，显示 Q 2 截距小于 0 ，未过 拟合，模型可靠，初步阐明“补肾活血”法治疗 OA 的科学性。

3.5 血浆潜在生物标志物的筛选

结合 OPLS-DA 模型的 VIP 值及 P ＜ 0.05 筛选出 OA 模型中 69 种差异表达代谢物。通过建立第 8 周和第 14 周各给药组样本的 S-plot 图，从淫羊藿组、川芎组、淫羊藿 - 川芎组和阳性药物组中筛选出 16 种差异表达代谢物，鉴定出 2 个潜在的生物标志物（表 1 和图 7 ）。 Lyso P C （ 16 ∶ 0 ）水平下调，该生物标志物仅受淫羊藿 - 川芎的调控，提示淫羊藿与川芎联合用药可能增强对 OA 代谢物的调控。

3.6 血浆的代谢通路分析

潜在的生物标志物属于甘油磷脂，如LysoPC 、PC ，与磷脂代谢通路相关，代谢途径可能为OA 的发病机制及淫羊藿- 川芎对OA 大鼠作用机制提供依据。

4 讨论

本研究采用UPLC-Q-TOF-MS 技术经时评价OA 造模过程中血浆内源性代谢物和MMP-13 水平的动态变化，直观揭示了OA 大鼠模型在代谢层面的发病轨迹，有助于为OA 的早期诊断提供依据。研究了淫羊藿- 川芎配伍对OA 大鼠代谢的调控作用，给药后各组代谢轮廓变化及回调，其结果与组织病理学和MMP-13 评价的药效作用一致，证明了其配伍治疗OA 的合理性。初步获得2 个生物标志物，分别为LysoPC （16 ∶0 ）和PC （18 ∶0/20 ∶4 ），与磷脂代谢通路相关。

课题组前期考察淫羊藿- 川芎不同配比（1 ∶1 、1 ∶2 、2 ∶1 ）对OA 大鼠的药效作用，病理组织学（HE 、SO 染色），MMP-13 、IL-1β 检测结果表明淫羊藿、川芎以1:1 配伍使对OA 大鼠的药效作用最佳。故选择淫羊藿- 川芎（1 ∶1 ）研究其配伍机制。

经时评价OA 造模过程中，第6 周OA 模型建立，第8 周OA 病变进一步加重。一般药物干预主要对早期和中期OA 有效。为体现OA 早期诊断、早期治疗的理念，本研究选择造模后6 周进行药物干预研究。随着造模时间的延长，模型组大鼠血浆中MMP-13 表达水平逐步升高，提示可将MMP-13 的表达水平作为评价OA 模型建立的一个评价指标。给药后，各给药组血浆中MMP-13 水平随给药时间的延长而逐步降低，提示MMP-13 水平可用于疗效评价。

本研究初步获得 2 个生物标志物，分别为LysoPC （16 ：0 ）和PC （18 ∶0/20 ∶4 ），LysoPCs 和PCs 是脂质双分子层的组成部分，可以调节跨膜信号传导[18] 。OA 与磷脂的结构和浓度的变化密切相关。PCs 的缺乏，特别是含有长链脂肪酸的不饱和PCs ，可能会增加摩擦并导致关节软骨损伤[19] 。PCs 可通过磷脂酶A2 （PLA2 ）水解生成LysoPCs ，而OA 软骨细胞中PLA2 的高表达可促进IL-1β 、TNF 等炎性因子的表达[20-21] 。Zhang 等[22] 研究发现，LysoPCs 与PCs 的比率可用于预测OA ，与正常人相比，OA 患者血浆中LysoPCs 水平较高，PCs 水平较低。氧化应激反应也可导致OA ，当活性氧存在时，可以促使PCs 转化为LysoPCs[23] 。OA 软骨细胞、滑膜成纤维细胞和脂肪细胞均可产生造成机械、生化刺激的活性氧和一氧化氮。中性粒细胞可以通过释放蛋白水解酶混合物以结合活性氧，导致软骨组织破坏[24-26] 。

本研究从代谢组学层面初步揭示淫羊藿 - 川芎对OA 大鼠的作用机制，证明其配伍合理性，动态监测OA 进展和淫羊藿- 川芎治疗情况的结果，为探索中医药治疗OA 代谢组学的潜在疗效和机制提供了更有效的途径。

利益冲突所有作者均声明不存在利益冲突

参考文献（略）

来 源：吴恩慧，章建华，张 琳，李清林，王 玲，尹 华．基于UPLC-Q-TOF-MS代谢组学技术研究淫羊藿-川芎对骨关节炎大鼠的调控作用 [J]. 中草药, 2023, 54(22): 7445-7453.

【福利时刻】科研服务（点击查看）：、、、、、。微信：1278317307。

【福利时刻】毕业季找工作点击链接，多家企业30W+年薪工作招聘任你选（赶紧点击进入注册，上传简历吧~ ） 。