黄芪桂枝五物汤治疗骨关节炎网络药理学及分子对接技术分析

doi:10.3969/j.issn.2095-4344.3129

中国组织工程研究 ›› 2021, Vol. 25 ›› Issue (14): 2224-2230.doi: 10.3969/j.issn.2095-4344.3129

• 组织构建相关数据分析 Date analysis of organization construction • 上一篇下一篇

黄芪桂枝五物汤治疗骨关节炎网络药理学及分子对接技术分析

韩杰1，章晓云1，2，陈日兰1，张驰1，廖子龙1，陈跃平1，陈锋1，廖建钊1

1广西中医药大学附属瑞康医院，广西壮族自治区南宁市 530011；2江西中医药大学，江西省南昌市 330004

收稿日期:2020-04-28 修回日期:2020-04-30 接受日期:2020-05-30 出版日期:2021-05-18 发布日期:2020-12-31
通讯作者: 陈日兰，主任医师，硕士生导师，广西中医药大学附属瑞康医院，广西壮族自治区南宁市 530011
作者简介:韩杰，男，1981年生，广西壮族自治区合浦县人，副主任医师，硕士生导师，主要从事中医药在骨伤疾病中的临床与基础研究。
基金资助:
国家自然科学基金(81960908)，项目负责人：陈日兰；广西中药壮瑶药制剂提升工程项目(GZZJ16-07)，项目负责人：陈日兰；广西中医药大学一流学科项目(2019XK026)，项目负责人：陈跃平；广西中医药大学一流学科项目(2019XK029)

Huangqi Guizhi Wuwu Decoction in the treatment of osteoarthritis: analysis based on network pharmacology and molecular docking technology

Han Jie1, Zhang Xiaoyun1, 2, Chen Rilan1, Zhang Chi1, Liao Zilong1, Chen Yueping1, Chen Feng1, Liao Jianzhao1

1Ruikang Hospital Affiliated to Guangxi University of Chinese Medicine, Nanning 530011, Guangxi Zhuang Autonomous Region, China; 2Jiangxi University of Traditional Chinese Medicine, Nanchang 330004, Jiangxi Province, China

Received:2020-04-28 Revised:2020-04-30 Accepted:2020-05-30 Online:2021-05-18 Published:2020-12-31
Contact: Chen Rilan, Chief physician, Master’s supervisor, Ruikang Hospital Affiliated to Guangxi University of Chinese Medicine, Nanning 530011, Guangxi Zhuang Autonomous Region, China
About author:Han Jie, Associate chief physician, Master’s supervisor, Ruikang Hospital Affiliated to Guangxi University of Chinese Medicine, Nanning 530011, Guangxi Zhuang Autonomous Region, China
Supported by:
the National Natural Science Foundation of China, No. 81960908 (to CRL); Guangxi Chinese Medicine/Zhuangyao Pharmaceutical Preparation Improvement Project, No. GZZJ16-07 (to CRL); First-Class Discipline Projects of Guangxi University of Chinese Medicine, No. 2019XK026 (to CYP) and 2019XK029

摘要/Abstract

摘要：

文题释义：
生物利用度：是反映所给药物进入人体循环的药量比例，它描述口服药物由胃肠道吸收，及经过肝脏而到达体循环血液中的药量占口服剂量的百分比，包括生物利用程度与生物利用速度。
分子对接：是通过受体的特征以及受体和药物分子之间的相互作用方式来进行药物设计的方法。主要研究分子间(如配体和受体)相互作用，并预测其结合模式和亲合力的一种理论模拟方法。近年来，分子对接方法已成为计算机辅助药物研究领域的一项重要技术。

背景：研究证明，黄芪桂枝五物汤可通过抑制血管新生减少软骨损伤，促进软骨修复，延缓病情发展，但其作用机制尚不明确。
目的：采用网络药理学方法及分子对接技术探讨黄芪桂枝五物汤治疗骨关节炎的分子机制。
方法：通过TCMSP数据库筛选黄芪桂枝五物汤的活性成分及对应靶点，同时检索OMIM、GeneCards、TTD数据库获取治疗骨关节炎的相关靶点，然后将药物预测靶点与疾病靶点进行交集得到黄芪桂枝五物汤治疗骨关节炎的靶点。接着，利用STRING数据库结合Cytoscape 软件构建黄芪桂枝五物汤治疗骨关节炎的“单味药-活性成分-作用靶点”网络和蛋白质相互作用网络，并使用 David 数据库对黄芪桂枝五物汤治疗骨关节炎的作用靶点进行GO功能富集分析和 KEGG 通路富集分析，最后运用PyMOL等软件对黄芪桂枝五物汤的关键活性成分与关键靶点进行分子对接验证。
结果与结论：①筛选得到63个黄芪桂枝五物汤的化合物和217个相关靶点，与治疗骨关节炎的靶点取交集后得到130个黄芪桂枝五物汤治疗骨关节炎的靶点，主要涉及AGE-RAGE信号通路、MAPK信号通路、核因子κB信号通路、肿瘤坏死因子信号通路、白细胞介素17信号通路等，可能通过对脂多糖的反应、对细菌来源分子的反应、对营养水平的反应、氧化应激反应等功能发挥治疗骨关节炎的作用；②蛋白互作分析提示JUN、AKT1、HSP90AA1、血管内皮生长因子A、白细胞介素6为蛋白互作网络中的核心靶点，共实现了5个关键活性成分与5个关键靶点的分子对接；③结果说明，黄芪桂枝五物汤通过多成分、多靶点、多通路调控的特点发挥治疗骨关节炎的作用。
https://orcid.org/0000-0003-4736-9131 (陈日兰)

