杜仲活性成分抗膝骨关节炎滑膜炎病变分子机制的网络药理学阐述

doi:10.3969/j.issn.2095-4344.3013

中国组织工程研究 ›› 2021, Vol. 25 ›› Issue (5): 765-771.doi: 10.3969/j.issn.2095-4344.3013

• 组织构建临床实践 clinical practice in tissue construction • 上一篇下一篇

杜仲活性成分抗膝骨关节炎滑膜炎病变分子机制的网络药理学阐述

黎永华1，冯强1，谭仁霆1，黄世福1，邱金龙1，尹恒2

1海南省儋州市中医医院骨伤科，海南省儋州市 571700；2南京中医药大学，江苏省南京市 210023

收稿日期:2020-02-14 修回日期:2020-02-24 接受日期:2020-03-25 出版日期:2021-02-18 发布日期:2020-12-01
通讯作者: 冯强，副主任中医师，海南省儋州市中医医院骨伤科，海南省儋州市 571700
作者简介:黎永华，男，1983年生，2007年长春中医药大学毕业，主治医师，主要从事中西医结合治疗骨伤科疾病的研究。
基金资助:
国家自然科学基金(81973878)

Molecular mechanism of Eucommia ulmoides active ingredients treating synovitis of knee osteoarthritis: an analysis based on network pharmacology

Li Yonghua1, Feng Qiang1, Tan Renting1, Huang Shifu1, Qiu Jinlong1, Yin Heng2

1Department of Orthopedics, Danzhou Traditional Chinese Medicine Hospital, Danzhou 571700, Hainan Province, China; 2Nanjing University of Chinese Medicine, Nanjing 210023, Jiangsu Province, China

Received:2020-02-14 Revised:2020-02-24 Accepted:2020-03-25 Online:2021-02-18 Published:2020-12-01
Contact: Feng Qiang, Associate chief physician, Department of Orthopedics, Danzhou Traditional Chinese Medicine Hospital, Danzhou 571700, Hainan Province, China
About author:Li Yonghua, Attending physician, Department of Orthopedics, Danzhou Traditional Chinese Medicine Hospital, Danzhou 571700, Hainan Province, China
Supported by:
the National Natural Science Foundation of China, No. 81973878

摘要/Abstract

摘要：

文题释义：
基因本体(gene ontology，GO)：是一个在生物信息学领域中广泛使用的本体，它涵盖生物学的3个方面：细胞组分、分子功能、生物过程。
京都基因与基因组百科全书(KEGG) ：是了解高级功能和生物系统(如细胞、生物和生态系统)，从分子水平信息，尤其是大型分子数据集生成的基因组测序和其他高通量实验技术的实用程序数据库资源，由日本京都大学生物信息学中心的Kanehisa实验室于1995年建立。KEGG是国际最常用的生物信息数据库之一，以“理解生物系统的高级功能和实用程序资源库”著称。

背景：研究显示杜仲有效成分存在抑制关节炎滑膜成纤维细胞生长的作用，但杜仲对膝骨关节炎病变过程的滑膜炎症病变作用的分子机制目前尚不明确。
目的：基于中药网络药理学理论与生物信息学研究方法，探讨中药杜仲治疗膝骨关节炎滑膜炎症病变的可能分子作用机制。
方法：利用中药系统药理学数据库与分析平台(TCMSP)筛选杜仲有效活性成分并预测其作用的疾病靶点。下载基因表达数据库(GEO)中关于膝骨关节炎滑膜炎症病变的基因芯片数据，筛选获得差异表达基因，与杜仲活性成分作用的靶点基因取交集获得目的靶点基因。利用CytoScape软件构建杜仲-活性成分-靶点基因调控网络；使用STRING网站对靶点基因进行蛋白互作网络构建分析，筛选其中关键基因；利用R语言Bioconductor程序包对靶点基因进行基因本体论(GO)与京都基因与基因组百科全书(KEGG)信号通路分析，筛选具有统计学差异(P < 0.05)的生物学过程或信号通路展开进一步分析。
结果与结论：①筛选获得杜仲的21个有效活性成分，包括丁子香萜、杜仲甙、杜仲脂素A、槲皮素等，作用于25个膝骨关节炎滑膜炎症病变靶点基因(7个上调基因、18个下调基因)；②靶点基因主要富集于细胞因子活性、细胞因子受体结合、激活转录因子结合等生物学过程；③杜仲有效活性成分可能通过调控白介素6、血管内皮生长因子A、细胞趋化因子8等关键基因的表达，进而调控靶点基因所在的丝裂原活化蛋白激酶信号通路、核因子κB信号通路、Toll样受体信号通路等发挥其治疗滑膜炎症病变的作用；④研究立足于对膝骨关节炎滑膜炎症病变中医病机的辨证论治，运用网络药理学与生物信息学理论方法探讨杜仲治疗该疾病的可能分子作用机制，可为进一步开展相关实验研究提供药效物质基础及思路方向参考。

