基于网络药理学探讨龟鹿二仙胶治疗骨关节炎的作用及机制

doi:10.3969/j.issn.2095-4344.2991

中国组织工程研究 ›› 2021, Vol. 25 ›› Issue (8): 1229-1236.doi: 10.3969/j.issn.2095-4344.2991

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

基于网络药理学探讨龟鹿二仙胶治疗骨关节炎的作用及机制

耿秋东1，葛海雅1，王和鸣1，2，3，李楠1， 2

1福建中医药大学中医学院，福建省福州市 350122；2中医骨伤及运动康复教育部重点实验室，福建省福州市 350122；3福建中医药大学附属康复医院，福建省福州市 350003

收稿日期:2020-05-14 修回日期:2020-05-19 接受日期:2020-06-01 出版日期:2021-03-18 发布日期:2020-12-12
通讯作者: 李楠，博士，教授，硕士生导师，福建中医药大学中医学院，福建省福州市 350122；中医骨伤及运动康复教育部重点实验室，福建省福州市 350122
作者简介:耿秋东，男，1991年生，河南省太康县人，汉族，福建中医药大学在读硕士，主要从事中医药防治骨病方向的研究。
基金资助:
国家自然科学基金面上项目(81973880)；福建省自然科学基金面上项目(2019J01349)

Role and mechanism of Guilu Erxianjiao in treatment of osteoarthritis based on network pharmacology

Geng Qiudong1, Ge Haiya1, Wang Heming1, 2, 3, Li Nan1, 2

1College of Traditional Chinese Medicine, Fujian University of Traditional Chinese Medicine, Fuzhou 350122, Fujian Province, China; 2Key Laboratory of Orthopedics & Traumatology and Rehabilitation of Traditional Chinese Medicine, Ministry of Education, Fuzhou 350122, Fujian Province, China; 3Rehabilitation Hospital Affiliated to Fujian University of Traditional Chinese Medicine, Fuzhou 350003, Fujian Province, China

Received:2020-05-14 Revised:2020-05-19 Accepted:2020-06-01 Online:2021-03-18 Published:2020-12-12
Contact: Li Nan, MD, Professor, Master’s supervisor, College of Traditional Chinese Medicine, Fujian University of Traditional Chinese Medicine, Fuzhou 350122, Fujian Province, China; Key Laboratory of Orthopedics & Traumatology and Rehabilitation of Traditional Chinese Medicine, Ministry of Education, Fuzhou 350122, Fujian Province, China
About author:Geng Qiudong, Master candidate, College of Traditional Chinese Medicine, Fujian University of Traditional Chinese Medicine, Fuzhou 350122, Fujian Province, China
Supported by:
the National Natural Science Foundation of China, No. 81973880; the Natural Science Foundation of Fujian Province, No. 2019J01349

摘要/Abstract

摘要：

文题释义：
网络药理学：在系统生物学、多向药理学及计算生物学等多学科技术基础上，发展出来的研究药物机制的新方法。
骨关节炎：一种退行性病变，系由于增龄、肥胖、劳损、创伤、关节先天性异常和关节畸形等诸多因素引起的关节软骨退化损伤、关节边缘和软骨下骨反应性增生，又称骨关节病、退行性关节炎、老年性关节炎和肥大性关节炎等。

背景：临床研究表明龟鹿二仙胶治疗骨关节炎效果显著，但其药理机制尚不明确。
目的：文章运用网络药理学研究方法，分析龟鹿二仙胶治疗骨关节炎的潜在靶点、作用网络及潜在机制。
方法：文章借助TCMSP，BATMAN-TCM和TCMID数据库获取龟鹿二仙胶的化学成分及对应靶点，OMIM，Gene Cards和Drug Bank数据库检索骨关节炎的已知疾病靶点；运用Perl语言将药物活性成分靶点与疾病靶点相映射得到龟鹿二仙胶治疗骨关节炎的作用靶点，运用Cytoscape3.7.2软件构建“药物-化合物-靶点”作用网络，运用STRING软件构建蛋白相互作用(PPI)网络，对所得PPI网络进行拓扑分析，筛选发挥治疗作用的关键成分与关键靶点；运用R语言对关键靶点进行GO富集分析和KEGG通路富集分析获知其潜在作用机制，构建龟鹿二仙胶“通路-靶点”，获知主要作用通路与靶点之间的关系。
结果与结论：①文章预测得到龟鹿二仙胶治疗骨关节炎的有效成分69个，靶点201个，关键成分有槲皮素、山奈酚及β-谷甾醇等6个，关键靶点有白细胞介素6、丝裂原活化蛋白激酶1和含半胱氨酸的天冬氨酸蛋白水解酶3等20个；②GO和KEGG分析结果显示，龟鹿二仙胶通过以上关键靶点，影响丝裂原活化蛋白激酶信号通路、白细胞介素17信号通路和肿瘤坏死因子信号通路等13个主要通路，集中于调控炎症、软骨代谢、细胞凋亡与增殖等功能来发挥防治骨关节炎的作用；③结果证实，以上关键靶点和通路构建了龟鹿二仙胶防治骨关节炎的作用网络，其多成分-多靶点-多途径的治疗特点为龟鹿二仙胶进一步的基础研究以及临床应用提供了新的思路与线索。

