李院魏　高驥　李雙印

【摘 要】目的：聯(lián)合GEO數(shù)據(jù)庫使用網(wǎng)絡(luò)藥理學(xué)探討四妙勇安湯治療類風(fēng)濕關(guān)節(jié)炎的主要成分和可能的作用機制。方法：通過GEO芯片數(shù)據(jù)挖掘OMIM數(shù)據(jù)庫和Genecards數(shù)據(jù)庫，檢索收集類風(fēng)濕關(guān)節(jié)炎靶點;通過TCMSP數(shù)據(jù)庫，以生物利用度 > 30%，類藥性 > 0.18為篩選條件，收集四妙勇安湯的活性成分及對應(yīng)靶點。將四妙勇安湯靶點與類風(fēng)濕關(guān)節(jié)炎靶點進行匹配獲取靶點基因，使用Cytoscape 3.7.2軟件構(gòu)建蛋白-蛋白相互作用網(wǎng)絡(luò)，并進行拓?fù)浞治龊Y選核心靶點。使用R軟件進行GO分析和KEGG通路分析。結(jié)果：四妙勇安湯有效成分共126個（如山柰酚、槲皮素、豆甾醇、β-谷甾醇、谷甾醇等），對應(yīng)靶點2044個;類風(fēng)濕關(guān)節(jié)炎靶點4360個;中藥-疾病靶點匹配獲得潛在靶點165個。GO分析生物過程獲得通路2329個，細(xì)胞組分獲得通路89個，分子功能獲得通路163個;KEGG通路分析共獲得通路160個。結(jié)論：本研究為更深層次探究四妙勇安湯及所含四妙勇安湯藥物成分方劑治療類風(fēng)濕關(guān)節(jié)炎作用機制奠定了基礎(chǔ)。

【關(guān)鍵詞】 類風(fēng)濕關(guān)節(jié)炎;四妙勇安湯;網(wǎng)絡(luò)藥理學(xué);信號通路;靶點;蛋白-蛋白相互作用網(wǎng)絡(luò)

On the Mechanism of Simiao Yong'an Tang（四妙勇安湯）in Treating Rheumatoid Arthritis Using Network

Pharmacology Combined with GEO Database

LI Yuan-wei，GAO Ji，LI Shuang-yin

【ABSTRACT】Objective：To explore the main components and possible mechanism of Simiao Yong'an Tang（四妙勇安湯）in the treatment of rheumatoid arthritis（RA）using network pharmacology

combined with GEO database.Methods：GEO chip database was used to mine the OMIM database and Genecards database，which were used to collect rheumatoid arthritis targets.Through TCMSP database，the bioactive components and corresponding targets of Simiao Yong'an Tang were collected under the screening conditions of bioavailability being over 30% and generic property over 0.18.The Simiao Yong'an Decoction target was matched with rheumatoid arthritis target to obtain target gene.The PPI network was constructed using software Cytoscape 3.7.2，and the core target was screened by topology analysis.Software R was used to analyze GO biological function and KEGG pathway enrichment.Results：There were 126 effective components in Simiao Yong'an Tang（such as kaempferol，quercetin，stigmasterol，β-sitosterol，sitosterol etc.）.There were 2044 corresponding targets，4360 rheumatoid arthritis targets and 165 potential targets obtained by Chinese medicine-disease target matching.Go analysis showed that there were 2329 pathways in biological process，89 pathways in cell component and 163 pathways in molecular function.A total of 160 pathways were obtained by KEGG pathway analysis.Conclusion：This study laid a foundation for further exploring the mechanism of Simiao Yong'an Tang in the treatment of rheumatoid arthritis.

【Keywords】 rheumatoid arthritis;

Simiao Yong'an Tang（四妙勇安湯）;network pharmacology;signal pathway;targets;PPI network

