張子函，李文鑫，羅香蓉，熊祖明，林益镕，張 森

（廣西醫(yī)科大學(xué)第一附屬醫(yī)院結(jié)直腸肛門外科，南寧 530021）

結(jié)直腸癌是最常見(jiàn)的惡性腫瘤之一，在2021年美國(guó)男性及女性的新發(fā)癌癥中均排在第3 位，在致死的惡性腫瘤中也排在第3 位[1]。早期結(jié)直腸癌主要通過(guò)手術(shù)來(lái)達(dá)到治療目的，但是一旦進(jìn)展到晚期則主要通過(guò)化療與靶向治療使疾病緩解。但部分患者對(duì)目前臨床廣泛使用的化療藥產(chǎn)生耐藥，迫切需要研發(fā)新的治療藥物。柚皮素是一種廣泛存在于柑橘類水果中的黃酮成分，先前的研究已經(jīng)證實(shí)，柚皮素具有一定的抗腫瘤作用[2-3]。但是柚皮素是否具有抗結(jié)直腸癌腫瘤的作用以及其可能的作用機(jī)制尚未見(jiàn)相關(guān)報(bào)道。本文運(yùn)用網(wǎng)絡(luò)藥理學(xué)的方法探究柚皮素作用于結(jié)直腸癌產(chǎn)生的作用及其機(jī)制，并通過(guò)實(shí)驗(yàn)進(jìn)一步驗(yàn)證柚皮素對(duì)結(jié)直腸癌細(xì)胞的抑制作用及核心基因表達(dá)的改變。

1 材料與方法

1.1 篩選柚皮素的作用靶點(diǎn)

柚皮素是在許多植物中都存在的化合物，在PubChem 數(shù)據(jù)庫(kù)（https://pubchem.ncbi.nlm.nih.gov/）中檢索柚皮素（naringenin），保存它的3D 格式并上傳至pharmmapper（http://www.lilab-ecust.cn/pharmmapper/）數(shù)據(jù)庫(kù)中得到柚皮素的主要藥效團(tuán)和其作用的靶點(diǎn)。同時(shí)利用SwissTargetPrediction（http://www.swisstargetprediction.ch/）數(shù)據(jù)庫(kù)獲得柚皮素的靶點(diǎn)基因，兩者取交集獲得最終作用靶點(diǎn)。

1.2 篩選結(jié)直腸癌的作用靶點(diǎn)

使用GeneCards 數(shù)據(jù)庫(kù)（https://www.genecards.org/），以“colorectal cancer”為關(guān)鍵詞檢索結(jié)直腸癌的作用靶點(diǎn)。同時(shí)使用DisGeNet 數(shù)據(jù)庫(kù)（https://www.disgenet.org/），以“colorectal cancer”為關(guān)鍵詞檢索結(jié)直腸癌的作用靶點(diǎn)。將兩個(gè)數(shù)據(jù)庫(kù)獲得的結(jié)直腸癌作用靶點(diǎn)合并后刪除重復(fù)項(xiàng)，獲得結(jié)直腸癌作用的最終靶點(diǎn)。

1.3 柚皮素治療結(jié)直腸癌的靶點(diǎn)

使用微生信在線網(wǎng)站（http://www.bioinformatics.com.cn/）中的jvenn 工具得到柚皮素作用靶點(diǎn)與結(jié)直腸癌靶點(diǎn)的交集，即柚皮素治療結(jié)直腸癌的作用靶點(diǎn)。

1.4 蛋白質(zhì)—蛋白質(zhì)相互作用（PPI）網(wǎng)絡(luò)的構(gòu)建

將獲得的柚皮素治療結(jié)直腸癌的作用靶點(diǎn)上傳至STRING（https://string-db.org/）數(shù)據(jù)庫(kù)，選擇物種為“homo sapiens”得到PPI相互作用網(wǎng)絡(luò)，下載其TSV格式。經(jīng)過(guò)處理后將數(shù)據(jù)導(dǎo)入Cytoscape中，并利用cytoHubba找到柚皮素作用于結(jié)直腸癌的核心基因。