https://orcid.org/0000-0001-6537-7464 (耿秋东)

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

关键词: 骨, 关节, 软骨, 网络药理学, 生物信息学, 靶点, 成分, 途径, 代谢, 细胞凋亡, 炎症, 中医药

Abstract: BACKGROUND: Clinical studies have shown that Guilu Erxianjiao has a remarkable effect on osteoarthritis, but its pharmacological mechanism is not yet clear.
OBJECTIVE: To analyze the potential targets, action networks and potential mechanisms of Guilu Erxianjiao for the treatment of osteoarthritis using network pharmacology.
METHODS: The chemical constituents and corresponding targets of Guilu Erxianjiao were obtained with the help of TCMSP, BATMAN-TCM and TCMID databases, and the known disease targets of osteoarthritis were searched in OMIM, GeneCards and Drug Bank databases. Perl language was used to map the targets of the active ingredients of the drug to the targets of the disease to obtain the targets of Guilu Erxianjiao for treating osteoarthritis. Cytoscape 3.7.2 software was used to build a “drug-compound-target” interaction network. STRING software was used to construct a protein-protein interaction network, followed by a topological analysis on the obtained protein-protein interaction network to screen the key components and key targets that play a therapeutic role. R language was used to perform gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Gnomes (KEGG) pathway enrichment analysis to understand its potential mechanism of action, construct the “pathway-target” of the Guilu Erxianjiao. The relationship between the main pathway and targets was then obtained.
RESULTS AND CONCLUSION: It is predicted that 69 active ingredients and 201 targets of Guilu Erxianjiao for treating osteoarthritis will be obtained. There are 6 key components such as quercetin, kaempferol and β-sitosterol, and 20 key targets such as interleukin-6, mitogen-activated protein kinase 1, and cysteine-containing aspartate proteolytic enzyme 3. The results of GO and KEGG analysis showed that, through the above key targets, the Guilu Erxianjiao affected 13 major pathways, including mitogen-activated protein kinase signaling pathway, interleukin-17 signaling pathway, and tumor necrosis factor signaling pathway, and played a role in preventing and treating osteoarthritis by the focus on regulating inflammation, cartilage metabolism, apoptosis and proliferation. In conclusion, the above key targets and pathways have constructed an action network for the prevention and treatment of osteoarthritis by Guilu Erxianjiao, and its multi-component-multi-target-multi-channel treatment features provide new ideas and clues for the further basic research and clinical application of Guilu Erxianjiao.

Key words: bone, joint, cartilage, network pharmacology, bioinformatics, target, ingredient, pathway, metabolism, cell apoptosis, inflammation, traditional Chinese medicine

中图分类号:

耿秋东, 葛海雅, 王和鸣, 李楠. 基于网络药理学探讨龟鹿二仙胶治疗骨关节炎的作用及机制[J]. 中国组织工程研究, 2021, 25(8): 1229-1236.

Geng Qiudong, Ge Haiya, Wang Heming, Li Nan. Role and mechanism of Guilu Erxianjiao in treatment of osteoarthritis based on network pharmacology[J]. Chinese Journal of Tissue Engineering Research, 2021, 25(8): 1229-1236.

图/表（结果） 13

2.1 活性成分筛选及靶点筛选共得到活性化合物95个，其中龟板胶19个，鹿角胶13个，人参22个，枸杞45个；共有对应疾病靶点3 493个，见表2。利用Perl语言将上述数据进行综合处理后得到2 163个药物对应靶点。

2.2 疾病靶点筛选结果将所得结果进行综合删除重复项后，得到骨关节炎相关靶点2 612个。
2.3 网络构建与结果分析得到201个共有靶点，对应龟鹿二仙胶活性成分69个。利用Cytoscape 3.7.2软件构建“中药-活性成分-疾病靶点”网络图见图1。靶点数排名前8位的成分为槲皮素(quercetin)、碳酸钙(calcium Carbonate)、山奈酚(kaempferol)、甘氨酸(glycine)、大豆黄素(glycitein)、β-谷甾醇(β-sitosterol)、丙氨酸(alanine)和原阿片碱(fumarine)，见表3，这些可能为龟鹿二仙胶治疗骨关节炎的有效成分。将映射所得的201个靶点利用String构建PPI网络，其产生相互作用的网络有201个节点，3 543条边，所得PPI网络导入Cytoscape 3.7.2软件，使用CytoNCA插件进行拓扑分析，得到20个靶点的网络，见图2。对这20个靶点进行分析，明确此20个靶点的基本功能信息，见表4。将这20个靶点所对应的成分进行网络构建得到图3所示网络，根据degree值大小筛选前8位的成分为槲皮素、山奈酚、人参皂苷Rh2、β-谷甾醇、大豆黄素、原阿片碱、碳酸钙和花生四烯酸盐，结合前面所得到的前8位成分取交集后得到槲皮素、山奈酚、β-谷甾醇、大豆黄素、原阿片碱和碳酸钙6个成分，提示这6个成分可能为龟鹿二仙胶治疗骨关节炎的关键成分。

