基于網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接技術(shù)研究甘草附子湯治療類風(fēng)濕性關(guān)節(jié)炎的作用機(jī)制*

藍(lán)紹航，龐宇舟，區(qū)佩琪，韋露秋，羅筱凡，呂 挺

1 廣西壯族自治區(qū)人民醫(yī)院，廣西 南寧 530022； 2 廣西中醫(yī)藥大學(xué)，廣西 南寧 530001

類風(fēng)濕性關(guān)節(jié)炎（rheumatoid arthritis，RA）是一種以滑膜的炎性反應(yīng)和異常增生為病理表現(xiàn)，影響關(guān)節(jié)和關(guān)節(jié)外結(jié)構(gòu)的慢性、全身性、自身免疫性疾?。?］，臨床以多關(guān)節(jié)炎、關(guān)節(jié)腫痛、軟骨和骨破壞，甚至關(guān)節(jié)畸形和功能喪失為主要表現(xiàn)特征。我國(guó)RA 的患病率約為0.42%，并呈現(xiàn)逐年上升的趨勢(shì)［2］。RA 至今病因未明，西醫(yī)認(rèn)為主要與遺傳、感染、內(nèi)分泌失調(diào)和藥物誘導(dǎo)等因素有關(guān)，多選用非甾體類抗炎藥、激素、抗風(fēng)濕藥、生物藥劑及小分子植入藥等治療，但這些藥物存在劑量大、給藥頻繁、毒副作用嚴(yán)重等缺點(diǎn)［3］。

RA 屬于中醫(yī)“痹證”范疇?！端貑?wèn)·痹論篇》曰：“風(fēng)、寒、濕三氣雜至，合而為痹也。”［4］《黃帝內(nèi)經(jīng)》首提“痹”之病名，明確提出痹證是風(fēng)、寒、濕三氣雜至，故多以祛風(fēng)、散寒、除濕為基本治法。甘草附子湯源自張仲景《傷寒雜病論》：“風(fēng)濕相摶，骨節(jié)疼煩，掣痛不得屈伸，近之則痛劇，汗出短氣，小便不利，惡風(fēng)不欲去衣，或身微腫者，甘草附子湯主之?！保?］方中桂枝溫通陽(yáng)氣，祛風(fēng)和營(yíng)；白術(shù)苦溫，健脾燥濕；附子辛熱扶陽(yáng)溫經(jīng)，長(zhǎng)于祛寒。附、桂、術(shù)三藥相伍，既除風(fēng)、寒、濕，又溫表陽(yáng)而固衛(wèi)氣，故被譽(yù)為治風(fēng)濕之圣藥［6］。又痹證急而劇，甘草既調(diào)和諸藥，又緩急止痛。全方合用，可有效抓住RA 病源，切中病機(jī)，扼其病勢(shì)，提高療效。近年來(lái)運(yùn)用甘草附子湯治療RA 也取得一定的研究進(jìn)展，廖常志運(yùn)用甘草附子湯治療RA 患者23 例，總有效率達(dá)100%［7］。王穎等［8］通過(guò)數(shù)據(jù)挖掘?qū)χ委烺A 常見(jiàn)藥物頻次總結(jié)后發(fā)現(xiàn)，桂枝和甘草分別排在出現(xiàn)頻次最高的第2 位和第4 位。研究發(fā)現(xiàn)，甘草附子湯能抑制T 淋巴細(xì)胞的增殖與分化，減少炎性細(xì)胞向滑膜組織的浸潤(rùn)，下調(diào)血清相關(guān)細(xì)胞因子白細(xì)胞介素2（interleukin-2，IL-2）、腫瘤壞死因子α（tumor necrosis factor，TNF-α）的表達(dá)水平，進(jìn)而起到抗炎和消除關(guān)節(jié)腫脹的作用［9］。此外，蔡悅等［10］研究發(fā)現(xiàn)，甘草附子湯可通過(guò)減少TNF-α、白細(xì)胞介素1β（interleukin-1β，IL-1β）等炎性細(xì)胞因子的產(chǎn)生，降低滑膜組織中基質(zhì)金屬蛋白酶1（matrix metalloproteinases-1）、MMP-3蛋白表達(dá)，從而抑制RA小鼠關(guān)節(jié)炎癥和骨侵蝕的發(fā)生。