中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

关键词: 骨关节炎, 网络药理学, 靶点, 分子对接, 黄芪, 桂枝, 蛋白, 因子

Abstract: BACKGROUND: Studies have shown that Huangqi Guizhi Wuwu Decoction can reduce cartilage damage, promote cartilage repair, and delay the disease development by inhibiting angiogenesis, but its mechanism of action is not yet clear.
OBJECTIVE: To investigate the molecular mechanism of Huangqi Guizhi Wuwu Decoction in the treatment of osteoarthritis based on network pharmacology and molecular docking technology.
METHODS: The active ingredients of Huangqi Guizhi Wuwu Decoction were screened by TCMSP database and corresponding targets were collected. The potential targets for treating osteoarthritis were obtained by OMIM, GeneCards, and TTD databases, and the intersected targets of Huangqi Guizhi Wuwu Decoction for treating osteoarthritis were obtained. The STRING database was used to construct the single drug-active ingredient-acting target network, followed by protein-protein interaction network analysis using Cytoscape software. David database was used for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Molecular docking of active ingredients and key targets was finally performed using PyMOL software.
RESULTS AND CONCLUSION: Sixty-three compounds of Huangqi Guizhi Wuwu Decoction and 217 related targets were screened, and 130 targets were obtained after the intersection, mainly involving AGE−RAGE signaling pathway, MAPK signaling pathway, nuclear factor-κB signaling pathway, tumor necrosis factor signaling pathway, and interleukin-17 signaling pathway. Treatment of osteoarthritis might be achieved by the response to lipopolysaccharide, molecules of bacterial origin, nutrient level, and oxidative stress. ESR1, PRKDC, HSPA8, EP300 and HSP90AA1 were the hub targets in the protein-protein interaction network. A total of five key active ingredients were docked with five key targets. To conclude, Huangqi Guizhi Wuwu Decoction plays a role in the treatment of osteoarthritis through multi-component, multi-target and multi-channel regulations.

Key words: osteoarthritis, network pharmacology, target, molecular docking, Radix Astragali, ramulus cinnamomi, protein, factor

中图分类号:

韩杰, 章晓云, 陈日兰, 张驰, 廖子龙, 陈跃平, 陈锋, 廖建钊. 黄芪桂枝五物汤治疗骨关节炎网络药理学及分子对接技术分析[J]. 中国组织工程研究, 2021, 25(14): 2224-2230.

Han Jie, Zhang Xiaoyun, Chen Rilan, Zhang Chi, Liao Zilong, Chen Yueping, Chen Feng, Liao Jianzhao . Huangqi Guizhi Wuwu Decoction in the treatment of osteoarthritis: analysis based on network pharmacology and molecular docking technology[J]. Chinese Journal of Tissue Engineering Research, 2021, 25(14): 2224-2230.

图/表（结果） 9

2.1 化合物筛选结果以口服生物利用度(OB)≥30%且药物相似性(DL)≥0.18作为筛选参数，共得到74个活性化学成分，黄芪的化合物有20个，桂枝的化合物有7个，白芍的化合物有13个，生姜的化合物有5个，大枣的化合物有29个。其中大枣、黄芪、白芍有1个交集化合物，大枣、生姜有1个交集化合物，大枣、桂枝、白芍、生姜有1个交集化合物，大枣、桂枝、白芍有1个交集化合物，大枣、黄芪有1个交集化合物，桂枝、白芍有1个交集化合物，黄芪、白芍有1个交集化合物，删除重复数据后剩余活性化学成分63个，活性成分的基本信息见表 1，由于筛选结果过多，因此仅列出共同作用的活性成分及每味药口服生物利用度(OB)值前5的活性成分。
2.2 预测黄芪桂枝五物汤治疗骨关节炎的靶点结果根据TCMSP数据库预测得到217个黄芪桂枝五物汤药效靶点。根据OMIM、GeneCards和TTD数据库得到 2 885个骨关节炎的治疗靶点。运用Venny平台将黄芪桂枝五物汤药效靶点与骨关节炎的治疗靶点取交集，得到130个黄芪桂枝五物汤治疗骨关节炎的靶点，如图1。

