基于網絡藥理學六君子湯治療胃輕癱的機制研究

彭美萍，郭思彤，藍 穎，劉曉霞，譚思濤

（廣西醫學科學院·廣西壯族自治區人民醫院，廣西 南寧 530021）

傳統中藥湯劑六君子湯源于《世醫得效方》，有益氣健脾、燥濕化痰之功效，主要用于脾胃不和、食欲不振，或反胃腹痛腹瀉、脾胃消化功能下降，即現代西醫理論中所說的慢性胃炎、胃潰瘍、胃輕癱等胃腸道疾病。近年來，也有很多高質量的多中心、雙盲、隨機試驗證實了六君子湯有助于恢復脾胃功能，減輕消化不良癥狀，改善胃排空情況[1]。 但是關于六君子湯發揮療效的關鍵靶點與分子機制鮮有報道。與研究化學藥的藥理特征以單一成分對應單一靶點相比，中藥具有多組分、多靶標和多途徑協同作用的特點。如何闡釋清楚其物質基礎及作用原理一直是中藥研究的一大難題，網絡藥理學可為中藥復方作用機制研究提供嶄新的思路與手段，這有助于揭示中藥復方的科學內涵、發掘復雜中藥的有效成分、預測作用靶點和新的適應證，為提升中醫藥創新研究和大數據挖掘提供理論支撐[2]。為此，本研究采用網絡藥理學方法，從六君子湯的物質基礎出發，著眼于整體角度分析六君子湯治療胃輕癱的分子機制。

1 材料與方法

1.1 六君子湯相關靶點篩選

利用中醫藥系統藥理學平臺（TCMSP 數據庫）[3]查找人參、半夏、白術、陳皮、茯苓、甘草六味中藥的化學組成成分，根據口服利用度（OB）≥30%且類藥性（DL）≥0.18 這兩項藥代動力學參數對活性化合物及其作用的蛋白質靶點進行初篩，結合化學專業數據庫[4],使用SwissADME（https://www.swissadme.ch）工具預測化合物的ADME（吸收、分布、代謝、排泄）性質并進行篩選，補充有效成分及其靶點，并參考相關文獻報道補充未預測到的活性化合物的已知靶點。匯總并刪除重復的靶點信息后，將化合物作用的蛋白質靶點全部在Uniprot 蛋白質數據庫（https://www.uniprot.org）進行規范。轉化為標準的蛋白靶點編號。

1.2 胃輕癱相關靶點篩選

以“Gastroparesis， Dysmotility，Esophageal Dysmotility，Gastric stasis，Gastrointestinal ulceration，Gastrointestinal Dysfunction，digestive disturbance” 作 為關鍵詞，收集OMIM 數據庫（https://www.omim.org）胃輕癱的潛在靶點信息。在DisGeNET數據庫（https://www.disgenet.org）中，選取與“Gastroparesis”聯系密切的相關疾病“Cyclical vomiting syndrome （disorder），Abnormality of the stomach，Abnormality of the gastric mucosa，Gastric ulcer，Functional Gastrointestinal Disorders”。根據查詢結果進行篩選，保留相關度較高的目標靶點作為胃輕癱的潛在靶點，合并兩個疾病數據庫靶點后，刪除重復值得到該疾病靶點。

1.3 六君子湯成分-胃輕癱靶點的蛋白質-蛋白質相互作用（PPI）網絡構建

為明晰六君子湯藥物相關靶點與胃輕癱靶點間的關聯，利用Excel 提取二者有交集的靶點并通過韋恩圖的方式來表達。用交集靶點通過STRING11.0 數據庫（https://string-db.org）構建PPI 網絡模型[5]，將生物種類設定為“Homosapiens”，最小互相作用閾值設定為“highest confidence ”（大于0.9），除此之外均為默認值，得到PPI網絡模型，通過CytoScape3.7.1 中的hubba 插件對PPI 網絡進行分析排序發現最佳的前20 個靶點，得到潛在的關鍵靶點基因，并通過分析其參與的生物學進程對其功能進行描述。

1.4 六君子湯成分-胃輕癱靶點功能與通路的富集分析

Metascape 平臺（https://metascape.org/gp/index.html）對基因擁有強大的分析功能，其覆蓋面也相當廣泛，整合了GO、KEGG、UniProt 和DrugBank 等多個權威的數據資源，可靠性強，通過該平臺可以對提交數據進行通路富集和生物過程注釋[6,7]。將六君子湯調節治療胃輕癱的靶點提交到Metascape 平臺進行數據分析，P值設置為小于0.01，除此之外均為默認值，分析其主要的生物學過程與代謝通路并進行富集分析。