目前對(duì)于甘草附子湯治療RA 的研究較多，但尚未對(duì)其作用機(jī)制進(jìn)行系統(tǒng)闡釋。網(wǎng)絡(luò)藥理學(xué)是基于系統(tǒng)生物學(xué)、藥理學(xué)、藥物化學(xué)和生物化學(xué)的一門(mén)綜合性的新興學(xué)科，為中藥發(fā)展提供新的研究思想，是藥物研究的新模式［11］，為明確藥物作用機(jī)制、新藥研發(fā)、復(fù)方用藥等提供理論依據(jù)和技術(shù)支撐［12］。同時(shí)，分子對(duì)接（molecular docking）也可揭示藥物治療疾病的分子機(jī)制，將藥理作用上升到分子水平上，為研究藥物機(jī)制奠定物質(zhì)基礎(chǔ)［13］。因此，本研究基于網(wǎng)絡(luò)藥理學(xué)結(jié)合分子對(duì)接技術(shù)，利用甘草附子湯中的化合物和受體分子虛擬結(jié)合，發(fā)現(xiàn)其作用靶點(diǎn)和結(jié)合位點(diǎn)，闡述甘草附子湯多成分、多靶點(diǎn)、多通路的協(xié)同作用機(jī)制，彌補(bǔ)片面的、不系統(tǒng)的研究，為后續(xù)深入研究提供相應(yīng)的理論基礎(chǔ)和科學(xué)依據(jù)。

1 資料與方法

1.1 藥物有效成分篩選及靶點(diǎn)預(yù)測(cè)在中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫(kù)與分析平臺(tái)（traditional Chinese medicine systems pharmacology database，and analysis platform，TCMSP）數(shù)據(jù)庫(kù)（http：//lsp.nwu.edu.cn/tcmsp.php）篩選甘草附子湯（甘草、附子、桂枝、白術(shù)）的全部化學(xué)成分，并以口服生物利用度（oral bioavailability，OB）≥30%、類藥性指數(shù)（drug-likeness，DL）≥0.18作為篩選標(biāo)準(zhǔn)，確定藥物的活性成分。然后在Related Targets一欄中找出其所對(duì)應(yīng)的作用靶點(diǎn)，將所有藥物靶點(diǎn)進(jìn)行合并整理，剔除重復(fù)的藥物靶點(diǎn)，在Uniport數(shù)據(jù)庫(kù)（https：//www.uniprot.org/）中輸入靶點(diǎn)名稱并限定物種為homo sapiens（human），對(duì)靶點(diǎn)進(jìn)行蛋白名稱標(biāo)準(zhǔn)化，獲得靶點(diǎn)對(duì)應(yīng)的標(biāo)準(zhǔn)基因名。

1.2 RA疾病靶點(diǎn)獲取以“rheumatoid arthritis”為檢索詞，在GeneCards（https：//www.genecards.org/）、OMIM（https：//omim.org/）、Pharm-Gkb（https：//www.pharmgkb.org/）、TTD（http：//db.idrblab.net/ttd/）、DrugBank（https：//www.drugbank.ca/）等數(shù)據(jù)庫(kù)獲取RA 相關(guān)靶點(diǎn)基因。對(duì)收集的所有RA 靶點(diǎn)進(jìn)行處理，刪除重復(fù)靶點(diǎn)，即為所需的RA 總靶點(diǎn)。再利用R 語(yǔ)言繪制RA 靶點(diǎn)基因交集Venny圖。

1.3 甘草附子湯與RA 交集靶點(diǎn)篩選將上述獲得的藥物靶點(diǎn)與疾病靶點(diǎn)基因通過(guò)R 語(yǔ)言繪制Venny 圖，對(duì)甘草附子湯活性成分靶點(diǎn)與RA 相關(guān)靶點(diǎn)基因進(jìn)行映射，得到藥物與疾病交集靶點(diǎn)。

1.4 構(gòu)建活性成分-核心靶點(diǎn)網(wǎng)絡(luò)將以上所得的甘草附子湯所有藥物的活性成分交集和RA 靶點(diǎn)基因交集數(shù)據(jù)導(dǎo)入Cytoscape 3.7.2 軟件，并對(duì)其參數(shù)進(jìn)行設(shè)置，得到活性成分-核心靶點(diǎn)相互作用網(wǎng)絡(luò)圖。