類風(fēng)濕關(guān)節(jié)炎（rheumatoid arthritis，RA）是一種以侵蝕關(guān)節(jié)為主要表現(xiàn)的全身性自身免疫性疾病。主要累及關(guān)節(jié)，也會出現(xiàn)類風(fēng)濕性結(jié)節(jié)、肺部受累血管炎，以及全身合并癥等關(guān)節(jié)外表現(xiàn)[1]。RA屬中醫(yī)學(xué)“尪痹”“歷節(jié)病”“痹病”等范疇。RA急性期紅腫熱痛與中醫(yī)的熱痹極其相似。中醫(yī)學(xué)認(rèn)為，毒熱之邪注于經(jīng)絡(luò)，流入關(guān)節(jié)，使氣機阻滯，血運不暢，導(dǎo)致紅腫熱痛。房定亞[2]提出運用四妙勇安湯治療RA急性期熱痹證。徐國山等[3]研究表明，在西藥基礎(chǔ)上聯(lián)合四妙勇安湯類方治療RA可增強療效并減少不良反應(yīng)。四妙勇安湯清熱活血，尤其是RA急性期瘀熱痹阻證，臨床療效顯著。

西醫(yī)臨床常用藥物包括非甾體抗炎藥、糖皮質(zhì)激素和改善病情抗風(fēng)濕藥（DMARDs）。傳統(tǒng)的DMARDs[如甲氨蝶呤、Janus激酶（JAK）抑制劑等]或生物來源的DMARDs（如腫瘤壞死因子抑制劑、共刺激調(diào)節(jié)劑、白細(xì)胞介素-6抑制劑等）[4]，雖然治療RA療效顯著，但存在一定的不良反應(yīng)，且JAK抑制劑等新型藥物價格高，患者難以承受。

四妙勇安湯聯(lián)合西藥治療RA效果顯著[3]，但其作用機制尚未明確。本文使用網(wǎng)絡(luò)藥理學(xué)方法挖掘現(xiàn)有數(shù)據(jù)庫對其作用機制進行探討，為四妙勇安湯及其類方治療RA的研究提供科學(xué)依據(jù)。

1 材料與方法

1.1 篩選RA對應(yīng)靶點 登錄美國國家生物技術(shù)信息中心（NCBI）選擇Gene Expression Omnibus（GEO）數(shù)據(jù)庫，以“rheumatoid arthritis”為關(guān)鍵詞搜索相關(guān)芯片，獲取編號為GSE93272的芯片數(shù)據(jù)原始文件和GPL570-55999的芯片基因注釋文件，此芯片數(shù)據(jù)中包含275個樣本，其中232個來自RA患者（試驗組），43個來自健康人群（對照組）。使用R軟件（https：//www.r-project.org/，version 3.6.2）對數(shù)據(jù)進行處理，以|logFC| > 0.5，P < 0.05為過濾條件獲取差異基因，即RA相關(guān)基因。按照P值篩選前20個差異基因繪制基因熱圖。登錄在線人類孟德爾遺傳數(shù)據(jù)庫（OMIM，https：//omim.org/）和Genecard（https：//www.genecards.org/），以“rheumatoid arthritis”為關(guān)鍵詞檢索獲取RA相關(guān)基因。最終整理匯總獲取RA有效靶點。

1.2 四妙勇安湯有效成分及對應(yīng)靶點的篩選 通過中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫及分析平臺（TCMSP，http：//tcmspw.com/tcmsp.php）對四妙勇安湯中的金銀花、玄參、當(dāng)歸和甘草的有效成分及靶點進行篩選，以生物利用度（OB）> 30%，類藥性（DL） > 0.18為篩選條件，獲得4味中藥的有效成分及對應(yīng)靶點。登錄蛋白質(zhì)數(shù)據(jù)庫（Uniprot Database，https：//www.uniprot.org/）對所獲得的基因進行標(biāo)準(zhǔn)化命名。

1.3 中藥復(fù)方調(diào)控網(wǎng)絡(luò)及PPI網(wǎng)絡(luò)的構(gòu)建 將獲取的RA靶點基因與四妙勇安湯有效靶點進行數(shù)據(jù)對比，預(yù)測中藥對RA有效的作用靶點，通過Cytoscape 3.7.2軟件進行可視化。利用Cytoscape中的Bisogenet插件進行蛋白-蛋白相互作用（PPI）網(wǎng)絡(luò)構(gòu)建，通過CytoNCA插件進行網(wǎng)絡(luò)拓?fù)浞治?，以度中心性（DC） > 61，介數(shù)中心性（BC） >?600進行過濾，得到PPI網(wǎng)絡(luò)。

