基于網(wǎng)絡(luò)藥理學(xué)探討四妙丸治療骨關(guān)節(jié)炎的作用機(jī)制

吳越 劉維

摘要 目的：采用網(wǎng)絡(luò)平臺(tái)藥理學(xué)的研究方法，對(duì)四妙丸（SMW）用于治療骨關(guān)節(jié)炎（OA）的藥理作用以及機(jī)制問(wèn)題進(jìn)行了研究和探討。方法：主要通過(guò)在TCMSP網(wǎng)絡(luò)平臺(tái)對(duì)四妙丸進(jìn)行活性藥物成分的篩選、在DrugBank數(shù)據(jù)庫(kù)的基礎(chǔ)上進(jìn)行了靶標(biāo)的預(yù)測(cè)，應(yīng)用Cytoscape 3.7.1軟件構(gòu)建了化合物-靶點(diǎn)數(shù)據(jù)分析網(wǎng)絡(luò);在TTD及DisGenet進(jìn)行數(shù)據(jù)庫(kù)的檢索，獲取了疾病的靶點(diǎn);通過(guò)構(gòu)建韋恩圖，取得了疾病與活性藥物作用靶點(diǎn)的交集，得到了四妙丸治療OA的作用靶點(diǎn);并對(duì)四妙丸進(jìn)行了GO功能注釋功能通路注釋和KEGG通路富集分析;構(gòu)建了四妙丸中藥-成分-靶點(diǎn)-信號(hào)通路網(wǎng)絡(luò)。結(jié)果：分析研究顯示共挖掘到了四妙丸中64種共同作用成分及163個(gè)潛在共同作用靶點(diǎn)，OA疾病989個(gè)成分及相關(guān)共同作用靶點(diǎn)，得到四妙丸-疾病共同靶點(diǎn)18個(gè)，共同信號(hào)作用的靶點(diǎn)主要富集于34個(gè)生物過(guò)程和26條信號(hào)通路上。結(jié)論：四妙丸中藥可能主要是通過(guò)調(diào)節(jié)PTGS2、TNF、AR等共同靶點(diǎn)，調(diào)控IL-17信號(hào)通路、TNF信號(hào)通路、松弛素、破骨細(xì)胞分化、松弛素等方式來(lái)抑制機(jī)體炎性反應(yīng)，進(jìn)一步調(diào)節(jié)機(jī)體免疫功能來(lái)治療骨關(guān)節(jié)炎。

關(guān)鍵詞 四妙丸;網(wǎng)絡(luò)藥理學(xué);骨關(guān)節(jié)炎;靶點(diǎn);信號(hào)通路

Mechanism of Simiao Pills in Treatment of Osteoarthritis Based on Network Pharmacology

WU Yue1，2，LIU Wei1，2

（1 First Teaching Hospital of Tianjin University of Traditional Chinese Medicine，Tianjin 300193，China; 2 Tianjin Key Laboratory of Traditional Chinese Medicine Formula and Syndrome Transformation Research，Tianjin 300193，China）

Abstract Objective：To study the mechanism of Simiao Pills in the treatment of osteoarthritis（OA），based on network pharmacology.Methods：By screening the active ingredients of Simiao Pills on TCMSP platform，the target in the database of drugbank was predicted.Cytoscape 3.7.1 software was used to construct the compound-target data analysis network; the database was searched in TTD and DisGenet to obtain the target of the disease; the Venn diagram was constructed to obtain the intersection of the target of the disease and the active drug.The target of Simiao Pill in the treatment of OA was obtained; GO function annotation functional pathway annotation and KEGG pathway enrichment analysis were performed on Simiao Pill; a network of Simiao Pill Chinese medicine-component-target-signaling pathway was constructed.Results：The research showed that 64 components and 163 potential targets in Simiao Pills，989 related targets in OA disease，and 18 common targets in Simiao Pills disease were obtained，which were mainly concentrated in 34 biological processes and 26 signal pathways.Conclusion：Simiao Pills may be used for the treatment of osteoarthritis by regulating PTGS2，TNF，AR and other targets，regulating IL-17 signaling pathway，TNF signaling pathway，relaxin，osteoclast differentiation and other pathways，so as to inhibit inflammatory response etc.

Keywords Simiao Pills; Network pharmacology; Osteoarthritis; Target; Signal pathway

中圖分類(lèi)號(hào)：R285文獻(xiàn)標(biāo)識(shí)碼：Adoi：10.3969/j.issn.1673-7202.2020.23.012