1.5 構(gòu)建蛋白質(zhì)相互作用（protein-protein interaction，PPl）網(wǎng)絡(luò)及核心基因預(yù)測(cè)將篩選得到的“Drug-Disease”中的靶點(diǎn)蛋白導(dǎo)入STRING 數(shù)據(jù)庫(kù)（https：//string-db.org/），限定研究物種為“Homo sapiens”，最低相互作用評(píng)分設(shè)置為“0.900”，同時(shí)設(shè)置隱藏游離點(diǎn)，其余參數(shù)保持默認(rèn)，進(jìn)行蛋白質(zhì)互作分析，得到PPI 網(wǎng)絡(luò)圖。將獲得數(shù)據(jù)導(dǎo)入Cytoscape 3.7.2 軟件，使用Network Analyzer插件對(duì)PPI網(wǎng)絡(luò)進(jìn)行拓?fù)浞治龊Y選關(guān)鍵基因。

1.6 GO 功能富集分析與KEGG 通路分析將“1.3”項(xiàng)下篩選得到的藥物疾病靶點(diǎn)交集通過(guò)導(dǎo)入R 語(yǔ)言軟件中進(jìn)行GO 功能富集分析和KEGG 通路分析。兩者皆以P＜0.01為篩選標(biāo)準(zhǔn)。

1.7 分子對(duì)接為驗(yàn)證篩選出的甘草附子湯治療RA 可能關(guān)鍵靶點(diǎn)的準(zhǔn)確性，將“1.1”項(xiàng)和“1.5”項(xiàng)中整理出的排名靠前的活性成分與藥物靶點(diǎn)進(jìn)行 分 子 對(duì) 接。在Pubchem（https：//pubchem.ncbi.nlm.nih.gov/）數(shù)據(jù)庫(kù)中下載活性成分的小分子配體2D 結(jié)構(gòu)，將2D 文件導(dǎo)入ChemBio3D Ultra 14.0軟件生成3D結(jié)構(gòu)，使用優(yōu)化模塊進(jìn)行分子結(jié)構(gòu)優(yōu)化，獲取最穩(wěn)定的分子構(gòu)象。從PDB（https：//www.rcsb.org/）數(shù)據(jù)庫(kù)下載靶點(diǎn)蛋白受體三維結(jié)構(gòu)，并將其導(dǎo)入PyMOLWin 軟件，輸入命令刪除靶蛋白的水分子和小分子配體，完成蛋白受體準(zhǔn)備。將蛋白受體文件導(dǎo)入AutoDock-Tools-1.5.6 軟件中，對(duì)蛋白分子加氫等處理，保存文件為pdbqt 格式，再將小分子配體文件導(dǎo)入AutoDockTools-1.5.6 軟件中，把文件轉(zhuǎn)換為pdbqt 格式。以靶點(diǎn)的原型配體所在位置的中心定義為活性對(duì)接口袋位點(diǎn)，口袋大小設(shè)為1.0 nm，其他參數(shù)根據(jù)需要適當(dāng)調(diào)整，確定活性口袋，使用Autodock vina軟件進(jìn)行對(duì)接。

2 結(jié)果

2.1 甘草附子湯有效活性成分篩選及靶點(diǎn)預(yù)測(cè)經(jīng)TCMSP 數(shù)據(jù)庫(kù)搜索，搜索關(guān)鍵詞為“甘草、桂枝、白術(shù)、附子”，以O(shè)B≥30%，DL≥0.18 為篩選條件，檢索到甘草附子湯4 味藥材的有效活性成分共有127 個(gè)，其中甘草92 個(gè)、桂枝7 個(gè)、白術(shù)7 個(gè)、附子21 個(gè)。整理后預(yù)測(cè)靶點(diǎn)為1895 個(gè)，經(jīng)Uniprot 數(shù)據(jù)庫(kù)標(biāo)準(zhǔn)化注釋并刪除重復(fù)靶點(diǎn)后得到133個(gè)作用靶點(diǎn)。見(jiàn)表1。

表1 “中藥-活性成分-疾病靶點(diǎn)”信息 個(gè)