1.4 GO和KEGG分析 通過R軟件中的Cluster Profiler[5]對四妙勇安湯和RA交集基因進行基因本體論（GO）[生物過程（BP）、分子功能（MF）、細(xì)胞組分（CC）]分析，京都基因和基因組百科全書（KEGG）通路分析，對P值進行FDR矯正，設(shè)置閾值校正P < 0.05，使用R軟件進行繪圖。通過Cytoscape軟件構(gòu)建KEGG與核心基因關(guān)系網(wǎng)絡(luò)，根據(jù)連接節(jié)點的數(shù)量多少設(shè)置每個節(jié)點的大小，做可視化處理。

2 結(jié) 果

2.1 RA對應(yīng)靶點 通過GEO芯片數(shù)據(jù)分析獲得對照組和RA差異基因，以|logFC| > 0.5，P < 0.05為過濾條件，共獲得238個差異基因。使用R軟件繪制熱圖，見圖1。OMIM數(shù)據(jù)庫獲取RA靶點48個。Genecard數(shù)據(jù)庫獲取RA靶點4157個。整理后共獲得4443個靶點基因，刪除重復(fù)基因83個，保留4360個靶點基因。

2.2 四妙勇安湯有效成分及對應(yīng)靶點 登錄TCMSP數(shù)據(jù)庫依次檢索4味中藥，以O(shè)B > 30%，DL > 0.18為檢索條件獲取金銀花有效成分23個，玄參有效成分9個，當(dāng)歸有效成分2個，甘草有效成分92個，其中山柰酚、槲皮素為金銀花、甘草的共有成分，豆甾醇為金銀花、當(dāng)歸的共有成分，β-谷甾醇為金銀花、當(dāng)歸、玄參的共有成分，谷甾醇為玄參、甘草的共有成分，及其相對應(yīng)的靶點共2355個。登錄蛋白質(zhì)數(shù)據(jù)庫對所獲得的靶點基因進行標(biāo)準(zhǔn)化命名，去除重復(fù)靶點和無法匹配的靶點，最后獲得2044個靶點。

2.3 中藥復(fù)方調(diào)控網(wǎng)絡(luò)的構(gòu)建及CytoNAC網(wǎng)絡(luò)拓?fù)浞治?將獲取的RA相關(guān)靶點和四妙勇安湯靶點進行匹配，獲得四妙勇安湯治療RA的潛在基因165個。通過Cytoscape軟件進行四妙勇安湯調(diào)控RA網(wǎng)絡(luò)構(gòu)建，見圖2。

通過Cytoscape軟件中的Bisogenet插件，選擇數(shù)據(jù)庫包括DIP、BIOGRID、HPRD、INTACT、MINI、BIND，進行PPI網(wǎng)絡(luò)構(gòu)建。使用Cytoscape軟件中的CytoNCA插件進行拓?fù)浞治?，獲得7644個節(jié)點。以DC > 61過濾獲得1549個節(jié)點;繼續(xù)以BC > 600依次過濾，最終獲得733個節(jié)點。依據(jù)BC值排名靠前的有神經(jīng)營養(yǎng)因子受體絡(luò)氨酸激酶1（NTRK1）、腫瘤蛋白質(zhì)P53（TP53）、Cullin 3（CUL3）、雌激素受體1（ESR1）、泛素C（UBC）、微染色體維護復(fù)合體組分（MCM2）、細(xì)胞周期蛋白依賴性激酶2（CDK2）、表皮生長因子受體（EGFR）等。

2.4 GO分析 使用R軟件對四妙勇安湯和RA的交集基因進行GO分析，包括BP、CC、MF，以校正P < 0.05為過濾條件，BP獲得通路2329個，CC獲得通路89個，MF獲得通路163個，篩選前20個做圖，見圖3。

2.5 KEGG通路分析 使用R軟件對四妙勇安湯與RA交集基因進行KEGG通路分析，設(shè)置校正P < 0.05過濾條件，獲得通路160個。使用R軟件對富集程度較高的前20個通路進行可視化，見圖4。通過Cytoscape軟件構(gòu)建KEGG與核心基因關(guān)系網(wǎng)絡(luò)，根據(jù)連接節(jié)點的數(shù)量多少設(shè)置每個節(jié)點的大小，做可視化處理，見圖5。其中DC較高的基因有絲裂原活化蛋白激酶1（MAPK1，17個）、MAPK3（17個）、絲氨酸/蘇氨酸蛋白激酶（AKT1，17個）、癌基因同源物A（RELA，15個）、螺旋環(huán)螺旋結(jié)構(gòu)域擴散激酶（CHUK，14個）、IκBκB（14個）、MAPK8（13個）、MAPK10（13個）、B淋巴細(xì)胞瘤-2相關(guān)X蛋白（BAX，13個）、RAF1（12個）、核轉(zhuǎn)錄因子-κB抑制蛋白（NFκBIA，12個）等。

