張 芳 沈 趣 吳 繁 王彬彬

結直腸癌發(fā)病率居所有癌癥的第三位，死亡率近年已躍居癌癥第二位[1-2]。結直腸癌病因尚不明確，但大量的研究證據(jù)表明結直腸癌的發(fā)生發(fā)展是有遺傳、環(huán)境和生活因素等多方面作用的結果，與結腸癌家族史、炎癥性腸病、紅肉和加工肉攝入、糖尿病、肥胖、吸煙、飲酒等危險因素相關[3]。結直腸癌屬中醫(yī)“癥瘕”“廟瘤”“積聚”“腸覃”“大瘕泄”“鎖肛痔”“腸風”“臟毒”等范疇[4-5]。中醫(yī)認為，結直腸癌發(fā)病多由先天稟賦不足、肝氣郁結、飲食失調、外感濕邪、脾腎虧虛等導致臟腑功能失調，濕熱內(nèi)生，蘊結腸道，積久而發(fā)病。臨床治療結直腸癌多從腸道濕熱論治，用藥以防風、白術使用最為頻繁，患者反饋獲益頗佳。白術健脾補中氣使脾氣自健，以補其虛；防風行氣活血化瘀使?jié)裥叙錾?，以治其實。目前防風-白術治療結直腸癌的分子機制和藥效物質基礎尚不明確，以致其在臨床上未被廣泛應用。

1 資料與方法

1.1 防風-白術主要化學成分篩選 本研究通過中藥系統(tǒng)藥理學數(shù)據(jù)庫（TCSMP，https//lsp.nwu.edu.cn/tcsmp.php）[6]檢索，以口服生物利用度≥30%，類藥性≥0.18 為閾值，篩選出防風-白術的主要化學成分。

1.2 防風-白術主要化學成分基因靶點預測 利用uniport（https://www.uniport.org/）數(shù)據(jù)庫將已篩選出的防風-白術主要化學成分對應的靶點與相應的人類已知基因靶點匹配[7]。

1.3 CRC 基因靶點預測 在GeneCards（https://www.genecard.org/）、OMIM（https://www.omim.org/）數(shù)據(jù)庫檢索結直腸癌相關基因，獲得結直腸癌相關基因。

1.4 網(wǎng)絡合并及分析 使用R 語言選取藥物主要成分基因靶點與疾病基因靶點交集，并繪制venn圖，并將防風-白術的活性成分與預測的基因交集關聯(lián)，很有可能得到防風-白術治療結直腸癌的活性成分及基因靶點，采用Cytoscape 3.7.1 軟件對所得結果進行可視化處理，并采用拓撲學分析各網(wǎng)絡節(jié)點的介度、緊密度以及連接度，遴選防風-白術治療結直腸癌的核心蛋白。利用String 數(shù)據(jù)庫（https://string-db.org/Version 11.0）預測靶標蛋白相互作用關系，把防風-白術預測作用靶標導入String 數(shù)據(jù)庫，獲得靶標蛋白相互作用關系。

1.5 防風-白術抗結直腸癌靶點生物功能分析及通路分析 應用R 語言對防風-白術治療結直腸癌的靶點進行GO（Gene Ontology）功能富集分析和KEGG（Kyoto Encyclopedia of Genes and Genomes）信號通路富集分析，探究防風-白術治療結直腸癌潛在靶點涉及的生物學過程及信號通路。

2 結果

2.1 防風-白術化學成分篩選及靶點預測 通過TCSMP 數(shù)據(jù)庫檢索，以口服生物利用度≥30%和類藥性≥0.18 為篩選條件并篩去無對應靶點的成分，得到防風-白術活性化合物紫花前胡素、珊瑚菜素、防風色酮醇、α-香樹精等共25 種。利用Uniport 數(shù)據(jù)庫將靶點名稱統(tǒng)一標準化處理后發(fā)現(xiàn)20 個活性成分具有對應靶點（見表1），共對應61 個靶點。

表1 防風-白術活性成分

2.2 防風-白術治療結直腸癌的候選靶點 通過GeneCards 數(shù)據(jù)庫檢索得到與結直腸癌相關靶點共9330 個，由于靶點數(shù)量大，選取其中relevance＞10 的靶點作為研究對象，得到密切關聯(lián)靶點1262 個。通過OMIM 數(shù)據(jù)庫檢索得到與結直腸癌相關基因341個，兩個數(shù)據(jù)庫共篩選得到基因靶點1603 個。通過R 語言對比分析結直腸癌與防風-白術有效活性成分對應靶標的共有部分，最終得到防風-白術活性成分治療結直腸癌的潛在靶點共35 個，并繪制成Venn圖。見表2，圖1。