2.3 药物-成分-靶点-疾病网络结果运用Cytoscape软件构建药物-成分-靶点-疾病网络，黄芪、桂枝、白芍、生姜和大枣共5种药物的37种化合物通过调控130个靶点治疗骨关节炎，如图2。

该网络中度值排名前5位的化合物分别是槲皮黄素、山柰酚、芒柄花黄素、异鼠李素、β-谷甾醇，分别能与97，37，19，17，15个靶点连接，对治疗骨关节炎具有重要的意义，其基本信息见表2。

2.4 GO富集分析与KEGG富集分析结果 GO富集分析所得到的GO条目共2 343个(P < 0.05)，其中生物过程(BP)条目 2 144个、细胞成分(CC)条目70个、分子功能(MF)条目129个；KEGG通路富集分析共得到146条通路。根据P排序，选取P值最显著的前10个条目以气泡图的形式展现，如图3。分析结果显示，黄芪桂枝五物汤治疗骨关节炎的生物过程主要涉及对脂多糖的反应、对细菌来源分子的反应、对营养水平的反应、氧化应激反应等；细胞成分主要涉及膜筏、膜微区、膜区等；分子功能主要涉及核受体活性、转录因子活性、细胞因子活性等；信号通路主要涉及AGE?RAGE信号通路、MAPK信号通路、核因子κB信号通路、肿瘤坏死因子信号通路、白细胞介素17信号通路等，如图3。

2.5 蛋白互作网络与蛋白模块借助String数据库和Cytoscape构建PPI网络图，见图4。图中共涉及节点65个、边125条，其中节点代表蛋白，边则代表蛋白的相互作用关系。度值排名前5的蛋白质为Jun原癌基因(Jun Proto-Oncogene，JUN)、AKT丝氨酸/苏氨酸激酶1(AKT Serine/Threonine Kinase 1，AKT1)、热休克蛋白90AA1(Heat Shock Protein 90 Alpha Family Class A Member 1，HSP90AA1)、血管内皮生长因子A(Vascular Endothelial Growth Factor A，VEGFA)、白细胞介素6(Interleukin 6，IL6)。这些度值较大的蛋白质在整个网络中起着关键的作用，同时这些蛋白所对应的靶点在黄芪桂枝五物汤治疗骨关节炎中具有重要的地位，可能是黄芪桂枝五物汤治疗骨关节炎的关键靶点，其基本信息见表3。

2.6 分子对接结果一般认为配体与受体结合越稳定结合能越低，发生作用的可能性越大，以结合能≤?5.0 kJ/mol作为标准，对预测的关键活性成分和关键靶点进行分子对接验证。分子对接结果显示黄芪桂枝五物汤的关键活性成分与关键靶点对接的亲和力远小于-5.0 kJ/mol，由此表明预测的黄芪桂枝五物汤关键活性成分与治疗骨关节炎的关键作用靶点具有较好的结合活性，证明该研究的预测较为可靠，结果见表4。选取亲和力最低的山柰酚和JUN进行分子对接展示，见图5。

参考文献

