陳江鋒 林芝嫻 郭 勇

癌性惡病質(zhì)被國際共識小組定義為一種多因素綜合征，其特征是骨骼肌質(zhì)量持續(xù)下降（伴或不伴脂肪減少），不能通過常規(guī)營養(yǎng)支持完全逆轉(zhuǎn)，導(dǎo)致進行性損害[1]。50%～80%的癌癥患者最終會發(fā)展為惡病質(zhì)，進一步加速了癌癥患者的死亡。因此，癌性惡病質(zhì)的治療一直是惡性腫瘤綜合治療中的重要組成部分。四君子湯被廣泛應(yīng)用于治療癌性惡病質(zhì)，其對疲勞、食欲不振、體重減輕等癥狀有較好的療效[2-3]。研究發(fā)現(xiàn)，四君子湯可以顯著降低血清中白細胞介素-1（interleukin 1，IL-1）、白細胞介素-6（interleukin 6，IL-6）、腫瘤壞死因子-α（tumor necrosis factor-α，TNF-α）、胰島素樣生長因子1（insulin-like growth factor 1，IGF-1）的水平，減緩癌性惡病質(zhì)骨骼肌萎縮的發(fā)生，改善癌性惡病質(zhì)患者的營養(yǎng)狀況，延緩疾病的進一步發(fā)展[4-6]。然而目前四君子湯對癌性惡病質(zhì)的具體作用機制尚不清楚，鑒于中醫(yī)治療的整體觀念和中藥方劑多化合物、多靶點、多通路協(xié)同作用的特點，本研究采用四個連續(xù)模塊（靶點預(yù)測、網(wǎng)絡(luò)構(gòu)建、分子對接、GO 和KEGG 富集分析）組成的綜合網(wǎng)絡(luò)藥理學(xué)方法探討四君子湯治療癌性惡病質(zhì)的潛在機制。

1 方 法

1.1 靶點預(yù)測

1.1.1 四君子湯的化合物 查詢中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫（TCMSP，http://lsp.nwu.edu.cn/tcmsp.php）和中醫(yī)藥綜合數(shù)據(jù)庫（TCMID，http://bidd.group/group2017/index.html），參照TCMSP 數(shù)據(jù)庫中常用的標(biāo)準(zhǔn)：口服生物利用度（OB）≥30%，藥物相似性（DL）≥0.18，Caco-2≥0，通過Pubchem（https://pubchem.ncbi.nlm.nih.gov/）搜索每個活性化合物的化學(xué)結(jié)構(gòu)，剔除沒有結(jié)構(gòu)信息的化合物。最終共搜索到105 個活性化合物。

1.1.2 四君子湯的復(fù)合靶點 Pharmmapper server（http://www.lilab-ecust.cn/pharmmapper/）是一個持續(xù)更新的綜合藥效基因匹配平臺。將105 種化合物的mol2 格式文件上傳至Pharmmapper server，generate confomers 設(shè)置為No，maximum generated conformations 設(shè)置為300。在高級選項中，選擇Druggable Pharmacophore Models（v2017，16159）作為靶點集，其他選項設(shè)置為默認，獲得預(yù)測的靶點。將搜索到的靶點導(dǎo)入UniProt 數(shù)據(jù)庫（http://stitch.embl.de），選擇的物種為Homo sapiens，去除重復(fù)、非人、非標(biāo)靶，通過檢索和轉(zhuǎn)化，最終得到四君子湯中活性化合物的450 個靶基因。

1.1.3 癌性惡病質(zhì)的靶點 通過DrugBank database（https://www.drugbank.ca），Online Mendelian Inheritance in Man database（OMIM，https://omim.org），Human gene database（GeneCards，https://www.genecards.org）和Therapeutic Target Database（TTD，https://db.idrblab.org/ttd/）搜集癌性惡病質(zhì)的相關(guān)靶點。以“Cancer cachexia”為關(guān)鍵詞進行搜索，刪除重復(fù)靶點，最終搜集到1189 個靶點。

1.1.4 蛋白相互作用網(wǎng)絡(luò)構(gòu)建 用Venny 2.1 軟件（http://bioinfogp.cnb.csic.es/tools/venny/）將四君子湯的活性化合物靶點與癌性惡病質(zhì)相關(guān)靶點進行匹配，選擇重疊靶點作為四君子湯治療癌性惡病質(zhì)的靶點，使用String（https://string-db.org/）對相關(guān)靶點進行處理，蛋白類型設(shè)置為Homo sapiens，minimum required interaction score 設(shè)置為“medium confidence”（>0.4），其他參數(shù)設(shè)置為默認，獲得最終的蛋白相互作用網(wǎng)絡(luò)。

1.2 網(wǎng)絡(luò)構(gòu)建

