姚茹冰，彭 浩，蔡孟成，李 霞 （海軍軍醫(yī)大學：. 中醫(yī)系, . 基礎醫(yī)學院，上海 200433）

類風濕關節(jié)炎（RA）以對稱性多關節(jié)滑膜炎為主要臨床表現(xiàn)，具有異質(zhì)性、系統(tǒng)性，呈慢性、進行性、侵襲性，如無恰當治療，則病情逐漸加重而形成永久性骨質(zhì)破壞并最終出現(xiàn)殘疾，甚至累及臟器和神經(jīng)系統(tǒng)而危及生命。通過早診斷、早治療，阻止關節(jié)炎癥病變發(fā)展成不可逆的骨破壞至關重要。目前常用的治療RA 藥物如非甾體抗炎藥、糖皮質(zhì)激素、改善病情的抗風濕藥、免疫抑制劑及生物制劑等，可在較短的時間改善患者的臨床癥狀，但長期使用存在不良反應增多、費用較高、單一靶向作用療效受限等問題，因此尋找有效、低毒、價廉及適合長期服用的RA 藥物，成為研究熱點。

RA 屬于中醫(yī)學“歷節(jié)風”、“頑痹”、“骨痹”、“鶴膝風”等范疇?！端貑枴け哉摗吩疲骸帮L寒濕三氣雜至，合而為痹”，治療以祛風散寒除濕，通絡止痛為主。青風藤為防己科植物青藤Sinomenium acutum（Thunb.）Rehd. et Wils.和毛青藤Sinomenium acutum（Thunb.）Rehd. et Wils. var. cinereum Rehd.et Wils.的干燥藤莖（《中國藥典》2020 年版)，其性味苦、辛、平，歸肝、脾經(jīng)，能有效祛風除濕，通絡止痛。《本草匯言》載“青風藤，散風寒濕痹之藥，能舒筋活血，正骨利髓?！薄侗静菥V目》認為有“治風濕流注，歷節(jié)鶴膝”之功?，F(xiàn)代研究顯示以青風藤為主要藥物的湯劑，可有效延緩疾病發(fā)展，改善患者臨床癥狀[1-2]；青風藤的有效成分青藤堿也具有較好的抗炎鎮(zhèn)痛和抑制自身抗體產(chǎn)生的作用，在治療RA 方面，具有較好的臨床療效和較低的副作用[3]。本文旨在通過對青風藤進行網(wǎng)絡藥理學分析，以期探索其物質(zhì)基礎及可能的作用機制，為進一步研究青風藤治療RA 提供理論依據(jù)。

1 材料與方法

1.1 青風藤活性成分的篩選

