王蘊(yùn)涵, 杜群, 李燕舞, 肖蘇婷, 喬楊
(廣州中醫(yī)藥大學(xué)科技創(chuàng)新中心,廣東廣州 510405)
急性肺損傷(acute lung injury)是由感染、創(chuàng)傷、有害氣體吸入、休克、中毒等多種因素所導(dǎo)致的以肺彌散功能障礙為特征的肺部疾病[1],在臨床中病死率較高,常危及生命[2]。目前,除了傳統(tǒng)的治療原發(fā)病、呼吸支持及使用糖皮質(zhì)激素外,國內(nèi)外關(guān)于其生理病理機(jī)制的研究也在不斷開展,旨在為急性肺損傷提供更好的治療方法。中醫(yī)學(xué)并無“急性肺損傷”這一疾病名稱,而是根據(jù)其發(fā)病時的癥狀、體征將其歸屬于“喘脫”“結(jié)胸”“暴喘”等病證的范疇[3]。臨床上急性肺損傷患者多為陽明燥實(shí),熱毒結(jié)于大腸,腑氣不通,形成陽明腑實(shí)之證,基于中醫(yī)藏象理論“肺與大腸相表里”“肺合大腸”,臨床上多采用肺腸同治之法[4]。大承氣湯出自《傷寒論》,是肺腸同治法臨床應(yīng)用的代表方劑,被廣泛應(yīng)用于急性肺損傷的治療,效果顯著[5-8];但其作用機(jī)制尚不清楚。因此,本研究通過網(wǎng)絡(luò)藥理學(xué)的方法初步探討大承氣湯治療急性肺損傷的主要活性成分及潛在靶點(diǎn),并加以體外實(shí)驗(yàn)驗(yàn)證,以期闡明其作用機(jī)理?,F(xiàn)將研究結(jié)果報道如下。
1.1 網(wǎng)絡(luò)藥理學(xué)預(yù)測
1.1.1 大承氣湯藥物化學(xué)成分收集及藥物靶點(diǎn)預(yù)測登錄中藥系統(tǒng)網(wǎng)絡(luò)藥理學(xué)數(shù)據(jù)庫選擇TCMSPTM(網(wǎng)址http:∕∕ibts.hkbu.edu.hk∕LSP∕tcmsp.php),設(shè)置搜索條目為Herb name。設(shè)置篩選條件為口服生物利用度(OB)≥30%、類藥性指數(shù)(DL)≥0.18,篩選大承氣湯的目標(biāo)靶點(diǎn)(targets information)。數(shù)據(jù)庫所獲得的藥物有效化學(xué)成分及靶點(diǎn)信息應(yīng)用Drugbank(https:∕∕www.drugbank.ca∕)使 其 格 式 標(biāo) 準(zhǔn)化。因TSMSP數(shù)據(jù)庫中尚未收錄芒硝的藥物成分及靶點(diǎn),登錄BATMAN數(shù)據(jù)庫(網(wǎng)址:http:∕∕bionet.ncpsb.org.cn),在搜索條目下以“mang xiao”為關(guān)鍵詞搜索,收集芒硝的藥物成分靶點(diǎn)。
1.1.2 疾病靶點(diǎn)來源登錄GeneCards人類基因數(shù)據(jù)庫[9](THE HUMAN GENE DATABASE),網(wǎng)址https:∕∕www.genecards.org∕。設(shè)置搜索條目為關(guān)鍵詞搜索人類孟德爾遺傳數(shù)據(jù)庫(Online Mendelian Inheritance in Man)[10],網(wǎng) 址https:∕∕www.ncbi.nlm.nih.gov∕omim。搜索條目下以“acute lung injury”為關(guān)鍵詞篩選,設(shè)置物種為人(Homo sapiens),獲取急性肺損傷的疾病靶點(diǎn)。對數(shù)據(jù)庫檢索到的疾病靶點(diǎn)進(jìn)行合并去重,作為最終的疾病靶點(diǎn)來源。
1.1.3 復(fù)方-疾病交集基因的獲取應(yīng)用R語言中的在線程序Draw Venn Diagrams對“1.1.1”項(xiàng)與“1.1.2”項(xiàng)中獲取到的大承氣湯藥物靶點(diǎn)與急性肺損傷疾病靶點(diǎn)進(jìn)行匹配,獲得交集基因靶點(diǎn),作為大承氣湯治療急性肺損傷的潛在作用靶點(diǎn)并繪制韋恩圖。
1.1.4 化學(xué)成分-靶點(diǎn)網(wǎng)絡(luò)圖的構(gòu)建將“1.1.3”項(xiàng)中獲取到的交集基因?qū)氲骄W(wǎng)絡(luò)可視化軟件Cytoscape 3.8.0,構(gòu)建化學(xué)成分-靶點(diǎn)作用網(wǎng)絡(luò)。對網(wǎng)絡(luò)進(jìn)行分析,以度值參數(shù)作為評價節(jié)點(diǎn)在網(wǎng)絡(luò)中重要性的標(biāo)準(zhǔn)。
1.1.5 大承氣湯治療急性肺損傷的蛋白相互作用(PPI)網(wǎng)絡(luò)構(gòu)建及核心蛋白篩選登錄String數(shù)據(jù)庫,網(wǎng) 址https:∕∕string-db.org∕。選 擇Multiple proteins,在搜索條目下輸入“1.1.3”項(xiàng)獲取到的交集基因,限定物種為人(Homo Sapiens),其余設(shè)置保持不變,構(gòu)建PPI功能網(wǎng)絡(luò),并繪制網(wǎng)絡(luò)圖。為篩選核心蛋白,將生成的PPI條目Node1、Node2以及Combined-scorce導(dǎo)入到Cytoscape3.8.0,應(yīng)用Network Analyse插件分析網(wǎng)絡(luò),選擇APP條目下的CytoNCA插件進(jìn)行網(wǎng)絡(luò)拓?fù)浞治觥?/p>
1.1.6 基因本體論(GO)和京都基因與基因組百科(KEGG)通路富集分析基于R語言對“1.1.3”項(xiàng)中獲取到的交集基因的細(xì)胞成分、分子功能、生物學(xué)過程以及信號通路進(jìn)行富集分析。首先將基因Symbol轉(zhuǎn)化為基因ID,選擇物種為人(Homo sapiens),設(shè)定閾值為P<0.05,輸出富集結(jié)果并繪制條形圖與氣泡圖。
1.2 體外實(shí)驗(yàn)驗(yàn)證
1.2.1 藥物、試劑與儀器木犀草素(luteolin)[高效液相色譜(HPLC)>98%,成都普菲德生物技術(shù)有限公司];蘆薈大黃素(aloe-emodin)(HPLC>97%,成都普菲德生物技術(shù)有限公司);柚皮素(naringenin)(HPLC>99%,成都普菲德生物技術(shù)有限公司)。脂多糖(lipopolysaccharide,LPS)(大腸桿菌055:B5)(美國Sigma-Aldrich公司);四甲基偶氮唑鹽(MTT)試劑盒、小鼠腫瘤壞死因子α(tumor necrosis factorα,TNF-α)酶聯(lián)免疫吸附分析(ELISA)試劑盒、小鼠白細(xì)胞介素6(interleukin 6,IL-6)ELISA試劑盒(碧云天生物技術(shù)公司);P53兔多克隆抗體(Protein-Tech Group公司);GAPDH兔多克隆抗體、辣根過氧化物酶(HRP)偶聯(lián)的兔免疫球蛋白(IgG)抗體(ABCAM公司)。酶標(biāo)儀(美國Thermo Fisher公司)。
1.2.2 細(xì)胞培養(yǎng)本實(shí)驗(yàn)選用小鼠單核巨噬細(xì)胞RAW264.7(購自上海中科院細(xì)胞庫)。在實(shí)驗(yàn)開始前先把RAW264.7細(xì)胞放置于37℃、體積分?jǐn)?shù)5%CO2的培養(yǎng)箱中,給予體積分?jǐn)?shù)10%胎牛血清(FBS)DMEM培養(yǎng)基處理以提供必要的生長條件。培養(yǎng)2~3 d達(dá)到90%聚集,再以1∶5的比例傳代培養(yǎng)用于實(shí)驗(yàn)。
1.2.3 MTT實(shí)驗(yàn)本研究采用MTT法分析評估木犀草素、蘆薈大黃素、柚皮素對RAW264.7細(xì)胞活力的影響,并由此確定合理的藥物濃度。首先將細(xì)胞以1×105個∕孔的密度接種在96孔板中,待細(xì)胞貼壁后,培養(yǎng)12 h,然后用含有不同濃度木犀草素、蘆薈大黃素、柚皮素(0、1.562 5、3.125、6.25、12.5、25、50、100μmol∕L)的培養(yǎng)基處理細(xì)胞24 h。接著,將20μL MTT(5 mg∕mL)添加到每個孔中,于37℃中溫育4 h后棄去培養(yǎng)基,于每個孔中加入100μL二甲基亞砜(DMSO)溶解已經(jīng)形成的甲瓚晶體,應(yīng)用酶標(biāo)儀,選擇最佳測量波長570 nm,測量每個孔的吸光度。
1.2.4 ELISA實(shí)驗(yàn)本研究采用ELISA法檢測細(xì)胞培養(yǎng)基上清TNF-α和IL-6水平。根據(jù)實(shí)驗(yàn)需求將RAW264.7細(xì)胞以3×106個∕孔的密度接種在6孔板中,給予與“1.2.3”項(xiàng)相同濃度、時間的藥物刺激,收集RAW264.7細(xì)胞培養(yǎng)基上清,根據(jù)說明書方法進(jìn)行操作。最后應(yīng)用酶標(biāo)儀,選擇最佳測量波長450 nm,測量每個孔的吸光度。
1.2.5 蛋白免疫印跡(Western Blot)實(shí)驗(yàn)同“1.2.4”項(xiàng)的細(xì)胞培養(yǎng)及刺激方式。去掉RAW264.7細(xì)胞培養(yǎng)基上清,使用預(yù)冷的PBS洗滌2次,加入結(jié)合蛋白酶抑制劑混合物的放射免疫沉淀分析(RIPA)緩沖液裂解處理后的RAW264.7細(xì)胞,提取總蛋白。通過8%~12%的十二烷基硫酸鈉-聚丙烯酰胺凝膠電泳(SDS-PAGE)分離蛋白,并將蛋白轉(zhuǎn)移至聚偏氟乙烯(PVDF)膜,用50 g∕L脫脂牛奶封閉。將PVDF膜與P53抗體4℃孵育過夜后,再與HRP標(biāo)記的IgG抗體孵育。應(yīng)用GE ImageQuant LAS4000mini(美國GEHealthcare公司)觀察目的條帶,并應(yīng)用Gel Doc XR系統(tǒng)(美國Bio-Rad Laboratories)采用光密度法定量,最后應(yīng)用Quantity One蛋白質(zhì)分析軟件(美國Bio Rad公司)進(jìn)一步分析。
2.1 復(fù)方化學(xué)成分及藥物靶點(diǎn)本次研究經(jīng)過合并去重,共收集大承氣湯藥物化學(xué)成分27個,預(yù)測靶點(diǎn)82個。
2.2 疾病靶點(diǎn)的獲取在GeneCards、OMIM數(shù)據(jù)庫中以“acute lung injury”為關(guān)鍵詞進(jìn)行檢索,通過R軟件的venn程輯包將2個數(shù)據(jù)庫中所得到的疾病靶點(diǎn)匯總并刪除重復(fù)。為保證數(shù)據(jù)的準(zhǔn)確性,GeneCards數(shù)據(jù)庫按照Relevance Score>10篩選,共獲得2 801個疾病靶點(diǎn),通過UniProt數(shù)據(jù)庫對疾病靶基因進(jìn)行校正,將其作為候選靶標(biāo)來源。
2.3 獲取復(fù)方-疾病交集基因基于R語言3.6.0軟件,在R語言中Bioconductor(https:∕∕www.bioconductor.org∕)網(wǎng)站安裝VennDiagram程輯包,將“1.1.1”項(xiàng)與“1.1.2”項(xiàng)中獲得的復(fù)方藥物靶點(diǎn)與疾病靶點(diǎn)數(shù)據(jù)導(dǎo)入到R軟件,應(yīng)用Draw Venn Diagrams在線程序獲取交集基因并繪制韋恩圖,見圖1。共獲得復(fù)方-疾病交集基因55個。
圖1 大承氣湯-急性肺損傷交集基因韋恩圖Figure 1 Venn diagram of intersection between Dachengqi Decoction-acute lung injury
2.4 構(gòu)建有效成分-靶點(diǎn)網(wǎng)絡(luò)圖運(yùn)用Cytoscape 3.8.0軟件中Import Network插件構(gòu)建復(fù)方有效化學(xué)成分-靶點(diǎn)作用網(wǎng)絡(luò),見圖2。網(wǎng)絡(luò)共有145個節(jié)點(diǎn),297條邊。核心靶標(biāo)顯示為方形,活性化學(xué)成分顯示為圓形,不同藥物顯示為不同顏色,其中,大黃顯示為紅色,芒硝顯示為綠色,枳實(shí)顯示為紫色,厚樸顯示為藍(lán)色。網(wǎng)絡(luò)中化合物與靶點(diǎn)的關(guān)系用每條邊來顯示。一個化學(xué)成分與多個交集基因靶點(diǎn)對應(yīng),且一個交集基因與靶點(diǎn)也同時對應(yīng)多個化學(xué)成分,化學(xué)成分又從屬于不同的中藥,這說明大承氣湯治療急性肺損傷的過程是多成分、多靶點(diǎn)的。
圖2 大承氣湯化學(xué)成分-靶點(diǎn)作用網(wǎng)絡(luò)Figure 2 Active ingredient-target network of Dachengqi Decoction
2.5 核心蛋白相互作用網(wǎng)絡(luò)構(gòu)建及核心蛋白篩選應(yīng)用String數(shù)據(jù)庫對交集基因初步構(gòu)建蛋白作用網(wǎng)絡(luò),見圖3。節(jié)點(diǎn)代表蛋白,邊代表蛋白與蛋白之間的作用關(guān)系。網(wǎng)絡(luò)連接程度越高說明蛋白之間關(guān)系越密切。將String數(shù)據(jù)庫得到的string_interactions.txt文件導(dǎo)入到Cytoscape 3.8.0軟件,應(yīng)用Network Analyse插件構(gòu)建網(wǎng)絡(luò),選擇CytoNCA插件進(jìn)行網(wǎng)絡(luò)拓?fù)浞治?,篩選網(wǎng)絡(luò)中的關(guān)鍵節(jié)點(diǎn)。核心蛋白篩選參數(shù):介度中心性(Betweenness Centrality),緊密中心性(Closeness Centrality),節(jié)度中心性(Degree Centrality),特征向量(Eigenvector Centrality),局 部 連 通 性(Local Average Connectivity-based method,LAC)。篩選條件:應(yīng)用R軟件的在線程輯包,對各交集基因上述核心蛋白篩選參數(shù)數(shù)值進(jìn)行測算,各組數(shù)值中大于本組中位值的基因予以保留,運(yùn)行R腳本,篩選核心蛋白,見表1,篩選流程見圖4。
圖3 大承氣湯與急性肺損傷的PPI網(wǎng)絡(luò)圖Figure 3 PPInetwork diagram of Dachengqi Decoction and acute lung injury
圖4 CytoNca核心蛋白基因篩選流程圖Figure 4 Screening flowchart of CytoNca core protein genes
表1 核心蛋白拓?fù)浞治鯰able 1 Topological analysis of core proteins
2.6 交集基因的GO功能富集分析基于R軟件3.6.0對大承氣湯治療急性肺損傷的55個作用靶點(diǎn)進(jìn)行GO功能富集分析。富集內(nèi)容包括細(xì)胞成分(celluar component,CC)、生物學(xué)過程(biological process,BP)、分子功能(molecular function,MF)。富集條件為P值=0.05且Q值=0.05。通過分析發(fā)現(xiàn),富集到了1 957個GO條目,細(xì)胞成分83個,生物學(xué)過程1 714個,分子功能160個。選取各排名前10位的富集條目,并根據(jù)每個項(xiàng)目的P值、Q值及富集在其上的基因數(shù)目繪制條形圖,見圖5。橫坐標(biāo)表示靶點(diǎn)數(shù),左邊表示BP、CC、MF,顏色表示P值,P值越小顏色越偏向紅色,P值越大則越偏向藍(lán)色。
圖5 大承氣湯治療急性肺損傷的GO功能富集分析Figure 5 GO function enrichment of Dachengqi Decoction in treating acute lung injury
2.7 交集基因的KEGG通路富集分析基于R語言對交集基因進(jìn)行KEGG通路富集分析。分析結(jié)果可見,55個復(fù)方疾病交集基因共富集到了126條通路(P<0.05)。選取排名前20位的富集通路,根據(jù)每個通路的P值、Q值及富集在其上的基因數(shù)目繪制氣泡圖,見圖6。橫坐標(biāo)表示富集到基因的數(shù)量,左邊表示通路名稱,顏色表示P值,P值越小顏色越偏向紅色,P值越大則越偏向藍(lán)色。主要涉及到P53信號通路。
圖6 大承氣湯與急性肺損傷交集基因KEGG通路富集圖Figure 6 KEGG pathway enrichment of Dachengqi Decoction in treating acute lung injury
2.8 木犀草素、蘆薈大黃素、柚皮素能降低LPS誘導(dǎo)的RAW264.7細(xì)胞中TNF-α、IL-6細(xì)胞因子水平MTT的數(shù)據(jù)結(jié)果表明,相較于正常對照組(Ctrl,藥物濃度為0μmol∕L),濃度小于100μmol∕L的柚皮素組和蘆薈大黃素組以及濃度小于12.5μmol∕L的木犀草素組對RAW264.7細(xì)胞沒有毒性,具體結(jié)果見圖7。藥效實(shí)驗(yàn)的結(jié)果表明,柚皮素和蘆薈大黃素從10μmol∕L濃度,木犀草素從2μmol∕L濃度開始,具有一定抑制TNF-α、IL-6細(xì)胞因子分泌的能力,且細(xì)胞因子的分泌水平呈劑量依賴性地降低,具體結(jié)果見圖8。
圖7 木犀草素、蘆薈大黃素、柚皮素對RAW264.7細(xì)胞活性的影響Figure 7 Effects of luteolin,aloe-emodin and naringenin on RAW264.7 cell activity
圖8 木犀草素、蘆薈大黃素、柚皮素對RAW264.7細(xì)胞TNF-α和IL-6分泌的影響Figure 8 Effects of luteolin,aloe-emodin and naringin on secretion of TNF-αand IL-6 in RAW264.7 cells
2.9 木犀草素、蘆薈大黃素、柚皮素降低LPS誘導(dǎo)的RAW264.7細(xì)胞中P53蛋白表達(dá)根據(jù)網(wǎng)絡(luò)藥理學(xué)預(yù)測結(jié)果,TP53基因在構(gòu)建的成分-靶點(diǎn)-通路網(wǎng)絡(luò)中,度值較大且控制著P53蛋白的表達(dá)。因此,本研究檢測P53蛋白表達(dá),對網(wǎng)絡(luò)藥理學(xué)預(yù)測靶點(diǎn)進(jìn)行驗(yàn)證。Western Blot實(shí)驗(yàn)結(jié)果顯示,與模型組比較,木犀草素、蘆薈大黃素、柚皮素的干預(yù)能顯著降低P53蛋白的表達(dá),且與給藥劑量呈現(xiàn)一定的量-效關(guān)系,具體結(jié)果見圖9。
本研究基于網(wǎng)絡(luò)藥理學(xué)的方法探討大承氣湯治療急性肺損傷的作用機(jī)制。根據(jù)OB和DL參數(shù),共篩選出大承氣湯化學(xué)成分27個,預(yù)測靶點(diǎn)經(jīng)合并去重收集82個。構(gòu)建有效成分-靶點(diǎn)網(wǎng)絡(luò)圖,結(jié)果表明,木犀草素(MOL000006 luteolin)、蘆薈大黃素(MOL000471 aloe-emodin)、柚皮素(MOL004328 naringenin)等是其主要的化學(xué)成分。既往的藥理學(xué)研究表明,木犀草素抗炎作用較強(qiáng),能抑制IL-6、TNF-α等炎癥細(xì)胞因子的表達(dá),并且能夠提高IL-10的表達(dá)水平來調(diào)控炎癥反應(yīng)[11]。周曉玲[12]研究發(fā)現(xiàn),蘆薈大黃素具有體外抗炎作用,能夠以劑量依賴的方式降低LPS刺激的RAW264.7細(xì)胞中NO、TNF-a、IL-1β的釋放。柚皮素具有較強(qiáng)的抗炎及抗氧化作用,薛傳優(yōu)等[13]研究發(fā)現(xiàn),柚皮素可以降低IL-1β、IL-6和TNF-α基因的表達(dá)水平,改善脂多糖誘導(dǎo)的H9c2心肌細(xì)胞的炎癥反應(yīng)和細(xì)胞凋亡。
經(jīng)過對55個交集基因蛋白進(jìn)行核心蛋白篩選,結(jié)果顯示TP53、JUN、MYC等是大承氣湯治療急性肺損傷起到重要作用的核心蛋白。Laurent等[14]發(fā)現(xiàn)TP53表達(dá)細(xì)胞的比例與炎癥強(qiáng)度相關(guān)。SONG等[15]研究發(fā)現(xiàn),c-Jun能夠激活NLRP3炎癥小體,促進(jìn)其磷酸化,介導(dǎo)炎癥反應(yīng)。黃杰中等[16]研究發(fā)現(xiàn),c-Myc與凋亡抑制因子呈正相關(guān),表明其可以顯著地抑制細(xì)胞凋亡。急性肺損傷的發(fā)病往往繼發(fā)于肺部疾患的炎癥反應(yīng),伴隨著肺部毛細(xì)血管的損傷,引起細(xì)胞凋亡??梢钥闯鯰P53、JUN、MYC均與急性肺損傷的發(fā)生發(fā)展具有一定的相關(guān)性。
交集基因KEGG通路富集分析顯示,其通路主要為P53信號通路。Wu等[17]的研究表明,通過抑制P53相關(guān)信號通路可以減輕類風(fēng)濕關(guān)節(jié)炎炎癥反應(yīng)。炎癥信號通路NF-κB蛋白表達(dá)水平亦可以影響凋亡因子P53的表達(dá)[18]。P53被認(rèn)為是NF-κB調(diào)控的細(xì)胞凋亡因子,可以直接或間接影響參與細(xì)胞呼吸的多種酶的表達(dá),并通過增加氧化磷酸化促進(jìn)活性氧的產(chǎn)生[19]。也有研究認(rèn)為,P53是調(diào)控細(xì)胞應(yīng)答多種胞外信號的重要轉(zhuǎn)錄因子,可以增加促炎因子的產(chǎn)生[20]。
本研究通過細(xì)胞實(shí)驗(yàn)證實(shí),木犀草素、蘆薈大黃素、柚皮素均具有減少巨噬細(xì)胞炎癥細(xì)胞因子TNF-α、IL-6釋放的作用。Western Blot結(jié)果顯示,3種活性成分均可劑量依賴性地降低P53蛋白的表達(dá),對機(jī)體的炎癥反應(yīng)具有一定的抑制作用。
綜上所述,本研究將網(wǎng)絡(luò)藥理學(xué)和生物信息學(xué)相結(jié)合,從網(wǎng)絡(luò)層面系統(tǒng)地分析預(yù)測了大承氣湯治療急性肺損傷的化合物-靶點(diǎn)-通路作用機(jī)制,使用細(xì)胞實(shí)驗(yàn)和分子生物學(xué)方法驗(yàn)證了其藥效與分子機(jī)制,這與中醫(yī)的整體觀念和辨證論治的理念相符合,為今后的相關(guān)研究和實(shí)驗(yàn)設(shè)計提供了參考和新思路。