2.4 GO富集分析将2.3中得到的20个靶点利用R语言软件借助Bioconductor数据包library(org.Hs.eg.db)对这20个靶点进行ID转换，然后运用R语言进行GO分析得到1 550个条目，其中生物学过程占1 512个，分子生物学功能占21个，细胞组分占17个，显示龟鹿二仙胶治疗骨关节炎主要与对氧化应激的反应、对脂多糖的反应、对细菌起源分子的反应、对活性氧的反应及细胞蛋白质定位建立的正调控等生物学过程以及生长因子受体结合、丝裂原活化蛋白激酶活性、受体配体活性、细胞因子受体结合、蛋白质丝氨酸/苏氨酸激酶活性等分子生物学功能密切相关，其过程涉及细胞内囊泡、膜筏、膜微区、内质网腔及核染色质等细胞组分。其主要生物学过程、分子生物学功能、细胞组分参与的基因见表5-7。

2.5 KEGG富集分析及“通路-关键靶点”关系网络构建KEGG分析结果见图3，共得到138条通路，将这138条通路信息利用Perl语言进行处理后，将所得数据导入Cytoscape 3.7.2软件构建“通路-关键靶点”关系网络见图4，根据degree值并排除宽泛通路后筛选出13条主要通路，具体主要通路和其涉及基因见表8。

2.6 龟鹿二仙胶治疗骨关节炎的化学成分及作用机制的文献验证通过检索共纳入文献17篇[8-11, 19-31]。分析所得文献，明晰龟鹿二仙胶的功用主治和方中单味药的功效，以及此方治疗骨关节炎的中医认识及作用机制，见图5。文献报道龟鹿二仙胶全方、拆方及其单药中主要化学成分槲皮素、山奈酚和甘氨酸等对骨关节炎均具有治疗作用[8-11, 19-26]；大豆黄素可能通过对成骨细胞的增殖作用，抑制破骨性骨吸收，防治骨关节炎的发生和进展[27-28]；人参皂苷能通过增高Bcl-2/Bax的比值、提高软骨细胞中Ⅱ型胶原基因的表达、抑制基质金属蛋白酶13的表达来改善骨关节炎患者症状[29-31]。其治疗骨关节炎的机制主要与对炎症介质、基质金属蛋白酶表达量和Ⅱ型胶原基因的调控相关，说明龟鹿二仙胶治疗骨关节炎可能与炎症、软骨代谢、细胞凋亡与增殖等通路有关。

参考文献