1.5 KEGG通路富集分析和GO富集分析

利用DAVID 在線數(shù)據(jù)庫(kù)（https://david.ncifcrf.gov/）對(duì)之前獲得的可能作用靶點(diǎn)進(jìn)行GO富集分析和KEGG富集分析。再使用R包ggpolt2分別繪制P值最小的20 個(gè)分子功能（molecular function，MF）、細(xì)胞組分（cellular component，CC）、生物過(guò)程（biological process，BP）以及20條KEGG通路的條形圖。

1.6 實(shí)驗(yàn)驗(yàn)證

1.6.1 CCK-8 檢測(cè)柚皮素對(duì)結(jié)直腸癌細(xì)胞的作用

柚皮素購(gòu)自上海源葉生物公司（貨號(hào)：B21596-20 mg）取對(duì)數(shù)生長(zhǎng)期的結(jié)直腸癌HCT116細(xì)胞，將細(xì)胞密度調(diào)整為5×104個(gè)/mL，接種于96 孔板中，每孔100 μL。在37 ℃、5%CO2的培養(yǎng)箱中培養(yǎng)24 h后，吸出培養(yǎng)基，每孔加入含有不同濃度柚皮素的新鮮培養(yǎng)基100 μL，濃度梯度為：0.1%DMSO、50 μmol/L、100 μmol/L、200 μmol/L、400 μmol/L、800 μmol/L。分別在37 ℃5% CO2的培養(yǎng)箱中培養(yǎng)24 h、48 h、72 h，每孔加入10 μL的CCK-8試劑混勻后在培養(yǎng)箱中放置90 min，用酶標(biāo)儀檢測(cè)450 nm波長(zhǎng)時(shí)的吸光度并計(jì)算柚皮素對(duì)結(jié)直腸癌細(xì)胞增殖的抑制作用并繪制折線圖。

1.6.2 實(shí)時(shí)熒光定量PCR（RT-qPCR）檢測(cè)核心基因的mRNA 表達(dá)水平 柚皮素作用48 h 后各孔加入600 μL Trizol 裂解后提取總RNA，之后逆轉(zhuǎn)錄成cDNA 進(jìn) 行PCR 擴(kuò)增。反 應(yīng)體系 為：cDNA 模 板2 μL，正向引物0.8 μL，反向引物0.8 μL，SYBR Premix Ex TaqII 10 μL，ddH2O 6.4 μL。反應(yīng)條件為：95 ℃，30 s，95 ℃5 s、60 ℃30 s重復(fù)45個(gè)循環(huán)。引物序列為：PIK3CA：上游：5’-CGGTGACTGTGTGGGACTTATTGAG-3’，下游：5’-TGTAGTGTGTGGCTGTTGAACTGC-3’；MET：上 游：ATCAGGAGGTGTTTGGAAAGAT-3’，下游：5’-CGACTGTATGTCAGCAGTATGA；ARK1C3-3’：上游：5’-AAGCTTTGGTCCACTTTTCATC-3’，下 游：5’-GGTCAACATAGTCCAATTGAGC-3’；HSD17B2：上游：5’-CTCGTTAGCCAGCAAGGACTTCTC-3’，下 游：5’-AGCAAGGCAGATCCACAAGTAAGC-3’；SNCA：上 游：5’-GACCAAAGAGCAAGTGACAAAT-3’，下游：5’-CAGAAGGCATTTCATAAGCCTC-3’；SRC：上游：5’-GGAGTGGTATTTTGGCAAGATC-3’，下 游：5’-GATCTTGTAGTGCTTCACGTTG-3’；ESR1：上游：5’-CCTCCTCATCCTCTCCCACATCAG-3’，下 游：5’-GCATCTCCAGCAGCAGGTCATAG-3’。實(shí)時(shí)熒光定量PCR的結(jié)果用相對(duì)定量法用2-ΔΔCT值來(lái)表示。

