Article Info

骨性关节炎（osteoarthritis，OA）是由于关节先天性异常、关节畸形等因素引起关节软骨退化损伤[1]，主要通过口服药物与手术缓解症状[2]，但治疗的相关不良反应多[3]。既往研究[4]证实，中药在治疗OA、改善生存质量、明显延长无病生存期等方面发挥积极作用。骨关节炎属中医学 “痹证” 或 “瘘证” 范畴[5]，正气不足是其发病内因，湿、热、寒、风是其发病外因[6]，治疗以补肝肾、盈气血为主。

丹芎活血片是蔡立民教授根据多年临床辨证经验提出的成方，其采用活血化瘀、固正扶阳的方法，正切中病机，研究显示其抗凝效果确切，能通过促进局部的血液循环改善关节功能，显著提高患者生存质量，但具体作用机制仍需进一步研究[7]，而采用单一的机制研究不能体现丹芎活血片治疗OA 的整体性与多向性。网络药理学将药物作用靶点和疾病靶点结合，建立疾病、靶点、药物及通路之间的相互作用网络，从而全面系统地分析药物的作用机制[8]。因此本研究对丹芎活血片进行网络药理学研究，旨在为后续基础与临床研究提供分子机制基础。

1 材料与方法

1.1 筛选丹芎活血片的成分与靶标

检索TCMID 数据库（http://tcmspw.com/tcmsp.php）和SymMap 数据库（https://www.symmap.org/），收集既往文献中当归[9]6 个成分、丹参[10]6 个成分、川芎[11]6 个成分、苏木[12]4 个成分、土鳖虫[13]5 个成分、血竭[14]6个成分、泽兰[15]5 个成分、路路通[16]5 个成分、粉防己[17]1 个成分、桑枝[18]6 个成分。以口服生物利用度（oral bioavailability，OB）≥30%，类药性（drug-likeness，DL）≥0.18 为条件，得到成分并查询对应的靶点。

1.2 筛选OA 的靶点

在Gene Cards（https://www.genecards.org/）、OMIM（https://omim.org/）、DisGeNET（https://www.disgenet.org/）等数据库检索 “Osteoarthritis”，获得疾病靶点。

1.3 “中药-成分-靶点” 网络的构建

用Venny 对药物与疾病靶点取交集获得共同靶点，并在Cytoscape 3.7.2 构建 “中药-成分-靶点” 网络。采用Cytoscape 分析网络。

1.4 蛋白质-蛋白质交互网络的构建

在STRING（https://string-db.org/）构建蛋白质-蛋白质交互（protein-protein interaction，PPI）网络，获得核心网络模块。

1.5 基因本体与京都基因与基因组百科全书途径富集分析

借助DAVID 数据库（https://DAVID.ncifcrf.gov/）进行基因本体（gene ontology，GO）、京都基因与基因组百科全书途径富集（Kyoto encyclopedia of genes and genomes pathway enrichment，KEGG）分析，用于预测药物治疗OA 的生物过程。

1.6 “成分-靶点-通路” 网络的构建

将成分、靶点和信号通路导入微生信网站（http://www.bioinformatics.com.cn）构建网络图，筛选排名前5位的成分进行分子对接验证。

1.7 有效成分-靶点分子对接

利用LEDOCK 分子对接软件，将有效成分与核心靶点进行分子对接验证。

2 结果

2.1 活性成分获取

筛选得到丹芎活血片的活性化合物50 个（表1），对应的靶点1128 个。

表1 丹芎活血片的成分信息表

2.2 疾病靶点

筛选得到丹芎活血片的疾病靶点774 个。用Venny 取药病靶点交集，得到50 个共同靶点，见图1。

图1 药物与OA 靶点交集的韦恩图

2.3 “药物-化合物-靶点” 网络图

该网络图度值排名前5 的有效成分为槲皮素，山柰酚，姜黄素，阿魏酸，藁本内酯（图2），这可能是治疗OA 的活性成分。

图2 “药物-化合物-靶点” 网络

2.4 PPI 结果

在阈值≥0.7 中，Cytoscape 3.7.2 筛选出STAT3、血管内皮生长因子A（vascular endothelial growth factor A，VEGFA）、白细胞介素-6（interleukin-6，IL-6）、肿瘤坏死因子（tumor necrosis factor，TNF）、JUN 有密切联系；在0.9 的高度阈值中，筛选显示STAT3、VEGFA、IL-6、TNF、JUN、MMP9 是核心成分（图3），以此作为分子对接的靶点。

图3 PPI 网络图