骨關(guān)節(jié)炎（Osteoarthritis，OA）主要是以慢性軟骨關(guān)節(jié)炎和軟骨炎的局部病灶及慢性關(guān)節(jié)肌肉病變、軟骨關(guān)節(jié)下骨肥厚和慢性關(guān)節(jié)局部邊緣骨贅等為主要臨床發(fā)病病理特征的慢性肢體關(guān)節(jié)疾病，常伴有局部肌肉疼痛或關(guān)節(jié)無(wú)力、韌帶松弛、關(guān)節(jié)對(duì)線不良、輕度滑膜炎和關(guān)節(jié)半月板變性。2010年中國(guó)50歲以上的骨關(guān)節(jié)炎人口患病率達(dá)到25.3%[1]，在60歲以上的骨關(guān)節(jié)炎人群中大約每年有5 000萬(wàn)人患有骨關(guān)節(jié)炎[2]，預(yù)計(jì)隨著亞洲人口年齡的進(jìn)一步增長(zhǎng)，亞洲骨關(guān)節(jié)炎患病率將進(jìn)一步的上升，有相關(guān)研究人員預(yù)測(cè)，65歲以上人口的比例在未來(lái)30年將逐步翻番，從2008年的6.8%進(jìn)一步上升甚至預(yù)計(jì)可能達(dá)到2040年的16.2%，使得社會(huì)的衛(wèi)生保健以及教育事業(yè)的發(fā)展面臨更嚴(yán)峻的挑戰(zhàn)[3]。本病主要屬于中醫(yī)“痹癥”“骨痹”范疇?，F(xiàn)代的研究結(jié)果表明，本病的發(fā)生與代謝異常、機(jī)械功能損傷、肥胖、及基因遺傳等一系列因素與有關(guān)[4-5]。我國(guó)現(xiàn)代醫(yī)學(xué)對(duì)于OA的診斷和治療主要目的是用于抗炎和調(diào)節(jié)骨代謝，常應(yīng)用非甾體抗炎藥、氨基葡萄糖、軟骨素、局部注射皮質(zhì)類(lèi)固醇或透明質(zhì)酸、外用制劑及手術(shù)等進(jìn)行骨關(guān)節(jié)的治療[6]。此外，關(guān)于骨關(guān)節(jié)炎的治療，現(xiàn)代中醫(yī)藥不僅有傳統(tǒng)辨證論治的特點(diǎn)，還有療效獨(dú)特且不良反應(yīng)較小的優(yōu)勢(shì)。但目前四妙丸對(duì)于OA的具體藥理作用和機(jī)制尚不明確，故嚴(yán)重限制了其臨床有效應(yīng)用和推廣。

四妙丸源自《成方便讀》，有清熱利濕的作用和功效，方中以黃柏為君，苦寒沉降，善除下焦?jié)駸帷７街幸陨n術(shù)燥濕健脾和胃為臣。薏苡仁祛濕清熱瀉濁為佐。牛膝活血，亦可通調(diào)經(jīng)絡(luò)，善引藥下行，為佐、使藥。在OA中證屬于濕熱痹阻，臨床中可根據(jù)病情靈活應(yīng)用[7-8]。本研究基于網(wǎng)絡(luò)藥理學(xué)方法探討四妙丸治療OA的主要有效化合物、作用靶點(diǎn)、作用通路，從而進(jìn)一步明確四妙丸治療OA的具體作用機(jī)制，為OA的臨床治療提供可能的思路。

1 材料與方法

1.1 四妙丸的主要成分及作用靶點(diǎn)的篩選 通過(guò)中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫(kù)（TCMSP）（http：//lsp.nwu.edu.cn/tcmsp.php）[9]檢索出四妙丸中蒼術(shù)、黃柏、牛膝、薏苡仁等4種中藥的靶蛋白化合物。對(duì)靶蛋白化合物成分進(jìn)行了篩選，其篩選的標(biāo)準(zhǔn)：中藥口服生物利用度（OB）篩選的閾值為OB≥30，化合物類(lèi)藥性（DL）篩選的閾值為DL≥0.18。通過(guò)DrugBank（https：//www.drugbank.ca/）[10]數(shù)據(jù)庫(kù)對(duì)化合物進(jìn)行靶點(diǎn)預(yù)測(cè)，剔除數(shù)據(jù)庫(kù)不能預(yù)測(cè)到靶點(diǎn)的化合物，在Uniprot（https：//www.uniprot.org/）[11]數(shù)據(jù)庫(kù)中將預(yù)測(cè)出的作用靶蛋白名轉(zhuǎn)換為作用靶點(diǎn)的簡(jiǎn)稱(chēng)，進(jìn)行標(biāo)準(zhǔn)化的處理。

1.2 建立四妙丸化合物-靶點(diǎn)網(wǎng)絡(luò) 對(duì)預(yù)測(cè)出的化合物名稱(chēng)進(jìn)行編碼以便于分析和作圖。將化合物及其靶點(diǎn)數(shù)據(jù)輸入至Cytoscape 3.7.1[12]中，從而建立了中藥化合物—靶點(diǎn)網(wǎng)絡(luò);然后使用network analyzer工具進(jìn)行網(wǎng)絡(luò)分析，通過(guò)自由度（Degree）判定節(jié)點(diǎn)的重要性，越高則代表越重要。

1.3 獲取骨關(guān)節(jié)炎疾病靶點(diǎn) 在Therapeutic Target Database（TTD）（https：//db.idrblab.org/ttd/）[13]以及DisGenet（http：//www.disgenet.org/web/DisGeNET/menu/search）[14]數(shù)據(jù)庫(kù)中檢索關(guān)鍵詞“Osteoarthritis”，來(lái)獲取骨關(guān)節(jié)炎疾病靶點(diǎn)。在Uniprot數(shù)據(jù)庫(kù)中得到靶點(diǎn)的簡(jiǎn)稱(chēng)，從而對(duì)其進(jìn)行了標(biāo)準(zhǔn)化，并構(gòu)建骨關(guān)節(jié)炎疾病靶點(diǎn)的數(shù)據(jù)集。