1.7 統(tǒng)計(jì)學(xué)方法

采用GraphPad Prism 軟件進(jìn)行繪圖與統(tǒng)計(jì)分析。計(jì)量資料以均數(shù)±標(biāo)準(zhǔn)差（±s）表示，組間比較采用t檢驗(yàn)，多組間比較采用單因素方差分析，組間兩兩比較采用LSD-t檢驗(yàn)，以P＜0.05為差異具有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 柚皮素的潛在作用靶點(diǎn)預(yù)測(cè)

將通過(guò)PubChem獲得的柚皮素3D結(jié)構(gòu)上傳至pharmmapper 后獲得柚皮素作用的靶點(diǎn)，共70 個(gè)。同時(shí)通過(guò)SwissTargetPrediction 獲得100 個(gè)柚皮素的作用靶點(diǎn)。刪除重復(fù)項(xiàng)后，總共獲得169 個(gè)柚皮素作用靶點(diǎn)。

2.2 結(jié)直腸癌的靶點(diǎn)預(yù)測(cè)及共同靶點(diǎn)篩選

在GenCards 數(shù)據(jù)庫(kù)中檢索得到結(jié)直腸癌的10 422 個(gè)相關(guān)靶點(diǎn)，在DisGeNet 數(shù)據(jù)庫(kù)中得到5 473 個(gè)結(jié)直腸癌相關(guān)靶點(diǎn)，刪除重復(fù)項(xiàng)后共獲得11 247 個(gè)結(jié)直腸癌相關(guān)靶點(diǎn)。通過(guò)jvenn 工具繪制維恩圖（圖1），得到141 個(gè)交集靶點(diǎn)，即柚皮素治療結(jié)直腸癌可能的作用靶點(diǎn)。

圖1 柚皮素作用靶點(diǎn)與結(jié)直腸癌相關(guān)基因交集

2.3 PPI網(wǎng)絡(luò)分析及核心靶點(diǎn)預(yù)測(cè)

在STRING數(shù)據(jù)庫(kù)中導(dǎo)入之前獲得的141個(gè)交集靶點(diǎn)，得到蛋白質(zhì)相互作用網(wǎng)絡(luò)（圖2）。之后在CytoScape 中獲取柚皮素作用于結(jié)直腸癌的核心基因，即SRC、CBL、PIK3CA、SYK、LCK、VEGFA、HCK、KDR、SNCA、AKR1C3、HSD17B1、HSD17B2、CYP3A4、CYP19A1、PIK3CB、MET、IGF1R、ESR1、INSR、VAV3（圖3）。

圖2 柚皮素作用于結(jié)直腸癌的核心基因

圖3 PPI網(wǎng)絡(luò)

2.4 GO生物功能分析和KEGG通路富集分析