[1] SILVERWOOD V, BLAGOJEVIC-BUCKNALL M, JINKS C, et al. Current evidence on risk factors for knee osteoarthritis in older adults: a systematic review and meta-analysis. Osteoarthritis Cartilage. 2015; 23(4):507-515.
[2] LANE NE, BRANDT K, HAWKER G, et al. OARSI-FDA initiative: defining the disease state of osteoarthritis. Osteoarthritis Cartilage. 2011; 19(5):478-482.
[3] 中医骨伤科临床诊疗指南·膝痹病(膝骨关节炎)[J].康复学报,2019,29(3):1-7.
[4] 林嘉辉,陈炳艺,李楠.龟鹿二仙胶治疗骨关节炎的研究概况[J].中华中医药杂志, 2017,32(12):5509-5512.
[5] 姜灵.龟鹿二仙胶汤治疗肝肾亏虚型膝关节骨性关节炎的临床疗效研究[D].大连:大连医科大学,2016.
[6] 刘俊锋.龟鹿二仙胶汤对膝骨关节炎早期干预的临床研究[D].成都:成都中医药大学, 2013.
[7] 张农.龟鹿二仙胶治疗老年增生性骨关节炎的临床观察[J].中医临床研究,2012,4(13): 79-80.
[8] 余添赐.龟鹿二仙胶不同提取部位对体外培养豚鼠关节软骨细胞活性的影响[D].福州:福建中医药大学,2014.
[9] 李楠,杜江,黄远鹏,等.龟鹿二仙胶及其拆方对SD大鼠膝骨关节炎的影响[J].中华中医药杂志,2013,28(6):1677-1680.
[10] 杜江,黄远鹏,李奕修,等.龟鹿二仙胶及其拆方对去势骨关节炎大鼠血清雌二醇浓度及膝关节软骨细胞Ⅱ型胶原表达的影响[J].中医正骨,2013,25(3):11-20.
[11] 杜江,黄远鹏,李奕修,等.龟鹿二仙胶及其拆方对去势骨关节炎大鼠关节软骨中Bax和Bcl-2的影响[J].中国中医骨伤科杂志, 2013,21(3):4-7.
[12] 刘艾林,杜冠华.网络药理学:药物发现的新思想[J].药学学报,2010,45(12):1472-1477.
[13] 许海玉,黄璐琦,卢鹏,等.基于体内ADME过程和网络药理学的中药现代研究思路[J].中国中药杂志,2012,37(2):142-145.
[14] Xu X, Zhang W, Huang C, et al. A novel chemometric method for the prediction of human oral bioavailability. Int J Mol Sci. 2012; 13(6):6964-6982.
[15] 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(1):15290.
[16] 侯磊,李晓宇,黄娜娜,等.四逆散干预焦虑症作用机制的网络药理学分析[J].中草药, 2019,50(21):5154-5161.
[17] SZKLARCZYK D, MORRIS J H, COOK H, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017; 45(D1):D362-D368.
[18] 张新雨,黄娜娜,孙凯滨,等.柴胡龙骨牡蛎汤治疗原发性高血压作用网络与机制预测[J].中草药,2019,50(21):5162-5169.
[19] KANZAKI N, SAITO K, MAEDA A, et al. Effect of a dietary supplement containing glucosamine hydrochloride, chondroitin sulfate and quercetin glycosides on symptomatic knee osteoarthritis: a randomized, double-blind, placebo-controlled study. J Sci Food Agric. 2012;92(4):862-869.
[20] LEYVA-LOPEZ N, GUTIERREZ-GRIJALVA E P, AMBRIZ-PEREZ D L, et al. Flavonoids as cytokine modulators: a possible therapy for inflammation-related diseases. Int J Mol Sci. 2016;17(6):921.
[21] QIU L, LUO Y, CHEN X. Quercetin attenuates mitochondrial dysfunction and biogenesis via upregulated AMPK/SIRT1 signaling pathway in OA rats. Biomed Pharmacother. 2018;103:1585-1591.
[22] BRITTI D, CRUPI R, IMPELLIZZERI D, et al. A novel composite formulation of palmitoylethanolamide and quercetin decreases inflammation and relieves pain in inflammatory and osteoarthritic pain models. BMC Vet Res. 2017;13(1): 229.
[23] HALEAGRAHARA N, HODGSON K, MIRANDA-HERNANDEZ S, et al. Flavonoid quercetin-methotrexate combination inhibits inflammatory mediators and matrix metalloproteinase expression, providing protection to joints in collagen-induced arthritis. Inflammopharmacology. 2018;26(5): 1219-1232.
[24] HU Y, GUI Z, ZHOU Y, et al. Quercetin alleviates rat osteoarthritis by inhibiting inflammation and apoptosis of chondrocytes, modulating synovial macrophages polarization to M2 macrophages. Free Radic Biol Med. 2019; 145:146-160.
[25] ZHUANG Z, YE G, HUANG B. Kaempferol Alleviates the Interleukin-1beta-Induced Inflammation in Rat Osteoarthritis Chondrocytes via Suppression of NF-kappaB. Med Sci Monit. 2017;23:3925-3931.
[26] de PAZ-LUGO P, LUPIANEZ JA, MELENDEZ-HEVIA E. High glycine concentration increases collagen synthesis by articular chondrocytes in vitro: acute glycine deficiency could be an important cause of osteoarthritis. Amino Acids. 2018,50(10):1357-1365.
[27] 王志平,禹博,刘明皓,等.大豆苷元对人成骨细胞经典的雌激素受体信号通路的调节作用[J].武警后勤学院学报(医学版), 2018,27(3):200-203.
[28] 王健,陶海荣.雌激素及相关化合物干预骨关节炎作用机制:关节保护与骨软骨细胞的修复[J].中国组织工程研究,2014, 18(33):5372-5376.
[29] 黎立,李晓旭,司裕,等.人参皂苷Rb2保护关节软骨细胞对抗氧化应激的作用机制研究[J].中国矫形外科杂志,2018,26(9): 845-849.
[30] WANG W, ZENG L, WANG ZM, et al. Ginsenoside Rb1 inhibits matrix metalloproteinase 13 through down-regulating Notch signaling pathway in osteoarthritis. Exp Biol Med (Maywood). 2015;240(12):1614-1621.
[31] 段超.人参皂苷Rg1防治兔膝关节骨性关节炎的实验研究[D].武汉:湖北中医药大学, 2016.
[32] 黄威,尹宗生.炎症与骨关节炎软骨退变[J].中国矫形外科杂志,2019,27(5):448-452.
[33] SOHN DH, SOKOLOVE J, SHARPE O, et al. Plasma proteins present in osteoarthritic synovial fluid can stimulate cytokine production via Toll-like receptor 4. Arthritis Res Ther. 2012;14(1):R7.
[34] CHOCKALINGAM PS, VARADARAJAN U, SHELDON R, et al. Involvement of protein kinase Czeta in interleukin-1beta induction of ADAMTS-4 and type 2 nitric oxide synthase via NF-kappaB signaling in primary human osteoarthritic chondrocytes. Arthritis Rheum. 2007;56(12):4074-4083.
[35] 陈喜德,魏波,杨柳菁,等.不同病期骨关节炎大鼠关节软骨中caspase-3与细胞凋亡的相关性[J].广东医学,2018,39(2):181-186.
[36] 冯顶丽,卓丽丹,芦笛,等.微小RNA调节间充质干细胞软骨分化机制的研究进展[J].国际口腔医学杂志,2018,45(6):640-645.
[37] WU Y H, LIU W, ZHANG L, et al. Effects of microRNA-24 targeting C-myc on apoptosis, proliferation, and cytokine expressions in chondrocytes of rats with osteoarthritis via MAPK signaling pathway. J Cell Biochem. 2018; 119(10):7944-7958.
[38] 邵进,张岩,王治,等.低氧诱导因子-1α参与骨发育及骨代谢调控的研究进展[J].中国骨质疏松杂志,2015,21(3):349-355.
[39] 田大川,李海乐,肖大伟,等.联合应用软骨再生支架与突变型HIF-1α修饰BMSCs分泌的外泌体对晚期软骨缺损修复的促进作用[J].吉林大学学报(医学版),2018,44(2): 216-222.
[40] 马笃军,彭力平,余阗,等.牛膝总皂苷对骨关节炎模型兔软骨修复及低氧诱导因子1信号通路的影响[J].中国组织工程研究,2019,23(27):4332-4337.
[41] 孙志江.自噬在晚期糖基化终产物诱导的软骨细胞凋亡中的作用[D].兰州:兰州大学, 2018.
[42] SORCI G, RIUZZI F, GIAMBANCO I, et al. RAGE in tissue homeostasis, repair and regeneration. Biochimica et Biophysica Acta (BBA). Mol Cell Res. 2013;1833(1):101-109.
[43] 徐小磊.TLR4/Akt/FoxO1轴在白藜芦醇抗肥胖相关骨性关节炎中的作用研究[D].沈阳:中国医科大学,2019.