https://orcid.org/0000-0001-5681-4757 (黎永华)

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

关键词: 杜仲, 骨, 膝骨关节炎, 滑膜炎症, 通路, 网络药理学, 生物信息学

Abstract: BACKGROUND: Studies have shown that the active ingredients of Eucommia ulmoides can inhibit the growth of synovial fibroblasts of arthritis, but the molecular mechanism underlying Eucommia ulmoides treating synovial inflammation of knee osteoarthritis is not yet clear.
OBJECTIVE: To investigate the possible molecular mechanism of Eucommia ulmoides ingredients treating synovitis of knee osteoarthritis based on the network pharmacology and bioinformatics.
METHODS: The active ingredients and targets of Eucommia ulmoides were screened by using the TCM systematic pharmacological analysis platform (TCMSP). The disease targets corresponding to synovitis of knee osteoarthritis were obtained from Gene Expression Omnibus (GEO). The protein-protein interaction network was constructed by STRING online database, analyzed and showed by the Cytoscape software. Gene ontology analysis and Kyoto encyclopedia of genes and genomes pathway enrichment analysis of the targets were conducted by R/Bioconductor. The enrichment results were obtained with a significant difference (P < 0.05).
RESULTS AND CONCLUSION: A total of 21 active ingredients of Eucommia ulmoides, including eugenol, eucommia glycoside, eucommia liposide A, quercetin, and 25 different expressed genes (7 up-regulated and 18 down-regulated differentially expressed genes) were obtained. The main biological processes of the differentially expressed genes were found different significantly in cytokine activity, cytokine receptor binding and activating transcription factor binding. The hub genes, such as interleukin-6, vascular endothelial growth factor A and chemotactic factor 8, may be regulated by the active ingredients to influence the signal of MAPK/NF-kappa B/Toll-like receptor signaling pathway. Based on the network pharmacology and bioinformatics approach, the study preliminarily validates the molecular mechanism of Eucommia ulmoides for treatment of synovitis of knee osteoarthritis, which may lay a foundation for further pharmacodynamics study and experiment.

Key words: Eucommia ulmoides, bone, knee osteoarthritis, synovitis, pathway, network pharmacology, bioinformatics

中图分类号:

黎永华, 冯强, 谭仁霆, 黄世福, 邱金龙, 尹恒. 杜仲活性成分抗膝骨关节炎滑膜炎病变分子机制的网络药理学阐述[J]. 中国组织工程研究, 2021, 25(5): 765-771.

Li Yonghua, Feng Qiang, Tan Renting, Huang Shifu, Qiu Jinlong, Yin Heng. Molecular mechanism of Eucommia ulmoides active ingredients treating synovitis of knee osteoarthritis: an analysis based on network pharmacology[J]. Chinese Journal of Tissue Engineering Research, 2021, 25(5): 765-771.

图/表（结果） 8

2.1 杜仲有效活性成分及作用靶点收集从TCMSP共检索获得杜仲活性成分147个，根据ADME药动学参数标准进行筛选，最终获得同时满足口服生物利用度(OB)>30%与类药性(DL)>0.18两个条件的杜仲有效活性成分28个，结果见表1。收集TCMSP提供的以上28个杜仲有效活性成分作用疾病靶点，对重复出现者进行剔除，共获得181个作用靶点基因。

2.2 滑膜炎症病变靶点基因收集借助R语言Limma程序包将GSE46750、GSE55235、GSE82107三个基因芯片获得的32个滑膜炎症组织样本基因数据与29个正常对照组滑膜组织样本基因数据进行对比分析，并根据P < 0.05与|FC|>2筛选显著差异表达的基因。结果获得表达上调基因171个，表达下调基因486个，对从基因芯片获得的全部基因根据校正后的P值及差异倍数取2的对数进行火山图绘制，展示上调与下调基因，见图1。

