0 引言 Introduction
骨关节炎是一种以关节疼痛、功能障碍为主要表现的骨科常见疾病，在中医属于“痹证”范畴[1]，中国历代医家认为“风、寒、湿”三邪是该病的主要病因[2]，三者可侵犯体表，阻碍经络气血运行，如《黄帝内经》中“风寒湿三气杂至，合而为痹”。现代医学发现骨关节炎发病机制为多种因素引起关节软骨退化及关节边缘炎性增生，西医主要以非类固醇类抗炎药、降钙素及二膦酸类药物进行消炎镇痛、改善关节功能为主，但上述药物的长期使用会对机体骨质、胃肠道及心血管造成损伤，具有一定的治疗局限性[3]。而中医中药具有药物毒性低、疗效确切的优势，但方剂中药种类较多，中药有效成分又过于复杂，这加大了通过药理实验研究复方作用机制的难度与经费消耗，限制了中药治疗骨关节炎的研究发展[4]。
网络药理学作为一门新学科被广泛运用于研究中药治疗疾病潜在的分子机制，它通过药物-基因-疾病多层次的网络中药机制分析，系统预测中药对致病基因潜在的药效靶点，科学阐述中药可能在疾病中的作用方式，为中药治疗疾病提供理论基础支持，这极大降低了科研难度与经费耗损，促进中医药事业的推广[5-6]。
五苓散源于《伤寒杂病论》，该方以泽泻为君，利水渗湿；以茯苓猪苓为臣，增强君药利水渗湿之力，以白术、桂枝为佐，既健脾利湿，又助阳化气，温经通络[7]。因此骨关节证属“风、寒、湿”三者都可在五苓散加减的基础上进行治疗，并且目前临床运用五苓散治疗膝关节炎疗效确切，但其药理作用、分子机制及药物活性物质尚未明确，因此通过网络药理学对五苓散治疗骨关节炎的作用机制进行系统、整体研究，为五苓散治疗骨关节炎提供一定理论基础支持。
1 资料和方法 Data and methods
1.1 活性成分的筛选 文章通过使用中药系统药理数据库和分析平台(TCMSP)(
1.2 活性成分靶点的预测 利用TCMSP数据库查询五苓散各个活性成分对应的靶蛋白，并通过蛋白数据库(UniProt)( Sapiens)检索，获取各个靶蛋白所对应的靶基因名，建立获取的靶基因名合集，删除重复基因，所得结果即为五苓散活性成分的靶点[9]。
1.3 骨关节炎治疗靶点的预测 以“骨关节炎”的英文“Osteoarthritis”为关键词分别在Gene Cards ( DigSee ( 病数据库检索相关靶点，整合3个数据库的查找结果并删除重复项，最后所得靶点即为骨关节炎的相关靶点。
1.4 潜在作用靶点的预测 将五苓散中活性成分的靶点和骨关节炎的相关靶点上传至Venny平台(
1.5 “单味药-活性成分-作用靶点”网络的构建 构建单位药活性成分及其作用靶点、靶点与骨关节炎的相互作用关系，并导入到Cytoscape 3.7.2 软件 ( /)， 通 过Cytoscape 软件中“merge”工具将五苓散各个中药的网络进行合并，形成五苓散治疗骨关节炎的“单味药-活性成分-作用靶点”网络图及“五苓散治疗骨关节炎的调控网络”[10-11]。
1.6 PPI网络的构建及网络拓扑分析Bisogenet (Version 3.0.0)是基于Cytoscape的一个新的基因网络构建、可视化和分析工具，搭载DIP，BIOGRID，HPRD，INTACT，MINT和BIND等分析蛋白质之间相互作用的数据库。为了进一步研究五苓散的作用机制，将五苓散的作用靶点导入到Bisogenet中，限定研究物种为“Homo Sapiens”获得蛋白互作关系，并运用Cytoscape中的CytoNCA插件，设置 degree 值 (DC)＞61 和 Betweenness值(BC)＞200，进行筛选富集程度高的蛋白互作关系，构建PPI核心网络图。
1.7 GO富集分析与KEGG富集分析通过DAVID数据库( 6.8)对潜在靶点进行GO功能富集分析与KEGG通路富集分析，分别研究五苓散治疗骨关节炎的主要生物功能与主要信号通路，其中P＜ 0.05代表富集结果显著，并通过R语言软件制作GO富集分析图及KEGG富集分析气泡图。
1.8 KEGG关系网络构建 先分别建立通路、基因相互作用关系和节点属性、基因相互作用关系，并导入Cytoscape 3.7.2软 件 ( /)中，选取前20个KEGG信号通路，构建五苓散治疗骨关节炎的“KEGG-基因”的关系网络图。
1.9 主要观察指标 ①主要的单味药-活性成分-作用关键化合物成分；②PPI网络及网络拓扑分析后五苓散关键蛋白基因；③筛选后关键基因在KEGG和GO富集的结果；④关键基因与关键化合物分子对接结果。
2 结果 Results
2.1 活性化合物的筛选结果 利用TCMSP数据库检索共收集得到五苓散的化学成分436个，其中茯苓、猪苓、泽泻、桂枝和白术分别为44，41，61，220和70个，根据设定的OB≥30%与DL≥0.18的筛选标准进行筛选，得到复方中共50个活性化学成分，猪苓、茯苓、白术、泽泻和桂枝的活性成分分别为11，15，7，10和7个，其中茯苓和猪苓有2个交集化合物，桂枝和猪苓有1个交集化合物，活性成分的基本信息见表1。
2.2 五苓散治疗骨关节炎潜在作用靶点预测 通过TCMSP数据库得到五苓散中活性成分靶点52个；通过检索Gene Cards，OMIM和DigSee数据库，共获得骨关节炎相关靶点2 854个，运用平台取得共有基因26个，即五苓散治疗骨关节炎的潜在作用靶点，见图1。
表1 ｜五苓散中潜在活性成分基本信息Table 1 ｜Basic information of potential active components inWulingpowder表柱：OB为口服生物利用度，DL为药物相似性来源 化学成分OB(%)DL白术 12-千里光酰基-8-反式白术三醇 62.40 0.22白术 14-乙酰基-12-千里光酰基-8-反式白术三醇 60.31 0.31白术 14-乙酰基-12-千里光酰基-8顺式白术三醇 63.37 0.30白术 α-香树脂醇 39.51 0.76白术 胆固醇衍生物 (3S，8S，9S，10R，13R，14S，17R)-10，13-dimethyl-17-[(2R，5S)-5-propan-2-yloctan-2-yl]-2，3，4，7，8，9，11，12，14，15，16，17-dodecahydro-1H-cyclopenta[a]phenanthren-3-ol 36.23 0.78白术 3β-乙酰氧基苍术酮54.07 0.22白术 8β-乙氧基苍朮内酯-Ⅲ 35.95 0.21茯苓 (2R)-2-[(3S，5R，10S，13R，14R，16R，17R)-3，16-dihydroxy-4，4，10，13，14-pentamethyl-2，3，5，6，12，15，16，17-octahydro-1H-cyclopenta[a]phenanthren-17-yl]-6-methylhept-5-enoic acid 30.93 0.81茯苓 栓菌酸38.71 0.80茯苓 7，9(11)-去氢茯苓酸 35.11 0.81茯苓、猪苓 啤酒甾醇 37.96 0.77茯苓 (2R)-2-[(3S，5R，10S，13R，14R，16R，17R)-3，16-dihydroxy-4，4，10，13，14-pentamethyl-2，3，5，6，12，15，16，17-octahydro-1H-cyclopenta[a]phenanthren-17-yl]-5-isopropyl-hex-5-enoic acid 31.07 0.82茯苓、猪苓 7，22-麦角甾二烯-3-醇43.51 0.72茯苓 麦角甾醇过氧化物 40.36 0.81茯苓 (2R)-2-[(5R，10S，13R，14R，16R，17R)-16-羟基-3-酮-4，4，10，13，14-五甲基-1，2，5，6，12，15，16，17-八氢环戊烯并[a]菲-17-基-5-异丙基六角-5-烯酸38.26 0.82茯苓 3β-羟基-24-甲基-8-羊毛甾烯-21-甲酸 38.70 0.81茯苓 茯苓酸 33.63 0.81茯苓 茯苓新酸 A 30.61 0.76茯苓 茯苓新酸 B 30.52 0.75茯苓 茯苓新酸 C 38.15 0.75茯苓 常春藤皂苷元 36.91 0.75茯苓 二氢齿孔酸 44.17 0.83桂枝 (-)-花旗松素 60.51 0.27桂枝 β-谷甾醇 36.91 0.75桂枝 谷甾醇 36.91 0.75桂枝 (+)-儿茶素 54.83 0.24桂枝 表儿茶素 48.96 0.24桂枝 花旗松素 57.84 0.27桂枝、猪苓 过氧麦角甾醇 44.39 0.82泽泻 谷甾醇 36.91 0.75泽泻 泽泻醇 B 34.47 0.82泽泻 泽泻醇B单乙酸酯 35.58 0.81泽泻 23-乙酰泽泻醇B 32.52 0.82泽泻 16β-甲氧基泽泻醇B单乙酸酯 32.43 0.77泽泻 泽泻醇 B 36.76 0.82泽泻 泽泻醇 C 32.70 0.82泽泻 泽泻醇C单乙酸酯 33.06 0.83泽泻 1-亚油酸甘油单酯 37.18 0.30泽泻 [(1S，3R)-1-[(2R)-3，3- 二甲基氧烷 -2-yl]-3-[(5R，8S，9S，10S，11S，14R)-11- 羟基 -4，4，8，10，14-五甲基 -3- 氧代 -1，2，5，6，7，9，11，12，15，16- 十氢环戊并[a]菲-17-yl]butyl]醋酸盐35.58 0.81猪苓 (22e，24r)-麦角甾-6-en-3β，5alpha，6β-三醇 30.20 0.76猪苓 (22e，24r)7，22-麦角甾-二烯-3-酮 44.88 0.72猪苓 麦角甾-7，22-二烯-3β-醇 43.51 0.72猪苓 5alpha，8alpha-表二氧基-(22e，24r)-麦角甾-6，22-二烯-3β-醇 44.39 0.82猪苓 7，22-麦角甾二烯-3-酮 44.88 0.72猪苓 7，22-麦角甾二烯-3-酮 46.18 0.72猪苓 多孔菌甾酮 E 45.71 0.85猪苓 多孔菌甾酮 G 33.43 0.81
2.3 “单味药-活性成分-作用靶点”网络 根据上述分析所得的潜在作用靶点的预测结果，运用Cytoscape软件中的“merge”功能构建“五苓散单味药-活性成分-作用靶点”网络，见图2。网络总共包括48个节点(22个活性成分节点和26个靶点节点)和60条边，不同颜色的椭圆形代表不同的药物，每条边表示单味药活性成分与作用靶点之间的关系。同时构建“五苓散治疗骨关节炎的调控网络”，见图3。该网络中度值排名前5位的化合物分别是β-谷甾醇(beta-sitosterol)、儿茶素(catechin)、花旗松素(taxifolin)、3β-乙酰氧基苍术酮(3β-acetoxyatractylone)、常春藤皂苷(hederagenin)分别能与15，7，6，5，4个靶点连接，对治疗骨关节炎具有重要的意义，其基本信息见表2。
2.4 PPI网络及网络拓扑分析结果 借助 Bisogenet(Version 3.0.0)中 的 DIP，BIOGRID，HPRD，INTACT，MINT，BIND等各大数据库和Cytoscape软件来处理共有蛋白基因，初步得到一个2 738个节点、63 635条边的PPI网络图。再运用Cytoscape中的CytoNCA插件，以DC＞61进行1次筛选，得到一个648个节点、26 153条边的 PPI网络图，再以 Betweenness值 (BC)＞4 000 进行 2 次筛选，得到一个12个节点、40条边的PPI网络图，其过程见图4。度值越高，在整个网络中起到的作用越大，度值排名前5的蛋白基因为神经营养酪氨酸受体激酶1(Neurotrophic ReceptorTyrosine Kinase 1，NTRK1)、 雌 激 素 受体 1(Estrogen Receptor 1，ESR1)、肿瘤蛋白p53(Tumor Protein p53，Tp53)、泛素连接酶CUL3蛋白(Cullin3，CUL3)微小染色体维持蛋白2(Minichromosome Maintenance Complex Component 2，MCM2)。这些度值较大的蛋白基因在整个网络中起着关键的作用。同时这些蛋白基因所对应的靶点可能是五苓散治疗骨关节炎的关键靶点。
表2 ｜关键化合物的基本信息Table 2 ｜The basic information of the key compounds化学成分 度值 来源 结构β-谷甾醇(beta-sitosterol) 15 桂枝儿茶素(catechin) 7 桂枝花旗松素(taxifolin) 6 桂枝3β-乙酰氧基苍术酮(3β-acetoxyatractylone) 5 白术常春藤皂苷(hederagenin) 4 茯苓
2.5 GO和KEGG富集分析结果 GO富集分析得到GO条目(P＜ 0.05)共709个，包括生物过程(BP)649个、细胞成分(CC)22个和分子功能(MF)38个，并依据P值排序选取每个模块排名前20位的条目，见图5。最终结果提示，生物学过程的相关条目与对辐射的反应(response to radiation )、对雌二醇的反应(response to estradiol)、对抗生素的反应(response to antibiotic)、对脂多糖的反应(response to lipopolysaccharide )、对细菌来源分子的反应(response to molecule of bacterial origin )等有关；细胞成分的相关条目主要包括膜筏(membrane raft)、膜微区(membrane microdomain)、膜区(membrane region)、线粒体外膜(mitochondrial outer membrane)、细胞器外膜(organelle outer membrane)等；分子功能的相关条目主要包括核受体活性 (nuclear receptor activity)、转录因子活性 (transcription factor activity)、半胱氨酸型内肽酶活性与细胞凋亡的关系 (cysteine-type endopeptidase activity involved in apoptotic process)、 类 固 醇激素受体活性(steroid hormone receptor activity)等。
KEGG通路富集分析共得到91条通路，根据P值排序，选取前20位的条目，见图6。关键蛋白模块的相关通路主要涉及AGE-RAGE信号通路(AGE-RAGE signaling pathway in diabetic complications)、雌激素信号通路(Estrogen signaling pathway)、白细胞介素17信号通路(IL-17 signaling pathway)、肿瘤坏死因子信号通路(TNF signaling pathway)和p53信号通路(p53 signaling pathway)等。
2.6 KEGG关系网络构建 结合五苓散成分-靶点网络、核心靶蛋白PPI网络，将上述可能参与到骨关节炎过程的前20条KEGG信号通路进行数据及属性文件编排，通过Cytoscape 3.7.2软件得到“KEGG关系网络图”，见图7。结果可以明显看出五苓散的各有效成分是通过多靶点、多信号通路共同起到治疗骨关节炎的作用，其中天冬氨酸蛋白水解酶3(CASP3)、Bax蛋白、BcL2蛋白、人源全长重组蛋白(RELA)和天冬氨酸蛋白水解酶8(CAPS8)为主要靶点蛋白。
2.7 分子对接结果 一般认为配体与受体结合越稳定结合能越低，发生作用的可能性越大，以结合能≤-5.0 kJ/mol作为标准，对预测的关键活性成分和关键靶点进行分子对接验证。分子对接结果显示五苓散的关键活性成分与关键靶点对接的亲和力远小于-5.0 kJ/mol，由此表明预测的五苓散的关键活性成分与治疗骨关节炎的关键作用靶点具有较好的结合活性，证明该研究的预测较为可靠，结果见表3。选取亲和力最低的儿茶素和NTRK1进行分子对接展示，见图8。
3 讨论 Discussion
文章通过对五苓散进行靶点筛选发现单味药-活性成分-作用靶点网络显示β-谷甾醇、儿茶素、花旗松素、3β-乙酰氧基苍术酮、常春藤皂苷能匹配到较多的靶点，可能是五苓散治疗骨关节炎的关键成分。3β-乙酰氧基苍术酮作为白术提取物的一种挥发油，虽然目前缺乏3β-乙酰氧基苍术酮相关研究报道，但是有研究发现另一种白术提取挥发油类物质苍术酮可以通过刺激PI3K/AKT通路影响下游mTOR信号因子[12]，调节细胞凋亡的发生，而作为同类物质的3β-乙酰氧基苍术酮可能也存在相关机制对细胞凋亡与炎症反应有一定影响。有实验研究发现儿茶素可以抑制炎症因子白细胞介素6与一氧化氮在大鼠细胞内的表达，阻碍炎症在骨关节内的发生，起到保护关节软骨细胞的作用[13]。有研究发现β-谷甾醇可通过刺激成骨细胞骨保护素mRNA表达促进成骨作用，同时下调破骨细胞分化因子的分泌抑制破骨作用达到护骨效果，并且作用卵巢颗粒细胞，上调雌二醇的表达[14]，而雌激素水平的下降是造成骨关节炎的关键原因之一已被多个研究证实[15]。然而花旗松素虽作为一种植物雌激素，却不能刺激雌激素受体通路作用于成骨细胞，但是它不仅可提升成骨细胞活性因子的表达，还可抑制肿瘤坏死因子α作用下核转录因子κB的核移位，阻止核转录因子κB对骨细胞的破坏作用[16]，并且随着对花旗松素的进一步研究，发现它抑制人骨髓巨噬细胞核转录因子κB配体诱导的破骨作用受体，抑制体内脂多糖诱导的骨丢失[17]。目前对常春藤皂苷的研究较少，但有研究发现它可减少骨表面破骨细胞数量和活性[18]，并认为它是通过减弱RANKL诱导的细胞内活性氧产生和降低MAPK信号通路(ERK和p38)激活，抑制下游c-Fos和NFATc1的诱导的骨丢失途径实现的，但是尚未有更多的研究对其进行佐证。
表3｜五苓散的关键活性成分与关键靶点对接的结合能Table 3 ｜Binding energies of the key active components and key targets ofWulingpowder?
蛋白互作网络表明ESR1，Tp53，CUL3和MCM2可能是五苓散治疗骨关节炎的关键靶点。骨关节炎主要由软骨细胞凋亡、炎性因子及相关水解蛋白破坏软骨细胞引起的，而上述4个蛋白对这一病因有着各自不同的关联。在雌激素内分泌系统中，雌激素受体1是信号转导途径中必不可少的介质，以往体外和体内实验都报道了激素受体1是对软骨细胞的潜在保护作用，其主要通过刺激生长因子如转化生长因子β、胰岛素样生长因子1的表达促进软骨细胞增殖，并抑制相关炎性细胞因子如白细胞介素1，6和肿瘤坏死因子α的分泌保护细胞不被破坏，它还可调节影响关键软骨形成相关蛋白基质金属蛋白酶(MMP)的表达[19]。另有研究认为骨关节炎的发生可能与激素受体1是基因的突变有关，激素受体1是基因单核苷酸多态性的变异是造成骨关节炎的因素之一，但通过循证医学方式发现此结论存在争议，需要更进一步的实验去证实[20]。而p53是肿瘤抑制基因的原型，最初发现被视作为与SV-40病毒的T抗原共沉淀的53-kD核磷蛋白，它对细胞周期，DNA修复，诱导细胞凋亡具有调控作用[21]。近年来有研究发现p53对炎性因子白细胞介素6具有抑制作用，可减轻软骨损伤，并抑制基质金属蛋白酶的活性，使其侵袭活性丧失，达到治疗骨关节炎的作用[22-23]。CUL1是Cullin基因组中成员之一，是泛素连接酶E3复合体的核心成分，CUL1蛋白所组装的复合体SCF(Skp-Cull-F box蛋白)是一个多亚基蛋白质复合物，具有典型的泛素连接酶活性，可介导许多细胞周期如癌细胞相关蛋白的凋亡降解过程[24]。有研究发现在骨肉瘤细胞中，当Cullin 1高表达时基质金属蛋白酶9的表达也随之升高[25]，反之则降低，这说明Cullin1对基质金属蛋白酶9有调控作用，而基质金属蛋白酶蛋白优势细胞降解蛋白，会加速细胞的凋亡，而这一结果也可能发生在骨关节炎软骨细胞中。MCM2是微小染色体维持蛋白家族成员之一，具有高度保守性，它在DNA复制过程中作为解旋酶对细胞增殖起着重要的作用，MCM蛋白表达的下调会引起细胞增殖能力的丢失，并且MCM mRNA 水平会跟随细胞周期变化而变化[26]。
图 1｜活性化合物-疾病交集基因韦恩图Figure 1 ｜Venn diagram of active compounds-disease intersecting genes图注：蓝色部分为骨关节炎疾病靶点，共2 828个，红色部分为五苓散活性成分靶点，共52个，两者交集部分为五苓散治疗骨关节炎的潜在作用靶点，共26个
图 2｜五苓散单味药-活性成分-作用靶点网络Figure 2 ｜Single drug-active components-targets network ofWulingpowder图注：图中黄色部分代表骨关节炎相关疾病靶点，紫蓝色代表桂枝活性成分，绿色代表茯苓活性成分，天蓝色代表猪苓成分，紫红色代表泽泻活性成分，灰色代表白术活性成分，红色代表上述5味药的共同活性成分。线条代表活性成分可作用于骨关节炎靶点，线越密集说明对活性成分对骨关节炎相关靶点作用越多，作用越突出
图 3｜五苓散治疗骨关节炎调控网络Figure 3 ｜Regulatory network ofWulingpowder in the treatment of osteoarthritis图注：图中黄色节点表示骨关节炎，蓝色节点表示与骨关节炎密切相关的作用靶点，红色节点表示五苓散，绿色节点表示与五苓散中的各种有效成分
图 4｜PPI 网络拓扑分析Figure 4 ｜Network topology analysis of protein-protein interaction network图注：图中红色的点代表关键基因；蓝色点代表关键基因拓扑分析的结果；第一张图黄色的点代表设置条件DC＞61后的结果；最后一张图代表DC＞61和 BC＞4 000筛选后的结果
图 5｜GO 功能富集分析Figure 5 ｜GO functional enrichment analysis图注：图中以富集分数为横坐标，P值大小以柱的颜色表示，颜色由红到蓝，P值由小变大，饼状图代表生物过程(BP)、细胞组成(CC)和分子功能(MF)各条目在所有GO富集分析条目的比例
图 6｜ KEGG 信号通路可视化结果Figure 6 ｜ Visualization results of KEGG signaling pathway图注：图中横坐标为该通路所富集的基因占人体总基因的比例，气泡大小与通路上富集的基因数呈正比，气泡颜色越红则代表P值越小，其富集程度越显著，可信度越高
图 7｜KEGG 关系网络图Figure 7 ｜ KEGG network diagram图注：内圈黄色长方形代表基因，外圈三角代表KEGG信号通路，如果KEGG信号通路和基因之间有连线则说明属于这个通路，图形越大，代表连接节点越多
图 8｜儿茶素和NTRK1的分子对接图Figure 8 ｜ Molecular docking diagram of catechin and NTRK1
从富集分析结果可知，在五苓散治疗骨关节炎的过程中主要涉及的信号通路为AGE-RAGE信号通路、白细胞介素17信号通路、肿瘤坏死因子信号通路、p53信号通路，这与软骨细胞的增殖修复、细胞凋亡与炎症机制的形成有紧密联系。晚期糖基化终末产物(AGEs)作为一种异质性生物大分子，是由葡萄糖、氨基酸及脂肪通过非酶糖基化反应形成的不可逆的终末期产物，被认为是加快人体衰老及关节慢性退化的原因之一。近些年的研究表明当晚期糖基化终末产物与其特异性受体(RAGE)结合会促进氧化应激及炎性反应[27]，导致机体产生病理改变(包括细胞的增生与坏死)，黄文舟等[28]的动物实验中发现晚期糖基化终末产物通过氧化应激诱导大鼠软骨细胞损伤，它主要上调了软骨细胞凋亡相关分子Bax/Bcl-2 和caspase-3水平的，并增加体内基质金属蛋白酶3，13表达及Ⅱ型胶原含量，而基质金属蛋白酶(MMPs)会引起胶原蛋白的分解，破坏胶原蛋白结构的稳定性，导致软骨弹性下降，增加关节磨损可能。而赵骏等[29]通过观察晚期糖基化终末产物对人关节软骨细胞炎性因子、基质降解因子表达及核因子κB活化水平的影响，发现当提升晚期糖基化终末产物浓度时，可诱导人关节软骨细胞肿瘤坏死因子α、基质金属蛋白酶13、白细胞介素1α，1β以及软骨细胞核中核因子κB p65表达增多，该类因子的表达升高也可破坏软骨细胞外基质，然而上述因子的表达也可通过上调核因子κB的表达实现，因此认为骨关节炎的发生可AGE/RAGE/核转录因子κB有关。
白细胞介素17作为一种由CD4+T细胞分化的T辅助17细胞产生的促炎性细胞因子，与骨关节炎发生有密切关联。有研究发现骨关节炎患者膝关节白细胞介素17含量明显高于正常膝关节，白细胞介素17与骨关节炎的疼痛水平和严重程度成正比[30-31]，BAI等[32]比较骨关节炎患者与正常人血清和滑膜液中白细胞介素17的水平与骨关节炎严重程度与上述研究结果一致。目前已证实白细胞介素17有5个受体(白细胞介素17RA-RD和SEF)和6个成员(白细胞介素17A-F)，当白细胞介素17与其受体结合会刺激成纤维细胞和上皮细胞释放白细胞介素1β、肿瘤坏死因子α、白细胞介素6，8，并与上述炎症因子协同作用，上调促进软骨基质降解的多种因子表达，如基质金属蛋白酶和一氧化氮等，加快软骨细胞损伤[33]。因此激活白细胞介素17信号通路可能会加快软骨细胞的凋亡，加重关节炎的严重程度。
骨关节炎的发生机制可能与炎症因子如C-反应蛋白、白细胞介素和肿瘤坏死因子家族成员在关节中高表达有关。其中肿瘤坏死因子α作为一种促炎细胞因子，可由软骨细胞，单核细胞，成骨细胞和滑膜组织产生的，并诱导软骨细胞和滑膜细胞自身产生其他细胞因子，参与炎症、免疫调节、增殖和凋亡等多种生理病理过程[34]。肿瘤坏死因子α的作用机制较为广泛，它可与其受体结合抑制细胞外基质的主要成分连接蛋白和Ⅱ型胶原的合成；它还可刺激促炎和促代谢介质的合成，如基质金属蛋白酶1，3，13，白细胞介素6，8等，这些介质具有引起关节炎症、诱导发热、降解软骨基质蛋白的能力，此外当阻断肿瘤坏死因子α受体可抑制软骨中一氧化氮的产生，而一氧化氮他是滑膜组织中潜在的凋亡介质[35-36]，这说明肿瘤坏死因子α信号通路引发骨关节炎的作用具有多机制特点，而寻找调控该信号通路治疗骨关节炎的发病是治疗新思路。
p53抑癌基因是细胞凋亡的重要调节因子，当细胞DNA受损时，p53蛋白在G1/S界面停止生长，直到损伤修复，如果损伤严重，p53提高在组织内的表达诱导细胞凋亡，有研究证实骨关节炎的发生与软骨细胞p53信号高表达有关[37]。但目前p53的细胞凋亡相关信号通路尚未明确，有学者认为p53作用机制主要为上调凋亡因子Bax表达[38]，加速细胞内线粒体凋亡，而p53信号表达的上调又受多种因素影响如细胞内活性氧、一氧化氮/p38/核转录因子κB信号、MAPK信号通路等。这说明调节p53在骨关节炎中的表达需要抑制上述因素相关蛋白的基础上，降低Bax的表达。
综上所述，文章利用网络药理学及分子对接技术对五苓散治疗骨关节炎的潜在活性成分、靶点、信号通路和作用机制等进行了分析，结果表明五苓散内的活性化合物β-谷甾醇、儿茶素、花旗松素和常春藤皂苷可能通过影响雌激素受体1、Tp53、CUL3和MCM2相关蛋白靶点的表达，调控AGE-RAGE、白细胞介素17、肿瘤坏死因子和p53信号通路对软骨细胞的增殖、修复、凋亡及细胞炎症反应的影响，发挥减缓软骨细胞凋亡、加强软骨细胞增殖修复及减轻炎症反应的作用，这为五苓散治疗骨关节炎的临床应用提供参考依据，同时为进一步深入研究和实验设计的优化提供新的思路和方法。
作者贡献：文章成文由第一作者完成。文章责任人为第二作者。第三作者参与审校及终审定稿。第四、五、六作者完成数据核对。
经费支持：该文章接受了“国家自然科学基金项目(，)、广西高校青年教师基础能力提升项目(2019KY0352)、广西中医药大学2019年校级科研课题(2019QN027)、广西自然科学基金课题(2015GXNSFAA)、广西中医药大学一流学科建设开放课题(2019XK029)、广西中管局自筹经费项目(GZZC)、广西卫生厅重点课题(S-05)、2016年全国名老中医黄有荣传承工作室建设项目、广西中医药大学岐黄工程高层次人才团队培育项目()”的资助。所有作者声明，经费支持没有影响文章观点和对研究数据客观结果的统计分析及其报道。
利益冲突：文章的全部作者声明，在课题研究和文章撰写过程，没有因其岗位角色影响文章观点和对数据结果的报道，不存在利益冲突。
文章查重：文章出版前已经过专业反剽窃文献检测系统进行3次查重。
文章外审：文章经小同行外审专家双盲外审，同行评议认为文章符合期刊发稿宗旨。
文章版权：文章出版前杂志已与全体作者授权人签署了版权相关协议。
开放获取声明：这是一篇开放获取文章,根据《知识共享许可协议》“署名-非商业性使用-相同方式共享4.0”条款,在合理引用的情况下,允许他人以非商业性目的基于原文内容编辑、调整和扩展,同时允许任何用户阅读、下载、拷贝、传递、打印、检索、超级链接该文献,并为之建立索引,用作软件的输入数据或其它任何合法用途。
4 参考文献 References
[1]贺自克,王上增,沈锦涛.基于Hedgehog通路探讨补肾活血方对实验性膝骨性关节炎兔退变软骨的保护作用[J].中国病理生理杂志,2020,36(4):700-706.
[2]钱凯,杜彦仪,韩隆胤,等.基于网络药理学探讨五苓散治疗类风湿关节炎的作用机制[J].中国实验方剂学杂志,2019,25(19):138-146.
[3]ABRAMOFF B, CALDERA FE. Osteoarthritis:pathology, diagnosis, and treatment options. Med Clin North Am. 2020;104(2):293-311.
[4]HAO DA C, XIAO PG. Network pharmacology:a Rosetta Stone for traditional Chinese medicine. Drug Dev Res. 2014;75(5):299-312.
[5]张彦琼,李梢.网络药理学与中医药现代研究的若干进展[J].中国药理学与毒理学杂志,2015,29(6):883-892.
[6]段贤春,黄石,彭代银,等.网络药理学在中药复方研究中的应用[J].中国药理学通报,2020,36(3):303-308.
[7]曲凡,杨锡燕.五苓散合黄连解毒汤加减治疗难治性高血压的临床效果[J].中国老年学杂志,2020,40(4):692-694.
[8]汝锦龙.中药系统药理学数据库和分析平台的构建和应用[D].咸阳:西北农林科技大学,2015.
[9]胡聪,黎红维,吴琳静,等.基于网络药理学的苍耳散治疗鼻炎作用机制研究[J/OL].中华中医药学刊:1-19[2020-06-14].
[10]浦冬青,刘政,梁佳玮,等.基于网络药理学及分子对接探讨热毒宁注射液抑制新型冠状病毒肺炎炎症风暴机制[J/OL].中药药理与临床:1-17[2020-06-14].
[11]廖姿,张家琼,付先芸,等.基于网络药理学的“大黄-桃仁”配伍治疗子宫腺肌病的作用机制研究[J/OL].中国中药杂志:1-10[2020-06-14].
[12]耿玮,梁巍,叶智斌,等.苍术酮对结直肠癌细胞HT29凋亡的机制[J].中成药,2018,40(4):937-940.
[13]孙东东,程海波,沈卫星,等.补肾活血方中5种黄酮类成分的体外抗炎活性[J].中国实验方剂学杂志,2015,21(7):137-141.
[14]曾莉萍,徐贤柱,饶华,等.杜仲叶β-谷甾醇对成骨细胞和卵巢颗粒细胞的影响[J].时珍国医国药,2012,23(5):1051-1053.
[15]柳海平,周明旺,李盛华,等.骨关节炎软骨细胞增殖相关信号通路研究进展[J].中国矫形外科杂志,2018,26(23):2163-2166.
[16]WANG YJ, ZHANG HQ, HAN HL, et al. Taxifolin enhances osteogenic differentiation of human bone marrow mesenchymal stem cells partially via NF-κB pathway. Biochem Biophys Res ;490(1):36-43.
[17]ZHANG HQ, WANG YJ, YANG GT, et al. Taxifolin Inhibits Receptor Activator of NF-κB Ligand-Induced Osteoclastogenesis of Human Bone Marrow-Derived Macrophages in vitro and Prevents Lipopolysaccharide-Induced Bone Loss in vivo. Pharmacology. 2019;103(1-2):101-109.
[18]TIAN K, SU Y, DING J, et al. Hederagenin protects mice against ovariectomy-induced bone loss by inhibiting RANKL-induced osteoclastogenesis and bone Sci. 2020;244:.
[19]李健,姚玉龙,张磊磊,等.雌激素在软骨中作用机制研究进展[J].中国矫形外科杂志,2016,24(23):2156-2159.
[20]MA H, WU W, YANG X, et al. Genetic effects of common polymorphisms in estrogen receptor alpha gene on osteoarthritis: a meta-analysis. Int J Clin Exp Med. 2015;8(8):-.
[21]HAYANO S, KOMATSU Y, PAN H, et BMP signaling in the neural crest inhibits nasal cartilage morphogenesis by inducing p53-mediated 2015;142(7):1357-1367.
[22]张国梁,王跃文,普日布苏荣,等.p53/miR-502-5p/TRAF2通路对骨关节炎软骨细胞损伤的影响[J].中国病理生理杂志,2018,34(7):1283-1290.
[23]庞琳娜,徐连那,董红宇.p53基因在类风湿关节炎骨代谢中的作用[J].中国骨质疏松杂志,2014,20(3):327-330, 337.
[24]ZHAO Y, XIONG X, SUN Y. Cullin-RING Ligase 5: Functional characterization and its role in human cancers. Semin Cancer Biol. 2020.Doi:10.1016/
[25]程庆华. Cullin-1调节骨肉瘤细胞增殖及转移机制研究[D].南京:南京医科大学,2018.
[26]郭会芹,赵焕,李楠,等.MCM2免疫细胞化学检测在宫颈癌筛查中的意义[J].中华肿瘤防治杂志,2009,16(16):1213-1216.
[27]刘畅,王娜,李玉坤,等.晚期糖基化终末产物对骨组织细胞代谢的影响[J].中国组织工程研究,2020,24(27):4376-4382.
[28]黄文舟,王丽丽,殷嫦嫦,等.晚期糖基化终末产物通过氧化应激诱导大鼠软骨细胞损伤[J].中国病理生理杂志,2016,32(11):2036-2042.
[29]赵峻,马翅,杨娜娜,等.晚期糖基化终末产物致软骨细胞炎性反应的潜在机制[J].中国组织工程研究,2018,22(24):3786-3791.
[30]LIU Y, PENG H, MENG Z, et al. Correlation of IL-17 level in synovia and severity of knee osteoarthritis. Med Sci Monit. 2015;21:1732-1736.
[31]CHEN B, DENG Y, TAN Y, et al. Association between severity of knee osteoarthritis and serum and synovial fluid interleukin 17 concentrations. J Int Med Res. 2014;42(1):138-144.
[32]BAI Y, GAO S, LIU Y, et al. Correlation between Interleukin-17 gene polymorphism and osteoarthritis susceptibility in Han Chinese population. BMC Med ;20(1):20.
[33]WANG Z, ZHENG C, ZHONG Y, et Can Induce Osteoarthritis in Rabbit Knee Joints Similar to Hulth’s Method. Biomed Res Int. 2017;2017:.
[34]KHOTIB J, UTAMI NW, GANI MA, et al.The change of proinflammatory cytokine tumor necrosis factor α level in the use of meloxicam in rat model of Basic Clin Physiol Pharmacol. 2019. Doi:10.1515/jbcpp-2019-0331.
[35]ABDEL GALIL SM, EZZELDIN N, FAWZY F,et al. The single-nucleotide polymorphism(SNP) of tumor necrosis factor α -308G/A gene is associated with early-onset primary knee osteoarthritis in an Egyptian female population. Clin Rheumatol. 2017;36(11):2525-2530.
[36]KOU S, WU Y. meta-analysis of tumor necrosis factor alpha -308 polymorphism and knee osteoarthritis risk. BMC Musculoskelet Disord. 2014;15:373.
[37]ZHU X, YANG S, LIN W, et al. Roles of cell cyle regulators cyclin d1, cdk4, and p53 in knee osteoarthritis. Genet Test Mol Biomarkers. 2016;20(9):529-534.
[38]周绪昌,邹军,吴伟.p53调控骨关节炎软骨细胞凋亡[J].中国生物化学与分子生物学报,2019,35(3):280-285.
文章来源：分子催化 网址: http://fzch.400nongye.com/lunwen/itemid-45865.shtml
上一篇： 暂无
下一篇： 有机化工论文_高效催化水解纤维素的固体酸特性及其影响因素