利用DAVID 數(shù)據(jù)庫(kù)對(duì)之前獲得的柚皮素可能作用于結(jié)直腸癌的141 個(gè)交集靶點(diǎn)進(jìn)行GO 生物功能分析和KEGG通路富集分析，獲得靶基因參與的信號(hào)通路和生物學(xué)過(guò)程。GO功能富集分析結(jié)果顯示，CC集中于細(xì)胞膜、血小板、肥大顆粒細(xì)胞、轉(zhuǎn)移酶復(fù)合物等（圖4A）。MF 主要集中于碳酸鹽脫水酶、蛋白酪氨酸激酶、水解酶等多種酶的活性，以及類固醇激素受體活性，類固醇結(jié)合，核受體、配體激活轉(zhuǎn)錄因子活性，激素結(jié)合等（圖4B）。BP 主要集中于蛋白質(zhì)自磷酸化，肽基酪氨酸的磷酸化及修飾，對(duì)淀粉樣蛋白-β的反應(yīng)，對(duì)類固醇激素的反應(yīng)，MAP 激酶活性的正相關(guān)性及活性，異生物代謝過(guò)程，生殖結(jié)構(gòu)與系統(tǒng)發(fā)育等（圖4C）。KEGG通路富集分析結(jié)果顯示，內(nèi)分泌抵抗（endocrine resistance）、EGFR 酪氨酸激酶抑制劑耐藥（EGFR tyrosine kinase inhibitor resistance）、PI3K-Akt 信號(hào)通路（PI3K-Akt signaling pathway）、癌癥中的蛋白多糖（proteoglycans in cancer）、前列腺癌（prostate cancer）、Ras 信號(hào)通路（ras signaling pathway）、花生四烯酸代謝（arachidonic acid metabolism）、Rap1 信號(hào)通路（Rap1 signaling pathway）、化學(xué)致癌作用—受體激活（Chemical carcinogenesis-receptor activation）、粘著力（focal adhesion），膀胱癌（bladder cancer）、脂質(zhì)和動(dòng)脈粥樣硬化（lipid and atherosclerosis）、催乳素信號(hào)通路（prolactin signaling pathway）、卵巢類固醇生成（ovarian steroidogenesis）、乳腺癌（breast cancer）、黑色素瘤（breast cancer）、非小細(xì)胞肺癌（Nonsmall cell lung cancer）、亞油酸代謝（linoleic acid metabolism）等是柚皮素的主要作用通路（圖4D）。由此可見(jiàn)，柚皮素對(duì)于結(jié)直腸癌可能具有一定的抑制作用。

圖4 GO分析和KEGG通路富集分析

2.5 實(shí)驗(yàn)驗(yàn)證

2.5.1 CCK-8 檢測(cè)柚皮素對(duì)結(jié)直腸癌細(xì)胞的抑制作用 在柚皮素濃度為50 μmol/L、100 μmol/L、200 μmol/L、400 μmol/L、800 μmol/L時(shí)，24 h的抑制率分別為6.12%、11.37%、26.16%、59.78%、81.43%；48 h 的抑制率分別為15.01%、30.37%、45.04%、85.50%、87.61%；72 h 的抑制率分別為14.35%、47.31%、69.28%、93.26%、91.03%。隨著柚皮素濃度的上升和作用時(shí)間的延長(zhǎng)，柚皮素對(duì)結(jié)直腸癌的抑制作用逐漸增強(qiáng)（P＜0.05），見(jiàn)圖5。

圖5 柚皮素對(duì)結(jié)直腸HCT116細(xì)胞增殖的抑制作用

2.5.2 RT-qPCR檢測(cè)核心基因的mRNA表達(dá)水平

給予200 μmol/L的柚皮素之后,與對(duì)照組相比，柚皮素處理過(guò)的細(xì)胞內(nèi)PIK3CA、MET、HSD17B2、SNCA、SRC、ESR1mRNA 表達(dá)水平降低（P＜0.05）；而ARK1C3mRNA 表達(dá)水平升高（P＜0.05），見(jiàn)圖6。

圖6 柚皮素對(duì)核心基因mRNA表達(dá)水平的影響

3 討論

有研究認(rèn)為陳皮、半夏、草豆蔻、甘草中的主要抗腫瘤成分為柚皮素[4-6]。柚皮素是天然二氫黃酮類化合物，研究顯示其可通過(guò)抑制細(xì)胞增殖而發(fā)揮抗腫瘤作用，如在HepG2肝癌細(xì)胞和人胎盤絨毛膜癌中[7-8]。同時(shí)，它也可通過(guò)誘導(dǎo)HepG2 肝癌細(xì)胞、胰腺癌等細(xì)胞凋亡而抑制腫瘤生長(zhǎng)[7，9]。柚皮素還可以抑制腫瘤細(xì)胞的侵襲遷移，如在神經(jīng)膠質(zhì)瘤細(xì)胞[2]、乳腺癌細(xì)胞[10]和前列腺癌細(xì)胞中[3]。但柚皮素在結(jié)直腸癌中發(fā)揮的抗腫瘤作用及其機(jī)制目前的研究相對(duì)較少。