2.3 杜仲有效活性成分调控滑膜炎症病变靶点基因收集将从TCMSP筛选后获得的28个杜仲有效活性成分作用的181个疾病靶点，分别与GEO基因芯片获得的滑膜炎症病变差异表达基因取交集，获得参与调控膝骨关节炎滑膜炎症病变的杜仲有效活性成分21个，上调表达基因7个，下调表达基因18个，以韦恩图(Venn Diagram)形式展示交集基因结果，见图2。

2.4 杜仲-活性成分-滑膜炎症靶点基因调控网络可视化将上述获得的21个作用于滑膜炎症病变的杜仲有效活性成分与差异表达基因数据共同导入至CytoScape软件中，构建杜仲-活性成分-滑膜炎症靶点基因调控网络，将分析结果以网络可视化图形呈现，网络图中共包含药物(杜仲)、活性成分以及调控滑膜炎症病变的上下调基因共47个节点，涉及77个相互作用关系，各节点间通过边线相连接，形成紧密的作用关系网络，见图3。

2.5 靶点基因互作网络关键基因筛选将杜仲有效活性成分作用滑膜炎症病变差异表达基因25个上传至STRING网站，设置物种为“Homo sapiens”(人类)，最小互作分数为0.4，运行靶点基因蛋白互作关系分析，将结果以TSV格式下载保存。通过R语言对靶点基因蛋白互作关系网络中靶点Degree值进行计算并排序，以条形图形式展现Degree值从高到低排名前20位的靶点基因及其Degree值，结果见图4。经计算，在差异基因互作网络中Degree值最高的是白细胞介素6(IL6)，其次为血管内皮生长因子A(VEGFA)、细胞趋化因子8(CXCL8)、FOS、JUN、基质金属蛋白酶9(MMP9)、PTGS2等，可能是杜仲活性成分调控膝骨关节炎滑膜炎症病变的关键基因。

2.6 靶点基因GO富集分析对杜仲有效活性成分作用滑膜炎症靶点基因进行整理，运行R语言Bioconductor clusterProfiler程序包对靶点基因进行GO富集分析[21]，筛选检验统计P < 0.05的基因产物注释过程，并按照P值从小到大进行排序，展示排名前20的生物学过程，见图5。GO富集分析结果显示，靶点基因主要富集于细胞因子活性(Cytokine activity)、细胞因子受体结合(Cytokine receptor binding)、超氧化物-生成NADPH氧化酶活性(Superoxide-generating NADPH oxidase activity)、激活转录因子结合(Activating transcription factor binding)、受体配体活性(Receptor ligand activity)、R-SMAD 结合(R-SMAD binding)、生长因子受体结合(Growth factor receptor binding)、血红素结合(Heme binding)等生物学过程。

2.7 靶点基因KEGG通路富集分析继续使用R语言Bioconductor clusterProfiler程序包对靶点基因进行KEGG通路富集分析，筛选出P < 0.05的富集信号通路，根据P 值从小到大排序，对排名前20的信号通路进行可视化展示，见图6。KEGG通路富集分析结果显示，靶点基因主要富集于丝裂原活化蛋白激酶信号通路(MAPK signaling pathway)、核因子κB信号通路(NF-kappa B signaling pathway)、Toll样受体信号通路(Toll-like receptor signaling pathway)、缺氧诱导因子1信号通路(HIF- 1 signaling pathway)、破骨细胞分化信号通路(Osteoclast differentiation)、细胞衰老信号通路(Cellular senescence)等通路。

使用R语言Bioconductor Pathview程序包对靶点基因富集的KEGG信号通路进行细节展示[22]，将P < 0.05的信号通路内容可视化，以便进一步了解靶点基因在各信号通路的各个阶段参与的信号转导过程。见图7，对靶点基因富集程度最高的MAPK信号通路细节展示，靶点基因CACNA1S、MYC、IL1B、FOS等分别富集于MAPK信号通路的不同环节，参与MAPK信号通路上下游的信号转导，调控基因转录、胞内外离子代谢等环节，并与肿瘤坏死因子信号通路(TNF signal pathway)、细胞增殖分化(Proliferation differentiation)、炎症(Inflammation)、细胞周期(Cell cycle)等产生级联信号转导沟通。

参考文献