3 討 論

RA是一種自身免疫性炎性風(fēng)濕病，影響許多組織和器官，引起慢性滑膜炎癥，最終導(dǎo)致關(guān)節(jié)破壞、慢性殘疾和預(yù)期壽命縮短[6]。近年來，國際新型藥物腫瘤壞死因子（TNF）抑制劑、JAK通路抑制劑等的研發(fā)和投入臨床，雖取得了良好的療效，但是由于藥物使用的局限性（如腫瘤、感染等患者不宜使用）和使用后的不良反應(yīng)，RA的臨床治療方法仍需不斷探索。中醫(yī)學(xué)認(rèn)為，RA的病因包括正虛與邪實兩個方面，正虛是發(fā)病的根本原因，主要為營衛(wèi)不和、氣血不足、脾胃虛弱、肝腎虧虛等內(nèi)在因素;邪實是發(fā)病的重要條件，多指風(fēng)、寒、濕、熱等六淫之邪[7]。RA急性期的炎癥反應(yīng)如紅腫熱痛、皮下結(jié)節(jié)、紅色皮疹、身熱、煩躁、口干、溲黃正是熱盛成毒，濕郁化濁，二邪交結(jié)，阻滯氣血所致。四妙勇安湯對于毒熱痹阻、灼傷血脈的病機具有極強的針對性，且該方免疫調(diào)節(jié)機制又為其用于治療RA提供了充分的理由[8]。

已有研究顯示，山柰酚是一種天然黃酮類化合物，具有抗氧化、抗癌、抗炎等多種治療作用[9]。由于其抗炎特性，可能被用于治療許多急慢性炎癥引起的疾病，包括椎間盤退變、結(jié)腸炎、絕經(jīng)后骨質(zhì)流失和急性肺損傷。此外，它對癌癥、肝損傷、肥胖和糖尿病均有益處，能抑制血管內(nèi)皮炎癥，保護腦神經(jīng)和心臟功能，并可能用于治療纖維增生性疾病，包括增生性瘢痕[10]。槲皮素是已知抗氧化活性的自由基清除和抗過敏特性的特點是刺激免疫系統(tǒng)，抗病毒活性，抑制組胺釋放，減少促炎細(xì)胞因子的產(chǎn)生，并抑制白細(xì)胞介素（IL）-4的產(chǎn)生;改善Th1/Th2平衡，抑制抗原特異性免疫球蛋白E抗體形成;對脂氧合酶、嗜酸性粒細(xì)胞和過氧化物酶等有抑制作用，對炎癥介質(zhì)有抑制作用[11-12]。

β-谷甾醇具有抗氧化劑、促血管生成、免疫調(diào)節(jié)、抗糖尿病、抗炎、抗癌和鎮(zhèn)痛等作用，且無明顯的不良反應(yīng)[13]。

在人類培養(yǎng)細(xì)胞中進行的全基因組基因敲除篩選顯示，CUL3對細(xì)胞活力有重要意義[14]，這表明CUL3甚至在內(nèi)皮細(xì)胞的維持中也起著至關(guān)重要的作用。另一篇文獻(xiàn)中也提到，CUL3的E3泛素連接酶（CRL）可以參與廣泛的細(xì)胞內(nèi)事件和調(diào)控機制，包括受體交叉信號傳遞、基因表達(dá)的表觀遺傳調(diào)控和細(xì)胞-細(xì)胞連接構(gòu)成[15]。雌激素受體1（ESR1）與骨質(zhì)疏松易感性有關(guān)[16-21]，這可能與RA后期導(dǎo)致關(guān)節(jié)骨質(zhì)破環(huán)有關(guān)。通過抑制CDK1、CDK2、CDK4、CDK7，可以阻止細(xì)胞周期的G1或G2期的細(xì)胞周期進程，減緩細(xì)胞炎癥，增殖和多種慢性疾病的發(fā)生[22]。