1.4 篩選重要靶點(diǎn)及建立蛋白相互作用網(wǎng)絡(luò) 運(yùn)用Venny2.1作圖平臺(tái)（http：//bioinfogp.cnb.csic.es/tools/venny/），通過(guò)映射四妙丸有效化合物的靶點(diǎn)和骨關(guān)節(jié)炎病靶來(lái)構(gòu)建韋恩圖。其共同靶點(diǎn)可能是四妙丸治療骨關(guān)節(jié)炎的重要靶點(diǎn)。在STRING平臺(tái)（https：//string-db.org/）[15]導(dǎo)入所獲取的共同靶點(diǎn)，構(gòu)建蛋白-蛋白相互作用網(wǎng)絡(luò)，生物種屬選擇“Homo sapiens”，網(wǎng)絡(luò)參數(shù)設(shè)置為“molecular action”并將最低相互作用的閾值參數(shù)設(shè)置為“medium confidence”（>0.400），其他參數(shù)的設(shè)置保持默認(rèn)。

1.5 重要靶點(diǎn)基因功能注釋和通路富集分析 應(yīng)用Cytoscape 3.7.1軟件中的ClueGo插件[16]對(duì)共同靶點(diǎn)作GO功能注釋和KEGG通路富集分析，取P≤0.05的結(jié)果進(jìn)行可視化分析，從而得到四妙丸治療骨關(guān)節(jié)炎可能的生物過(guò)程和信號(hào)通路。最后使用merge功能構(gòu)建四妙丸化合物-靶點(diǎn)-信號(hào)通路網(wǎng)絡(luò)。

2 結(jié)果

2.1 建立四妙丸化合物-靶點(diǎn)網(wǎng)絡(luò) 本研究得到四妙丸四味中藥的86種化合物，以及909個(gè)可能的相關(guān)靶點(diǎn)，靶點(diǎn)去重后有163個(gè);其中有22個(gè)化合物沒(méi)有找到相關(guān)靶點(diǎn)，故未參與分析。最終在四妙丸中得到64種化合物，蒼術(shù)中有4種，黃柏中有25種，牛膝中有17種，薏苡仁中有18種。見(jiàn)表1。該網(wǎng)絡(luò)中含有227個(gè)節(jié)點(diǎn)和908個(gè)相互關(guān)系（見(jiàn)圖1）?！白杂啥取迸琶叭幕衔锓謩e是槲皮素（Quercetin）（自由度值為75）、山柰酚（Kaempferol）（自由度值為44）、異延胡索單酚堿（Isocorypalmine）（自由度值為36）。

2.2 獲取共同靶點(diǎn)及蛋白相互作用網(wǎng)絡(luò)構(gòu)建 將獲取的骨關(guān)節(jié)炎989個(gè)疾病靶點(diǎn)與163個(gè)四妙丸化合物的作用靶點(diǎn)錄入韋恩圖在線繪制平臺(tái)，二者取完交集后獲得18個(gè)疾病-化合物共同靶點(diǎn)，該共同靶點(diǎn)在四妙丸治療骨關(guān)節(jié)炎中具有重要意義。見(jiàn)圖2。蛋白-蛋白相互作用網(wǎng)絡(luò)共獲得18個(gè)相互作用節(jié)點(diǎn)和108條相互作用關(guān)系。見(jiàn)圖3。