[1] CAI W, LI H, ZHANG Y, et al. Identification of key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis. Peer J.2020; 8(8390):1-18
[2] MATHIESSEN A, CONAGHAN PG. Synovitis in osteoarthritis: current understanding with therapeutic implications. Arthritis Res Ther. 2017;19(1):18.
[3] MIAO CG, YANG YY, HE X, et al. Wnt signaling pathway in rheumatoid arthritis, with special emphasis on the different roles in synovial inflammation and bone remodeling. Cell Signal. 2013;25(10):2069-2078.
[4] LI N, XU Q, LIU Q, et al. Leonurine attenuates fibroblast-like synoviocyte-mediated synovial inflammation and joint destruction in rheumatoid arthritis. Rheumatology (Oxford, England), 2017; 56(8): 1417-1427.
[5] 陈卫衡,刘献祥,童培建,等.膝骨关节炎中医诊疗专家共识(2015年版) [J].中医正骨, 2015, 27(07): 4-5.
[6] 梁子聪,易艾晶.基于文献数据挖掘调节膝骨关节炎炎症指标的中药内服方药用药规律研究[J]. 黔南民族医专学报, 2018, 31(2): 90-94.
[7] 刘磊, 张舒, 周悦悦. 复方杜仲健骨颗粒联合硫酸氨基葡萄糖治疗膝骨关节炎的临床研究[J]. 现代药物与临床, 2019, 34(11): 3343-6.
[8] 杨亚旭.蜂毒素联合杜仲黄酮治疗类风湿关节炎的实验研究[D].扬州:扬州大学, 2017.
[9] 张彦琼,李梢.网络药理学与中医药现代研究的若干进展[J].中国药理学与毒理学杂志, 2015, 29(6): 883-92.
[10] LIU H, WANG J, ZHOU W, et al. Systems approaches and polypharmacology for drug discovery from herbal medicines: an example using licorice.J Ethnopharmacol. 2013;146(3): 773-793.
[11] 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.
[12] 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.
[13] WANG S, WANG H, LU Y. Tianfoshen oral liquid: a CFDA approved clinical traditional Chinese medicine, normalizes major cellular pathways disordered during colorectal carcinogenesis.Oncotarget.2017;8(9): 14549-14569.
[14] RITCHIE ME, PHIPSON B, WU D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
[15] SHANNON P, MARKIEL A, OZIER O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498-2504.
[16] 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.
[17] SZKLARCZYK D, GABLE A L, LYON D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.Nucleic Acids Res. 2019; 47(D1):D607-D613.
[18] DE SANTIAGO I, CARROLL T. Analysis of ChIP-seq Data in R/Bioconductor.Methods Mol Biol. 2018;1689:195-226.
[19] DU J, YUAN Z, MA Z, et al. KEGG-PATH: Kyoto encyclopedia of genes and genomes-based pathway analysis using a path analysis model. Mol Biosyst. 2014;10(9):2441-2447.
[20] PRESS OU. Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. 2017;45(D1):D331-D338.
[21] YU G, WANG L G, HAN Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters]. OMICS. 2012;16(5):284-287.
[22] LUO W,BROUWER C.Pathview: an R/Bioconductor package for pathway-based data integration and visualization.Bioinformatics. 2013;29(14):1830-1831.
[23] 吕苏梅, 张瑞丽.中老年膝骨关节炎的流行病学研究进展[J].中国老年学杂志, 2016, 36(16):4133-5.
[24] BENITO MJ, VEALE DJ, FITZGERALD O, et al. Synovial tissue inflammation in early and late osteoarthritis. Ann Rheum Dis. 2005;64(9): 1263-1267.
[25] O’NEILL TW, PARKES MJ, MARICAR N, et al. Synovial tissue volume: a treatment target in knee osteoarthritis (OA). Ann Rheum Dis. 2016;75(1):84-90.
[26] ZHAI Y, GAO GD, XU SY. Basic research progress of knee osteoarthritis. Zhongguo Gu Shang. 2012;25(1):83-87.
[27] 陈达,陈后煌,邵翔,等.乌头汤抑制骨关节炎炎症反应的作用机制探讨 [J]. 风湿病与关节炎, 2016, 5(8): 62-66.
[28] CENGIC T, TRKULJA V, PAVELIC SK, et al. Association of TGFB1 29C/T and IL6 -572G/C polymorphisms with developmental hip dysplasia: a case-control study in adults with severe osteoarthritis. Int Orthop. 2015;39(4): 793-798.
[29] FUJIHASHI K, KONO Y, KIYONO H. Effects of IL6 on B cells in mucosal immune response and inflammation. Res Immunol. 1992;143(7): 744-749.
[30] FRANCKHAUSER S, ELIAS I, ROTTER SOPASAKIS V, et al. Overexpression of Il6 leads to hyperinsulinaemia, liver inflammation and reduced body weight in mice. Diabetologia. 2008;51(7): 1306-1316.
[31] PEARSON MJ, HERNDLER-BRANDSTETTER D, TARIQ MA, et al. IL-6 secretion in osteoarthritis patients is mediated by chondrocyte-synovial fibroblast cross-talk and is enhanced by obesity. Sci Rep. 2017;7(1):3451.
[32] ZUPAN J, VRTACNIK P, COR A, et al. VEGF-A is associated with early degenerative changes in cartilage and subchondral bone. Growth Factors. 2018;36(5-6):263-273.
[33] HAMILTON JL, NAGAO M, LEVINE BR, et al. Targeting VEGF and Its Receptors for the Treatment of Osteoarthritis and Associated Pain.J Bone Miner Res. 2016;31(5):911-924.
[34] ERDEM H, PAY S, MUSABAK U, et al. Synovial angiostatic non-ELR CXC chemokines in inflammatory arthritides: does CXCL4 designate chronicity of synovitis? . Rheumatol Int. 2007;27(10):969-973.
[35] YANG CC, LIN CY, WANG HS, et al. Matrix metalloproteases and tissue inhibitors of metalloproteinases in medial plica and pannus-like tissue contribute to knee osteoarthritis progression.PLoS One. 2013;8(11):e79662.
[36] LIN J, WU G, ZHAO Z, et al. Bioinformatics analysis to identify key genes and pathways influencing synovial inflammation in osteoarthritis. Mol Med Rep. 2018;18(6):5594-5602.
[37] WANG D, QIAO J, ZHAO X, et al. Thymoquinone Inhibits IL-1beta-Induced Inflammation in Human Osteoarthritis Chondrocytes by Suppressing NF-kappaB and MAPKs Signaling Pathway.Inflammation. 2015;38(6):2235-2241.
[38] HASEEB A, HAQQI T M. Immunopathogenesis of osteoarthritis. Clin Immunol. 2013;146(3): 185-196.
[39] WU L, HUANG X, LI L, et al. Insights on biology and pathology of HIF-1alpha/-2alpha, TGFbeta/BMP, Wnt/beta-catenin, and NF-kappaB pathways in osteoarthritis.Curr Pharm Des. 2012;18(22):3293-3312.
[40] ZHANG FJ, LUO W, LEI GH. Role of HIF-1alpha and HIF-2alpha in osteoarthritis. Joint bone spine.2015; 82(3): 144-147.
[41] FERNANDEZ-TORRES J, MARTINEZ-NAVA GA, GUTIERREZ-RUIZ MC, et al. Role of HIF-1alpha signaling pathway in osteoarthritis: a systematic review. Rev Bras Reumatol Engl Ed. 2017;57(2):162-173.
[42] NAIR A, KANDA V, BUSH-JOSEPH C, et al. Synovial fluid from patients with early osteoarthritis modulates fibroblast-like synoviocyte responses to toll-like receptor 4 and toll-like receptor 2 ligands via soluble CD14. Arthritis Rheum. 2012;64(7):2268-6877.
[43] STANNUS O, JONES G, CICUTTINI F, et al. Circulating levels of IL-6 and TNF-alpha are associated with knee radiographic osteoarthritis and knee cartilage loss in older adults. Osteoarthritis Cartilage. 2010; 18(11):1441-1447.

