基于網絡藥理學和生物信息學的長春新堿治療乳腺癌分子機制研究

蔣正立 陳 靜 夏愛曉 戴丹平

乳腺癌是發生在乳腺腺上皮組織的惡性腫瘤。根據腫瘤的分期和患者身體狀況，可采用手術、化療、內分泌治療、靶向治療及中醫藥輔助治療等多種手段，而化療方案中常用藥物有蒽環類、紫杉醇類以及長春新堿等。長春新堿（vincristine，VCR）治療乳腺癌具有較好的臨床療效[1-3]。本文基于網絡藥理學和生物信息學方法，通過相關數據庫篩選VCR 治療乳腺癌的靶點，分析其生物學功能和作用信號通路，從而探討VCR 治療乳腺癌的分子作用機制。

1 資料與方法

1.1 長春新堿相關靶點預測 VCR 是夾竹桃科植物長春花中提取出的生物堿，以“Vincristine”及其結構為關鍵詞，分別從Drugbank、Swiss Target Prediction、Stitch 數據庫預測篩選VCR 相關靶點。

1.2 GEO 基因芯片數據分析乳腺癌相關靶點 從GEO（https://www.ncbi.nlm.nih.gov/geo/）數據庫下載編號GSE29044 的基因芯片數據原始文件，采用Geo2r在線工具進行差異基因分析，以|logFC|≥2 和P<0.05為篩選條件進行篩選，繪制前50 個差異基因熱圖。

1.3 乳腺癌相關靶點獲取 以“Breast cancer”為關鍵詞，分別通過CTD（http://ctdbase.org/）和OMIM（http://www.omim.org/）數據庫查找乳腺癌相關靶標。并與GEO 基因芯片數據分析得到乳腺癌相關靶點合并刪除重復，即為乳腺癌疾病靶點。

1.4 VCR 治療乳腺癌相關靶點PPI 網絡建立及核心靶點獲取 將預測得到VCR 靶點與乳腺癌靶點導入Venny 在線工具取交集，即為VCR 治療乳腺癌相關靶點，導入String 數據庫，建立VCR 治療乳腺癌相關靶點PPI 網絡，將結果導入Cytoscape3.7.1 軟件，對網絡進行拓撲學分析，篩選核心靶點。

1.5 GO 功能和KEGG 通路富集 為進一步闡明VCR 治療乳腺癌相關機制，將VCR 治療乳腺癌潛在靶點通過DAVID（https://david.ncifcrf.gov/）數據庫進行GO 功能和KEGG 通路富集分析，并用Prism 軟件和在線繪圖站Ehbio 將結果繪制成柱狀圖與氣泡圖。

1.6 分子對接 采用分子對接方法研究受體蛋白EGFR、TP53、AKT1 與ACR 配體小分子之間的結合模式，以結合能≤-5.0 kJ/mol 為標準，評價VCR 與相關靶點結合作用。

2 結果

2.1 長春新堿作用靶點 從Drugbank（https://go.drugbank.com/）數據庫中得到15個靶點，Swiss Target Prediction 數據庫中獲得10 個靶點，Stitch 數據庫中獲得100 個靶點，刪去重復靶點，共得到124個VCR 作用靶點。

2.2 乳腺癌相關基因分析與檢索 篩選出319 個明顯影響和改變的差異表達基因，通過CTD 和OMIM數據庫共篩選得到718 個疾病靶點。兩者合并，刪除重復，共得到1003 個乳腺癌靶點。

2.3 長春新堿治療乳腺癌靶標分析共得到26 個（映射率為2.4%）VCR 治療乳腺癌潛在靶標（圖2），包括表皮生長因子受體（EGFR）、腫瘤抑制蛋白p53（TP53）、RAC-α 絲氨酸/蘇氨酸-蛋白激酶（AKT1）、周期蛋白A2（CCNA2）、90kDa 熱休克蛋白αA1（HSP90AA1）、雌激素受體基因（ESR1）、雷帕霉素靶蛋白（MTOR）、B 淋巴細胞瘤-2 基因（BCL2）、非受體酪氨酸激酶（SRC）、多藥耐藥性蛋白（ABCB1）等。

2.4 VCR 治療乳腺癌相關靶點PPI 網絡建立 通過String 數據庫建立PPI 網絡，包含26 個節點和133 條相互作用關系（見圖2），進一步拓撲學分析，發現其中點度中心性（DC）、接近中心性（CC）和中介中心性（BC）大于均值（點度中心性=10.23，中介中心性=15.85，接近中心性=0.025）的關鍵節點有10 個，占總節點數的38.5%，其中度值前三為：表皮生長因子受體（EGFR）、腫瘤抑制蛋白p53（TP53）、RAC-α絲氨酸/蘇氨酸-蛋白激酶（AKT1），提示它們在整個網絡中發揮較為重要的作用。