2.2 藥物靶點(diǎn)與疾病靶點(diǎn)交集經(jīng)分析得到藥物潛在活性成分靶點(diǎn)133個(gè)，RA疾病基因3750個(gè)，見(jiàn)圖1；二者得出的交集靶點(diǎn)基因?yàn)?33個(gè)，見(jiàn)圖2。

圖1 疾病基因合并Venny圖

圖2 藥物靶點(diǎn)和疾病靶點(diǎn)Venny圖

2.3 活性成分-核心靶點(diǎn)相互作用網(wǎng)絡(luò)將“2.2”得到的133 個(gè)交集靶點(diǎn)運(yùn)用Perl 相關(guān)腳本進(jìn)行文件處理，所得數(shù)據(jù)導(dǎo)入Cytoscape 3.7.2軟件進(jìn)行可視化網(wǎng)絡(luò)構(gòu)建，得到活性成分-核心靶點(diǎn)相互作用網(wǎng)絡(luò)圖，見(jiàn)圖3。圖3中圓形節(jié)點(diǎn)代表藥物的有效成分，綠、紅、橙、藍(lán)色分別代表甘草、白術(shù)、桂枝、附子的有效成分，多個(gè)顏色代表多個(gè)藥物共同的有效成分，灰色線條代表邊，黃色三角形節(jié)點(diǎn)代表中藥和疾病的交集靶點(diǎn)。圖中共有246個(gè)節(jié)點(diǎn)，其中來(lái)源于化合物的有99個(gè)，來(lái)源于基因的有147 個(gè)；經(jīng)統(tǒng)計(jì)degree 較高的有效成分為槲皮素（quercetin）、山柰酚（kaempferol）、柚皮素（naringenin）、芒柄花素（formononetin）等。見(jiàn)表2。

圖3 活性成分-核心靶點(diǎn)相互作用網(wǎng)絡(luò)圖

表2 度值大于20的前4位化學(xué)成分

2.4 PPl 網(wǎng)絡(luò)及核心基因預(yù)測(cè)將“2.2”項(xiàng)得到的133 個(gè)交集靶點(diǎn)導(dǎo)入STRING 數(shù)據(jù)庫(kù)，設(shè)定物種為“Homo sapiens”，得出string_interactions文件，并以tsv 格式保存。網(wǎng)絡(luò)圖共有132 個(gè)節(jié)點(diǎn)、435條邊，平均每個(gè)節(jié)點(diǎn)的度值為6.59，平均局部聚類系數(shù)為0.489，PPI enrichment p-value值為＜1.0e-16，鄰接節(jié)點(diǎn)數(shù)目越多，說(shuō)明這些蛋白在網(wǎng)絡(luò)中的作用越顯著。見(jiàn)圖4。將tsv 格式文件數(shù)據(jù)導(dǎo)入Cytoscape 3.7.2 軟件中，使用Network Analyzer 插件對(duì)PPI 網(wǎng)絡(luò)進(jìn)行拓?fù)浞治龊Y選關(guān)鍵基因，處于核心位置的靶點(diǎn)有JAK2、JUN、IL-10、TNF、CSF2、CCL2、JAK3、IL-1B、JAK1、MAPK14、IL-2、IL-6等，提示這些基因在PPI中發(fā)揮了關(guān)鍵作用。見(jiàn)圖5。

圖4 交集基因PPI網(wǎng)絡(luò)圖

2.5 GO 富集分析和KEGG 通路富集分析經(jīng)GO富集分析共得到1912 條富集結(jié)果，其中生物過(guò)程（biological process，BP）1792 個(gè)，主要涉及對(duì)細(xì)胞因子產(chǎn)生的正調(diào)控、對(duì)細(xì)胞黏附調(diào)節(jié)、對(duì)炎癥反應(yīng)的調(diào)節(jié)、對(duì)細(xì)菌源分子的反應(yīng)、對(duì)脂多糖的反應(yīng)等生物過(guò)程的調(diào)節(jié)；細(xì)胞組分（cellular component，CC）28 個(gè)，主要涉及質(zhì)膜外側(cè)、膜筏、膜微區(qū)、膜區(qū)等細(xì)胞組分；分子功能（molecular function，MF）89個(gè)，主要涉及細(xì)胞因子受體結(jié)合、細(xì)胞因子活性、受體配體活性、信號(hào)受體激活劑活性等。提示甘草附子湯在治療RA 過(guò)程中通過(guò)參與調(diào)控多種生物學(xué)過(guò)程而發(fā)揮作用。見(jiàn)圖6。