中藥系統(tǒng)藥理學數(shù)據(jù)庫和分析平臺(TCMSP）是國內(nèi)相對完善的中藥成分靶點數(shù)據(jù)庫，收錄了中藥502 味，相關靶點3 339 個，同時對每個成分的藥動學參數(shù)進行了詳細的描述。通過TCMSP 獲取青風藤的化學成分。再根據(jù)化合物口服生物利用度（OB）與類藥性指數(shù)（DL）篩選活性成分，其中OB 閾值設為OB≥30%，DL 閾值設為DL≥0.18[4-5]。下載活性成分的結構，將活性成分的結構保存為mol2 格式。

1.2 青風藤潛在作用靶點的預測

DRAR-CPI（ https://cpi.bio-x.cn/drar/） 是 基 于DOCK 分子對接程序開發(fā)的反向分子對接服務器，可對接藥物分子和蛋白的三維結構和特征，并通過蛋白-配體相互作用能的親和性以函數(shù)打分并排序。對接分數(shù)（Z-score）代表藥物分子與蛋白的相互作用強度。服務器研發(fā)者認為，Z-score ＜-1 的蛋白質(zhì)是良好的潛在靶點。通過DRAR-CPI 服務器篩選出青風藤有效成分的蛋白靶點后，利用UniProt數(shù)據(jù)庫對靶點處理，得到標準名稱，獲取基因靶點。

1.3 RA 相關靶點的篩選

通過對Geencards 數(shù)據(jù)庫、OMIM 數(shù)據(jù)庫進行檢索，得到RA 的相關靶點。

1.4 藥物-化合物-疾病-靶點網(wǎng)絡構建及可視化分析

為明確青風藤與RA 潛在靶點間的相互作用，利用Venn 在線軟件將藥物與疾病靶點取交集，通過Venn 圖的形式展現(xiàn)藥物靶點與疾病靶點間存在的交集靶點，運用Cytoscape 軟件進行可視化分析。

1.5 靶蛋白相互作用（PPI）網(wǎng)絡的構建

將篩選得到的藥物-疾病共同靶點上傳至String 在線平臺數(shù)據(jù)庫，構建蛋白質(zhì)相互作用網(wǎng)絡模型，獲得PPI 網(wǎng)絡。

1.6 功能富集分析與通路分析

利用clusterProfiler 程序包（3.6 版）對藥物-疾病共同靶點進行GO 分析和KEGG 通路分析，設定閾值P＜0.05，并按照涉及的靶點數(shù)目多少進行排序，篩選GO 富集分析中分子功能、細胞組成和生物過程的前10 個條目以及KEGG 富集的前10 個信號通路，進一步篩選與RA 相關的信號通路，探討青風藤抗RA 的生物學過程及信號通路。

2 結果

2.1 青風藤有效化學成分的篩選

從TCMSP 共獲得青風藤16 個化學成分，以OB＞30%，DL＞0.18 為篩選條件，共篩選出有效活性成分6 個（表1）。

表1 青風藤有效成分

2.2 青風藤潛在作用靶點的預測

登陸DRAR-CPI 服務器并上傳青風藤有效活性成分結構的mol2 格式文件，篩選Z-score ＜-1的蛋白靶點的PDB ID 導入Uniprot 數(shù)據(jù)庫，獲取6 個有效活性成分的潛在基因靶點合計383 個。其中，存在多個有效成分對應同一基因靶點的情況，經(jīng)去重處理后，共獲得青風藤潛在基因靶點176 個。

2.3 RA 靶點的篩選

通過Geencards 數(shù)據(jù)庫得到RA 相關靶點4329 個，通過篩選設置relevance score≥15，最終獲得305 個相關靶點，OMIM 數(shù)據(jù)庫得到與RA 相關靶點13 個，篩重后獲得疾病相關靶點305 個。

2.4 “中藥-化合物-靶點-疾病”的網(wǎng)絡構建及可視化分析

利用Venn 2.1 軟件將“青風藤”活性成分對應的176 個靶點與RA 對應的305 個靶點進行交集，獲得15 個共有靶點（圖1）。將15 個共有靶點導入Cytoscape 3.7.2 軟件，構建“青風藤-化合物-靶點-RA”的網(wǎng)絡圖并進行可視化分析（圖1）。該網(wǎng)絡包含了23 個節(jié)點，100 條邊，其中，三角形代表藥物；燕尾形代表疾??；菱形代表活性成分；橢圓形代表靶點；連接化合物與靶點的邊表示它們之間有作用關系。依據(jù)網(wǎng)絡拓撲學性質(zhì)可知，靶點度值(degree)是拓撲結構中的重要參數(shù)，其表示網(wǎng)絡中和節(jié)點相連的路線的條數(shù)，度值高的節(jié)點很可能在整個網(wǎng)絡中起到橋梁作用，有較多節(jié)點的化合物或藥物靶點在整個網(wǎng)絡中可能起到關鍵的作用。本研究篩選度值較大的節(jié)點進行分析，排名前3 位的化合物分別為青藤堿（sinomenine）、拉茲馬寧堿（16-epi-isositsirikine）及烏心石環(huán)氧內(nèi)酯（michelenolide），其分別能與8、6 及6 個靶點蛋白發(fā)生作用。從靶點的角度看，排名前4 位的是IL-10、B2M、NR3C1 及S100A9，分別能與4、4、3 及3 個化合物發(fā)生相互作用，體現(xiàn)了中藥多靶點的特性。

圖1 “青風藤-化合物-靶點-RA”可視化網(wǎng)絡圖

2.5 靶點PPI 網(wǎng)絡分析

將Venn 圖獲得的15 個藥物-疾病共同靶點上傳至String 在線數(shù)據(jù)庫平臺得到PPI 網(wǎng)絡（圖2），同時得到網(wǎng)絡中相關節(jié)點的度值。在此網(wǎng)絡中包含節(jié)點15 個，邊55 條，度值均值7.33。根據(jù)“度值＞均值”篩選出關鍵靶點7 個，包括：IL-10、IL-4、INS、MAPK8、ELANE、MAPK1、MAPK14。度值最高的是IL-10，為14；其次為IL-4，度值為13；INS度值為12。度值大的靶點提示在網(wǎng)絡調(diào)控中起著關鍵作用，度值大的靶點很可能是青風藤治療RA 的關鍵靶點。

2.6 青風藤治療RA 潛在靶點的GO 生物功能注釋

圖2 靶蛋白PPI 網(wǎng)絡圖

GO 是Gene Ontology 的簡稱，主要包括生物過程（BP）、細胞組成（CC）及分子功能（MF）3 個部分。本研究GO 功能富集分析顯示500 個生物過程（BP）、18 個細胞組成（CC）、28 個分子功能（MF），分別對BP、CC 及MF 前10 個條目進行可視化分析，得到圖3，其中縱坐標表示富集條目，橫坐標表示富集個數(shù)。圖4 顯示共同作用靶點在BP 中，主要富集條目為對細菌起源分子的反應、細胞對生物刺激的反應，調(diào)節(jié)DNA 結合轉(zhuǎn)錄因子活性、積極調(diào)節(jié)細胞因子的產(chǎn)生、調(diào)節(jié)炎癥反應，以及對中性粒細胞活化的調(diào)節(jié)、對脂多糖的反應等；共同作用靶點在CC 中主要富集條目為分泌顆粒內(nèi)腔、細胞質(zhì)囊泡腔、內(nèi)體腔、液泡腔等；共同作用靶點在MF 中主要富集條目為MAP 激酶活性、Toll樣受體結合、細胞因子受體結合、蛋白質(zhì)絲氨酸/蘇氨酸激酶活性、磷酸酶結合、細胞因子活性等。

2.7 青風藤治療RA 潛在靶點的KEGG 通路富集分析

KEGG 通路富集分析顯示77 條通路，對前10 條通路進行可視化分析，得到圖4，其中縱坐標表示富集通路的名稱，橫坐標表示富集基因個數(shù)。前10 條通路顯示，靶點主要富集于FcεRI 信號通路、百日咳、利什曼病、結核、IL-17 信號通路、恰加斯?。绹F蟲?。?、C 型凝集素受體信號通路、T 細胞受體信號通路、弓形蟲病、FoxO 信號通路等，多為感染、炎癥和免疫相關通路。其中涉及到的靶點基因，出現(xiàn)頻次較高的是MAPK1、MAPK14、MAPK8、IL-10 及IL-4。

圖3 靶點的GO 分析（含BP、CC 及MF）柱狀圖

圖4 靶點KEGG 通路富集分析圖

3 討論

痹病，是人體正氣不足或臟腑功能失調(diào)，風寒濕熱等邪為患，痰濁瘀血滯留，引起經(jīng)脈氣血不通，出現(xiàn)肢體關節(jié)疼痛腫脹及活動不利等為特征的疾病。中藥青風藤功能祛風除濕，通絡止痛，在痹病治療中應用較多，療效可靠。由于中藥是一個多成分多靶標并且可以調(diào)節(jié)機體內(nèi)多個生物學途徑的復雜混合物系統(tǒng)，中藥中的化學成分可以作用于機體內(nèi)多個靶標從而發(fā)揮協(xié)同作用來治療疾病，給中藥作用機制的研究帶來了一定的挑戰(zhàn)。近年來，網(wǎng)絡藥理學發(fā)展十分迅速，已成為在藥理學研究中一個很有價值的工具，現(xiàn)今也逐漸應用到中醫(yī)藥研究領域，對單位中藥及中藥復方中活性藥物成分的預測和分析、闡述復方的分子作用機制有重要意義[6-9]。

本研究通過TCMSP 篩選青風藤有效活性成分6 個，分別是β-谷固醇、拉茲馬寧堿、千金藤啶堿、青藤堿、烏心石環(huán)氧內(nèi)酯、甘露聚糖。通過DRAR-CPI 服務器篩選出青風藤靶點176 個。使用GeenCards 及OMIM 數(shù)據(jù)庫得到與RA 相關靶點305 個。藥物靶點與疾病靶點進行Venn 圖繪制，獲得共同相關的15 個靶點。將這15 個共同靶點上傳至String 在線數(shù)據(jù)庫平臺得到PPI 網(wǎng)絡，根據(jù)“度值＞均值”篩選出關鍵靶點7 個，包括：IL-10、IL-4、INS、MAPK8、ELANE、MAPK1、MAPK14。GO 功能富集分析顯示500 個BP、18 個CC、28 個MF。KEGG 通路富集分析顯示77 條通路，排名前10 的信號通路集中在感染、炎癥和免疫等相關通路。

已有的研究提示，β-谷甾醇是植物甾醇的重要組成成分，具有抗炎、調(diào)節(jié)人體甾體激素、抗氧化及抗腫瘤等作用[10]。Kripa 等研究發(fā)現(xiàn)，含有β-谷甾醇的白術醇提物可明顯降低佐劑型關節(jié)炎大鼠血清和組織中的TNF-α、C-反應蛋白和IL-2 的含量，證明含有β-谷甾醇的白術具有顯著的抗炎作用[11]。拉茲馬寧堿可調(diào)節(jié)COX-2 表達[12]。千金藤啶堿不僅具有鎮(zhèn)痛及對多巴胺D1 激動、D2 拮抗、D3 拮抗等多重作用，對心臟也具有較強的生物活性[13]。青藤堿具有免疫抑制、抗炎、鎮(zhèn)痛、抗心律失常、降血壓等多種藥理作用，已經(jīng)被證明在治療多種疾病中有藥理學活性[12]。

本研究中藥物與疾病共同的15 個靶點，主要參與對細菌起源分子的反應、細胞對生物刺激的反應、調(diào)節(jié)DNA 結合轉(zhuǎn)錄因子活性、調(diào)節(jié)細胞因子的產(chǎn)生、調(diào)節(jié)炎癥反應、對中性粒細胞活化的調(diào)節(jié)、對脂多糖的反應等生物學過程。主要參與FcεRI 信號通路、百日咳、利什曼病、結核、IL-17 信號通路、恰加斯病（美國錐蟲?。?、C 型凝集素受體信號通路、T 細胞受體信號通路、弓形蟲病、FoxO 信號通路等感染、炎癥和免疫相關通路。這些通路中涉及次數(shù)較多的靶基因依次為MAPK1、MAPK14、MAPK8、IL-10 及IL-4?？傊?，提示青風藤具有多成分、多靶點、多通路、多途徑的作用特點，通過不同靶點和通路發(fā)揮藥效以治療RA。

RA 特征性的病理改變包括關節(jié)滑膜組織炎性細胞浸潤、血管翳形成及對關節(jié)骨和軟骨以及周圍組織造成破壞與侵蝕，在免疫介導的RA 的炎癥、增生和骨破壞三大病理過程中，MAPK 都起著重要的調(diào)節(jié)作用。IL-10、IL-4 作為抑炎和抗炎因子，參與了RA 的炎癥過程，高水平的細胞因子不斷刺激T 細胞、破骨細胞、滑膜細胞，導致細胞內(nèi)的蛋白激酶異常激活和過度的信號傳導，從而調(diào)控相關基因的過量表達。研究表明，對膠原誘導的關節(jié)炎大鼠給予MAPK 信號轉(zhuǎn)導通路抑制劑后，在抑制滑膜炎癥、骨質(zhì)破壞和關節(jié)軟骨的破壞方面，與未給信號通路抑制劑組相比有顯著性差異[14]。以MAPK信號通路作為治療靶點，然后針對其信號傳遞中的各個階段抑制其活化，已經(jīng)成為治療RA 的一個熱點。研究顯示，用青藤堿與針刺聯(lián)合治療大鼠膠原性關節(jié)炎，關節(jié)炎評分降低，癥狀好轉(zhuǎn)。同時，核因子κB 和磷酸化的p38 絲裂原活化蛋白激酶（MAPK）蛋白表達被抑制，提示可能通過核因子κB 和MAPK 信號通路發(fā)揮作用[15]。

本研究通過中藥活性成分的篩選、靶點預測、網(wǎng)絡構建與分析，從“藥物-化合物-靶點-疾病”的關聯(lián)性對青風藤治療RA 的相關作用機制進行分析，找出了關鍵靶點和通路，MAPK 信號通路也可作為后續(xù)研究的參考。為中藥青風藤治療RA 的臨床和實驗研究提供了部分依據(jù)和理論支持。同時，本研究尚存在一些局限性，如數(shù)據(jù)庫中數(shù)據(jù)可能不完全，獲取小分子化合物及作用靶標數(shù)量有限，不能完全揭示中藥藥理作用；依托網(wǎng)絡構建的數(shù)據(jù)來源中存在對某些研究熱點傾向性，可能對本研究結果產(chǎn)生影響，但仍需進一步做驗證實驗，以準確揭示中藥抗RA 有效成分及作用機制。