引言

骨关节炎是最常见的关节疾病，是导致残疾的主要原因。个体遗传易感性、衰老、肥胖和关节排列不良等几个危险因素与骨关节炎密切相关[1-2]。在中医学中并无骨关节炎的病名，但据其临床表现可归属于中医学骨痹范畴，《中医骨伤科临床诊疗指南·膝痹病(膝骨关节炎)》将其分为肝肾亏虚、气血虚弱、寒湿痹阻、湿热痹阻及气滞血瘀共5个证型，其治疗原则多为滋补肝肾、补气养血、温经散寒、清热除湿和活血化瘀等[3]。目前骨关节炎的诊断仅是基于临床症状和影像学标准，病因至今未明，为骨伤科治疗难题之一，因此，亟需探明有效的骨关节炎防治策略。因中医药治疗骨关节炎的不良反应较小，临床疗效肯定，近年来研究者对中药复方治疗骨关节炎作用机制的探究也越来越多。龟鹿二仙胶最早见于《医便》，由龟板胶、鹿角胶、人参和枸杞组成，诸药合用，共奏补益肝肾、填精益髓及助阳益气之功[4]。临床研究表明龟鹿二仙胶方治疗骨关节炎效果显著，能明显改善患者的临床症状[5-7]。此外，龟鹿二仙胶防治骨关节炎的作用机制也取得了一定的进展，研究表明龟鹿二仙胶能抑制软骨细胞凋亡、延缓软骨外基质降解及减少炎症因子表达[8-11]，但由于龟鹿二仙胶复方成分较复杂，其作用于骨关节炎的具体机制仍缺乏有力证据证明。
网络药理学集合了数据挖掘、靶点预测和高通量数据整合等多种研究方法，它把中药成分靶点间的相互作用和联系系统的展现出来，还能整体分析中药复方多成分的作用机制，其系统性和整体性与中医药辨证论治、整体观念的原则具有很大程度的一致性，深入推动了中药复方药效分子机制的研究[12-13]。文章基于网络药理学，通过数据挖掘与分析，筛选和系统预测龟鹿二仙胶治疗骨关节炎的潜在作用靶点和信号通路，为龟鹿二仙胶临床精准定位提供药理学依据，为进一步的基础研究以及临床应用提供了新的思路与线索。中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

材料方法