1.2.1 網(wǎng)絡(luò)構(gòu)建的方法 Cytoscape 是一個開放性平臺，用于可視化復(fù)雜網(wǎng)絡(luò)并能將其與任何屬性的數(shù)據(jù)進行處理。將化合物和蛋白相互作用網(wǎng)絡(luò)的數(shù)據(jù)輸入到Cytoscape 中，并構(gòu)建不同的網(wǎng)絡(luò):復(fù)合靶點網(wǎng)絡(luò)和化合物-靶點網(wǎng)絡(luò)。

1.2.2 拓撲分析 通過拓撲參數(shù)分析，對復(fù)合靶點網(wǎng)絡(luò)進行分析篩選。通過度中心性（degree centrality，DC）、接近中心性（closeness centrality，CC）、介度中心性（betweeness centrality，BC）三個參數(shù)對重要靶點進行評價。這三個參數(shù)的大小代表了網(wǎng)絡(luò)中節(jié)點的重要性。重要節(jié)點表示網(wǎng)絡(luò)中輸出值高。篩選標(biāo)準(zhǔn)為DC、CC、BC 大于1 倍中位數(shù)。最后，篩選出重要靶點。

1.3 分子對接 通過分子對接驗證四君子湯中化合物與重要靶點的結(jié)合活性，篩選出潛在的活性成分。我們之前搜集到的105 個化合物的小分子結(jié)構(gòu)都是從Pubchem 數(shù)據(jù)庫中獲取的，并保存為mol2 格式文件。所有具有重要靶點蛋白晶體結(jié)構(gòu)的蛋白配體均從RCSB 蛋白數(shù)據(jù)庫（https://www.rcsb.org/）中獲得。SystemsDock（http://systemsdock.unit.oist.jp/iddp/home/index）是一個基于網(wǎng)絡(luò)藥物學(xué)預(yù)測和分析的網(wǎng)絡(luò)服務(wù)器，使用高精度的對接模擬和分子通路圖來充分表達配體的選擇性，揭示配體在復(fù)雜分子網(wǎng)絡(luò)中的作用。因此，我們將這兩種數(shù)據(jù)導(dǎo)入SystemsDock 中進行分子對接，并以對接得分作為評價標(biāo)準(zhǔn)。

1.4 GO 和KEGG 富集分析 將用于治療癌性惡病質(zhì)的四君子湯的相關(guān)靶點導(dǎo)入到DAVID 數(shù)據(jù)庫中（DAVID，http://david.abcc.ncifcrf.gov/）。所有的靶基因都被導(dǎo)入到基因列表中，物種被定義為人類。通過上述數(shù)據(jù)庫檢索和轉(zhuǎn)換，進行GO 生物過程和KEGG代謝通路富集分析。P<0.05 為顯著的生物過程。使用在線作圖網(wǎng)站Omishare Tools（http://www.omicshare.com/tools/index.php/）將結(jié)果可視化。

2 結(jié)果

2.1 靶點預(yù)測 經(jīng)過篩選和刪除，在TCMSP 和TCMID 數(shù)據(jù)庫中搜集到105 個四君子湯的化合物。其中黨參16 個，白術(shù)2 個，茯苓13 個，甘草74 個。通過Pharmmapper server 對105 個化合物的靶點進行搜索，經(jīng)UniProt 數(shù)據(jù)庫處理后得到450 個潛在化合物靶點。通過檢索DrugBank、OMIM、GeneCards 和TTD 數(shù)據(jù)庫，共獲得1189 個癌性惡病質(zhì)的相關(guān)靶點。最后，通過Venny 2.1 軟件獲得藥物化合物與疾病的52 個交集靶點（見圖1）。

圖1 四君子湯作用靶點和疾病靶點的Venny 圖

2.2 網(wǎng)絡(luò)構(gòu)建 將蛋白相互作用網(wǎng)絡(luò)的數(shù)據(jù)導(dǎo)入Cytoscape 軟件，以確定靶蛋白之間的相互作用網(wǎng)絡(luò)。圖中的節(jié)點表示靶蛋白，邊表示靶蛋白之間的相關(guān)性。剔除相關(guān)性較低的節(jié)點，包括50 個節(jié)點和140條邊。拓撲分析的主要參數(shù)為度、接近中心性和介度中心性。初次篩選度為3.5，介度中心性為0.007233，接近中心性為0.404986，得到17 個節(jié)點、54 條邊。第二次篩選度是6，介度中心性為0.030593，接近中心性為0.445491，最后得到了8 個節(jié)點和25 條邊（見圖2）。在蛋白相互作用網(wǎng)絡(luò)中，根據(jù)度的大小，INS（25），ALB（25），VEGFA（16），CCND1（15），CTNNB1（14），CYCS（13），AR（12）和ACTA1（8）被認為是重要的靶點（見圖3）。以這8 個重要靶點為基礎(chǔ)，構(gòu)建化合物-重要靶點網(wǎng)絡(luò)（見圖4），這個網(wǎng)絡(luò)包括96 種化合物，其中ZINC00519174（PubChem CID：928837）、Gancaonin H（PubChem CID：5481949）、Frutinone A（PubChem CID：441965）和DFV（Pub－Chem CID：114829）與5 個重要靶點均相關(guān)。