2.3 GO功能注釋 本研究通過(guò)對(duì)18個(gè)共同靶點(diǎn)進(jìn)行GO功能注釋?zhuān)瑥亩玫紾O分析網(wǎng)絡(luò)圖。見(jiàn)圖4。其共同靶點(diǎn)主要富集在參與蛋白質(zhì)去乙?；恼{(diào)控（Regulation of Protein Deacetylation），細(xì)胞遷移參與新生血管生成（Cell Migration Involved in Sprouting Angiogenesis），急性炎性反應(yīng)的調(diào)節(jié)（Regulation of Acute Inflammatory Response），一氧化氮代謝過(guò)程（Nitric Oxide Metabolic Process），血管平滑肌細(xì)胞增殖（Vascular Smooth Muscle Cell Proliferation），膠原分解的過(guò)程（Collagen Catabolic Process），等34個(gè)生物過(guò)程。這表明四妙丸對(duì)骨關(guān)節(jié)炎的治療可能通過(guò)細(xì)胞因子分泌調(diào)節(jié)、抗氧化應(yīng)激、蛋白降解、血管生成、組織重塑等。

2.4 KEGG通路富集分析 本研究對(duì)18個(gè)共同靶點(diǎn)進(jìn)行KEGG通路分析，然后在omicshare（https：//www.omicshare.com/）在線平臺(tái)上對(duì)結(jié)果進(jìn)行了可視化分析，得到KEGG分析氣泡圖。見(jiàn)圖5。共同靶點(diǎn)主要富集于TNF信號(hào)通路、破骨細(xì)胞分化、IL-17信號(hào)通路、弓形蟲(chóng)病、恰加斯?。绹?guó)錐蟲(chóng)病）、利什曼病、流體剪切應(yīng)力與動(dòng)脈粥樣硬化、糖尿病并發(fā)癥的AGE-RAGE信號(hào)通路、松弛素信號(hào)通路、FcεRI信號(hào)通路等26條通路。

2.5 化合物-靶點(diǎn)-通路網(wǎng)絡(luò)的構(gòu)建 四妙丸中藥化合物-共同靶點(diǎn)-信號(hào)通路網(wǎng)絡(luò)關(guān)系圖包含88個(gè)節(jié)點(diǎn)，235個(gè)相互作用關(guān)系。見(jiàn)圖6。根據(jù)自由度高低可得出較為重要的化合物、靶點(diǎn)及通路。見(jiàn)表2。

3 討論

四妙丸有清熱利濕的藥理作用和治療功效，能夠有效治療濕熱痹阻證骨關(guān)節(jié)炎。根據(jù)現(xiàn)代基礎(chǔ)研究表明，蒼術(shù)具有抗炎、抗過(guò)敏等作用[17]，黃柏具有調(diào)節(jié)機(jī)體自身免疫、抗炎、抗細(xì)胞凋亡等作用[18-19]，牛膝具有抗炎、擴(kuò)張下肢血管、抗骨質(zhì)疏松等作用[20]，薏苡仁具有抗凝血、抗炎等作用[21]。

本研究發(fā)現(xiàn)在中藥-主要成分-靶點(diǎn)-信號(hào)通路網(wǎng)絡(luò)中槲皮素（Quercetin）、漢黃芩素（Wogonin）、山柰酚（Kaempferol）、巴馬汀（Palmatine）、黃連堿（Coptisine）、黃芩素（Baicalein）等的自由度最高。其中槲皮素不僅對(duì)關(guān)節(jié)軟骨基質(zhì)具有優(yōu)先的親和力，還表明它可以通過(guò)抑制和減少活性氧種類(lèi)（ROS）和下游因子（如MMP）來(lái)有效保護(hù)關(guān)節(jié)軟骨。由于槲皮素本身是一種很強(qiáng)的抗氧化劑，它可以有效地通過(guò)抑制ROS的產(chǎn)生和NOX4的表達(dá)來(lái)提供潛在的軟骨保護(hù)作用[22]。另外還有一些研究結(jié)果表明槲皮素具有與塞來(lái)昔布相當(dāng)?shù)目寡趸饔?，可增加SOD的含量，下調(diào)軟骨基質(zhì)金屬蛋白酶13（MMP-13）的表達(dá)并可以通過(guò)上調(diào)軟骨細(xì)胞中TIMP-1的基因表達(dá)，減弱機(jī)體對(duì)氧化應(yīng)激的反應(yīng)，減輕關(guān)節(jié)軟骨細(xì)胞外基質(zhì)降解，從而有效改善骨關(guān)節(jié)炎表現(xiàn)[23]。另外還有一些研究結(jié)果指出漢黃芩素有效地抑制了兔關(guān)節(jié)軟骨基質(zhì)細(xì)胞中潛在的IL-1β誘導(dǎo)的基質(zhì)金屬蛋白酶1（MMP-1）、MMP-13和adamts-4的基因表達(dá)。同時(shí)它還恢復(fù)了IL-1β抑制的Ⅱ型膠原的基因表達(dá)。因此通過(guò)調(diào)節(jié)骨關(guān)節(jié)炎關(guān)節(jié)軟骨破壞過(guò)程中不同蛋白酶的基因表達(dá)，以及在轉(zhuǎn)錄水平上促進(jìn)II型膠原的基因轉(zhuǎn)錄表達(dá)，可有效支持漢黃芩素對(duì)軟骨的保護(hù)作用[24]。漢黃芩素可通過(guò)插入機(jī)制較好地表現(xiàn)出與軟骨細(xì)胞的DNA結(jié)合能力，其在細(xì)胞核內(nèi)可以通過(guò)抑制IL-1β介導(dǎo)的ROS誘導(dǎo)、DNA的碎片化、外源性和內(nèi)在性誘導(dǎo)的凋亡途徑，體現(xiàn)出對(duì)骨關(guān)節(jié)炎軟骨細(xì)胞的保護(hù)作用[25]。研究表明，山柰酚能顯著地抑制炎性反應(yīng)遞質(zhì)如誘導(dǎo)型一氧化氮合酶和環(huán)氧化酶2的白細(xì)胞介素（IL）1β誘導(dǎo)蛋白表達(dá)。并且山柰酚還抑制了由IL-1β誘導(dǎo)的常見(jiàn)基質(zhì)降解酶，從而消除了膠原II的降解。山柰酚的抗炎作用是通過(guò)抑制絲裂原活化蛋白激酶相關(guān)的細(xì)胞外信號(hào)調(diào)節(jié)激酶和p38信號(hào)通路來(lái)介導(dǎo)的[26]。山柰酚的可以顯著地阻止脂多糖引發(fā)的軟骨細(xì)胞凋亡和炎性反應(yīng)，其能顯著下調(diào)MIR-146a，抑制了去整合素的表達(dá)，從而充分發(fā)揮對(duì)軟骨的保護(hù)作用[27]。有研究已經(jīng)證實(shí)了巴馬汀對(duì)IL-1β刺激誘導(dǎo)的兔軟骨細(xì)胞炎性反應(yīng)具有保護(hù)作用，其對(duì)MMP的抑制作用可能與Wnt/β-連環(huán)蛋白和Hedgehog信號(hào)通路的抑制有關(guān)[28]。研究表明，黃連素通過(guò)抑制NF-kB信號(hào)通路從而抑制IL-1β誘導(dǎo)的炎性反應(yīng)，其能顯著抑制一氧化氮和前列腺素E2的產(chǎn)生，抑制IL-1β誘導(dǎo)的人OA軟骨細(xì)胞中誘導(dǎo)型一氧化氮合酶和環(huán)氧化酶2的表達(dá)。同時(shí)也顯著地抑制了MMP-3和MMP-13在IL-1β刺激的人OA軟骨細(xì)胞中的表達(dá)，并且顯著抑制了人OA軟骨細(xì)胞中IL-1β誘導(dǎo)的NF-kB活化[29]。黃芩素的抗軟骨細(xì)胞凋亡作用主要是由于抑制一氧化氮的產(chǎn)生和阻斷Caspase下游的信號(hào)通路，其作用能顯著地降低關(guān)節(jié)軟骨細(xì)胞中MMP-3和MMP-13的分泌。黃芩素對(duì)小鼠關(guān)節(jié)軟骨外植體的糖胺聚糖和Ⅱ型膠原沉積的恢復(fù)進(jìn)一步證實(shí)了其抗分解代謝的重要作用[30]。故四妙丸治療骨關(guān)節(jié)炎的過(guò)程中較為重要的成分可能是槲皮素（Quercetin）、漢黃芩素（Wogonin）、山柰酚（Kaempferol）、巴馬?。≒almatine）、黃連堿（Coptisine）、黃芩素（Baicalein）等。