1.1 数据库、软件及分析平台作者于2020年2月至5月运用生物信息学和数据挖掘技术检索以下10个数据库的文献，见表1。

1.2 分析方法
1.2.1 龟鹿二仙胶活性成分的获取及靶点筛选利用TCMSP，BATMAN-TCM和TCMID数据库收集龟鹿二仙胶中龟板胶、鹿角胶、人参及枸杞的化学成分，基于口服生物利用度(OB)≥30%和类药性(DL)≥0.18进行筛选[14-15]，同时查阅文献对未纳入数据库但已报道具生物活性和药理作用的也一并纳入，最终确定其活性成分。通过 TCMSP和Drug Bank数据库查询活性成分相对应的靶标蛋白，并利用 Unitprot 数据库查询靶标蛋白对应的基因名称，剔除非人类的靶点蛋白[16]，最终整理得到龟鹿二仙胶的活性成分及其对应靶点。
1.2.2 骨关节炎靶点筛选预测利用GeneCards和OMIM数据库，以 “osteoarthritis”为关键词，对骨关节炎相关靶点进行检索及筛选，将所有筛选结果合并去重后得到骨关节炎的相关靶标信息。
1.2.3 网络构建与分析利用Perl语言将所得药物活性成分靶点与骨关节炎疾病靶点相映射，得到龟鹿二仙胶治疗骨关节炎的预测靶点。将交集得到的靶点导入String数据库[17]，构建蛋白质-蛋白质相互作用网络，并用Cytoscape 3.7.2软件中的CytoNCA插件对所得的蛋白质相互作用网络数据进行分析，借助度中心性(DC)和介度中心性(BC)指标进行筛选，首先筛选大于2倍中位数的度中心性值的节点，筛选完成后根据初筛结果继续筛选大于每个指标所对应的中位数的节点[18]，得到龟鹿二仙胶治疗骨关节炎的关键靶点。
1.2.4 基因本体论(gene ontology，GO)富集分析首先，使用R语言(版本3.6.2)软件将龟鹿二仙胶治疗骨关节炎的靶点进行ID转换，为下一步做基础；其次，在R语言中调用生物类的程序包，如clusterProfiler，bioconductor等对关键靶点进行GO富集分析(P < 0.05)，得到GO富集分析结果，其中包括龟鹿二仙胶治疗骨关节炎靶点的生物学过程(biological process，BP)、分子生物学功能(molecular function，MF)和细胞组分(cellular component，CC)等具体数据，通过对基因的富集分析预测龟鹿二仙胶治疗骨关节炎的可能机制。
1.2.5 京都基因和基因组百科全书(Kyoto Encyclopedia of Genes and Gnomes，KEGG)富集分析及KEGG关系网络的构建使用R语言对关键靶点进行KEGG富集分析(P < 0.05)，得到KEGG富集分析具体通路信息，将得到的结果利用Perl语言进行处理后，运用Cytoscape 3.7.2构建“通路-关键靶点”关系网络，分析龟鹿二仙胶治疗骨关节炎的“通路-关键靶点”之间的关系。
1.2.6 龟鹿二仙胶治疗骨关节炎的化学成分及作用机制的文献验证通过使用检索词结合不同数据库的检索方式，检索以下3个数据库：知网、万方和PubMed数据库。中文检索词包括：“骨关节炎”“龟鹿二仙胶”“中医”“槲皮素”“山奈酚”等；英文检索词包括：“Osteoarthritis，Guilu Erxianjiao，quercetin，kaempferol”等，挖掘分析龟鹿二仙胶治疗骨关节炎的相关成分与机制，为网络药理学所分析的结果提供文献验证。
1.3 主要观察指标关键靶点的GO和KEGG富集结果；关键靶点在龟鹿二仙胶治疗骨关节炎相关通路中的分布情况。

讨论