圖1 結直腸癌與防風-白術相關靶點韋恩圖

表2 防風-白術抗結直腸癌預測靶標

2.3 防風-白術活性成分抗結直腸癌的作用靶標分析 將篩選的活性化合物及靶點導入Cytoscape

3.7.1 軟件構建“疾病-藥物-成分-靶標”網(wǎng)絡圖（見圖2），其中，三角形節(jié)點代表防風-白術有效活性成分，圓形節(jié)點代表作用靶標，菱形節(jié)點代表結直腸癌，六邊形節(jié)點代表防風-白術，邊線代表藥物、疾病、活性成分與靶標之間的相互聯(lián)系。漢黃芩素（Wogonin）對應的靶點最多，共有29 個靶點與其直接相連，前列腺素內(nèi)過氧化物合酶2（PTGS2）與有效化合物及結直腸癌間聯(lián)系邊線最多，提示漢黃芩素可作為防風-白術抗結直腸癌的關鍵藥理成分，通過影響PTGS2 發(fā)揮抗結直腸癌作用。從圖2 中可以看出，防風的多個活性成分可作用于一個靶點，也可以作用于不同靶點，同時，35 個作用靶點涉及多條信號通路，且多條信號通路中又可通過共有靶點相連，體現(xiàn)了防風-白術多成分、多靶點、多通路協(xié)同抗結直腸癌的作用機制。

圖2 “疾病-藥物-成分-靶標”網(wǎng)絡圖

2.4 防風-白術與結直腸癌相關靶標PPI 網(wǎng)絡的構建與分析 在String 數(shù)據(jù)庫平臺構建共有靶點的蛋白互作PPI 網(wǎng)絡（見圖3），將數(shù)據(jù)導入Cytoscape 3.7.1 軟件，利用軟件分析和處理（見圖4），網(wǎng)絡中節(jié)點大小及顏色深淺反應Degree 值的大小。PPI 網(wǎng)絡前十的核心基因分別是CASP3、AKT1、JUN、CCND1、IL-6、ESR1、PTGS2、RELA、MAPK14、CDKN1A。

圖3 String 平臺PPI 網(wǎng)絡互作圖

圖4 靶標互作網(wǎng)絡圖

2.5 GO 功能富集分析及KEGG 通路富集分析 采用Conbioconductor 數(shù)據(jù)庫及R 語言對防風-白術、CRC 的交集靶點進行GO 功能富集分析，設定標準為P 值＜0.05，GO 分析得到具有顯著意義的65 條GO 功能條目，均為分子功能，前20 條GO 分析結果如圖5 所示，主要涉及結合DNA 結合轉錄因子（DNA-binding transcription factor binding）、核受體活性（nuclear receptor activity）、配體激活的轉錄因子活性（ligand-activated transcription factor activity）等過程。KEGG 分析共得到139 條具有顯著意義的通路，前20 條與結直腸癌相關度較高的通路如圖6所示，圓點越大顏色越紅，表明該通路富集的靶點越多，P 值也越小，排名靠前的信號通路主要有人巨細胞病毒感染（Human cytomegalovirus infection）、糖尿病并發(fā)癥中的AGE-RAGE 信號通路（AGE-RAGE signaling pathway in diabetic complications）、乙型肝炎（Hepatitis B）、Kaposi 肉瘤相關皰疹病毒感染（Kaposi sarcoma-associated herpesvirus infection）等。

圖5 核心靶點的GO 分析

圖6 核心靶點的KEGG 分析

3 討論

本研究基于網(wǎng)絡藥理學探討防風-白術治療結直腸癌的潛在作用機制，構建疾病-藥物-成分-標靶網(wǎng)絡，該網(wǎng)絡表明了防風-白術治療結直腸癌具有多成分、多靶點、多通路協(xié)同作用的特點，該網(wǎng)絡提示抗結直腸癌的核心成分為漢黃芩素、β 谷甾醇、3β-羥基蒼術酮。漢黃芩素能夠通過下調綠膽素還原酶A（BLVARA），達到抑制結直腸癌細胞系HT-29 及SW620 細胞的增殖，抑制腫瘤細胞的遷移和侵襲能力，促進腫瘤細胞凋亡，抑制腫瘤上皮-間質轉化（EMT）進程，從而發(fā)揮抗結直腸癌作用[8]。Β-谷甾醇可以抑制腫瘤細胞增殖、誘導期凋亡，抑制其遷移和侵襲，通過Akt/GSK-3β 信號通路抑制EMT[9-10]，并且能夠調節(jié)巨噬細胞的功能，降低IgG、IgG1 等特異性抗體水平，抑制炎癥細胞因子產(chǎn)生[11]，從而減輕腫瘤病人的炎癥反應。

