基于生物信息学与孟德尔随机化的活络消灵丸治疗骨关节炎分子机制解析

《Chinese Journal of Analytical Chemistry》：Decoding the molecular mechanism of Huoluoxiaoling pill in the treatment of osteoarthritis: a study based on bioinformatics and Mendelian randomization

本研究针对骨关节炎治疗中药物机制不明、靶点复杂等问题，通过整合网络药理学、孟德尔随机化与单细胞转录组分析，系统阐释了活络消灵丸通过“豆甾醇-TNF信号通路-CXCL10-巨噬细胞”轴调控免疫微环境的治疗新机制，为中药复方精准治疗骨关节炎提供理论依据。

　　

骨关节炎（Osteoarthritis, OA）作为一种常见的退行性关节疾病，严重影响患者生活质量，其病理特征包括关节软骨破坏、滑膜炎症及骨赘形成。当前治疗方法如非甾体抗炎药易引发胃肠道不良反应，而手术干预成本高昂且适用人群有限。传统中药活络消灵丸（Huoluoxiaoling Pill, HLXLP）在临床实践中展现出良好疗效，但其多成分、多靶点的特性使得作用机制尚未明确。

为系统解析HLXLP治疗OA的分子机制，研究团队整合生物信息学与孟德尔随机化（Mendelian Randomization, MR）技术，通过网络药理学筛选药物活性成分与靶点，结合差异表达分析、加权基因共表达网络分析（WGCNA）及单细胞转录组测序（scRNA-seq）等技术，从多层次揭示HLXLP的核心作用通路。该研究发表于《Chinese Journal of Analytical Chemistry》。

研究主要采用以下关键技术方法：

1. 1.
   通过网络药理学数据库（TCMSP、HERB、ETCM）筛选HLXLP活性成分及靶点；
2. 2.
   基于GEO数据库的GSE55235数据集进行差异表达基因分析；
3. 3.
   利用WGCNA识别OA相关基因模块，并通过MCODE插件提取关键簇基因；
4. 4.
   结合免疫浸润分析与scRNA-seq（GSE176308数据集）明确细胞类型与靶点分布；
5. 5.
   采用孟德尔随机化验证基因与OA的因果关系；
6. 6.
   通过分子对接与分子动力学模拟验证核心成分与靶点的结合稳定性。

3.1 潜在靶点获取

从HLXLP中鉴定出245个药物靶点，构建“药物-成分-靶点”网络，涵盖344个节点和1418条边。

3.2 OA差异表达基因分析

从GSE55235数据集筛选出716个差异表达基因（DEGs），其中上调基因382个、下调基因334个。基因集富集分析（GSEA）显示，TNF信号通路、细胞因子-细胞因子受体相互作用等通路显著富集。WGCNA分析进一步识别出与OA高度相关的turquoise模块（相关性0.95），其与DEGs交集获得297个疾病核心基因。

3.3 关键基因筛选与富集分析

药物靶点与疾病基因交集得到18个共同基因，通过PPI网络和MCODE分析提取23个关键基因。GO富集分析表明这些基因主要参与单核细胞分化等过程，KEGG通路富集提示其与IL-17信号通路、TNF信号通路密切相关。

3.4 免疫浸润与ROC分析

免疫浸润结果显示CXCL10与M1型巨噬细胞呈强正相关。ROC曲线分析表明CXCL8、CXCL10、CCL2等11个基因对OA具有较高诊断价值（AUC>0.5）。

3.5 单细胞转录组分析

scRNA-seq鉴定出3种细胞类型（软骨细胞、巨噬细胞、成纤维细胞）和8个细胞亚群，通路活性评分显示HLXLP主要作用于软骨细胞和巨噬细胞。

3.6 孟德尔随机化验证

MR分析证实CXCL10（OR_IVW=1.002）和CCL2（OR_IVW=1.014）为OA风险因子，CXCL8（OR_IVW=0.997）为保护因子，反向MR未发现反向因果关系。

3.8 分子对接与动力学模拟

核心基因（CXCL10、CXCL8、CCL2）与HLXLP活性成分（豆甾醇、槲皮素等）的分子对接显示豆甾醇与CXCL10结合最稳定（Vina评分=-9.6 kcal/mol）。分子动力学模拟表明CXCL10-豆甾醇复合物RMSD值稳定在1.9 ?，氢键数量持续为1，证实结合构象高度稳定。

讨论与结论

本研究首次通过多组学整合策略揭示HLXLP治疗OA的核心机制为“豆甾醇-TNF信号通路-CXCL10-巨噬细胞”轴。CXCL10通过趋化M1型巨噬细胞加剧炎症，而豆甾醇通过稳定结合CXCL10阻断其与CXCR3受体互作，从而抑制巨噬细胞活化。TNF信号通路作为上游调控因子，进一步强化了CXCL10的炎症促进作用。分子动力学模拟从原子层面验证了豆甾醇-CXCL10结合的可行性，为中药成分的靶向治疗提供结构生物学依据。

该研究不仅阐明了HLXLP的多靶点作用网络，还通过MR技术规避了传统观察性研究的混杂偏倚，为中药复方机制研究提供了新范式。未来可通过纳米载体递送豆甾醇等活性成分，结合智能康复设备实现OA的精准防治。