2 結果

2.1 六君子湯活性成分靶點的獲取

初步提取的中藥化學成分結果如下：人參190 種、半夏116 種、白術55 種、陳皮63 種、茯苓34 種、甘草280 種，經藥物ADME 篩選后共獲得白術7 種活性成分，陳皮5 種活性成分，茯苓15 種活性成分，甘草92 種活性成分，繼續篩選刪除無對應靶點成分。各成分初步提取的作用靶點結果如下：人參116 個、半夏97 個、陳皮66 個、茯苓23 個、甘草230 個，合并除去相同的數值后共得到靶點267個。六君子湯成分-靶點網絡圖的構建，核心靶點的提取，主要產生藥效的活性成分的提取都可以通過CytoScape3.7.1 中的“NetworkAnalyzer”分析有效成分及靶點的網絡拓撲參數得到并將結果可視化，見表1 與圖1。

圖1 六君子湯成分-靶點網絡圖

表1 六君子湯有效成分篩選結果

2.2 胃輕癱相關靶點的獲取

由OMIM 及DisGeNET 數據庫根據疾病相關性進行篩選，分別獲得胃輕癱靶點1989、328 個。結合DRUGBANK 數據庫補充相關靶點，合并去除相同的數值，最終得到1902 個胃輕癱相關靶點。

2.3 六君子湯成分-胃輕癱靶點PPI 網絡的構建

將篩選的六君子湯活性成分靶點與胃輕癱靶點取交集，繪制韋恩圖，得到六君子湯成分-胃輕癱共同靶點68 個，見圖2。六君子湯靶點PPI 網絡可以通過STRING11.0 平臺構建，見圖3。通過CytoScape3.7.1 中的hubba 插件對PPI 網絡進一步分析得到潛在的關鍵靶點基因，見圖4。分析得出AKT1 在網絡中的連接度、節點介度 、緊密度分別為51、0.0552、0.7791，推測AKT1 為六君子湯治療胃輕癱的最主要靶點。TNF、

圖2 六君子湯成分-胃輕癱靶點韋恩圖

圖3 六君子湯-胃輕癱靶點PPI 網絡

圖4 六君子湯-胃輕癱靶點PPI 網絡中的Module

IL-6、IL-1β、VEGFA、MAPK3、PPARG、PTGS2、CASP3、MMP9 亦為相對重要的靶點。

由于PPI 網絡中蛋白與蛋白之間的關系是相互作用的，因此可視化結果多以無向圖表示。PPI 復雜網絡中存在部分連線密度較高的模塊，其內部網絡是PPI 網絡的核心子網，是潛在的具有生物學意義的集合[8]。因此，為更精確地分析六君子湯治療胃輕癱的作用機制，在得到六君子湯PPI 網絡后，要應用Metascape 數據平臺中的MCODE 分析功能，通過分子復合物檢測算法對相互作用關系進行分析聚類（見圖4）。分別保留PPI 網絡與聚類模塊中排名前三的生物學進程對其進行功能描述，見表2。

表2 聚類模塊評分前3

2.4 靶點功能與通路的富集分析

把六君子湯治療胃輕癱相關靶點通過Metascape 數據平臺進行信號通路分析，并對結果可視化，見圖5。由結果可見六君子湯主要參與的生物學過程（GO BP）包括對脂多糖的應答、對細胞遷移的正調控、對無機物質的反應、對損傷的應答、凋亡信號通路。參與的通路（KEGG）主要有糖尿病并發癥中的AGE-RAGE 信號通路、流體剪切應力與動脈粥樣硬化、VEGF 信號通路、TNF 信號通路、麻疹等。相關靶點治療胃輕癱的功能（GO MF）主要富集于細胞因子活性、血紅素結合、蛋白質同源二聚化活性、蛋白酶結合、磷酸酶結合等，以及細胞組分（GO CC）的富集，見圖6。

圖5 靶點-通路網絡

圖6 六君子湯主要成分潛在靶點的富集分析

3 討論