[1] Li D, Li G, Chen Y, et al. Astragaloside IV protects ATDC5 cells from lipopolysaccharide-caused damage through regulating miR-203/MyD88.Pharm Biol. 2020;58(1):89-97.
[2] Zhao J, Huang H, Liang G, et al. Effects and safety of the combination of platelet-rich plasma (PRP) and hyaluronic acid (HA) in the treatment of knee osteoarthritis: a systematic review and meta-analysis. BMC Musculoskelet Disord. 2020; 21(1):224.
[3] Liu T. Super-aging and social security for the most elderly in China. Z Gerontol Geriatr. 2018;51(1): 105-112.
[4] Mol MF, Runhaar J, Bos PK, et al. Effectiveness of intramuscular gluteal glucocorticoid injection versus intra-articular glucocorticoid injection in knee osteoarthritis: design of a multicenter randomized, 24 weeks comparative parallel-group trial. BMC Musculoskelet Disord. 2020;21(1):225.
[5] 陈日兰,邓凯烽,韦星成,等.电针对膝关节骨性关节炎患者疼痛改善及关节功能影响的荟萃分析[J].中国组织工程研究,2020,24(21): 3438-3444.
[6] 邓凯烽,尚鑫阳,朱圣旺,等.穴位注射改善膝骨关节炎患者疼痛及关节功能的Meta分析[J].中国组织工程研究,2020,24(18):2879-2887.
[7] 王骥.中医药疗法治疗膝关节骨性关节炎研究进展[J].中国临床医生杂志,2015,43(3): 36-39.
[8] 王清.黄芪桂枝五物汤结合关节镜治疗膝骨性关节炎的疗效观察[D]. 昆明:云南中医学院, 2017.
[9] 赵乐,李艳彦,王永辉,等.黄芪桂枝五物汤对骨关节炎大鼠血管新生的作用[J].中国实验方剂学杂志,2019,25(3):87-93.
[10] Li J, Zhao P, Li Y, et al. Systems pharmacology-based dissection of mechanisms of Chinese medicinal formula Bufei Yishen as an effective treatment for chronic obstructive pulmonary disease. Sci Rep. 2015;5:15290.
[11] Ru J, Li P, Wang J, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13.
[12] 赵乐,李艳彦,王永辉,等.黄芪桂枝五物汤对阳虚寒凝型骨关节炎大鼠免疫相关细胞因子的影响[J].中国实验方剂学杂志,2017, 23(7):160-166.
[13] 汤晓晨.黄芪甲苷干预人膝骨关节炎退变关节软骨细胞的效应及机制研究[D].南京:南京中医药大学, 2015.
[14] 胡鹏飞.芍药苷在骨性关节炎中的作用及其机制研究[D].杭州:浙江大学,2017.
[15] Wei B, Zhang Y, Tang L, et al. Protective effects of quercetin against inflammation and oxidative stress in a rabbit model of knee osteoarthritis.Drug Dev Res. 2019;80(3):360-367.
[16] Permatasari DA, Karliana D, Iskandarsyah I, et al. Quercetin prevent proteoglycan destruction by inhibits matrix metalloproteinase-9, matrix metalloproteinase-13, a disintegrin and metalloproteinase with thrombospondin motifs-5 expressions on osteoarthritis model rats. J Adv Pharm Technol Res. 2019;10(1):2-8.
[17] Zhuang Z, Ye G, Huang B. Kaempferol Alleviates the Interleukin-1β-Induced Inflammation in Rat Osteoarthritis Chondrocytes via Suppression of NF-κB.Med Sci Monit. 2017;23:3925-3931.
[18] Cho IA, Kim TH, Lim H, et al. Formononetin Antagonizes the Interleukin-1β-Induced Catabolic Effects Through Suppressing Inflammation in Primary Rat Chondrocytes. Inflammation. 2019; 42(4):1426-1440.
[19] Zhou F, Mei J, Yuan K, et al. Isorhamnetin attenuates osteoarthritis by inhibiting osteoclastogenesis and protecting chondrocytes through modulating reactive oxygen species homeostasis. J Cell Mol Med. 2019;23(6): 4395-4407.
[20] Cheng BC, Fu XQ, Guo H, et al. The genus Rosa and arthritis: Overview on pharmacological perspectives. Pharmacol Res.2016;114:219-234.
[21] Millerand M, Berenbaum F, Jacques C. Danger signals and inflammaging in osteoarthritis.Clin Exp Rheumatol. 2019;120(5):48-56.
[22] Rosenberg JH, Rai V, Dilisio MF, et al. Damage-associated molecular patterns in the pathogenesis of osteoarthritis:potentially novel therapeutic targets.Mol Cell Biochem. 2017; 434(1-2):171-179.
[23] Wei Y, Jin Z, Zhang H, et al. The Transient Receptor Potential Channel, Vanilloid 5, Induces Chondrocyte Apoptosis via Ca2+ CaMKII-Dependent MAPK and Akt/ mTOR Pathways in a Rat Osteoarthritis Model. Cell Physiol Biochem. 2018;51(5):2309-2323.
[24] Ran J, Ma C, Xu K, et al. Schisandrin B ameliorated chondrocytes inflammation and osteoarthritis via suppression of NF-κB and MAPK signal pathways.Drug Des Devel Ther. 2018;12: 1195-1204.
[25] Zhou F, Mei J, Han X, et al. Kinsenoside attenuates osteoarthritis by repolarizing macrophages through inactivating NF-κB/MAPK signaling and protecting chondrocytes. Acta Pharm Sin B. 2019;9(5):973-985.
[26] Yan S, Ding H, Peng J, et al. Down-regulation of protease-activated receptor 2 ameliorated osteoarthritis in rats through regulation of MAPK/NF-κB signaling pathway in vivo and in vitro. Biosci Rep. 2020;40(3):BSR20192620.
[27] Chow YY, Chin KY. The Role of Inflammation in the Pathogenesis of Osteoarthritis. Mediators Inflamm. 2020;2020:8293921.
[28] Zhang G, Gu M, Xu Y, et al. A comprehensive analysis on the effects of 1,25(OH)2D3 on primary chondrocytes cultured from patients with osteoarthritis. Gene. 2020;730:144322.
[29] Gao N, Ding L, Pang J, et al. Metabonomic-Transcriptome Integration Analysis on Osteoarthritis and Rheumatoid Arthritis.Int J Genomics. 2020;2020:5925126.
[30] Gao X, Sun Y, Li X. Identification of key gene modules and transcription factors for human osteoarthritis by weighted gene co-expression network analysis.Exp Ther Med. 2019;18(4): 2479-2490.
[31] Qiong J, Xia Z, Jing L, et al. Synovial mesenchymal stem cells effectively alleviate osteoarthritis through promoting the proliferation and differentiation of meniscus chondrocytes.Eur Rev Med Pharmacol Sci. 2020;24(4):1645-1655.
[32] Xue H, Tu Y, Ma T, et al. Lactoferrin Inhibits IL-1β-Induced Chondrocyte Apoptosis Through AKT1-Induced CREB1 Activation.Cell Physiol Biochem. 2015;36(6):2456-2465.
[33] Zhao X, Wang T, Cai B, et al. MicroRNA-495 enhances chondrocyte apoptosis, senescence and promotes the progression of osteoarthritis by targeting AKT1.Am J Transl Res. 2019;11(4): 2232-2244.
[34] Siebelt M, Jahr H, Groen HC, et al. Hsp90 inhibition protects against biomechanically induced osteoarthritis in rats. Arthritis Rheum. 2013;65(8):2102-2112.
[35] Boehm AK, Seth M, Mayr KG, et al. Hsp90 mediates insulin-like growth factor 1 and interleukin-1beta signaling in an age-dependent manner in equine articular chondrocytes. Arthritis Rheum. 2007;56(7): 2335-2343.
[36] Guan M, Zhu Y, Liao B, et al.Low-intensity pulsed ultrasound inhibits VEGFA expression in chondrocytes and protects against cartilage degeneration in experimental osteoarthritis.FEBS Open Bio. 2020;10(3):434-443.
[37] Takano S, Uchida K, Shoji S, et al. Vascular Endothelial Growth Factor Is Regulated by the Canonical and Noncanonical Transforming Growth Factor-β Pathway in Synovial Fibroblasts Derived from Osteoarthritis Patients.Biomed Res Int. 2019;2019:6959056.
[38] Sakata S, Kunimatsu R, Tsuka Y, et al. High-Frequency Near-Infrared Diode Laser Irradiation Attenuates IL-1β-Induced Expression of Inflammatory Cytokines and Matrix Metalloproteinases in Human Primary Chondrocytes. J Clin Med. 2020;9(3):881.
[39] Li D, Li G, Chen Y, et al. Astragaloside IV protects ATDC5 cells from lipopolysaccharide-caused damage through regulating miR-203/MyD88. Pharm Biol. 2020;58(1):89-97.