引言

膝骨关节炎是骨科临床常见的慢性退行性疾病，世界卫生组织流行病学研究显示，该病在全球范围内60岁以上人群中男性发病率约为9.6%，女性约为18%，并有25%的致残率[1]。长久以来研究者多认为膝骨关节炎的发病以关节软骨退化、软骨下骨增生、骨赘形成等为主要病理改变[2]，但对病变过程中的非感染性关节滑膜炎症改变却未足够重视。近来有研究显示滑膜炎症与调节关节软骨代谢、细胞因子与蛋白酶的理化作用等密切相关[3- 4]，但目前对关节滑膜炎症发病的分子机制研究仍较少。
根据患者膝关节肿胀、疼痛、积液、活动受限等常见表现，中医将膝骨关节炎归类为“痹症”“骨痹”范畴，其病机总属肝肾亏虚[5]；滑膜炎症病变为关节炎病变发展的一部分，基本病机亦属肝肾亏虚。中药杜仲具有补肝肾、强筋骨之功，是中医骨伤科治疗膝骨关节炎的常见药物，被广泛应用于中药方剂或中成药配伍中[6-7]。研究显示杜仲有效成分存在抑制关节炎滑膜成纤维细胞生长的作用[8]，具有治疗滑膜炎症病变的前景，但杜仲对膝骨关节炎病变过程的滑膜炎症病变作用的分子机制目前尚不明确。
获益于生物信息学与大数据分析的发展，中药网络药理学分析也飞速发展，成为一个探索预测中药治疗疾病潜在作用机制的重要工具，为实施具体实验提供了方向[9]。此次研究基于网络药理学分析与生物信息学数据挖掘理论方法，探索杜仲活性成分对膝骨关节炎滑膜炎症病变影响的可能分子机制，以期为膝骨关节炎的治疗提供新的靶点与方向。
中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