圖2 蛋白相互作用網(wǎng)絡(luò)的拓撲篩選

圖3 四君子湯治療癌性惡病質(zhì)的潛在靶點相互作用網(wǎng)絡(luò)

圖4 化合物-重要靶點網(wǎng)絡(luò)

2.3 分子對接

2.3.1 結(jié)合活性驗證 選取8 個重要靶點（INS、ALB、VEGFA、CCND1、CTNNB1、CYCS、AR 和 ACTA1），導(dǎo)入SystemsDock，與105 個四君子湯的化合物進行對接分析。超過7.0 分的化合物有168 個（20%），5.0～7.0 分的化合物有502 個（59.8%），4.25～5.0 分的化合物有68 個（8.1%）。共有102 種化合物（12.1%）的對接分?jǐn)?shù)小于4.25。對接分?jǐn)?shù)大于7.0 說明該化合物對靶點具有較強的結(jié)合活性，大于5.0 說明該化合物對靶點具有較好的結(jié)合活性，大于4.25說明該化合物對靶點具有一定的結(jié)合活性。因此，四君子湯的活性化合物與8 個重要靶點具有良好的結(jié)合活性。

2.3.2 核心化合物的篩選 對接結(jié)果表明，與8 個靶點對接得分最高的化合物為ergosterol peroxide（PubChem CID：5351516，和CYCS），glyasperin B（PubChem CID：480784，和VEGFA），dehydroeburicoic acid（PubChem CID：15250826，和 ALB），glyasperin F（PubChem CID：392442，和INS），ergosterol peroxide（PubChem CID：5351516，和AR），taraxerol（PubChem CID：92097，和ACTA1），7，9（11）-dehydropachymic acid（PubChem CID：15226717，和CCND1）和 5' -prenylbutein（PubChem CID：11267805，和CTNNB1）。使用PyMOL 軟件對這8 對靶點化合物進行優(yōu)化，對接的特異性蛋白配體相互作用如圖5 所示。8 個靶點得分均在6 分以上的化合物為7，9（11）-dehydropachymic acid（PubChem CID：15226717），taraxerol（PubChem CID：92097），Cerevisterol（PubChem CID：10181133），hederagenin（Pub-Chem CID：73299），polyporenic acid C（PubChem CID：9805290），16 alpha-hydroxydehydrotrametenolic acid（PubChem CID：10743008），pachymic acid（Pub-Chem CID：5484385）和poricoic acid C（PubChem CID：56668247）。滿足上述兩個條件的化合物[7，9（11）-dehydropachymic acid 和taraxerol]被認為是可能通過作用于重要靶點而影響癌性惡病質(zhì)的核心化合物。

圖5 蛋白質(zhì)-配體的對接模擬

2.4 GO 和KEGG 富集分析

2.4.1 GO 富集分析 用DAVID 工具對四君子湯治療癌性惡病質(zhì)的52 個靶點進行GO 富集分析，進一步明確交集靶點的可能作用。生物過程（biological process）、分子功能（molecular function）和細胞組分（cellular component）的比對條目見圖6。

圖6 四君子湯治療癌性惡病質(zhì)潛在靶點的生物學(xué)功能分析

2.4.2 KEGG 通路富集分析 通過KEGG 通路富集分析篩選出四君子湯對癌性惡病質(zhì)中36 條生物通路有顯著影響（見表1）。其作用機制主要涉及pathways in cancer，HTLV-I infection，viral carcinogenesis，human cytomegalovirus infection，PI3K-Akt signaling pathway，human papillomavirus infection 和thyroid hormone synthesis 等（見圖7）。對36 條信號通路進行分析和整合，繪制出最終的通路圖（見圖8）。

表1 四君子湯靶點的通路分析

圖7 四君子湯治療癌性惡病質(zhì)的KEGG 富集關(guān)鍵靶點氣泡圖

圖8 四君子湯治療癌性惡病質(zhì)的潛在機制

3 討論

在臨床實踐中，四君子湯對癌性惡病質(zhì)的治療有顯著作用，可以減輕化療引起的毒性或肌肉萎縮[5，7-9]。本文從網(wǎng)絡(luò)藥理學(xué)的角度探討了四君子湯治療癌性惡病質(zhì)的潛在機制。