在對六君子湯治療胃輕癱作用機制的網絡藥理學研究中對“成分-成分靶點網絡”“胃輕癱相關蛋白相互作用網絡”“成分-胃輕癱靶點蛋白相互作用網絡”進行綜合分析得到，六君子湯治療胃輕癱的10 個潛在核心靶點AKT1、TNF、IL-6、IL-1β、VEGFA、MAPK3、PPARG、PTGS2、CASP3、MMP9。

AKT1 即蛋白激酶B，可通過一系列下游底物的絲氨酸或蘇氨酸磷酸化的介導參與調控細胞的代謝、增殖、細胞存活、生長和血管生成等各種進程。AKT 調控的通路是人類癌癥中最常見的激活信號通路。一半以上的人類腫瘤含有過度活化的AKT[9-10]。VEGFA 即血管內皮生長因子A，是調節脂肪組織中血管生成的關鍵因子。VEGFA 可誘導內皮細胞增殖，促進細胞遷移，抑制細胞凋亡，誘導血管通透性的改變[11]，也會造成血管系統功能紊亂，生成高密度和高通透性的微血管促進腫瘤擴張。腫瘤組織血管生成主要依賴于VEGFA 的應答。

KEGG 通路富集分析顯示，六君子湯治療胃輕癱所調控的核心靶點主要參與AGE-RAGE 信號通路、VEGF信號通路、流體剪切應力與動脈粥樣硬化。GO 富集分析顯示，六君子湯所調控的核心靶點顯著富集對脂多糖的應答、對細胞遷移的正調控、對無機物的應答等生物過程。通過對關鍵通路及靶點的綜合分析，發現PI3KAKT 通路作為子網絡在上述多個關鍵通路中都有很高的參與度，該途徑是AKT 激活的關鍵節點。PI3K-AKT 信號通路被不同種類的細胞刺激或遭遇毒性損傷而激活，參與調節細胞促進生長和分裂，創建代謝反應、輔助繁殖等基本功能[12]。當酪氨酸激酶（RTK）或G 蛋白偶聯受體（GPCR）與生長因子結合，會分別激活PI3K 的兩個異源二聚體，催化細胞膜上PIP2 轉化為PIP3，生成物PIP3 又成為下游通路的第二信使幫助激活AKT，進而通過磷酸化下游底物參與凋亡、蛋白質合成、代謝和細胞周期來控制關鍵的細胞過程[13]。而且該通路為關鍵通路AGE-RAGE 信號通路的下游，通過PI3K-AKT 通路調節NF-κB 的表達可以對TNF、IL-6、IL-1β 等炎癥因子產生雙向調節的作用。

通過復雜網絡分析藥理特性的研究對于準確識別中藥成分的復雜機制具有重要意義。KEGG 通路富集分析表明靶點PTGS2 參與如花生四烯酸、TNF、NF-κB 和VEGF信號通路等炎癥途徑。PTGS2作為一種炎癥介質，是一種限速酶，可催化花生四烯酸產生前列腺素，發揮促炎癥作用并破壞內環境穩態，從而促進白細胞分泌并加重炎癥反應[14]。許多炎癥性疾病的發病機制與花生四烯酸代謝途徑相關。PTGS2 是前列腺素的誘導酶，不僅參與炎癥反應前列腺素的生物合成，而且是該途徑的關鍵酶，這可以作為后續驗證PTGS2 作為關鍵蛋白的有力基礎。

研究結果表明六君子湯作為一個復方中藥制劑含有多種有效成分可以通過多個靶點干預疾病的進程，而同一靶點可在不同的生物學過程及信號通路發揮作用，體現了六君子湯多通路、多靶點聯合作用的復雜性。本文借助網絡藥理思維，構建六君子湯的“組分-靶點-疾病途徑”關聯網絡并預測其潛在的作用機制。使六君子湯治療胃輕癱更好地應用于臨床，為后續評價六君子湯臨床療效指標的挖掘提供了重要線索，也為深入研究六君子湯可能存在的作用機制提供了新的思路和技術手段。但由于網絡藥理學方法僅僅是單純的數據分析，有自身的局限性，加上中藥復方制劑成分較多，在熬制時可能會產生復雜的化學反應，所以熬制后六君子湯的化學成分變化以及活性成分在體內的代謝物并不在本次研究分析范圍內。本研究主要利用數據庫提取數據分析，結合生物信息學，對六君子湯做了初步的功能和調控作用的注釋和預測，后續將在此數據挖掘的基礎上進一步對六君子湯的有效活性成分進行更深入的篩選與通過動物實驗進行驗證，從而進一步明晰六君子湯的主要調控靶點。