GO和KEGG通路分析富集結(jié)果顯示，氧化應(yīng)激過程與炎癥和炎癥反應(yīng)的聯(lián)系密切[23]。已有研究表明，TNF和IL-17等炎癥通路與RA發(fā)病密切相關(guān)[18-19]。RA急性期的炎癥反應(yīng)明顯，四妙勇安湯可能通過影響氧化應(yīng)激生物過程干預(yù)急性期的炎癥反應(yīng)，緩解癥狀。氧化還原反應(yīng)（ROS）分子還參與信號轉(zhuǎn)導(dǎo)途徑和依賴氧化還原調(diào)控一系列生物活性，類如骨髓來源的抑制細(xì)胞（MDSCs），這些異源髓樣細(xì)胞受ROS控制，但它們也使用ROS實現(xiàn)抑制功能[24]。諸如慢性炎癥、感染和癌癥之類的病理狀況會誘導(dǎo)MDSCs，MDSCs由未成熟髓細(xì)胞的異質(zhì)群體組成[25-26]。流體剪切應(yīng)力和動脈粥樣硬化通路，有文獻(xiàn)表明，生理層流切變可以促進細(xì)胞在流動方向上的伸長和定向，從而抑制增殖，刺激抗炎基因的表達(dá)并抑制炎性途徑的表達(dá)[27]。現(xiàn)有臨床試驗表明，TNF、IL-6，以及可能的粒細(xì)胞-單核細(xì)胞集落刺激因子對RA的發(fā)病過程至關(guān)重要[28]。

綜上所述，四妙勇安湯可能是通過多種活性成分，如山柰酚、槲皮素、β-谷甾醇等，調(diào)控RA相關(guān)基因CUL3、ESAR1、CDK，通過氧化應(yīng)激過程、ROS等生物過程，進一步調(diào)控多種信號通路，如年齡-憤怒信號通路、前列腺癌、流體剪切應(yīng)力和動脈粥樣硬化通路、TNF、IL-17信號、JAK-STAT信號傳送途徑等通路抑制炎癥反應(yīng)和干預(yù)免疫調(diào)節(jié)，綜合作用機體治療RA。本文通過網(wǎng)絡(luò)藥理學(xué)對四妙勇安湯治療RA的作用機制進行探討，可為臨床應(yīng)用四妙勇安湯及其類方治療RA提供理論依據(jù)。

參考文獻(xiàn)

[1] SMOLEN JS，ALETAHA D，MCINNES IB.Rheumatoid arthritis[J].Lancet，2016，388（10055）：2023-2038.

[2] 房定亞，劉云峰.四妙勇安湯治療類風(fēng)濕關(guān)節(jié)炎[J].中醫(yī)函授通訊，1984（3）：134-135.

[3] 徐國山，俞曉芳，孔美珠，等.清熱活血法治療瘀熱型類風(fēng)濕關(guān)節(jié)炎臨床觀察[J].中華中醫(yī)藥雜志，2018，33（3）：1167-1170.

[4] BURMESTER GR，POPE JE.Novel treatment strategies in rheumatoid arthritis[J].Lancet，2017，389（10086）：2338-2348.

[5] YU G，WANG LG，HAN Y，et al.clusterProfiler：an R package for comparing biological themes among gene clusters[J].Omics，2012，16（5）：284-287.

[6] SONG X，LIN Q.Genomics，transcriptomics and proteomics to elucidate the pathogenesis of rheumatoid arthritis[J].Rheumatol Int，2017，37（8）：1257-1265.

[7] 吳閔，姚曉玲，姚血明，等.類風(fēng)濕關(guān)節(jié)炎中醫(yī)證候分型研究進展[J].風(fēng)濕病與關(guān)節(jié)炎，2018，7（11）：71-74.

[8] 李亞慧，劉珊.周彩云教授運用四妙勇安湯治療風(fēng)濕病驗案舉隅[J].風(fēng)濕病與關(guān)節(jié)炎，2017，6（3）：41-43.

[9] IMRAN M，RAUF A，SHAH ZA，et al.Chemo-preventive and therapeutic effect of the dietary flavonoid kaempferol：a comprehensive review[J].Phytother Res，2019，33（2）：263-275.