在中藥-化合物-靶點(diǎn)-信號(hào)通路網(wǎng)絡(luò)中自由度最高的前6位靶點(diǎn)為前列腺素G/H合成酶2（PTGS2），絲裂原活化蛋白激酶1（MAPK1），絲裂原活化蛋白激酶14（MAPK14），絲裂原活化蛋白激酶8（MAPK8），腫瘤壞死因子（TNF），雄激素受體（AR）等。PTGS2是炎性反應(yīng)性PG生物蛋白合成的一個(gè)重要關(guān)鍵酶，其可以作為炎性反應(yīng)生物蛋白的標(biāo)志物[31]。PTGS2在OA軟骨和滑膜組織中過(guò)度的表達(dá)，損害蛋白聚糖的合成并增強(qiáng)了軟骨中對(duì)軟骨聚集蛋白聚糖和II型膠原的降解[32]。有更多的研究結(jié)果表明，阿片類(lèi)藥物誘導(dǎo)的對(duì)OA抗炎作用可能是由它們對(duì)MPGES和PTGS2表達(dá)的影響介導(dǎo)的[33]。TNF-α可與PGRN相互作用，誘導(dǎo)對(duì)軟骨細(xì)胞降解蛋白酶的表達(dá)和軟骨的丟失[34]。TNF-α可以通過(guò)激活NF-κB信號(hào)通路，上調(diào)各種MMP和ADAMT的水平，從而參與介導(dǎo)OA軟骨的降解[35]。研究結(jié)果表明，TNF/TNFR的配體寡聚誘導(dǎo)了銜接蛋白的募集，如Fas相關(guān)死亡結(jié)構(gòu)域蛋白和Tradd，這些銜接蛋白通過(guò)同源死亡結(jié)構(gòu)域相互作用與受體胞質(zhì)尾部結(jié)合，Procaspase-8通過(guò)同源死亡效應(yīng)器結(jié)構(gòu)域相互作用被銜接蛋白招募到復(fù)合物中，形成死亡誘導(dǎo)信號(hào)復(fù)合物，啟動(dòng)半胱天冬酶進(jìn)而激活細(xì)胞并執(zhí)行半胱天冬酶，其進(jìn)而切割受體死亡底物，最終可能導(dǎo)致受體細(xì)胞凋亡[36]，細(xì)胞凋亡被認(rèn)為是導(dǎo)致OA軟骨老化的一個(gè)關(guān)鍵因素[37]。睪酮的缺乏可誘導(dǎo)去勢(shì)家兔細(xì)胞凋亡，減少軟骨細(xì)胞的增殖和數(shù)量[38]。MA對(duì)雄性小鼠進(jìn)行了體內(nèi)研究（睪丸切除術(shù)與對(duì)照組），并確定完整小鼠表現(xiàn)出更嚴(yán)重的OA。此外，去勢(shì)雄性小鼠補(bǔ)充了二氫睪酮后，與正常小鼠比較，出現(xiàn)了嚴(yán)重的OA[39]。MAPK1、MAPK14、MAPK8卻少見(jiàn)相關(guān)報(bào)道。故四妙丸治療骨關(guān)節(jié)炎的作用可能主要是通過(guò)調(diào)節(jié)PTGS2、TNF、AR等靶點(diǎn)來(lái)實(shí)現(xiàn)。