圖6 交集基因GO條形圖

經(jīng)KEGG 通路富集分析共得到131條通路富集結(jié)果，主要涉及細(xì)胞因子-細(xì)胞因子受體相互作用通路、恰加斯病、類風(fēng)濕性關(guān)節(jié)炎、新型冠狀病毒肺炎、肺結(jié)核、人巨細(xì)胞病毒感染、TNF信號(hào)通路）、破骨細(xì)胞分化、人T細(xì)胞白血病病毒1感染等。提示甘草附子湯可能通過(guò)以上通路發(fā)揮對(duì)RA的治療作用。見(jiàn)圖7。

圖7 交集基因KEGG條形圖

2.6 潛在活性成分-核心靶點(diǎn)分子對(duì)接分析根據(jù)甘草附子湯“活性成分-核心靶點(diǎn)”網(wǎng)絡(luò)中靶點(diǎn)與成分連接的度值，選取活性化合物中degree＞20 的槲皮素、山柰酚、柚皮素、芒柄花素4 個(gè)主要有效成分，分別與經(jīng)拓?fù)浞治龊Y選后的關(guān)鍵核心靶點(diǎn)IL-6、JUN、MAPK14 進(jìn)行分子對(duì)接。分子之間結(jié)合能量較低的構(gòu)象則較為穩(wěn)定，篩選對(duì)接結(jié)果中結(jié)合能量低的分子構(gòu)象即為最優(yōu)對(duì)接結(jié)果，結(jié)果表明槲皮素、山柰酚、柚皮素、芒柄花素與IL-6、JUN、MAPK14 的結(jié)合能均≤-6.0 kcal/mol，展現(xiàn)出良好的結(jié)合力。見(jiàn)圖8，表3。

圖8 潛在活性成分-核心靶點(diǎn)分子對(duì)接分析

表3 槲皮素、山柰酚、柚皮素和芒柄花素與RA潛在靶點(diǎn)結(jié)合能量 kcal/mol

3 討論

甘草附子湯臨床應(yīng)用廣泛，為深入研究其治療RA 的作用機(jī)制，本研究采用網(wǎng)絡(luò)藥理學(xué)的方法，共篩選出127 個(gè)活性成分及133 個(gè)靶點(diǎn)，其中槲皮素、山柰酚、柚皮素、芒柄花素等活性成分對(duì)應(yīng)靶點(diǎn)較多，是甘草附子湯治療RA 的主要成分。槲皮素具有抗炎、抗氧化、調(diào)節(jié)免疫等作用［14-15］，可促進(jìn)RA成纖維樣滑膜細(xì)胞的凋亡［16］，降低IL-1β誘導(dǎo)的前列腺素E2（prostaglandin E2，PGE2）產(chǎn)物量與MMP-1、MMP-3、MMP-13的含量表達(dá)，抑制NF-κB激活，減弱炎癥環(huán)境下軟骨細(xì)胞內(nèi)MMP-13 的產(chǎn)生、基質(zhì)降解和細(xì)胞凋亡，從而起到保護(hù)軟骨的作用［17-18］。山柰酚具有顯著的抗炎作用，可通過(guò)抑制絲裂原活化蛋白激酶（MAPK）、磷脂酰肌醇-3-激酶（PI3K）等信號(hào)通路的表達(dá)，調(diào)控炎癥反應(yīng)過(guò)程，減少炎癥因子產(chǎn)生，從而有效抑制炎癥的發(fā)生發(fā)展［19］。柚皮素、芒柄花素等活性成分也有抗炎、抗氧化、調(diào)節(jié)免疫等作用［14-15］，但關(guān)于其治療RA 的作用目前鮮見(jiàn)報(bào)道，這可作為將來(lái)研究中藥活性成分治療RA 的新方向。上述成分在經(jīng)過(guò)分子對(duì)接后進(jìn)一步驗(yàn)證其可信度，為下一步臨床研究治療RA提供了新依據(jù)。

