类风湿关节炎(RA)作为一种常见的关节疾病,其患病率在不同地区呈现出较大差异。据研究显示,全球范围内该病的总体发病率为 0.5%~1.0%
[1]。在我国大陆地区,其发病率为0.42%,且存在明显的性别差异,男性患病率远低于女性。同时,RA可发生于任何年龄段,并且患病率会随着年龄的增长而增加,发病率与年龄呈正相关
[2]。Safiri等
[3]研究了1990~2017年195个国家的RA情况,从性别、年龄、社会人口指数等多维度分析了患病率、发病率及伤残调整生命年的水平与变化趋势,发现年龄标准化的患病率和发病率呈上升态势,加拿大、巴拉圭、危地马拉等国尤为明显。目前,RA已成为一项重大的全球公共卫生挑战。部分患者病情进展迅速,在发病后的2年内,关节周围会出现不可逆性损伤,这使得该病具有较高的致残率和死亡率,给患者家庭乃至整个社会都带来了沉重的负担
[4, 5]。
串联质量标签(TMT)技术是当前应用最广泛的定量方法之一
[6]。该技术在多个研究领域具有重要应用价值,其可以用于筛选和验证候选药物靶点,快速识别不同表达的蛋白质,发现潜在的治疗靶点
[7];可以用于药物作用机制的研究,通过比较药物处理组和对照组的蛋白质表达差异,揭示药物的作用机理;可以帮助评估药物的安全性和药效,为药物研发提供重要依据
[9]。
中药具有整体调节、多靶点作用的独特优势,在RA的干预机制研究领域日益深入,为临床治疗提供了新思路。豨莶丸(XXW)作为一味经典的中药制剂,仅由单味豨莶草加工制成,其药用价值在历代医籍中均有记载,具有显著的祛风除湿、通痹止痛、健体强身之效
[9, 10]。本研究基于TMT定量蛋白质组学分析豨莶丸治疗RA的分子机制。将围绕“模型复制→分子机制解析→多学科交叉验证”的逻辑主线,通过体内实验,结合免疫组化、免疫荧光等技术,探讨豨莶丸治疗RA的信号通路及核心靶点,为临床应用豨莶丸防治RA提供科学依据,进一步促进豨莶丸及其活性成分的研发应用,是一项十分有意义的创新性尝试。
1 材料和方法
1.1 实验材料
1.1.1 实验动物
健康SD大鼠48只(雄性,180±20 g),由贵州中医药大学实验动物研究所提供《SYXK(黔)2021-0005》,经贵州中医药大学动物伦理委员会批准(伦理批号:20250523003)。
1.1.2 主要试剂与药物
甲酸(Thermo Fisher Scientific);BCA定量试剂盒;SDS-PAGE蛋白上样缓冲液(碧云天);Trypsin (Promega);二硫苏糖醇 (DTT, Sigma);TMT 6plex Isobaric Mass Tag Labeling kit (Thermo Fisher Scientific)。肿瘤坏死因子-α(TNF-α)、白细胞介素-6(IL-6)、白细胞介素10(IL-10)ELISA试剂盒购自深圳子科生物科技有限公司。豨莶丸(江西民济药业有限公司,9 g/丸)
1.1.3 主要仪器
1.2 方法
1.2.1 复制模型及给药
48只大鼠适应环境后随机分为6组(豨莶丸高、中、低剂量组、雷公藤多苷片组、正常对照组和模型对照组,8只/组),按照文献复制大鼠胶原诱导关节炎(CIA) 模型
[11, 12]。造模成功2周后开始药物干预,豨莶丸低、中、高剂量换算后灌胃给药分别为200、400、800 mg/kg
[13],雷公藤多苷片(LGTDGP)9 mg/kg,1次/d。正常对照组和模型对照组灌胃等体积 1% 羧甲基纤维素钠,连续灌胃给药3周。
1.2.2 踝关节肿胀度测量
给药第 0 天及第 21 天,用游标卡尺测量每组大鼠右侧踝关节周长,并称量大鼠体质量。踝关节的肿胀度:周长(cm)/体质量(g)的比值。
1.2.3 病理组织学分析
实验结束后,腹腔注射麻醉大鼠,腹主动脉取血。取大鼠踝关节浸入4%多聚甲醛中固定24 h,固定完成后,依次进行脱钙处理,待组织硬度适宜后进行石蜡包埋,制作切片。在光学显微镜下观察踝关节病理组织学变化。
1.2.4 血清相关细胞因子水平检测
腹主动脉取血,静止30 min,4 ℃离心机4000 r/min离心30 min,分离出上层血清,置-80 ℃冰箱保存备用。使用ELISA试剂盒检测IL-10、IL-6、TNF-α 细胞因子的含量。TNF-α的检测范围为10 pg/mL~320 pg/mL,IL-6的检测范围为5 pg/mL~160 pg/mL, IL-10检测范围为12.5 pg/mL~400 pg/mL。
1.2.5 蛋白定性定量分析
取适量样本溶解后,精确称取100 mg进行胰蛋白酶消化处理。按照TMT标记试剂盒说明书的操作步骤进行标记
[14, 15]。标记完成并混合后的肽段,采用 Agilent 1260 Infinity II HPLC系统进行分级。每份样品通过纳升流速 Easy nLC 系统实现分离,所用缓冲液配置如下:A液为含0.1%甲酸的水溶液,B 液为含 0.1% 甲酸的 80% 乙腈水溶液。分离梯度设置具体为:0~3 min,B 液占比维持 6%;3~45 min,B 液占比从 6% 线性提升至 28%;45~50 min,B 液占比从 28% 线性提升至 38%;55~60 min,B 液占比保持 100%。实验所用数据库为 Rattus_ norvegicus_10116。
1.2.6 核心靶点筛选
通过STRING(
https://string-db.org/)和Cytoscape(
https://cytoscape.org/)联合分析差异蛋白相互作用,基于STRING数据库检索差异蛋白间潜在的相互作用,并构建蛋白质相互作用网络,明确差异蛋白之间的相关性。利用Cytoscape输出核心靶点
[16]。
1.2.7 核心靶点ROC分析
从标准化表达矩阵中获取CCL5、STAT1、GZMB和IL7R基因的表达情况,绘制ROC曲线并进行分析,计算ROC曲线下面积(AUC),评估CCL5、STAT1、GZMB和IL7R基因用于RA诊断的价值。
1.2.8 免疫组化检测大鼠踝关节CCL5和STAT1表达
免疫组化检测大鼠踝关节CCL5和STAT1的表达。每个样品封闭20 min后分别加稀释好的一抗,4 ℃过夜孵育,PBS 冲洗,加入稀释的二抗。阳性细胞染色程度评分:未染色(阴性)为 0 分,淡黄色(弱阳性)为1分、棕黄色(阳性)为2分、棕褐色(强阳性)为3分。阳性细胞计分:无阳性细胞为0分,1%~10%阳性细胞为1分、11%~50%阳性细胞为2分,51%~80%阳性细胞为3分,>80%阳性细胞为4分。免疫组化染色结果为阳性细胞染色程度评分与阳性细胞计分的乘积。
1.2.9 免疫荧光检测大鼠踝关节CCL5和STAT1蛋白阳性表达
踝关节组织用4%多聚甲醛进行固定,将一抗CCL5和STAT1按照一定比例稀释(1∶200),4 ℃ 孵育过夜。二抗室温孵育50 min,显色、复染、封片等。选择合适的荧光通道观察 CCL5和STAT1的染色情况,随机选取5个不重叠的高倍视野进行拍照。采用 Image J 软件对所拍图像进行分析,计算 CCL5和STAT1 的荧光强度(积分光密度值,IOD)和平均光密度(AOD),以此量化分析CCL5和STAT1的表达水平。
1.2.10 核心靶点与免疫细胞的相关性分析
使用 R 语言(R4.3.0),基于 ggplot2、tidyverse 及 ggExtra 包,绘制CCL5、STAT1、GZMB和IL17R与免疫细胞的相关性可视化图。输入免疫细胞浸润结果文件、训练数据集风险文件,去除正常对照组样本
[17]。横坐标代表相关的免疫特征基因,纵坐标代表免疫细胞类型,模块中颜色深浅代表该免疫特征基因与免疫细胞间的相关性的强弱。
1.3 统计学分析
数据分析采用 SPSS 统计软件。符合正态分布的数据以均数±标准差呈现,多组间的比较采用单因素方差分析,当P<0.05时,认为差异具有统计学意义。
2 结果
2.1 大鼠踝关节肿胀度及相关细胞因子检测
表2结果显示,与正常对照组比较,模型组大鼠踝关节肿胀度增加(
P<0.01)。与模型组比较,雷公藤多苷片(LGTDGP)、豨莶丸高剂量组(XXW-H)、中剂量组(XXW-M)大鼠踝关节肿胀度降低(
P<0.01,
P<0.05);与正常对照组比较,模型组大鼠血清中TNF-
α、IL-6的含量明显升高(
P<0.01),IL-10含量显著降低(
P<0.01);与模型组比较,LGTDGP、XXW-H和XXW-M的TNF-α、IL-6的含量降低(
P<0.01,
P<0.05),IL-10含量增加(
P<0.05,
P<0.01)。
2.2 大鼠踝关节组织病理学分析
正常大鼠踝关节结构正常,无炎性细胞浸润,滑膜细胞排列整齐,软骨表面光滑,关节腔内未见渗出液(
图1);与正常对照组比较,模型组大鼠关节破坏明显,周围大量中性粒细胞浸润,滑膜增生肥厚,纤维组织增生,软骨及骨骼损伤(
图1)。LGTDGP、XXW-H和XXW-M病理改变明显(
图1),炎细胞减少,血管增生减轻,豨莶丸低剂量组病理改变不明显(
图1)。
2.3 差异表达蛋白及核心靶点分析
基于TMT鉴定结果,选择蛋白丰度改变倍数>1.2或<0.8,
P<0.05的蛋白为差异蛋白(
图2A)。与正常对照组比较,模型组鉴定出346个差异蛋白;与模型组比较,给予XXW-H干预后,鉴定出287个差异蛋白。模型组和XXW-H组共同交集出160个差异蛋白,其中98个差异蛋白高表达,62个差异蛋白低表达(
图2B)。通过对160个差异蛋白进行相互作用分析,发现趋化因子(CCL5),信号转导与转录激活因子1(STAT1),颗粒酶 B(GZMB)和白介素7受体(IL7R)为核心靶点(
图2C)。
2.4 通路分析
160个差异蛋白主要参与56条信号通路(
P≤ 0.01),输出前10条信号通路,其中与RA相关的通路主要为辅助性T 细胞17 (Th17) 细胞分化、破骨细胞分化、Th1和Th2细胞分化等(
表3)。
2.5 差异表达基因ROC分析
从标准化表达矩阵中获取CCL5、STAT1、GZMB和IL7R基因的表达情况,绘制ROC曲线并进行分析。CCL5、STAT1的ROC曲线下面积大于0.9;GZMB和IL7R曲线下面积大于0.8(
图3)。CCL5、STAT1蛋白的表达采用免疫组化和免疫荧光进一步验证分析。
2.6 免疫组化检测大鼠踝关节中 CCL5和STAT1蛋白表达
结果显示(
表4,图
4、
5),与空白对照组比较,模型组大鼠踝关节中CCL5、STAT1的水平升高(
P<0.01);与模型组比较,经XXW-H和LGTDGP干预后,CCL5和STAT1的表达降低(
P<0.01)。
2.7 免疫荧光检测踝关节中CCL5、STAT1蛋白表达
与空白对照组比较,模型组大鼠踝关节中CCL5、STAT1的表达水平升高(
P<0.01);与模型组比较,经XXW-H和LGTDGP干预后,CCL5、STAT1的表达降低(
P<0.01,
表5,图
6、
7)。
2.8 差异表达基因与免疫细胞浸润之间的关系
基于 R 语言绘制CCL5、STAT1、GZMB和IL7R与免疫细胞的热图,纵坐标代表免疫细胞,横坐标代表CCL5、STAT1、GZMB和IL7R基因,模块中颜色代表CCL5、STAT1、GZMB和IL7R基因与免疫细胞间的相关性,发现CCL5,STAT1、GZMB和IL7R基因与多种免疫细胞相关(
图8)。4个核心靶点(CCL5,STAT1、GZMB、IL7R)能同时与中央记忆 CD4 T 细胞、效应记忆 CD4 T 细胞、骨髓源性抑制细胞、调节性T细胞、效应记忆 CD8 T 细胞、活化的 CD8T 细胞、活化的 CD4 T 细胞呈正相关(
P<0.05)。
3 讨论
豨莶丸作为治疗RA的潜在中药,遵循中医药辨证施治的核心原则,充分彰显了中医药在 RA 治疗中整体调节、个体化干预的独特与显著优势。目前以豨莶草为成分开发的中成药豨莶通栓丸(胶囊)、豨莶风湿片(胶囊)、豨桐丸(豨莶草与臭梧桐配伍)、豨红通络口服液已经应用到临床上。豨莶草的活性成分(Phynova)对风湿性疼痛、关节和肌肉疼痛以及轻微运动损伤有益,已经获英国药品及保健品管理署(MHRA)批准发售,且现有研究表明豨莶草具有抗炎、缓解痹痛、免疫抑制等作用
[17,18]。本研究发现豨莶丸通过调节炎症因子发挥抗RA作用,TNF-
α主要涉及系统性炎症,具有广泛的生物学效应,由活化的单核细胞、巨噬细胞和 T 淋巴细胞产生,主要功能就是调节免疫细胞
[19]。目前以 TNF-
α为靶点的生物药在自身免疫疾病中发挥作用。TNF-
α 是重要的炎症介质,被认为是 RA 发生、发展的首要动员细胞因子,能启动并维持 RA 炎症反应,许多药物都以 TNF-
α为靶向目标进行研究;IL-6导致 T细胞的增殖,循环加重RA的炎症反应过程
[20];IL-10具有免疫抑制及调控作用,已有研究发现IL-10在RA的治疗发挥作用
[21]。本研究发现豨莶丸能降低大鼠血清中TNF-
α和IL-6的含量,升高IL-10的含量发挥抗炎作用,对RA大鼠具有保护作用。
在蛋白组学研究中,核心靶点的发现和研究具有重要意义,能更深入地理解细胞内蛋白质之间的相互作用和调控机制。本研究发现160个差异蛋白主要通过 Th17 细胞分化、破骨细胞分化、Th1和Th2细胞分化信号通路参与RA的调节,Th17细胞在RA的发生、发展中发挥着重要作用。Th17 细胞主要分泌白细胞介素- 17(IL-17)、白细胞介素- 22(IL-22)等促炎因子,参与 RA 的发生
[22]。通过核心靶点的筛选本研究发现豨莶丸主要调节CCL5、STAT1、GZMB和IL7R的表达治疗RA,其中CCL5、STAT1的ROC曲线下面积大于0.9,表明其在诊断 RA 方面具有较高价值。CCL5是典型的炎症源性趋化因子,具有介导受体极化和定向树突状细胞、嗜碱性粒细胞等作用。CCL5 通过影响细胞免疫、体液免疫以及炎症过程在自身免疫疾病发挥着重要作用,如RA、系统性血管炎、自身免疫性溶血性贫血等
[23, 24];IL7是主要的促炎细胞因子,由辅助性 T 细胞(Th17)及先天性免疫细胞等分泌,在感染、自身免疫病、炎症等疾病中发挥主要作用
[25];RA滑膜组织中 STAT1/STAT3 被激活,导致抗凋亡分子表达上调,从而抑制细胞凋亡过程
[26];GZMB 是一种丝氨酸蛋白酶,主要参与机体的病理过程。目前研究发现 GZMB 在RA中高表达,RA的发生进展与GZMB的表达水平密切相关
[27, 28]。免疫细胞浸润也发现4个核心靶点(CCL5、STAT1、GZMB、IL7R)能同时与中央记忆 CD4 T 细胞、效应记忆 CD4 T 细胞、骨髓源性抑制细胞、效应记忆 CD8 T 细胞、活化的 CD8T 细胞、活化的 CD4 T 细胞呈正相关,进一步促进RA的发生、发展
[29, 30]。
综上所述,豨莶丸通过降低CIA大鼠血清中IL-6和TNF-α含量、增加IL-10含量发挥抗炎作用,起到对CIA大鼠关节的保护作用。在 ROC曲线分析中,CCL5、STAT1 的AUC值大于0.9,GZMB、IL7R的AUC值大于0.8,上述 4 个核心靶点作为RA的诊断标志物具有潜在价值。通过免疫组化和免疫荧光进一步验证CCL5、STAT1 在CIA大鼠踝关节的表达情况。发现CCL5、STAT1 在CIA大鼠踝关节中高表达,经豨莶丸治疗后表达下降。本研究已从动物实验水平为豨莶丸RA提供了理论依据,后续将进一步从体外实验水平展开深入分析。通过体内与体外研究的结合,有望弥补当前豨莶丸治疗 RA 作用机制研究中的不足,为其临床应用提供更全面的理论支持。
贵州中医药大学国家与省级科技创新人才团队培育项目(贵中医TD合字[2022]004号)