根據(jù)富集分析結(jié)果，四妙丸治療骨關(guān)節(jié)炎主要涉及的通路可能是IL-17信號(hào)通路、流體剪切應(yīng)力與動(dòng)脈粥樣硬化、糖尿病并發(fā)癥的AGE-RAGE信號(hào)通路、松弛素信號(hào)通路、TNF信號(hào)通路、破骨細(xì)胞分化等。IL-17主要影響關(guān)節(jié)中的軟骨細(xì)胞和成纖維細(xì)胞樣滑膜細(xì)胞（FLS），其細(xì)胞表面表達(dá)IL-17受體[40]。IL-17已被證明能有效抑制關(guān)節(jié)中的軟骨細(xì)胞合成蛋白聚糖，并且能促進(jìn)軟骨細(xì)胞基質(zhì)金屬蛋白酶的產(chǎn)生[41-43]。此外，IL-17影響其他對(duì)軟骨有一定負(fù)面影響的細(xì)胞因子和化合物的分泌，如IL-1、TNF、IL-6、NO和PGE2[44-46]?；颊叩难搴突褐蠭L-17水平明顯升高，與OA病變的影像學(xué)和臨床表現(xiàn)表現(xiàn)呈正相關(guān)[47]。TNF-α是導(dǎo)致關(guān)節(jié)軟骨退變的關(guān)鍵炎性反應(yīng)因子之一，參與關(guān)節(jié)滑膜炎癥的反應(yīng)、軟骨細(xì)胞凋亡的啟動(dòng)、軟骨損傷和修復(fù)之間的不平衡以及最終加速關(guān)節(jié)軟骨退化[48]。松弛素已被證明能在各種情況下引起抗炎作用。在大鼠和豚鼠的炎性反應(yīng)模型中，松弛素通過(guò)一氧化氮的產(chǎn)生介導(dǎo)抑制肥大細(xì)胞釋放組胺[49]。松弛素還抑制中性粒細(xì)胞[50]和肥大細(xì)胞[51]等炎性反應(yīng)細(xì)胞流入受損器官。研究結(jié)果表明，CD14-細(xì)胞促進(jìn)破骨細(xì)胞的形成不僅發(fā)生在RA中，而且也發(fā)生在OA中，一種罕見(jiàn)的侵蝕性炎癥性O(shè)A可發(fā)生在手指間關(guān)節(jié)[52]。而流體剪切應(yīng)力與動(dòng)脈粥樣硬化，糖尿病并發(fā)癥的AGE-RAGE信號(hào)通路相關(guān)報(bào)道較少。此外研究發(fā)現(xiàn)四妙丸治療OA主要涉及蛋白質(zhì)去乙酰化的調(diào)控，細(xì)胞遷移參與新生血管生成，急性炎性反應(yīng)的調(diào)節(jié)，一氧化氮代謝過(guò)程，血管平滑肌細(xì)胞增殖，膠原分解過(guò)程等生物過(guò)程。由此可推測(cè)，四妙丸可能通過(guò)抑制炎性反應(yīng)、調(diào)節(jié)免疫功能來(lái)治療骨關(guān)節(jié)炎。

綜上所述，本研究通過(guò)網(wǎng)絡(luò)藥理學(xué)挖掘到四妙丸中64種化合物及163個(gè)靶點(diǎn)，骨關(guān)節(jié)炎疾病靶點(diǎn)989個(gè)，其中四妙丸-骨關(guān)節(jié)炎疾病共同靶點(diǎn)18個(gè)，對(duì)共同靶點(diǎn)進(jìn)行GO功能注釋及KEGG信號(hào)通路富集分析，得出四妙丸可能是通過(guò)調(diào)節(jié)PTGS2、TNF、AR等靶點(diǎn)，調(diào)控IL-17信號(hào)通路、TNF信號(hào)通路、松弛素、破骨細(xì)胞分化等通路，從而抑制炎性反應(yīng)、調(diào)節(jié)機(jī)體免疫功能來(lái)治療骨關(guān)節(jié)炎。由于各種平臺(tái)數(shù)據(jù)庫(kù)、相關(guān)軟件算法的限制，本研究的預(yù)測(cè)結(jié)果可為四妙丸的基礎(chǔ)研究提供線索，為了提高結(jié)果的可靠性，其具體作用機(jī)制還需進(jìn)一步的實(shí)驗(yàn)來(lái)驗(yàn)證。

參考文獻(xiàn)

[1]Tang X，Wang S，Zhan S，et al.The Prevalence of Symptomatic Knee Osteoarthritis in China：Results From the China Health and Retirement Longitudinal Study[J].Arthritis & Rheumatology，2016，68（3）：648-653.

[2]李寧華.中老年人群骨關(guān)節(jié)炎的流行病學(xué)特征[J].中國(guó)臨床康復(fù)，2005，9（38）：133-135.

[3]Fransen M，Bridgett L，March L，et al.The epidemiology of osteoarthritis in Asia[J].Int J Rheum Dis，2011，14（2）：113-121.

[4]Ding C，Martel-Pelletier J，Pelletier J P，et al.Meniscal tear as an osteoarthritis risk factor in a largely non-osteoarthritic cohort：a cross-sectional study[J].J Rheumatol，2007，34（4）：776-784.

[5]Hunter D J，Zhang Y Q，Niu J B，et al.The association of meniscal pathologic changes with cartilage loss in symptomatic knee osteoarthritis[J].Arthritis Rheum，2006，54（3）：795-801.

[6]Schmidt T W.Approach to Osteoarthritis Management for the Primary Care Provider[J].Primary Care：Clinics in Office Practice，2018，45（2）：361-378.