通過(guò)PPI 網(wǎng)絡(luò)拓?fù)浞治龊蠛Y選到JAK2、JUN、IL-10、TNF、CSF2、CCL2、JAK3、IL1B、JAK1、MAPK14、IL-2、IL-6 等12 個(gè)核心靶點(diǎn)，其中以IL-6 最為關(guān)鍵。IL-6 是活動(dòng)期RA 患者血清及關(guān)節(jié)液中含量最高的細(xì)胞因子之一，其功能涉及多種與炎癥相關(guān)的疾病狀態(tài)，包括對(duì)糖尿病和全身性RA 的可疑性［20-21］。同時(shí)，破骨細(xì)胞活性與關(guān)節(jié)破壞程度密切相關(guān)，IL-6 作為一種炎性因子，可引起人體內(nèi)破骨細(xì)胞與成骨細(xì)胞之間的平衡失調(diào)，是造成RA發(fā)生發(fā)展的主因之一［22］。高濃度的IL-6 可刺激骨吸收的調(diào)節(jié)作用，進(jìn)一步刺激破骨細(xì)胞的前體細(xì)胞，激發(fā)破骨細(xì)胞活性，導(dǎo)致骨質(zhì)破壞，通過(guò)多通路、多靶點(diǎn)參與RA 的炎癥反應(yīng)過(guò)程［23-24］。抑制IL-6的表達(dá)已成為治療RA的新方法［25］，托珠單抗（tocilizum）是目前廣泛應(yīng)用于臨床治療RA 的強(qiáng)效藥，可通過(guò)結(jié)合IL-6 跨膜及溶解性受體來(lái)抑制IL-6 介導(dǎo)的信號(hào)轉(zhuǎn)導(dǎo)，從而有效地緩解癥狀，預(yù)防關(guān)節(jié)破壞［26-27］。絲裂原活化蛋白激酶（MAKP）可參與細(xì)胞分化、增殖、生長(zhǎng)、免疫調(diào)控等多種活動(dòng)［28］，有研究發(fā)現(xiàn)，microRNA-137 通過(guò)MAPK 信號(hào)通路和IL-6/COX-II 活化在PMS 誘導(dǎo)的RA 急性加重中起重要作用［29］。轉(zhuǎn)錄因子AP-1（JUN）由多種蛋白復(fù)合物組成，廣泛參與調(diào)節(jié)細(xì)胞生命，如細(xì)胞分化、增殖、凋亡等過(guò)程，與外界環(huán)境壓力、細(xì)菌及病毒感染、多種疾病和病理?yè)p傷的發(fā)生發(fā)展密切相關(guān)，如神經(jīng)退行性疾病、RA、腫瘤等［30-31］。TNF靶點(diǎn)也非常重要，TNF 參與機(jī)體的免疫調(diào)節(jié)和炎癥反應(yīng)［32］，而TNF-α 在維護(hù)骨、軟骨組織的完整性及促進(jìn)損傷修復(fù)方面發(fā)揮重要作用［33］。相關(guān)研究表明，其余靶點(diǎn)基因也與RA 密切相關(guān)，在此不再贅述。由此推測(cè)，甘草附子湯中的中藥有效成分可能是通過(guò)上述靶點(diǎn)發(fā)揮治療RA的作用。