中医学认为骨关节炎属“痹证”和“骨痹”范畴。骨痹的病因多考虑为外感风邪，遇寒或受热，加之湿邪偏盛，诸邪裹结，合而为痹。而素体骨髓亏虚，肌腠不能抵御邪气，外邪乘隙入侵，是导致骨痹发生的内在因素。中医治疗骨痹讲究辨证论治，旨在通过中草药的配伍，以达到通络止痛、补虚泻实之功效。龟鹿二仙胶方中的龟板胶入肝、肾、肺三经，有滋阴潜阳和益肾健骨之用；鹿角胶归经肝肾，能温补肝肾和益精养血；人参入脾、肺、心经，有大补元气之功；枸杞入肝、肾经，能治肝肾亏虚、腰膝酸软之症。四药合用，能对肾元虚损和精血不足所致的筋骨形体失养从根本上进行治疗。
网络药理学能够从中药成分、靶点与疾病间相互作用的整体性和系统性出发进行研究，特别适宜研究中药“多成分-多靶点-多途径”的作用关系，探明中药的药效物质基础及具体作用机制。文章基于网络药理学方法分析龟鹿二仙胶治疗骨关节炎“化学成分-靶点-通路”之间的关系，探讨了其潜在作用机制。文章通过富集分析得出了龟鹿二仙胶治疗骨关节炎的主要成分有69个，其中槲皮素、山奈酚和β-谷甾醇等成分可能为龟鹿二仙胶治疗骨关节炎的关键成分；主要治疗靶点有201个，其中的白细胞介素6、丝裂原活化蛋白激酶1和含半胱氨酸的天冬氨酸蛋白水解酶3等20个靶点为核心靶点。
在骨关节炎发展过程中，炎症、软骨代谢、细胞凋亡与增殖几种因素是密切关联的，炎症可以导致软骨组织细胞外基质的降解，也可以通过线粒体依赖途径、内质网应激途径以及死亡受体介导的细胞信号转导途径引起软骨细胞凋亡[32]。在龟鹿二仙胶治疗骨关节炎所针对的核心靶点中，白细胞介素1β由关节中的软骨细胞、成骨细胞、脂肪细胞和滑膜细胞、浸润滑膜的免疫细胞等产生[33]，能强烈诱导蛋白水解酶(如基质金属蛋白酶和聚集蛋白聚糖酶)的表达和释放，并抑制细胞外基质的表达，包括Ⅱ型胶原蛋白和聚集蛋白聚糖[34]，它还与其他细胞因子如白细胞介素6、白细胞介素8、单核细胞趋化蛋白1和趋化因子2起协同作用，进一步增强炎症反应。此外，丝氨酸/苏氨酸蛋白激酶能调控新陈代谢、细胞增殖、细胞存活、生长和血管生成等过程，这些功能是丝氨酸、苏氨酸单独或者联合通过磷酸化一系列下游底物来介导的，如miR-543能够提升蛋白激酶和B淋巴细胞瘤-2基因表达水平，使软骨细胞的增殖得到提升，有效对抗白细胞介素1β诱导引发的软骨细胞损伤[35-36]。分析以上结果可知，龟鹿二仙胶治疗骨关节炎的作用与调控炎症、软骨代谢、细胞凋亡与增殖等因素密切相关。
KEGG富集分析结果显示，龟鹿二仙胶防治骨关节炎可能主要通过调控丝裂原活化蛋白激酶、晚期糖基化终末产物-糖基化终末产物受体、叉头框转录因子O和低氧诱导因子1等13条信号通路发挥作用。其中丝裂原活化蛋白激酶信号通路与骨关节炎软骨损伤具有重要联系，它的激活能调控炎症反应、软骨代谢、细胞凋亡与增殖等多种生物学反应。研究证实，上调miR-24下调C-myc可以通过失活丝裂原活化蛋白激酶信号通路来抑制软骨细胞凋亡并促进软骨细胞增殖，进而阻止骨关节炎的发生和发展[37]。
正常的关节软骨是一个低氧组织，作为对低氧反应的主要调节因子，低氧诱导因子1α能显著调节骨的发育和重建，在软骨细胞生长和分化过程中具有关键作用，软骨缺乏低氧诱导因子1α能导致大量软骨细胞死亡[38-39]。马笃军等[40]用牛膝总皂苷干预兔骨关节炎模型，发现低氧诱导因子1信号通路在实验兔骨关节炎模型软骨修复过程中发挥了较大作用。有研究发现晚期糖基化终末产物如戊糖素在关节软骨中的含量随着年龄增长而增加，而戊糖素含量与关节软骨的胶原网络脆性具有较高联系性；此外，软骨与另一些晚期糖基化终末产物如苏阿糖等发生糖化作用后，胶原网络的刚性会发生增长，且增长幅度与晚期糖基化终末产物含量呈高相关性；胶原网络两种特性的增加则会提升其易感性，关节软骨更易发生机械损伤[41]。糖基化终末产物受体与晚期糖基化终末产物结合可以引发一系列病变，机制十分复杂。研究表明，糖基化终末产物受体与衰老、慢性炎症及细胞应激有密切联系，在细胞稳态及损伤后修复等过程中都有参与[42]。
目前，中医药防治骨关节炎的临床研究和动物实验很多，但是明确针对晚期糖基化终末产物−糖基化终末产物受体信号通路的研究仍然较少，晚期糖基化终末产物作为衰老导致骨关节炎发生的一个重要因素，中医药通过调控此信号通路防治骨关节炎的具体作用机制仍待进一步研究。此外，肥胖因素已成为骨关节炎患者大量增加的一个重要原因之一，肥胖导致的慢性轻度炎症与骨关节炎进程密切相关。
Toll样受体是一类非特异性免疫受体，其家族成员Toll样受体4在肥胖相关的骨关节炎进程中有至关重要的作用，阻断Toll样受体4信号通路即可起到治疗肥胖相关骨关节炎的作用[43]，但是目前对中医药防治肥胖相关骨关节炎的研究相对较少，因此可以进一步开展相关研究。
文章运用网络药理学的研究技术预测了龟鹿二仙胶治疗骨关节炎的潜在活性成分、关键靶点及其作用途径网络，为更深入地探讨其作用机制以及更广泛地应用于临床提供了一定依据。另一方面以文献挖掘得出的信息进行比对，得出对此次龟鹿二仙胶治疗骨关节炎的作用靶点、作用网络具有一定的预测性和准确性，但是在预测结果中仍有部分活性成分、靶点及作用网络的临床及基础验证不足，在后期可对这些活性成分、靶点及作用网络进行进一步验证，为临床的应用服务提供更充足的证据。中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