2.5 GO 和KEGG

本研究中，GO 结果显示，分子功能富集于聚合酶Ⅱ启动子转录阳性、细胞质、蛋白质结合等（P<0.05）（图4，封四）。KEGG 结果显示，在癌症、乙型肝炎、核因子κB （nuclear factor κB，NF -κB）、PI3K -Akt 和JAK-STAT 信号通路中富集（图5）。

图4 GO分析图（见内文第5页）

图5 KEGG 分析图

2.6 “成分-靶点-通路” 网络图

将成分、靶点、前10 条通路使用Cytoscape 3.7.2构建 “成分-靶点-通路” 网络图。前5 个成分为阿魏酸、芥子酸、异阿魏酸、姜黄素、山柰酚（图6），以此作为分子对接的成分。

图6 “成分-靶点-通路” 网络图

2.7 分子对接

以6 个核心靶点与5 个有效成分进行分子对接（图7），结果显示，STAT3、VEGFA、IL-6、TNF、JUN、MMP9 的结合能均小于-5 kJ/mol，说明有较好的结合活性。

图7 药物和靶点的分子对接过程图

3 讨论

本研究通过网络分析，得到的槲皮素、山柰酚、姜黄素、阿魏酸已有文献证明对OA 有疗效，并通过PPI获得STAT3、MMP、IL-6、TNF、JUN。其中IL-6 严重破坏关节软骨组织平衡，而TNF-α 可以阻断软骨胶原的合成，导致软骨基质降解和破坏[19]。槲皮素可降低IL-6、IL-1β 和TNF-α 的分泌水平，并在蛋白水平下

调Nlrp3、caspase-1、proIL-1β 和IL-1β 的表达[20]，表现其抗氧化性和抗炎活性，从而治疗OA[21]。山柰酚可抑制bFGF-FGFR3-RSK2 信号抑制类风湿性关节炎成纤维样滑膜细胞的增殖和迁移，抑制活性T 细胞介质的炎症细胞因子（如IL-17、IL-21 和TNF-α）的释放[22]，而Tanideh 等[23]在利用山柰酚治疗前交叉韧带横切术（anterior cruciate ligament transection，ACLT）诱导的大鼠的实验中，发现大鼠体内的MMP-13、MMP-3、TNF-α、IL-1 明显降低。姜黄素同样可以充当抗炎药，抑制软骨细胞中NF-κB 介导的IL-1β/TNF-α 分解代谢信号通路[24]，由此可见，丹芎活血片中的活性成分能从多方面抑制关节病变的过程，从而发挥治疗OA 的作用。

KEGG 结果显示，调节PI3K/AKT 信号通路可以治疗OA。OA 患者中miR-146a-5P 表达增加，促进软骨细胞凋亡和炎症反应[25]。在OA 大鼠模型中，山柰酚抑制了miR-146a 和Decorin 在软骨组织中的表达，从而激活PI3K/Akt/mTOR 信号通路，与KEGG 的结果一致[26]。NF-κB 信号通路也与治疗OA 相关，槲皮素通过抑制TNF-α 和NF-κB 之间的正反馈调节，从而起到治疗OA 的作用[27]。为了客观展示药物治疗OA的分子机制，本研究结合了药物与疾病的映射靶点和信号通路的结果，绘制了科学假说图（图8，封四），为全面深入地了解作用机制奠定了基础。

图8 丹芎活血片治疗OA 的分子机制科学假说图（见内文第8 页）

本研究借助网络药理学技术，预测槲皮素、山柰酚、姜黄素、阿魏酸通过靶向PI3K/Akt、NF-κB 信号通路中的TNF、NF-κB、AKT 等分子起作用，从分子水平方向揭示了丹芎活血片治疗OA 的潜在作用机制，为后期药物研发及疾病治疗提供一定理论依据。

[参考文献]

[1]屈红林，刘瑞莲.老年退行性骨关节病的相关研究进展[J].中国老年学杂志，2014，34（6）：1720-1722.

[2]Martel-Pelletier J，Barr AJ，Cicuttini FM，et al.Osteoarthritis[J].Nat Rev Dis Primers，2016，2：16072.

[3]Hermann W，Lambova S，Muller-Ladner U.Current Treatment Options for Osteoarthritis[J].Curr Rheumatol Rev，2018，14（2）：108-116.

[4]Hou PW，Fu PK，Hsu HC，et al.Traditional Chinese medicine in patients with osteoarthritis of the knee[J].J Tradit Complement Med，2015，5（4）：182-196.