[10] REN J，LU Y，QIAN Y，et al.Recent progress regarding kaempferol for the treatment of various diseases[J].Exp Ther Med，2019，18（4）：2759-2776.

[11] MLCEK J，JURIKOVA T，SKROVANKOVA S，et al.Quercetin and its anti-allergic immune response[J].Molecules，2016，21（5）：623.

[12] LI Y，YAO J，HAN C，et al.Quercetin，inflammation and immunity[J].Nutrients，2016，8（3）：167.

[13] BIN SAYEED MS，KARIM SMR，SHARMIN T，et al.

Critical analysis on characterization，systemic effect，and therapeutic potential of beta-sitosterol：a plant-derived orphan phytosterol[J].Medicines（Basel，Switzerland），2016，3（4）：29.

[14] SINGER JD，GURIAN-WEST M，CLURMAN B，et al.

Cullin-3 targets cyclin E for ubiquitination and controls S phase in mammalian cells[J].Genes Dev，1999，13（18）：2375-2387.

[15] SAKAUE T，MAEKAWA M，NAKAYAMA H，et al.Prospect of divergent roles for the CUL3 system in vascular endothelial cell function and angiogenesis[J].J Biochem，2017，162（4）：237-245.

[16] SANO M，INOUE S，HOSOI T，et al.Association of?estrogen receptor dinucleotide repeat polymorphism with osteoporosis[J].Biochem Biophys Res Commun，1995，217（1）：378-383.

[17] KOBAYASHI S，INOUE S，HOSOI T，et al.Association of bone mineral density with polymorphism of the estrogen receptor gene[J].J Bone Miner Res，1996，11（3）：306-311.

[18] 馬武開，陸道敏，姚血明，等.NF-κB、IL-17與膠原誘導(dǎo)性關(guān)節(jié)炎大鼠關(guān)節(jié)滑膜細(xì)胞及血清表達(dá)的關(guān)系[J].風(fēng)濕病與關(guān)節(jié)炎，2015，4（2）：26-29.

[19] ALAM J，JANTAN I，BUKHARI SNA.Rheumatoid arthritis：recent advances on its etiology，role of cytokines and pharmacotherapy[J].Biomedicine Pharmacother，2017，92（3）：615-633.

[20] WILLING M，SOWERS M，ARON D，et al.Bone mineral density and its change in white women：estrogen and vitamin D receptor genotypes and their interaction[J].J Bone Miner Res，1998，13（4）：695-705.

[21] ALBAGHA OM，MCGUIGAN FE，REID DM，et al.Estrogen receptor alpha gene polymorphisms and bone mineral density：haplotype analysis in women from the United Kingdom[J].J Bone Miner Res，2001，16（1）：128-134.

[22] SRIKUMAR T，PADMANABHAN J.Potential use of flavopiridol in treatment of chronic diseases[J].Adv Exp Med Biol，2016，929（1）：209-228.

[23] SIES H.Oxidative stress：a concept in redox biology and medicine[J].Redox Biol，2015，4（1）：180-183.

[24] OHL K，TENBROCK K.Reactive oxygen species as regulators of mdsc-mediated immune suppression[J].Front Immunol，2018，30（9）：2499.

[25] CORZO CA，COTTER MJ，CHENG P，et al.Mechanism regulating reactive oxygen species in tumor-induced myeloid-derived suppressor cells[J].J Immunol，2009，182（9）：5693-5701.

[26] GABRILOVICH DI，NAGARAJ S.Myeloid-derived suppressor cells as regulators of the immune system[J].Nat Rev Immunol，2009，9（3）：162-174.

[27] BAEYENS N，BANDYOPADHYAY C，COON BG，et al.

Endothelial fluid shear stress sensing in vascular health and disease[J].J Clin Invest，2016，126（3）：821-828.

[28] FELDMANN M，MAINI SRN.Role of cytokines in rheumatoid arthritis：an education in pathophysiology and therapeutics[J].Immunol Rev，2008，223（1）：7-19.

收稿日期：2021-02-25;修回日期：2021-03-28

作者單位：鄭州市第三人民醫(yī)院，河南 鄭州 450000

通信作者：李院魏 河南省鄭州市管城區(qū)南順城街136號，964564470@qq.com，（0371）66905254