通過(guò)網(wǎng)絡(luò)藥理學(xué)分析顯示，PIK3CA、SRC、HSD17B2、SNCA、CYP19A1、PIK3CB、INSR、MET、ESR1、ARK1C3等為柚皮素治療結(jié)直腸癌的核心基因。其中，SRC在超過(guò)80%的結(jié)直腸癌患者中異常過(guò)度表達(dá)和激活，柚皮素作用于SRC基因就可能發(fā)揮其抗結(jié)直腸腫瘤作用[11]。MET 也稱為肝細(xì)胞生長(zhǎng)因子受體，是一種酪氨酸激酶，被認(rèn)為是結(jié)直腸癌早期侵襲和轉(zhuǎn)移的重要因素，可作為重要的預(yù)后指標(biāo)[12]。研究顯示，激活MET可以誘導(dǎo)結(jié)直腸癌細(xì)胞的侵襲和轉(zhuǎn)移[12]，同時(shí)，MET 也參與了癌細(xì)胞的免疫逃逸[13]。提示柚皮素可能通過(guò)作用于MET 基因抑制結(jié)直腸癌的轉(zhuǎn)移，改善患者預(yù)后。另有研究顯示，ESR1 的表達(dá)與結(jié)直腸癌細(xì)胞對(duì)五氟尿嘧啶的敏感性呈負(fù)相關(guān)[14]，柚皮素降低ESR1的表達(dá)可能會(huì)增高結(jié)直腸癌細(xì)胞對(duì)于化療藥物的敏感性從而達(dá)到更好的治療效果。SNCA基因被證實(shí)與胃腸道腫瘤存在一定的關(guān)系[15]，柚皮素作用于這一基因可能對(duì)結(jié)直腸癌產(chǎn)生一定的影響。

本研究通路富集結(jié)果顯示，PI3K-Akt信號(hào)通路是柚皮素抗結(jié)直腸癌可能發(fā)揮作用的通路，PIK3CA在結(jié)直腸癌組織中的表達(dá)高于正常腸上皮組織[16]，PI3K-Akt信號(hào)通路與結(jié)直腸癌相關(guān)，先前的研究證實(shí)了PI3K-Akt 信號(hào)通路可以作為結(jié)直腸癌的治療靶點(diǎn)[17-18]。同時(shí)，柚皮素可作用于催乳素信號(hào)通路，卵巢類固醇生成信號(hào)通路，乳腺癌信號(hào)通路，而這些通路都與體內(nèi)雌激素相關(guān)。柚皮素可能通過(guò)作用于雌激素受體而發(fā)揮抗結(jié)直腸癌的作用。而花生四烯酸代謝、亞油酸代謝通路是體內(nèi)的脂代謝通路，柚皮素也可能通過(guò)改變脂肪組織而發(fā)揮抗結(jié)直腸癌的作用。

CCK-8實(shí)驗(yàn)顯示，柚皮素對(duì)結(jié)直腸癌細(xì)胞的生長(zhǎng)具有抑制作用。同時(shí)RT-qPCR實(shí)驗(yàn)顯示，結(jié)直腸癌細(xì)胞經(jīng)過(guò)柚皮素處理后核心基因PIK3CA、SRC、HSD17B2、MET、ESR1和ARK1C3的表達(dá)量會(huì)發(fā)生顯著性改變，進(jìn)一步表明柚皮素通過(guò)靶向這些基因產(chǎn)生抗結(jié)直腸癌作用。

本研究通過(guò)網(wǎng)絡(luò)藥理學(xué)的方法，預(yù)測(cè)了柚皮素的抗結(jié)直腸癌作用及其可能的作用機(jī)制，通過(guò)實(shí)驗(yàn)驗(yàn)證了柚皮素對(duì)結(jié)直腸癌細(xì)胞具有抑制作用并使核心基因在結(jié)直腸癌細(xì)胞中的表達(dá)情況改變，可為后續(xù)柚皮素治療結(jié)直腸癌的進(jìn)一步研究提供理論依據(jù)。