基于网络药理学探讨龟鹿二仙胶治疗骨关节炎的作用及机制

Role and mechanism of Guilu Erxianjiao in treatment of osteoarthritis based on network pharmacology

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 13

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	徐峰, 康辉, 魏坦军, 席金涛. 椎弓根螺钉不同固定方法治疗胸腰椎骨折的生物力学分析[J]. 中国组织工程研究, 2021, 25(9): 1313-1317.
[2]	姜勇, 罗翼, 丁永利, 周勇, 闵理, 唐凡, 张闻力, 段宏, 屠重棋. 骶骨精准切除影响骨盆稳定性的冯米斯应力特征及临床验证[J]. 中国组织工程研究, 2021, 25(9): 1318-1323.
[3]	张同同, 王中华, 文杰, 宋玉鑫, 刘林. 3D打印模型在颈椎肿瘤手术切除与重建中的应用[J]. 中国组织工程研究, 2021, 25(9): 1335-1339.
[4]	张宇, 田少奇, 曾国波, 胡川. 初次下肢全关节置换后发生心肌梗死的危险因素[J]. 中国组织工程研究, 2021, 25(9): 1340-1345.
[5]	李大地, 朱梁, 郑力, 赵凤朝. 全膝关节置换疗效及患者满意度与人格特征的关联性[J]. 中国组织工程研究, 2021, 25(9): 1346-1350.
[6]	韦玮, 李剑, 黄林海, 兰敏东, 卢显威, 黄绍东. 全膝或全髋关节置换后老年人首次活动时跌倒恐惧的影响因素[J]. 中国组织工程研究, 2021, 25(9): 1351-1355.
[7]	王金军, 邓增发, 刘康, 何智勇, 余新平, 梁建基, 李晨, 郭洲洋. 全膝关节置换静脉滴注氨甲环酸联合含氨甲环酸鸡尾酒局部应用的止血效果及安全性[J]. 中国组织工程研究, 2021, 25(9): 1356-1361.
[8]	肖国庆, 刘选泽, 严钰皓, 钟喜红. 后交叉韧带替代型假体全膝关节置换术后屈曲受限的影响因素[J]. 中国组织工程研究, 2021, 25(9): 1362-1367.
[9]	彭智浩, 冯宗权, 邹勇根, 牛国庆, 吴峰. 活动平台单髁置换后下肢力线与外侧间室骨关节炎进展的关系[J]. 中国组织工程研究, 2021, 25(9): 1368-1374.
[10]	黄泽晓, 杨妹, 林诗炜, 何和与. 血清n-3多不饱和脂肪酸水平与全膝关节置换早期股四头肌肌力变化的相关性[J]. 中国组织工程研究, 2021, 25(9): 1375-1380.
[11]	张冲, 刘志昂, 姚帅辉, 高军胜, 姜岩, 张陆. 局部应用氨甲环酸减少老年股骨颈骨折全髋关节置换后引流的安全和有效性[J]. 中国组织工程研究, 2021, 25(9): 1381-1386.
[12]	姚汝斌, 王仕永, 杨开舜. 微创经椎间孔椎间融合治疗单节段腰椎管狭窄症对腰椎-骨盆平衡的改善作用[J]. 中国组织工程研究, 2021, 25(9): 1387-1392.
[13]	王海莹, 吕冰, 李辉, 王顺义. 减压植骨融合内固定治疗退变性腰椎滑脱：基于脊柱骨盆参数预测患者的功能预后[J]. 中国组织工程研究, 2021, 25(9): 1393-1397.
[14]	吕振, 白金柱. 基于McKenzie技术的腰椎运动链训练应用于腰椎间孔镜术后分期康复的前瞻性研究[J]. 中国组织工程研究, 2021, 25(9): 1398-1403.
[15]	陈心敏, 李文标, 熊凯凯, 熊晓燕, 郑利钦, 李木生, 郑永泽, 林梓凌. 钉道强化股骨近端防旋髓内钉治疗老年A3.3型股骨转子间骨折：最佳骨水泥量有限元分析[J]. 中国组织工程研究, 2021, 25(9): 1404-1409.