3.1 調(diào)節(jié)代謝 惡性腫瘤可通過自身代謝產(chǎn)生的有毒物質(zhì)或通過改變機體代謝而引起惡病質(zhì)，主要表現(xiàn)為葡萄糖代謝，脂質(zhì)代謝和蛋白質(zhì)代謝三個方面的異常。四君子湯通過其化合物如glycitein（Pub－Chem CID：5317750），isorhamnetin（PubChem CID：5281654），quercetin（PubChem CID：5280343），11-hydroxyrankinidine（PubChem CID：5318332）和AC1LCTJ1（PubChem CID：637112），可刺激INS 與其受體結(jié)合，調(diào)控PI3K-Akt 信號通路，上調(diào)Akt2 的表達，影響下游靶點，抑制PCK2 的表達，最終抑制胰島素抵抗[10-11]，調(diào)節(jié)糖代謝，抑制癌性惡病質(zhì)的發(fā)展。在四君子湯中，影響FABP3 的活性成分包括5 alpha-stigmastan-3，6-dione（PubChem CID：13992092），Ergosta -7，22 e -dien -3beta -ol（PubChem CID：5283628），poricoic acid C（PubChem CID：56668247），glyasperin C（PubChem CID：480859），gancaonin A（PubChem CID：5317478），和licocoumarone（Pub－Chem CID：503731），這些活性成分通過影響FABP3來調(diào)控PPAR 信號通路，進而下調(diào)下游靶基因轉(zhuǎn)錄表達水平，抑制脂質(zhì)分解和脂肪酸氧化，從而降低機體能量消耗[12]。此外，四君子湯中的Taraxerol（Pub－Chem CID：92097）、perlolyrine（PubChem CID：160179）、Stigmasterol（PubChem CID：5280794）等成分可能通過調(diào)控mTOR 下游信號通路磷酸化4EBP1，增加eIF4E 和eIF4G 的結(jié)合活性，最終促進骨骼肌蛋白合成[13]?；衔锶?1-hydroxyrankinidine（PubChem CID：5318332）和AC1LCTJ1（PubChem CID：637112）作用于INS，Quercetin 作用于SKP1，進一步調(diào)控PI3K/Akt/mTOR 通路[14]，從而抑制蛋白分解。

3.2 抑制腫瘤 惡性腫瘤可通過促進腫瘤細胞增殖、抑制腫瘤細胞凋亡、促進腫瘤血管生成等途徑導(dǎo)致癌性惡病質(zhì)的發(fā)生發(fā)展。四君子湯的活性成分licochalcone B（PubChem CID：5318999）、isorhamnetin（PubChem CID：5281654）、calycosin（PubChem CID：5280448）等可能通過影響Wnt 和JAK-STAT 信號通路，抑制CCND1、CTNNB1、STAT5A 的表達[15-17]，調(diào)控細胞周期，抑制細胞增殖，發(fā)揮其抗腫瘤作用。7，9（11）-dehydropachymic acid（PubChem CID：15226717）、poriferasta-7，22 e-dien-3beta-ol（Pub-Chem CID：5283663）、luteolin（PubChem CID：5280445）等活性成分可能作用于IL6ST 和AR[18-19]，調(diào)節(jié)信號轉(zhuǎn)導(dǎo)，最終抑制腫瘤細胞增殖。在四君子湯中，licochalcone B（PubChem CID：5318999）作 用 于STAT5A 和CYCS[20-22]，抑制JAK-STAT 等信號通路下游表達，促進腫瘤細胞凋亡。licochalcone B（Pub－Chem CID：5318999）對STAT5A 和CYCS 的調(diào)控可能是四君子湯對腫瘤細胞凋亡調(diào)控的關(guān)鍵部分。此外，四君子湯可能通過calycosin（PubChem CID：5280448）作用于VEGFA，調(diào)控PI3K-Akt、HIF-1、VEGF 信號通路下游基因的表達，最終抑制腫瘤血管生成[23-24]。

4 結(jié)論

運用網(wǎng)絡(luò)藥理學(xué)和分子對接方法預(yù)測四君子湯治療癌性惡病質(zhì)的化合物、靶點和通路。其潛在機制可能主要集中在調(diào)節(jié)代謝和抑制腫瘤兩個方面，調(diào)節(jié)代謝包括調(diào)節(jié)葡萄糖、脂質(zhì)和蛋白質(zhì)的代謝，抑制腫瘤包括防止細胞增殖，促進腫瘤細胞凋亡，抑制血管生成。特異性信號通路分析表明，四君子湯主要通過PI3K-Akt 和Wnt 信號通路發(fā)揮治療癌性惡病質(zhì)的作用。因此，四君子湯通過多個化合物、多個靶點、多個通路的直接或間接協(xié)同作用，能夠有效地對抗癌性惡病質(zhì)。然而，本研究僅是基于數(shù)據(jù)分析，尚需要進一步的體內(nèi)和體外實驗來驗證。