GO 富集分析主要涉及對(duì)細(xì)胞因子產(chǎn)生的正調(diào)控、對(duì)細(xì)胞黏附調(diào)節(jié)、對(duì)炎癥反應(yīng)的調(diào)節(jié)、對(duì)細(xì)菌源分子的反應(yīng)、對(duì)脂多糖的反應(yīng)等生物過(guò)程的調(diào)節(jié)。其中，細(xì)胞因子是RA 炎癥及關(guān)節(jié)損傷的重要介質(zhì)，是免疫原、絲裂原或其他刺激劑誘導(dǎo)多種細(xì)胞產(chǎn)生的低分子量可溶性蛋白質(zhì)，具有調(diào)節(jié)固有免疫和適應(yīng)性免疫、血細(xì)胞生成、細(xì)胞生長(zhǎng)及損傷組織修復(fù)等多種功能，細(xì)胞因子網(wǎng)絡(luò)的不平衡可能誘導(dǎo)和維持炎癥及相應(yīng)的組織損傷，在RA 發(fā)病機(jī)制中有重要作用，與RA 最密切相關(guān)的細(xì)胞因子包括TNF-α、IL-6、IL-2 等［34-36］。脂多糖可刺激內(nèi)皮細(xì)胞誘導(dǎo)轉(zhuǎn)錄因子核因子κB（nuclear factor kappa-B，NF-κB）活性升高，并上調(diào)其他炎癥因子（如TNF-α、IL-1β）的表達(dá)，導(dǎo)致炎癥的持續(xù)發(fā)生以及骨質(zhì)的持續(xù)破壞并與其發(fā)生協(xié)同作用，產(chǎn)生炎癥性疾?。?7］。由此可知，在RA 的病理過(guò)程中，炎癥反應(yīng)是重要環(huán)節(jié)，如何使用中藥緩解RA的炎癥反應(yīng)是治療的重要內(nèi)容。

KEGG 通路富集分析結(jié)果顯示，甘草附子湯作用RA 的靶點(diǎn)主要涉及細(xì)胞因子-細(xì)胞因子受體相互作用通路、恰加斯病、類風(fēng)濕性關(guān)節(jié)炎、新型冠狀病毒肺炎、肺結(jié)核、人巨細(xì)胞病毒感染、TNF 信號(hào)通路、破骨細(xì)胞分化通路、人T 細(xì)胞白血病病毒1 感染等。以上通路大多數(shù)涉及炎癥反應(yīng)和骨破壞，而炎癥反應(yīng)和骨破壞是RA 的主要臨床表現(xiàn)，骨平衡被破壞是導(dǎo)致骨侵蝕的主因，炎癥反應(yīng)直接或間接誘導(dǎo)破骨細(xì)胞分化，并使其功能減少，導(dǎo)致RA骨破壞［38］。TNF信號(hào)通路是炎癥反應(yīng)的重要通路，可以激活MAPK/NF-κβ等加劇炎癥反應(yīng)［39］。良好的TNF 抑制劑臨床療效也驗(yàn)證了該因子及其通路在RA 中的重要地位［40］。除上述信號(hào)通路外，尚有多條信號(hào)通路參與甘草附子湯治療RA 的過(guò)程，這些信號(hào)通路為我們進(jìn)一步研究甘草附子湯與RA 的關(guān)系提供一定依據(jù)和思路。由此預(yù)測(cè)，甘草附子湯主要是通過(guò)調(diào)節(jié)上述通路下調(diào)炎性因子表達(dá)，調(diào)節(jié)骨代謝，干預(yù)骨破壞，調(diào)節(jié)免疫功能等發(fā)揮治療RA的作用，從而改善患者的生活質(zhì)量。

為驗(yàn)證網(wǎng)絡(luò)預(yù)測(cè)的準(zhǔn)確性，利用分子對(duì)接技術(shù)，選取甘草附子湯中4 個(gè)核心活性成分（槲皮素、山柰酚、柚皮素、芒柄花素）與關(guān)鍵靶點(diǎn)IL-6、JUN、MAPK14 結(jié)合，結(jié)果提示結(jié)合能量均≤-6.0 kcal/mol，匹配能量較低，說(shuō)明小分子與大分子結(jié)合的構(gòu)象穩(wěn)定，結(jié)合性好，展現(xiàn)出良好的親和力，也驗(yàn)證了本研究成果。

綜上所述，RA 是一種病理機(jī)制非常復(fù)雜的疾病，甘草附子湯治療RA 的作用機(jī)制是一個(gè)多基因、多生物過(guò)程及多信號(hào)通路共同參與過(guò)程，符合中藥治療疾病多成分、多靶點(diǎn)的特征。本研究通過(guò)網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接技術(shù)，初步預(yù)測(cè)甘草附子湯對(duì)RA 的主要有效成分、作用靶點(diǎn)、相關(guān)通路及協(xié)同作用的分子機(jī)制，為后續(xù)進(jìn)一步開(kāi)展機(jī)制研究、藥物開(kāi)發(fā)及臨床試驗(yàn)等奠定基礎(chǔ)，也提供了理論依據(jù)、線索與思路。