引言

骨关节炎是一种退化性疾病，其特征在于关节软骨进行性破坏和软骨下骨代偿性增生，往往会导致疼痛、僵硬和整个关节活动度下降[1-2]。由于人口老龄化和肥胖病的流行，骨关节炎的发病率正在上升[3]，研究表明每年1 000人中有40.2例患有骨关节炎[4]，将给社会带来巨大的经济负担。然而，目前尚未完全阐明骨关节炎的发病机制，且缺乏有效延缓骨关节炎进展的治疗手段。目前临床上治疗骨关节炎的药物均为对症治疗，并不能有效阻止骨关节炎病情的发展[5-6]。随着现代医学对临床疾病的治疗模式的改变，中医药在骨关节炎治疗过程中的重要性已逐渐被大家所认可[7]。因此寻求一种安全有效的中药方剂在治疗骨关节炎的过程中具有非常重要的临床意义。
黄芪桂枝五物汤源于张仲景所著《金匮要略》中，具有益气温经、活血通痹的功效。临床已证实黄芪桂枝五物汤加减治疗骨关节炎具有很好的疗效[8]，同时现代药理学基础实验也证明了黄芪桂枝五物汤可通过抑制血管新生减少软骨损伤，促进软骨修复，延缓病情发展[9]。虽然临床与基础实验都已证实黄芪桂枝五物汤可有效改善骨关节炎，但其作用机制尚不明确。由于黄芪桂枝五物汤的多成分、多靶点作用机制较复杂，实验研究难以系统分析药效化合物成分与靶点蛋白的联系。因此作者通过网络药理学方法分析黄芪桂枝五物汤作用于骨关节炎的靶点，并采用生物信息学技术分析靶点的功能和相关信号通路，为黄芪桂枝五物汤治疗骨关节炎的作用机制提供理论基础。
中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

材料方法