2.5 GO 功能和KEGG 通路富集 采用DAVID 數據庫對VCR 治療乳腺癌潛在靶點進行GO 富集分析，共得到140 個顯著條目（P<0.05），如圖3 所示，包括生物過程95 個，主要涉及蛋白質自磷酸化、磷脂酰肌醇介導的信號傳導等；細胞組成11 個，主要涉及質膜、磷脂酰肌醇3-激酶復合物、核等；分子功能34 個，主要涉及ATP 結合、一氧化氮合酶調節劑活性、激酶活性等。KEGG 富集結果共得到65 條顯著通路（P<0.05），排名前20 條通路如圖4 所示，主要涉及磷脂酰肌醇3 激酶/蛋白激酶B（PI3K-Akt）、表皮生長因子受體（ErbB）、缺氧誘導因子-1（HIF-1）等信號通路。

2.6 分子對接 采用分子對接的方法研究受體蛋白EGFR、TP53、AKT1 與ACR 配體小分子之間的結合模式。EGFR 與與VCR 配體小分子結合能為-7.66 kJ/mol，TP53 與VCR 配體小分子結合能為-7.89 kJ/mol，AKT1 與與VCR 配體小分子結合能為-5.07 kJ/mol，由此可見VCR 與三個靶蛋白結合能均小于-5 kJ/mol，表明具有較好的結合能力。見表1。

表1 長春新堿與核心靶點分子對接結果

3 討論

近年來，從基礎研究到臨床觀察多項證據顯示，VCR 對于治療乳腺癌顯示了一定效果，不良反應未見升高，可改善患者的生存質量[4-5]，特別是VCR新型給藥系統的研究開發，如靶向制劑的研制，可降低VCR 毒副作用，提高治療效果[6-8]。

本 文 從Drugbank、Swiss Target Prediction 及Stitch 數據庫共得到124 個成分靶點，通過GEO 基因芯片數據、CTD 和OMIM 數據庫共得到1003 個乳腺癌靶點，因此得到26 個VCR 治療乳腺癌潛在靶標，建立PPI 網絡進行拓撲學分析，獲得關鍵節點10個，分別是EGFR、TP53、AKT1 等。同時，對前三個核心治療靶點與VCR 進行分子對接，提示VCR 可能通過這些靶點起到治療乳腺癌的作用。EGFR 靶點對腫瘤生長、發展以及腫瘤干細胞的維持有著非常重要作用，EGFR 表達與乳腺癌患者腫瘤大小、淋巴結轉移、病理學分型相關[9-10]，在三陰型乳腺癌中的陽性表達率最高[11]，高表達也提示乳腺癌組織的惡性程度更高、侵襲性更強。TP53 能調節細胞周期和避免細胞癌變發生，VCR 作用于乳腺癌p53 靶點，mRNA 和蛋白水平升高[12]，p53 參與乳腺癌的生物進程與信號通路，使其作為治療乳腺癌的生物標志物和治療靶點成為可能[13-14]。AKT1 參與多種生物學過程，在乳腺癌浸潤進展過程中發揮重要作用[15]，雖未見文獻報道VCR 與此關鍵基因相關性，基于分子對接結果，很可能參與VCR 的體內作用效應機制。采用DAVID 數據庫對VCR 治療乳腺癌潛在靶標進行GO 功能和KEGG 富集分析，得到140 個顯著條目和65 條顯著通路，主要涉及PI3K-Akt、ErbB、HIF-1 等信號通路。PI3K-Akt 信號通路涉及細胞增殖、分化、凋亡等多種功能的調節，過度激活的PI3K-Akt 信號通路在乳腺癌的發生、發展中發揮重要作用，為治療的新興靶點。多個研究提示，PI3K-Akt 通路為治療乳腺癌的潛在靶標[16-17]。以上的關鍵靶點如EGFR、AKT1 等可通過PI3K-Akt 通路，抑制細胞凋亡，VCR 可通過抑制PI3K-Akt 通路，可提高實體腫瘤對VCR 的敏感性[18-19]。ErbB 也是乳腺癌的一個重要分子標志物和治療靶點，其異常激活可引起細胞生長增殖失控、惡性轉化及腫瘤浸潤轉移[20]，乳腺癌c-ErbB-2 陽性表達率較高，與淋巴結轉移和腫瘤標志物水平相關[21]。HIF-1 也參與腫瘤細胞的增殖、對抗凋亡和耐藥等重要過程[22]，是VCR 治療乳腺癌潛在靶點的信號通路，尚需實驗驗證。

綜上所述，本研究篩選了VCR 作用靶點和乳腺癌靶點，得到VCR 治療乳腺癌的靶點，分析了VCR治療乳腺癌的生物學功能和作用信號通路，為進一步探討VCR 抗乳腺癌作用分子機制及靶點提供思路。