[5]刘力生，王文，姚崇华.2009年基层版《中国高血压防治指南》（摘录）[J].柳州医学，2011，24（2）：104-117.

[6]刘淑刚，王金榜.现代中医对膝骨关节炎的认识[J].现代中西医结合杂志，2013，22（13）：1473-1475.

[7]黄祖辉，蔡立民，谭志超，等.活血灵片在关节置换术后预防深静脉血栓形成的临床疗效[J].内蒙古中医药，2020，39（9）：41-42.

[8]权琦，陈平，黄圆圆，等.基于网络药理学研究皂角刺治疗乳腺癌的作用机制[J].中国当代医药，2021，28（9）：4-8.

[9]王圣鑫，余爱明，闫向丽，等.基于顶空固相微萃取气相质谱联用技术分析当归川芎药对配伍前后挥发性组分变化规律[J].中华中医药学刊，2020，38（5）：118-124.

[10]王云龙，房岐，郑超.丹参化学成分、药理作用及质量控制研究进展[J].中国药业，2020，29（15）：6-10.

[11]李芊，吴效科.川芎化学成分及药理作用研究新进展[J].化学工程师，2020，34（1）：62-64.

[12]邓成杰，刘爽，徐晓云，等.苏木化学成分及药理作用的研究进展[J].中国现代中药，2020，22（5）：810-826.

[13]王立娜，王颖，朱明珠，等.土鳖虫的活性成分及药理研究进展[J].化工时刊，2017，31（6）：34-36.

[14]林忆龙，文亦磊，黄洪，等.龙血竭化学成分及药理作用研究进展[J].中国民族民间医药，2020，29（6）：50-55.

[15]李余钊，文琰章，杨新洲，等.紫茎泽兰醋酸乙酯部位化学成分研究[J].中草药，2020，51（4）：932-936.

[16]马志良，赛桑杰，多杰.柱前衍生HPLC/FLD-APCI/MS 法同时测定藏药悬钩木中6 种三萜酸类成分的含量[J].中国药房，2019，30（16）：2243-2247.

[17]康凯，张静娇，刘斯亮，等.固相萃取技术-HPLC 法同时测定滑膜炎片/颗粒中防己诺林碱和粉防己碱的含量[J].中国药物评价，2020，37（3）：196-200.

[18]冯馨慧，刘志明，王海英.桑枝活性成分的提取及其抑菌和抗氧化作用[J].江苏农业科学，2020，48（7）：217-221.

[19]Saha S，Kundu J，Verma RJ，et al.Albumin coated polymer nanoparticles loaded with plant extract derived quercetin for modulation of inflammation[J].Materialia，2020，9（3）：7-11.

[20]刘俊.槲皮素对兔骨关节炎模型关节液中IL-1β、TNF-α水平的影响[D].桂林：桂林医学院，2017.

[21]Qiu L，Luo Y，Chen X.Quercetin attenuates mitochondrial dysfunction and biogenesis via upregulated AMPK/SIRT1 signaling pathway in OA rats[J].Biomed Pharmacother ，2018，103：1585-1591.

[22]Estakhri F，Panjehshahin M R，Tanideh N，et al.The effect of kaempferol and apigenin on allogenic synovial membranederived stem cells therapy in knee osteoarthritic male rats[J].Knee，2020，27（3）：817-832.

[23]Schulze-Tanzil G，Mobasheri A，Sendzik J，et al.Effects of curcumin （diferuloylmethane） on nuclear factor κB signaling in interleukin-1β-stimulated chondrocytes[J].Ann N Y Acad Sci，2004，1030：578-586.

[24]Shakibaei M，John T，Schulze -Tanzil G，et al.Suppression of NF -κB activation by curcumin leads to inhibition of expression of cyclo-oxygenase-2 and matrix metalloproteinase -9 in human articular chondrocytes：Implications for the treatment of osteoarthritis[J].Biochem Pharmacol，2007，73（9）：1434-1445.

[25]Zhao G，Gu W.Effects of miR -146a -5p on chondrocyte interleukin-1β-induced inflammation and apoptosis involving thioredoxin interacting protein regulation[J].J Int Med Res，2020，48（11）：1220769102.

[26]Jiang R，Hao P，Yu G，et al.Kaempferol protects chondrogenic ATDC5 cells against inflammatory injury triggered by lipopolysaccharide through down-regulating miR-146a[J].Int Immunopharmacol，2019，69：373-381.

[27]茅华杰.槲皮素对骨性关节炎的作用及其机制研究[D].杭州：浙江大学，2019.