將藥物-疾病交集靶點進行蛋白互作網(wǎng)絡的分析，利用Cytoscape 的cytohubba 插件MCC 算法計算出PPI 網(wǎng)絡前十位核心基因，與CRC 相關度較高的靶點有CASP3、AKT1、JUN、CCND1、IL-6、ESR1、PTGS2 等。故推測上述靶點與防風-白術治療結直腸癌的作用機制密切相關。CASP3 通過切割各種基本細胞蛋白在凋亡通路中起著關鍵作用，在結直腸癌、卵巢癌和乳腺癌等幾種類型癌癥中，CASP3mRNA 水平明顯低于正常人[12]，通過調節(jié)CASP3 水平治療結直腸癌[13]。AKT1 是PI3K/Akt 信號通路中的核心靶標，調控細胞代謝、細胞周期、細胞增殖、細胞凋亡和血管生成等方面[14]。異常激活的Akt 可通過抑制蛋白水解酶Caspase9 活性促使結直腸癌細胞逃避細胞凋亡，同時激活下游mTOR 標靶引發(fā)一系列反應進而導致結直腸癌細胞增殖、血管生成和侵襲等[15]。JUN在細胞增殖、凋亡、應激反應及腫瘤的發(fā)生發(fā)展中發(fā)揮中要作用[16]。CCND1 基因編碼的細胞周期蛋白，具有促進細胞增殖和基因轉錄等功能，可在多種癌癥中成高表達，阻斷CCND1 編碼細胞周期蛋白的治療策略將有助于干擾多種癌癥特異性新高通路達到治療癌癥的目的[17]。IL-6 持續(xù)可激活各種炎癥/免疫/凋亡基因，如Bcl-2 基因家族和VEGF，這些基因導致腫瘤細胞入侵、血管生成和轉移[18]。ESR1 在結直腸癌組織中的表達水平顯著高于正常的結直腸組織[19]，給予結直腸癌患者5-FU 治療后可顯著降低患者癌組織中ESR1 mRNA 的表達水平及無病生存期（DFS）[20]。PTGS2 參與大腸癌發(fā)展的早期階段，最初影響腫瘤的基質微環(huán)境，然后影響上皮細胞，并且PTGS2 表達與結直腸癌中端粒酶活性之間成負相關關系[21]。

GO 功能富集分析共得到65 個富集項，且均為分子功能，防風-白術治療結直腸癌可能與核受體活性、配體激活的轉錄因子活性、結合DNA 結合轉錄因子、蛋白質絲氨酸/蘇氨酸活性、類固醇結合等條目密切相關。促進膽汁酸核受體活性可以明顯抑制直腸腫瘤細胞增殖能力[22-23]?？寡住⒖咕?、抗脂質過氧化、KEGG 通路富集提示，防風-白術抗結直腸癌的信號通路主要集中在惡性腫瘤、CRC 相關通路等，與AGE-RAGE 信號通路、乙型肝炎、人巨細胞病毒感染、IL-17、EB 病毒感染等信號通路密切相關。AGERAGE 信號通路激活可促進直腸癌HCT116 細胞的增殖、遷移、侵襲和克隆形成能力，并維持腫瘤細胞干性，其機制可能與AGE/RAGE/通路介導的EMT 有關[24]。結直腸癌患者血清HBs Ag 表達可能與結直腸癌的肝轉移及浸潤深度有關[25]，慢性HBV 感染與結直腸癌風險增加顯著相關[26]。人巨細胞病毒編碼的US28 基因使結直腸癌患者有更大的化療耐藥性及更高的侵襲能力，并加速細胞周期進展[27]。IL-17 在結直腸癌患者的血清和組織中表達上調，在結直腸癌的轉移和預后中有重要作用，可以抑制結直腸癌的發(fā)生[28]。EB 病毒誘導基因3 可能介導結直腸癌細胞對5-FU 的耐藥過程[29]。

綜上所述，本研究應用網(wǎng)絡藥理學方法預測了防風-白術治療結直腸癌的作用靶點及作用機制，闡明防風-白術通過多成分-多靶點-多通路治療結直腸癌，為進一步實驗研究和臨床應用提供了新思路、指明了新方向。但是，基于網(wǎng)絡藥理學研究藥物治療疾病的靶點仍存在局限性，其雖然能較為系統(tǒng)地預測活性成分、作用靶點與作用通路，但預測結果不能完全證明其正確性，其作用機制需要開展體外和體內(nèi)實驗進行驗證。