材料方法

1.1 杜仲有效活性成分筛选及作用靶点预测采用中药系统药理学数据库及分析平台(TCMSP，http://tcmspw.com/tcmsp.php)对杜仲活性成分进行收集，根据口服中药在人体内的药动学特性(ADME)参数对杜仲活性成分进行筛选[10- 11]。以药物的口服生物利用度(oral bioavailability，OB)大于30%、类药性(drug like，DL)大于0.18为标准筛选获得杜仲有效活性成分[12-13]，并将TCMSP提供的杜仲有效活性成分作用疾病靶点信息收集整理。
1.2 滑膜炎症靶点基因挖掘筛选在美国国立生物技术信息中心基因表达数据库(GENE EXPRESSION OMNIBUS，GEO)以“Knee Osteoarthritis”(膝骨关节炎)、“Synovitis”(滑膜炎)为关键词检索相关基因芯片，获取编号为GSE46750、GSE55235、GSE82107的基因芯片原始数据及对应的基因注释文件。GSE46750芯片包含12个正常患者滑膜组织与12个非感染性炎症滑膜组织，GSE55235芯片包含10个正常关节滑膜组织与10个非感染性关节炎滑膜组织，GSE82107芯片包含7个正常对照组滑膜组织与10个非感染性炎症滑膜组织。
运用R语言对获取的3个基因芯片原始数据进行多矩阵数据归一与背景校正，运行R语言Limma程序包对基因芯片进行差异基因分析[14]，以检验统计P < 0.05、差异倍数(fold change，FC)绝对值>2为筛选条件，获得骨关节炎滑膜炎症病变的差异表达靶点基因。差异基因的FC>2则为上调基因，FC<2则为下调基因。
1.3 杜仲有效活性成分调控滑膜炎症靶标基因获取将收集的杜仲有效活性成分作用疾病靶点基因，分别与滑膜炎症病变的上调基因与下调基因取交集，获得杜仲有效活性成分作用于膝骨关节炎滑膜炎症病变的上调基因与下调基因。
1.4 杜仲-活性成分-靶点基因调控网络构建采用生物信息网络可视化分析软件CytoScape software 3.6.1对杜仲、有效活性成分及滑膜炎症差异表达靶点基因进行互作调控网络的构建[15]。将获得的有效活性成分、上调下调基因分类整理导入CytoScape软件中，以“节点”(Node)表示杜仲、活性成分及靶点基因，以“边线”(Edge)表示不同节点间的相互作用关系[16]，节点、边线互相连接，构成可视化网络。
1.5 靶点基因蛋白互作网络构建与关键基因筛选采用基因蛋白互作在线分析工具STRING(https://string-db.org/)对获得的杜仲有效活性成分作用滑膜炎症靶点基因进行蛋白质-蛋白质互作(protein-protein interaction，PPI)网络构建[17]。STRING网站对网络中靶点基因进行连接度(Degree)计算，Degree值越大，代表该靶点在网络中与更多基因发生互作关系，提示其重要性可能越大。根据Degree值大小进行靶点基因重要性排序，筛选出其中的关键基因。
1.6 靶点基因信号通路富集与基因本体论分析使用R语言生物信息学分析工具包Bioconductor对获得的杜仲有效活性成分作用滑膜炎症靶点基因进行京都基因与基因组百科全书(Kyoto Encyclopedia of genes and genomes，KEGG)通路富集分析与基因本体论(gene ontology，GO)分析[18-20]，预测靶点基因主要参与的生物学功能注释过程与分子信号通路，以检验统计P < 0.05为差异有统计学意义。

讨论

膝骨关节炎是导致中老年人运动能力丧失的常见疾病，发病率高、致残率高，在全球老龄化趋势下给国家社会带来了沉重的经济负担[23]。膝关节非感染性滑膜炎症病变常被认为是膝骨关节炎过程中关节软骨退化、骨赘增生等结构变化带来的继发反应[24]。O’NEILL等[25]研究发现，中重度膝骨关节炎患者滑膜存在病理改变，不单纯是软骨退化、骨赘增生带来的继发反应，滑膜炎症病变对关节软骨、骨质会带来进一步破坏，导致膝骨关节炎患者关节肿痛等症状的加重，但目前滑膜炎症病变发病机制尚不明确，从分子机制探索其发病与治疗有着重要的理论价值与治疗导向作用[25]。
中医理论一般将膝骨关节炎归于“痹症”“骨痹”范畴，其基本病机为肝肾亏虚，多因机体阳虚血瘀、精血不足、风寒湿邪外侵等原因，致气血凝涩、骨骼失养、筋脉痹阻发而为病[26-27]。杜仲为杜仲科植物杜仲的树皮，性甘微辛温，入肺肝肾经，功擅补肝肾，强筋骨，是中医骨伤科治疗骨关节炎的常用药物。此次研究立足中医辨证论治原则，结合网络药理学与生物信息学理论方法，对中药杜仲有效活性成分调控膝骨关节炎滑膜炎症病变基因产物相关分子机制展开探索。
此次研究基于中药系统药理学数据库及分析平台TCMSP共筛选获得杜仲中的丁子香萜、杜仲甙、杜仲脂素A、槲皮素等21个有效活性成分，直接作用于GEO基因芯片筛选获得的关节滑膜炎症病变的7个上调表达基因和18个下调表达基因，并根据药物-有效活性成分-差异表达基因之间的联系绘制出复杂调控网络图。
从图3网络图中可见，杜仲有效活性成分对滑膜炎症病变差异表达靶点基因调控作用存在以下几个特点：①山柰酚、槲皮素、β-胡萝卜素等活性成分均同时与多个靶点基因存在作用关系，提示同一活性成分对滑膜炎症病变的调控作用可能通过不同靶点基因发挥不同的作用；②靶点基因PTGS2、基质金属蛋白酶1、JUN等同时受到多个杜仲活性成分的作用，提示不同活性成分对同一靶点基因的调控可能存在协同作用；③部分活性成分同时作用于上调基因与下调基因，如槲皮素同时作用于22个靶点基因，其中包括6个下调基因，16个上调基因，提示其对关节滑膜炎症病变的作用可能是双向调控的。中药杜仲治疗膝骨关节炎滑膜炎症病变存在多成分、多靶点、协同调控、双向调控等特点，机制复杂，与中医辨证论治、整体观念思想相契合。
膝骨关节炎中滑膜炎症病变发病机制复杂，从GEO数据库提供的相关基因芯片测序结果获得了数以千计的基因与其生物学过程相关，根据统计学差异筛选出676个表达差异明显的基因，其中上调基因171个，表达下调基因486个，下调基因数量明显高于上调基因，占差异表达基因总数的71.89%。进一步与杜仲有效活性成分直接作用的疾病靶点基因取交集，获得其参与调控的7个上调基因与18个下调基因。靶点基因之间存在复杂的互作关系，利用STRING网站对所有靶点基因进行蛋白互作网络构建与分析，根据靶点基因之间连接度Degree值大小排序，初步认识靶点基因之间的联系与重要性。
在Degree值较高的前几位关键靶点基因中，分别包括白细胞介素6、血管内皮生长因子A、JUN、PTGS2等下调表达基因，细胞趋化因子8、FOS、基质金属蛋白酶9等上调表达基因。白细胞介素6是骨关节炎发病过程的重要分子[28]，研究证明白细胞介素6表达水平与骨骼肌肉系统非感染性炎症反应存在明显相关性，白细胞介素6的高表达将增加炎症组织中B淋巴细胞反应区域从而调动机体免疫系统，并影响关节软骨基质的退化[29-31]。血管内皮生长因子A通过刺激血管内皮细胞有丝分裂活动而促进血管新生，这对局部营养物质交换及炎症反应扩散是有利的[32]，血管内皮生长因子A适度的基因表达有助于滑膜细胞在缺氧环境中的修复，但过高的表达量反而导致滑膜细胞异常增殖及滑膜炎症的加剧[33]。细胞趋化因子8表达水平与滑膜炎症中血管新生及中性粒细胞游走密切相关[34]，CAI等[1]通过将正常关节滑膜组织与膝骨关节炎患者滑膜组织进行免疫浸润分析，发现关节炎患者中细胞趋化因子8表达水平明显升高，可能与细胞趋化因子8激活记忆B淋巴细胞、CD4+T淋巴细胞等免疫细胞的活动有关。基质金属蛋白酶9的作用主要是细胞外基质的降解和变革，高表达的基质金属蛋白酶9与滑膜炎症病变加重存在相关性，基质金属蛋白酶9是治疗滑膜炎症的重要靶点[35-36]。此次研究中筛选获得的多个杜仲有效活性成分调控的滑膜炎症靶点基因已被证明是对膝骨关节炎发病产生重要作用的靶点，是治疗滑膜炎症、减轻关节炎症状的重要靶点，杜仲活性成分可能通过对上述关键靶点基因的表达双向调控，从而调控机体的生物学过程、信号通路转导等环节发挥治疗作用。
靶点基因的KEGG信号通路富集分析提示杜仲活性成分调控滑膜炎症病变可能主要从炎症相关、细胞周期调控等通路发挥作用。
(1)细胞周期相关通路：①MAPK信号通路：丝裂原活化蛋白激酶是真核生物中非常保守的丝氨酸蛋白激酶，其下游信号通路与细胞增殖、分化，运动或凋亡有着密切联系，滑膜炎症产生的炎症因子激活MAPK信号通路，可导致基质金属蛋白酶升高而破坏关节软骨细胞[37]，也可能通过调控基因转录、胞内外离子代谢等环节影响滑膜细胞的增殖、分化与凋亡，同时MAPK信号通路与TNF信号通路、P53信号通路、炎症信号通路等均存在通路之间的沟通，调控作用复杂；②NF-kappa B信号通路：NF-kappa B信号通路激活诱导滑膜细胞产生大量趋化因子、细胞因子及基质金属蛋白酶，促进关节软骨分解破坏，加重对滑膜组织的刺激作用[38-39]；③缺氧诱导因子1信号通路：缺氧诱导因子1是细胞缺氧状态下广泛参与外环境酸碱调控、细胞能量代谢、血管新生与细胞周期等生理进程的重要分子，介导软骨细胞、滑膜细胞在缺氧状态下的损伤修复[40- 41]。
(2)炎症相关通路：①Toll样受体信号通路：骨关节炎患者机体CD14免疫因子可激活滑膜成纤维样细胞的Toll样受体(TLR2、TLR4)，产生炎症反应[42]；②肿瘤坏死因子信号通路：肿瘤坏死因子具有明确的炎症介导作用，肿瘤坏死因子α可诱导白细胞介素6产生并进一步提高肿瘤坏死因子α活性，激活分解软骨、滑膜的蛋白酶，介导炎症发生[43]。
结论：综上所述，研究基于对膝骨关节炎滑膜炎症病变的中医肝肾亏虚病机辨证，利用生物信息学及网络药理学方法探索补肝肾、强筋骨中药杜仲有效活性成分调控滑膜炎症病变的可能分子机制。研究筛选获得21个参与调控膝骨关节炎滑膜炎症病变的杜仲活性成分与25个靶点基因，其调控作用可能通过抑制局部炎症、调控软骨细胞与滑膜细胞的细胞周期活动发挥。膝骨关节炎及其滑膜炎症病变机制复杂，中药杜仲活性成分具有治疗前景，后续仍需要设计严谨的具体实验对药物成分作用、靶点基因作用进一步验证，探索各活性成分在相关通路中对靶点基因的表达干预作用与调控作用，并最终服务于临床疾病治疗中，促进中医药现代化发展。中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

杜仲活性成分抗膝骨关节炎滑膜炎病变分子机制的网络药理学阐述

Molecular mechanism of Eucommia ulmoides active ingredients treating synovitis of knee osteoarthritis: an analysis based on network pharmacology

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 8

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

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