1.1 设计构建药物-成分-靶点-疾病网络，进行GO 功能富集分析和 KEGG 通路富集分析。
1.2 时间及地点研究于2020年2至3月在广西中医药大学附属瑞康医院完成。
1.3 材料
1.3.1 化合物筛选与靶点预测以药物口服生物利用度(oral bioavailability，OB)≥30%，药物相似性(drug likeness，DL)≥0.18作为筛选参数[10]，在TCMSP数据库(http://tcmspw.com/tcmsp.php)检索黄芪、桂枝、白芍、生姜、大枣的化合物成分及相关靶点[11]。通过UniProt数据库(https://www.uniprot.org/)筛选人源靶点并获取基因名。
1.3.2 骨关节炎治疗靶点预测通过OMIM数据库(https://omim.org/)、GeneCards数据库(http://www.genecards/, Version4.12)和TTD数据库(https://db.idrblab.org/ttd/)检索骨关节炎的治疗靶点，综合上述3个数据库的靶点，去除重复值，得到骨关节炎的治疗靶点。
1.3.3 黄芪桂枝五物汤治疗骨关节炎的靶点预测运用Venny平台(版本2.1，http://bioinfogp.cnb.csic.es/tools/venny/)将黄芪桂枝五物汤的化合物调控的靶点与骨关节炎的治疗靶点取交集，得到黄芪桂枝五物汤治疗骨关节炎的靶点。
1.4 方法
1.4.1 构建药物-成分-靶点-疾病网络先分别建立每种单味药的活性成分、活性成分作用靶点、靶点与骨关节炎相互作用关系，并导入到 Cytoscape3.7.2 软件(http://www. cytoscape.org/)中，再利用该软件的“merge”工具将5种单味药的网络进行合并，构建黄芪桂枝五物汤治疗骨关节炎的药物-成分-靶点-疾病网络图。
1.4.2 GO富集分析与KEGG富集分析运用DAVID数据库(版本6.8，https://david.ncifcrf.gov/)做GO富集分析与KEGG富集分析，FDR < 0.05代表富集结果显著[8]。最后利用R语言绘制GO富集分析气泡图及KEGG富集分析气泡图。
1.4.3 构建蛋白互作网络通过STRING数据库(http://string-db.org/cgi/input.pl)分析黄芪桂枝五物汤治疗骨关节炎的靶点之间的相互作用关系[9]，限定研究物种为“Homo Sapiens”获得蛋白互作关系，设置连接评分(Combined score)大于0.98，再导出蛋白质相互作用关系的数据文件，选取文件中的 node1、node2 和 Combined score三列导入到
Cytoscape中，最后利用Cytoscape软件的“NetworkAnalyzer”工具，设置节点大小与颜色反映度值大小，边的粗细反映Combined score的高低，构建 PPI网络图。
1.4.4 化合物与靶点分子对接为了验证预测的关键活性成分和关键靶点之间的相互作用关系，对二者进行分子对接。先用ChemOffice软件制作药物关键化合物的3D结构，再从PDB数据库下载核心蛋白基因的3D结构，利用PyMOL软件对蛋白质进行去水、去磷酸根等操作，利用Auto Dock1.5.6软件将化合物及核心蛋白基因pdb格式转为pdbqt格式并寻找活性口袋，最后运行Vina进行对接。根据化合物与靶点的结合能来评估两者的结合活性，目前对于活性成分的靶点筛选尚无统一标准，根据文献报道，选取结合能≤−5.0 kJ/mol的活性成分作为黄芪桂枝五物汤治疗骨关节炎靶点的筛选依据[10]，以评价网络分析预测的可靠性。
1.5 主要观察指标 ①化合物筛选结果；②预测黄芪桂枝五物汤治疗骨关节炎的靶点结果；③GO富集分析与KEGG富集分析结果；④蛋白互作网络；⑤分子对接结果。

讨论

中医将骨关节炎归为“骨痹”的范畴，并认为其病机为人体素虚，风、寒、湿邪乘虚而入，闭阻经络，导致气血运行不畅。目前临床针对骨关节炎的治疗，主要以非类固醇抗炎药缓解患处的疼痛，无法延缓疾病的进展。为改善患肢功能，患者往往在骨关节炎终末期接受关节置换。与西医的治疗理念相比，中医讲究的是整体观念与辨证论治，其治疗骨关节炎具有独特的优势。在骨关节炎早期，治以益气补阳，能够扶正气以抵御寒邪，阻止骨关节炎的进展。黄芪桂枝五物汤益气温经、活血通痹，使营卫调和，气血畅行充足，寒邪外散，则骨痹可除。方中黄芪为君，益气补卫，其得桂枝，固表而不留邪，桂枝得黄芪，散邪而不伤正，益气而振奋卫阳，桂、芍相伍，共调营卫和表里，两药为臣。生姜辛温，疏散风邪，以助桂枝之力；大枣甘温，养血益气，以资黄芪、芍药之功；与生姜为伍，又能和营卫，调诸药，以为佐使。实验研究证实其能够通过调节白细胞介素4修复阳虚寒凝型骨关节炎模型大鼠软骨细胞[12]，方中主要成分黄芪甲苷和芍药苷可抑制炎症因子白细胞介素1β的表达，从而抑制软骨细胞凋亡，促进关节软骨修复[13-14]。但由于黄芪桂枝五物汤治疗骨关节炎的具体分子机制尚未明确，此次研究通过网络药理学及分子对接技术从微观角度进行阐述分析。
此次研究通过对黄芪桂枝五物汤进行初步筛选共预测出黄芪桂枝五物汤中217个活性成分靶点和2 885个骨关节炎的治疗靶点，利用Venny平台进行映射取交集后共获得130个黄芪桂枝五物汤治疗骨关节炎的潜在作用靶点。单味药-活性成分-作用靶点网络显示槲皮黄素、山柰酚、芒柄花黄素、异鼠李素、β-谷甾醇能匹配到较多的靶点，可能是黄芪桂枝五物汤治疗骨关节炎的关键有效成分。研究表明槲皮黄素通过上调超氧化物歧化酶和基质金属蛋白酶抑制剂1，下调基质金属蛋白酶13，从而减弱氧化应激反应并抑制软骨细胞外基质的降解以改善骨关节炎病情[15]，还可通过抑制肿瘤坏死因子α和基质金属蛋白酶9的表达来预防炎症和蛋白聚糖降解以减缓关节软骨退变[16]。山萘酚通过抑制核因子κB途径，减少白细胞介素1β的分泌，可以预防、停止或延迟骨关节炎病情发展[17]。芒柄花黄素可抑制白细胞介素1β诱导软骨细胞及关节软骨分解代谢，改善骨关节炎状况，可能是治疗和预防骨关节炎的潜在有效的候选药物[18]。异鼠李素可抑制破骨细胞活性和软骨细胞凋亡从而避免软骨损伤，从而达到治疗骨关节炎的效果[19]。β-谷甾醇可降低炎性细胞因子的产生，并降低氧化应激，进而抑制炎症过程，可用于临床骨关节炎的防治[20]。作者认为黄芪桂枝五物汤通过有效成分槲皮黄素、山柰酚、芒柄花黄素、异鼠李素、β-谷甾醇的作用，以减弱氧化应激反应，抑制炎症的发生，软骨细胞的凋亡，从而发挥治疗骨关节炎的作用。
GO富集结果显示黄芪桂枝五物汤治疗骨关节炎的130个潜在作用靶点基因富集生物学过程主要涉及对脂多糖的反应、对细菌来源分子的反应、对营养水平的反应、氧化应激反应等，均与炎症的形成与发展有着密切的联系。KEGG通路富集结果显示黄芪桂枝五物汤治疗骨关节炎的信号通路主要有AGE−RAGE信号通路、MAPK信号通路、核因子κB信号通路、肿瘤坏死因子信号通路、白细胞介素17信号通路等。研究表明AGE−RAGE可激活促炎因子的分泌，导致关节出现炎症反应，同时可充当炎症细胞因子来激活先天免疫细胞，进一步导致骨关节炎病情的发展[21-22]。MAPK信号通路参与软骨细胞凋亡过程[23]，诸多中药提取物可通过MAPK信号通路发挥作用，改善骨关节炎病情[24-25]。通过体内和体外MAPK/NF-κB信号通路的调控下调PAR2可有效延缓骨关节炎病情的发展[26]。CHOW等[27]研究表明骨关节炎发病时，关节炎症可刺激白细胞介素1β的释放，进而诱导MAPK、PI3K/AKT、核因子κB等信号通路信号传导途径的激活，从而导致更多炎症因子的释放，研究中还指出白细胞介素17诱导滑膜成纤维细胞和软骨细胞释放白细胞介素6、白细胞介素8和肿瘤坏死因子α，导致炎症和软骨分解。相关研究KEGG分析也表明肿瘤坏死因子信号途径、白细胞介素17信号途径、核因子κB信号途径是骨关节炎疾病发生发展过程中重要的信号通路[28-29]。作者认为黄芪桂枝五物汤主要通过影响炎症相关及软骨细胞凋亡相关通路，调控炎症因子和软骨细胞凋亡水平调控骨关节炎病情的发展和转归。
蛋白互作网络表明JUN、AKT1、HSP90AA1、血管内皮生长因子A和白细胞介素6可能是四妙散治疗GA的关键靶点。研究表明FoxO，MAPK和PI3K/Akt信号通路对骨关节炎的发病机制至关重要，JUN是这些通路激活的关键基因，其可激活相关通路传导调节软骨细胞[30-31]。AKT1对软骨细胞凋亡起着非常重要的作用，其水平的高低直接影响了骨关节炎病情的发展[32-33]。研究表明Hsp90在软骨细胞应激反应中起着关键作用，其参与软骨代谢的多种途径，抑制Hsp90可保护或改善软骨退变情况，阻止骨关节炎病情发展[34-35]。减轻白细胞介素1β诱导的小鼠原代软骨细胞中血管内皮生长因子A的表达可有效减缓骨关节炎的进展[36]，血管内皮生长因子与骨关节炎的严重程度和疼痛密切相关[37]。另有研究表明骨关节炎的严重程度还与炎性因子有关，炎性细胞因子白细胞介素6、白细胞介素8和肿瘤坏死因子α水平的高低均会影响骨关节炎病情的发展[38-39]。
综上所述，此次研究初步预测了黄芪桂枝五物汤治疗骨关节炎的化合物-靶点作用机制，分析了靶点的功能和相关信号通路，并完成了部分蛋白与化合物的分子对接，初步验证了部分研究结果。但是此次研究收集的药效靶点可能不够全面，相关数据库的更新进度可能滞后，预测得到的靶点未展开体外和体内实验验证，因此存在一定的局限性，关于黄芪桂枝五物汤治疗骨关节炎的作用机制验证将成为下一步研究的重点。
中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

黄芪桂枝五物汤治疗骨关节炎网络药理学及分子对接技术分析

Huangqi Guizhi Wuwu Decoction in the treatment of osteoarthritis: analysis based on network pharmacology and molecular docking technology

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 9

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	林清凡, 解一新, 陈婉清, 叶振忠, 陈幼芳. 人胎盘源间充质干细胞条件培养液可上调缺氧状态下BeWo细胞活力和紧密连接因子的表达[J]. 中国组织工程研究, 2021, 25(在线): 4970-4975.
[2]	彭智浩, 冯宗权, 邹勇根, 牛国庆, 吴峰. 活动平台单髁置换后下肢力线与外侧间室骨关节炎进展的关系[J]. 中国组织工程研究, 2021, 25(9): 1368-1374.
[3]	袁家威, 张海涛, 揭珂, 曹厚然, 曾意荣. 基于网络药理学研究桃红四物汤治疗假体周围感染的潜在靶点和机制[J]. 中国组织工程研究, 2021, 25(9): 1428-1433.
[4]	黄登承, 王志科, 曹学伟. 体外冲击波疗法治疗中老年膝骨关节炎短期疗效对比的荟萃分析[J]. 中国组织工程研究, 2021, 25(9): 1471-1476.
[5]	吴训, 孟娟红, 张建运, 王亮. 浓缩生长因子修复兔髁突全层软骨损伤[J]. 中国组织工程研究, 2021, 25(8): 1166-1171.
[6]	申晋波, 张林. 急性力竭运动导致大鼠跟腱微损伤的超微结构变化及机制[J]. 中国组织工程研究, 2021, 25(8): 1190-1195.
[7]	李静, 谢建山, 崔慧林, 曹锡梅, 杨艳萍, 李海荣. 不同毛色小鼠背部皮肤中甘油二酯激酶ζ和蛋白激酶CβⅡ的表达和定位[J]. 中国组织工程研究, 2021, 25(8): 1196-1200.
[8]	柴乐, 吕建兰, 胡劲涛, 胡华辉, 许庆军, 余进伟, 全仁夫. 诱导急性脊髓损伤模型大鼠炎症反应信号通路的变化[J]. 中国组织工程研究, 2021, 25(8): 1218-1223.
[9]	耿秋东, 葛海雅, 王和鸣, 李楠. 基于网络药理学探讨龟鹿二仙胶治疗骨关节炎的作用及机制[J]. 中国组织工程研究, 2021, 25(8): 1229-1236.
[10]	刘翔翔, 黄云梅, 陈文列, 林如辉, 卢小冬, 李钻芳, 许亚晔, 黄美雅, 李西海. 早期骨关节炎模型大鼠半月板白区细胞的超微结构变化[J]. 中国组织工程研究, 2021, 25(8): 1237-1242.
[11]	谭婧玉, 刘海文. 成骨细胞特异因子2基因家族的全基因组鉴定、分类及进化分析[J]. 中国组织工程研究, 2021, 25(8): 1243-1248.
[12]	李中峰, 陈明海, 凡一诺, 魏秋实, 何伟, 陈镇秋. 右归饮治疗激素性股骨头坏死作用机制的网络药理学分析[J]. 中国组织工程研究, 2021, 25(8): 1256-1263.
[13]	顾霞, 赵敏, 王平义, 李一梅, 李文华. 低氧诱导因子1α与低氧相关疾病信号通路的关系[J]. 中国组织工程研究, 2021, 25(8): 1284-1289.
[14]	姬志祥, 蓝常贡. 尿酸盐转运蛋白在痛风中的多态性和治疗相关性[J]. 中国组织工程研究, 2021, 25(8): 1290-1298.
[15]	刘聪, 刘肃. miR-17-5p调控低氧诱导因子1α介导脂肪细胞分化及血管生成的分子机制[J]. 中国组织工程研究, 2021, 25(7): 1069-1074.