[7]肖玉蓮，馬合肖，潘樂(lè)意.四妙丸加中藥熏洗治療膝骨性關(guān)節(jié)炎急性發(fā)作臨床觀察[J].中醫(yī)學(xué)報(bào)，2011，26（6）：735-736.

[8]蔣淑蓮，陳保安，吳民祥，等.四妙丸治療膝骨性關(guān)節(jié)炎倂積液臨床療效觀察[J].內(nèi)蒙古中醫(yī)藥，2014，33（35）：24-25.

[9]Ru J，Li P，Wang J，et al.TCMSP：a database of systems pharmacology for drug discovery from herbal medicines[J].J Cheminform，2014，6：13.

[10]Wishart D S，Knox C，Guo A C，et al.DrugBank：a knowledgebase for drugs，drug actions and drug targets[J].Nucleic Acids Res，2008，36（Database issue）：D901-D906.

[11]UniProt：the universal protein knowledgebase[J].Nucleic Acids Res，2017，45（D1）：D158-D169.

[12]Franz M，Lopes C T，Huck G，et al.Cytoscape.js：a graph theory library for visualisation and analysis[J].Bioinformatics，2016，32（2）：309-311.

[13]Liu X，Zhu F，Ma X，et al.The Therapeutic Target Database：an internet resource for the primary targets of approved，clinical trial and experimental drugs[J].Expert Opin Ther Targets，2011，15（8）：903-912.

[14]Pinero J，Queralt-Rosinach N，Bravo A，et al.DisGeNET：a discovery platform for the dynamical exploration of human diseases and their genes[J].Database（Oxford），2015，2015：v28.

[15]von Mering C，Huynen M，Jaeggi D，et al.STRING：a database of predicted functional associations between proteins[J].Nucleic Acids Res，2003，31（1）：258-261.

[16]Bindea G，Mlecnik B，Hackl H，et al.ClueGO：a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks[J].Bioinformatics，2009，25（8）：1091-1093.

[17]Jun X，F(xiàn)u P，Lei Y，et al.Pharmacological effects of medicinal components of Atractylodes lancea（Thunb.）DC[J].Chin Med，2018，13：59.

[18]張冠英，董瑞娟，廉蓮.川黃柏、關(guān)黃柏的化學(xué)成分及藥理活性研究進(jìn)展[J].沈陽(yáng)藥科大學(xué)學(xué)報(bào)，2012，29（10）：812-821.

[19]李先寬，馮杉，鄭艷超，等.黃柏與關(guān)黃柏的化學(xué)成分及生物活性研究進(jìn)展[J].藥物評(píng)價(jià)研究，2019，42（5）：1033-1037.

[20]田碩，苗明三.牛膝的化學(xué)、藥理及應(yīng)用特點(diǎn)探討[J].中醫(yī)學(xué)報(bào)，2014，29（8）：1186-1188.

[21]吳巖，原永芳.薏苡仁的化學(xué)成分和藥理活性研究進(jìn)展[J].華西藥學(xué)雜志，2010，25（1）：111-113.

[22]Mok S，F(xiàn)u S，Cheuk Y，et al.Intra-Articular Delivery of Quercetin Using Thermosensitive Hydrogel Attenuate Cartilage Degradation in an Osteoarthritis Rat Model[J].CARTILAGE，2018：764875239.

[23]Wei B，Zhang Y，Tang L，et al.Protective effects of quercetin against inflammation and oxidative stress in a rabbit model of knee osteoarthritis[J].Drug Development Research，2019，80（3）：360-367.

[24]Park J S，Lee H J，Lee D Y，et al.Chondroprotective Effects of Wogonin in Experimental Models of Osteoarthritis in vitro and in vivo[J].Biomol Ther（Seoul），2015，23（5）：442-448.

[25]Khan N M，Ahmad I，Ansari M Y，et al.Wogonin，a natural flavonoid，intercalates with genomic DNA and exhibits protective effects in IL-1beta stimulated osteoarthritis chondrocytes[J].Chem Biol Interact，2017，274：13-23.

[26]Huang X，Pan Q，Mao Z，et al.Kaempferol inhibits interleukin-1β stimulated matrix metalloproteinases by suppressing the MAPK-associated ERK and P38 signaling pathways[J].Molecular Medicine Reports，2018，18（3）：2697-2704.

[27]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].International Immunopharmacology，2019，69：373-381.

[28]Zhou X，Lin X，Xiong Y，et al.Chondroprotective effects of palmatine on osteoarthritis in vivo and in vitro：A possible mechanism of inhibiting the Wnt/β-catenin and Hedgehog signaling pathways[J].International Immunopharmacology，2016，34：129-138.

[29]Zhou K，Hu L，Liao W，et al.Coptisine Prevented IL-β-Induced Expression of Inflammatory Mediators in Chondrocytes[J].Inflammation，2016，39（4）：1558-1565.

[30]Zhang X，Zhu Y，Chen X，et al.Baicalein ameliorates inflammatory-related apoptotic and catabolic phenotypes in human chondrocytes[J].International Immunopharmacology，2014，21（2）：301-308.

[31]Benelli R，Vene R，F(xiàn)errari N.Prostaglandin-endoperoxide synthase 2（cyclooxygenase-2），a complex target for colorectal cancer prevention and therapy[J].Transl Res，2018，196：42-61.

[32]Attur M，Al-Mussawir H E，Patel J，et al.Prostaglandin E2 exerts catabolic effects in osteoarthritis cartilage：evidence for signaling via the EP4 receptor[J].J Immunol，2008，181（7）：5082-5088.

[33]Schwarzbach S V，Melo C F，Xavier P，et al.Morphine，but not methadone，inhibits microsomal prostaglandin E synthase-1 and prostaglandin-endoperoxide synthase 2 in lipopolysaccharide-stimulated horse synoviocytes[J].Biochimie，2019，160：28-33.

[34]Zhao Y，Liu B，Tian Q，et al.Progranulin protects against osteoarthritis through interacting with TNF-α and β-Catenin signalling[J].Annals of the Rheumatic Diseases，2015，74（12）：2244-2253.

[35]Kapoor M，Martel-Pelletier J，Lajeunesse D，et al.Role of proinflammatory cytokines in the pathophysiology of osteoarthritis[J].Nat Rev Rheumatol，2011，7（1）：33-42.

[36]Aizawa T，Kon T，Einhorn T A，et al.Induction of apoptosis in chondrocytes by tumor necrosis factor-alpha[J].J Orthop Res，2001，19（5）：785-796.

[37]Pelletier J P，Jovanovic D V，Lascau-Coman V，et al.Selective inhibition of inducible nitric oxide synthase reduces progression of experimental osteoarthritis in vivo：possible link with the reduction in chondrocyte apoptosis and caspase 3 level[J].Arthritis Rheum，2000，43（6）：1290-1299.

[38]Irie T，Aizawa T，Kokubun S.The role of sex hormones in the kinetics of chondrocytes in the growth plate.A study in the rabbit[J].J Bone Joint Surg Br，2005，87（9）：1278-1284.

[39]Ma H L，Blanchet T J，Peluso D，et al.Osteoarthritis severity is sex dependent in a surgical mouse model[J].Osteoarthritis Cartilage，2007，15（6）：695-700.

[40]Honorati M C，Neri S，Cattini L，et al.Interleukin-17，a regulator of angiogenic factor release by synovial fibroblasts[J].Osteoarthritis Cartilage，2006，14（4）：345-352.

[41]Benderdour M，Tardif G，Pelletier J P，et al.Interleukin 17（IL-17）induces collagenase-3 production in human osteoarthritic chondrocytes via AP-1 dependent activation：differential activation of AP-1 members by IL-17 and IL-1beta[J].J Rheumatol，2002，29（6）：1262-1272.

[42]Martel-Pelletier J，Mineau F，Jovanovic D，et al.Mitogen-activated protein kinase and nuclear factor kappaB together regulate interleukin-17-induced nitric oxide production in human osteoarthritic chondrocytes：possible role of transactivating factor mitogen-activated protein kinase-activated proten kinase（MAPKAPK）[J].Arthritis Rheum，1999，42（11）：2399-2409.

[43]Lubberts E，Joosten L A，van de Loo F A，et al.Reduction of interleukin-17-induced inhibition of chondrocyte proteoglycan synthesis in intact murine articular cartilage by interleukin-4[J].Arthritis Rheum，2000，43（6）：1300-1306.

[44]Honorati M C，Bovara M，Cattini L，et al.Contribution of interleukin 17 to human cartilage degradation and synovial inflammation in osteoarthritis[J].Osteoarthritis Cartilage，2002，10（10）：799-807.

[45]Legrand A，F(xiàn)ermor B，F(xiàn)ink C，et al.Interleukin-1，tumor necrosis factor alpha，and interleukin-17 synergistically up-regulate nitric oxide and prostaglandin E2 production in explants of human osteoarthritic knee menisci[J].Arthritis Rheum，2001，44（9）：2078-2083.

[46]Attur M G，Patel R N，Abramson S B，et al.Interleukin-17 up-regulation of nitric oxide production in human osteoarthritis cartilage[J].Arthritis Rheum，1997，40（6）：1050-1053.

[47]Chen B，Deng Y，Tan Y，et al.Association between severity of knee osteoarthritis and serum and synovial fluid interleukin 17 concentrations[J].J Int Med Res，2014，42（1）：138-144.

[48]Urech D M，F(xiàn)eige U，Ewert S，et al.Anti-inflammatory and cartilage-protecting effects of an intra-articularly injected anti-TNF{alpha} single-chain Fv antibody（ESBA105）designed for local therapeutic use[J].Ann Rheum Dis，2010，69（2）：443-449.

[49]Bani D，Ballati L，Masini E，et al.Relaxin counteracts asthma-like reaction induced by inhaled antigen in sensitized guinea pigs[J].Endocrinology，1997，138（5）：1909-1915.

[50]Masini E，Nistri S，Vannacci A，et al.Relaxin inhibits the activation of human neutrophils：involvement of the nitric oxide pathway[J].Endocrinology，2004，145（3）：1106-1112.

[51]Masini E，Bani D，Bigazzi M，et al.Effects of relaxin on mast cells.In vitro and in vivo studies in rats and guinea pigs[J].J Clin Invest，1994，94（5）：1974-1980.

[52]Danks L，Sabokbar A，Gundle R，et al.Synovial macrophage-osteoclast differentiation in inflammatory arthritis[J].Ann Rheum Dis，2002，61（10）：916-921.

（2019-12-04收稿 責(zé)任編輯：蒼寧）