劉敏，李羿南，張青*

·信息技術(shù)與中醫(yī)藥·

基于網(wǎng)絡(luò)藥理學(xué)聯(lián)合GEO數(shù)據(jù)庫及分子對接技術(shù)探討化瘀丸治療乳腺癌的作用機(jī)制

劉敏1，李羿南2，張青1*

1.首都醫(yī)科大學(xué)附屬北京中醫(yī)醫(yī)院，北京 100010；2.北京中醫(yī)藥大學(xué)，北京 100029

通過網(wǎng)絡(luò)藥理學(xué)聯(lián)合GEO數(shù)據(jù)庫差異基因分析及分子對接技術(shù)研究化瘀丸治療乳腺癌的有效成分及作用機(jī)制。利用TCMSP數(shù)據(jù)庫檢索化瘀丸中12味中藥信息，獲得化瘀丸的活性成分和分子靶點(diǎn)。通過GEO數(shù)據(jù)庫檢索“breast cancer”篩選后獲得GSE139038數(shù)據(jù)，通過R軟件分析得出差異基因。通過DrugBank、OMIM、TTD和CTD數(shù)據(jù)庫篩選出乳腺癌的相關(guān)靶點(diǎn)，利用UniProt數(shù)據(jù)庫對靶點(diǎn)進(jìn)行基因注釋，通過Cytoscape3.2.1軟件構(gòu)建化瘀丸與乳腺癌的成分-靶點(diǎn)網(wǎng)絡(luò)圖，使用插件BisoGenet和CytoNCA繪制核心靶點(diǎn)拓?fù)渚W(wǎng)絡(luò)，并進(jìn)行藥物與疾病交集靶點(diǎn)的GO功能與KEGG通路富集分析，最后將網(wǎng)絡(luò)中度值較高的成分與核心靶點(diǎn)進(jìn)行分子對接。從GEO數(shù)據(jù)庫下載分析得到乳腺癌上調(diào)和下調(diào)最明顯的差異基因40個。獲得化瘀丸主要活性成分7個，包括山柰酚（kaempferol）、谷甾醇（sitosterol）、異鼠李素（isorhamnetin）、常春藤皂苷元（hederagenin）、槲皮素（quercetin）、β-谷甾醇（beta-sitosterol）和豆甾醇（stigmasterol）等，作用于JUN、FOS、MYC、E2F1、E2F2等32個基因靶點(diǎn)。GO功能與KEGG通路富集分析表明，化瘀丸治療乳腺癌主要通過調(diào)節(jié)氧化應(yīng)激、激活轉(zhuǎn)錄因子結(jié)合、趨化因子受體結(jié)合等生物反應(yīng)，作用于病毒感染通路、細(xì)胞衰老信號通路、IL-17信號通路、肝癌通路等發(fā)揮治療作用?；鐾柚委熑橄侔┑幕钚猿煞譃殚纹に氐?，主要涉及病毒感染、IL-17及其他與癌癥相關(guān)的通路。

化瘀丸；乳腺癌；網(wǎng)絡(luò)藥理學(xué)；生物信息學(xué)；分子對接；作用機(jī)制

乳腺癌是女性最常見的惡性腫瘤之一，據(jù)國際癌癥研究中心GLOBOCAN對全球185個國家/地區(qū)的36種癌癥的發(fā)病率、死亡率及癌癥發(fā)展趨勢的統(tǒng)計數(shù)據(jù)，2020年新發(fā)乳腺癌2 261 419例，占總體癌癥發(fā)病的11.7%[1]。目前，西醫(yī)治療乳腺癌的主要方法是手術(shù)、化療、靶向治療、免疫治療、放射治療等，中醫(yī)藥治療乳腺癌效果顯著且不良反應(yīng)較小。采用中醫(yī)藥防治乳腺癌成為研究熱點(diǎn)。國家級名老中醫(yī)郁仁存教授以益氣活血、化瘀解毒立法，創(chuàng)立化瘀丸，臨床用于治療乳腺癌等腫瘤效果顯著[2]。本研究通過網(wǎng)絡(luò)藥理學(xué)聯(lián)合基因芯片數(shù)據(jù)庫GEO（Gene Expression Omnibus），獲取正常乳腺與乳腺癌標(biāo)本的基因芯片，分析兩者存在的差異表達(dá)基因，為乳腺癌腫瘤標(biāo)記物的進(jìn)一步研究提供參考，并通過分子對接技術(shù)評價受體和藥物分子親和關(guān)系。將網(wǎng)絡(luò)藥理學(xué)、GEO數(shù)據(jù)庫和分子對接技術(shù)相結(jié)合，對化瘀丸的主要活性成分進(jìn)行篩選，分析活性成分的作用靶點(diǎn)和生物通路，深入探索化瘀丸治療乳腺癌的分子機(jī)制，為乳腺癌患者的精準(zhǔn)治療提供依據(jù)。

1 資料與方法

1.1 化瘀丸活性成分收集與篩選

化瘀丸由砂仁、桃仁、赤芍、黃芪、拳參、白芷、莪術(shù)、郁金、王不留行、延胡索、紅參、西紅花12味組成。通過TCMSP（https://tcmsp-e.com/）檢索化瘀丸的活性成分，獲取分子信息，包括活性成分的分子名稱、口服生物利用度（OB）、類藥性（DL）等。以O(shè)B≥30%和DL≥0.18[3]為篩選條件，得到活性成分。

1.2 化瘀丸作用靶點(diǎn)與乳腺癌靶點(diǎn)預(yù)測

通過TCMSP獲得活性成分對應(yīng)靶點(diǎn)信息，利用UniProt數(shù)據(jù)庫（http://www.uniprot.org/）進(jìn)行校正，獲得對應(yīng)的基因名稱。將靶點(diǎn)基因的全稱轉(zhuǎn)化為簡稱，并刪除無對應(yīng)基因的靶點(diǎn)。利用DrugBank數(shù)據(jù)庫（https://www.drugbank.ca）、OMIM數(shù)據(jù)庫（https://www.omim.org/）、CTD數(shù)據(jù)庫（http://ctdbase.org/）和TTD數(shù)據(jù)庫（http://bidd.nus.edu.sg/BIDD-Databases/TTD/TTD.asp），以“breast cancer”為關(guān)鍵詞檢索，得到疾病相關(guān)的基因及靶點(diǎn)信息。將藥物靶點(diǎn)與疾病靶點(diǎn)取交集，獲得化瘀丸與乳腺癌的交集靶點(diǎn)，即化瘀丸治療乳腺癌的潛在作用靶點(diǎn)。

1.3 乳腺癌GEO芯片差異分析

在GEO數(shù)據(jù)庫（https://www.ncbi.nlm.nih.gov/geo/）檢索“breast cancer”，篩選并下載GSE139038芯片，使用R軟件設(shè)定差異倍數(shù)（fold change，F(xiàn)C）絕對值＞1且校正后值（adjPvalue）＜0.05作為顯著差異基因的篩選條件，得出差異分析結(jié)果，并繪制熱圖：使用R語言中的packages、pheatmap、Limma程序包對數(shù)據(jù)進(jìn)行二次分析，結(jié)合值和FC值進(jìn)行篩選，選取上調(diào)、下調(diào)最明顯的前20個基因，并定義分組（C為正常對照組，T為乳腺癌組）顏色、高度，繪制熱圖。火山圖繪制：設(shè)置軸與軸最大值，將基因統(tǒng)一設(shè)置為無差異基因，用黑色表示，再通過、logFC進(jìn)行篩選。當(dāng)＜0.05且logFC＞1，表示該基因上調(diào)，用紅色表示；＜0.05且logFC＜-1，表示該基因下調(diào)，用綠色表示。

1.4 中藥-活性成分-靶點(diǎn)-疾病網(wǎng)絡(luò)構(gòu)建

將化瘀丸治療乳腺癌的潛在作用靶點(diǎn)即交集基因輸入Cytoscape3.2.1軟件，構(gòu)建中藥-活性成分-靶點(diǎn)-疾病網(wǎng)絡(luò)。用節(jié)點(diǎn)（node）表示分子、靶蛋白，用邊（edge）表示活性成分與靶點(diǎn)之間的關(guān)系。采用Network Analyzer插件進(jìn)行網(wǎng)絡(luò)特征分析，以明確化瘀丸中重要的活性成分及靶點(diǎn)。

1.5 蛋白互作網(wǎng)絡(luò)構(gòu)建與分析

利用Cytoscape3.2.1軟件的BisoGenet和CytoNCA插件對蛋白互作（PPI）網(wǎng)絡(luò)節(jié)點(diǎn)進(jìn)行分析，設(shè)置度中心性（DC）＞61、介度中心性（BC）＞600[4]，篩選重要基因。獲得的相關(guān)活性成分可能直接或間接調(diào)控這些靶點(diǎn)發(fā)揮重要作用。

1.6 GO功能富集分析

使用R軟件對化瘀丸治療乳腺癌的基因進(jìn)行GO富集分析。根據(jù)富集的數(shù)目和顯著程度（＜0.05，＜0.05），選取前20條富集條目，分別進(jìn)行生物過程（biological process，BP）、細(xì)胞成分（cellular component，CC）、分子功能（molecular function，MF）的富集，繪制柱狀圖。

1.7 KEGG富集分析

使用R軟件對化瘀丸治療乳腺癌的基因進(jìn)行KEGG富集分析，使用DOSE、clusterProfiler、enrichplot、colorspace、stringi、ggplot2程序包進(jìn)行分析。根據(jù)基因富集的數(shù)目和顯著程度（＜0.05，＜0.05），取前20條KEGG通路，將得到的KEGG通路導(dǎo)入Cytoscape3.2.1軟件，根據(jù)degree值大小，繪制通路與基因網(wǎng)絡(luò)圖。

1.8 分子對接

從RSCB PDB數(shù)據(jù)庫（https://www.rcsb.org/）下載核心靶點(diǎn)及PPI網(wǎng)絡(luò)中核心蛋白的3D結(jié)構(gòu)PDB格式文件，在PubChem中下載配體，對接軟件為AutoDock Tools1.5.6，選擇默認(rèn)設(shè)置，設(shè)置Grid Box為整個蛋白分子，運(yùn)行Local Search Parameters進(jìn)行分子對接，對接成功后，選擇分?jǐn)?shù)最高的配體-受體構(gòu)象繪制對接示意圖。

2 結(jié)果

2.1 化瘀丸活性成分

通過TCMSP數(shù)據(jù)庫，以O(shè)B≥30%、DL≥0.18為條件篩選，獲得化瘀丸共計136個化合物，其中砂仁9個，桃仁19個，黃芪17個，拳參6個，白芷20個，莪術(shù)1個，郁金3個，王不留行4個，延胡索49個，紅參3個，西紅花5個。

2.2 化瘀丸作用靶點(diǎn)與乳腺癌疾病靶點(diǎn)信息

將TCMSP數(shù)據(jù)庫預(yù)測的化瘀丸成分靶點(diǎn)通過UniProt數(shù)據(jù)庫進(jìn)行基因注釋后得到注釋基因，通過DrugBank、OMIM、CTD、TTD數(shù)據(jù)庫檢索后得到相關(guān)基因，取交集，共獲得40個交集靶點(diǎn)，分別是PTGS2、DPP4、ADRB2、JUN、NR3C2、ADH1C、PPARG、MMP9、TOP2A、E2F2、IGHG1、OLR1、CAV1、ABCG2、IGF2、ADRA2A、MMP1、MYC、CXCL11、GSTA1、MAOA、STAT1、F3、CXCL2、GSTA2、CHEK1、GSTM2、CCNB1、CXCL10、CAT、CCNA2、AKR1C3、THBD、SPP1、GSR、ADH1B、FOS、SERPINE1、E2F1、AKR1C1。

2.3 乳腺癌GEO芯片差異分析

GEO數(shù)據(jù)庫下載GSE139038芯片數(shù)據(jù)，獲取相關(guān)的正常乳腺組織和乳腺癌標(biāo)本的芯片，分析二者存在的差異表達(dá)基因。樣本中正常對照組24個樣本，乳腺癌組41個樣本。利用R軟件和相關(guān)的軟件包對GSE139038進(jìn)行過濾機(jī)標(biāo)準(zhǔn)化，篩選差異基因。共得到乳腺癌上調(diào)和下調(diào)最明顯的差異基因40個，其中上調(diào)10個，下調(diào)30個。正常對照組與乳腺癌組差異基因火山圖及基因熱圖見本文OSID碼。

2.4 化瘀丸治療乳腺癌的活性成分及靶點(diǎn)網(wǎng)絡(luò)

將化瘀丸與乳腺癌差異基因取交集后，導(dǎo)入Cytoscape3.2.1軟件，根據(jù)不同靶點(diǎn)與有效成分對應(yīng)關(guān)系和數(shù)目進(jìn)行分類展示。得到潛在基因靶點(diǎn)40個，對應(yīng)有效成分100個。其中山柰酚（kaempferol）、谷甾醇（sitosterol）、異鼠李素（isorhamnetin）、常春藤皂苷元（hederagenin）、槲皮素（quercetin）、β-谷甾醇（beta-sitosterol）和豆甾醇（stigmasterol）為共有有效成分。在7個共有有效成分中，槲皮素對應(yīng)靶點(diǎn)最多，因此槲皮素可能是化瘀丸治療乳腺癌的關(guān)鍵活性成分，見圖1。

2.5 化瘀丸治療乳腺癌基因的蛋白互作網(wǎng)絡(luò)

將化瘀丸與乳腺癌差異基因取交集后，通過Cytoscape3.2.1軟件繪制PPI網(wǎng)絡(luò)圖，藍(lán)色代表化瘀丸與乳腺癌差異基因取交集后的所有靶點(diǎn)，黃色代表核心基因，紅色代表目標(biāo)基因（見本文OSID碼）。借助BisoGenet和CytoNCA插件，分別設(shè)置度中心性（DC）＞61、介度中心性（BC）＞600，構(gòu)建拓?fù)渚W(wǎng)絡(luò)，篩選出目標(biāo)基因MYC，見圖2。

注：黃色為靶點(diǎn)基因，紅色為化瘀丸共同有效成分，藍(lán)色為元胡有效成分，綠色為白芷有效成分，深藍(lán)為王不留行有效成分，紫色為砂仁有效成分，玫紅為拳參有效成分，橘黃為桃仁有效成分，淡黃為黃芪有效成分，淡藍(lán)為郁金有效成分，深綠為紅參有效成分

圖2 化瘀丸治療乳腺癌基因拓?fù)渚W(wǎng)絡(luò)（DC＞61，BC＞600）

2.6 GO功能富集分析

利用R軟件，將化瘀丸治療乳腺癌的基因進(jìn)行GO富集分析。根據(jù)富集顯著程度，繪制BP（前20條）、CC（3條）、MF（前20條）柱狀圖。BP主要涉及機(jī)械刺激反應(yīng)、對含嘌呤化合物的反應(yīng)、抗氧化反應(yīng)、成纖維細(xì)胞增殖的調(diào)節(jié)、脂肪酸反應(yīng)等（見圖3）；CC主要涉及RNA聚合酶Ⅱ轉(zhuǎn)錄因子復(fù)合物、核轉(zhuǎn)錄因子復(fù)合體和轉(zhuǎn)錄因子復(fù)合物（見圖4）；MF主要涉及抗氧化活性、過氧化物酶的活性、氧化還原酶活性（以過氧化物為受體）、趨化因子受體結(jié)合、谷胱甘肽轉(zhuǎn)移酶活性、激活轉(zhuǎn)錄因子結(jié)合等（見圖5）。

圖3 化瘀丸治療乳腺癌靶點(diǎn)GO富集分析BP條目（前20條）

圖4 化瘀丸治療乳腺癌靶點(diǎn)GO富集分析CC條目

2.7 KEGG通路富集分析及基因通路網(wǎng)絡(luò)圖

利用R軟件對化瘀丸治療乳腺癌的基因進(jìn)行富集分析，共富集到42條通路。將KEGG通路導(dǎo)入Cytoscape3.2.1軟件，取前20條通路。根據(jù)Degree值對通路的富集程度，繪制通路基因網(wǎng)絡(luò)圖（見 圖6）。其中富集靶點(diǎn)較多的通路有：液體剪切應(yīng)力和動脈粥樣硬化通路（Fluid shear stress and atherosclerosis）、乙肝通路（Hepatitis B）和卡波西肉瘤相關(guān)的皰疹病毒感染通路（Kaposi sarcoma-associated herpesvirus infection）均富集8個靶點(diǎn)，IL-17信號通路（IL-17 signaling pathway）、細(xì)胞衰老信號通路（Cellular senescence）、肝癌通路（Hepatocellular carcinoma）、愛潑斯坦-巴爾病毒感染（Epstein-Barr virus infection）、人T細(xì)胞白血病病毒1型感染（Human T-cell leukemia virus 1 infection）通路均富集7個靶點(diǎn)。KEGG通路對應(yīng)的前5位靶點(diǎn)基因分別是JUN、FOS、E2F1、E2F2、MYC。

圖5 化瘀丸治療乳腺癌靶點(diǎn)GO富集分析MF條目（前20條）

圖6 化瘀丸治療乳腺癌KEGG通路與基因網(wǎng)絡(luò)

2.8 化瘀丸藥效成分槲皮素與5種蛋白的分子對接

將化瘀丸活性成分槲皮素與圖6中富集前5位的核心靶點(diǎn)蛋白使用AutoDock軟件進(jìn)行分子對接，對接結(jié)果見圖7。分子與靶點(diǎn)蛋白的最低結(jié)合能均為0，說明配體和受體均可以自發(fā)結(jié)合。

圖7 化瘀丸有效成分槲皮素與核心靶點(diǎn)分子對接模式圖

3 討論

乳腺癌是中國女性最常見的惡性腫瘤，發(fā)病率位居中國女性惡性腫瘤首位，發(fā)病年齡日趨年輕化[5]。國家級名老中醫(yī)郁仁存教授認(rèn)為痰、毒、瘀作為病理產(chǎn)物貫穿乳腺癌發(fā)生和發(fā)展的整個過程，針對痰瘀毒結(jié)的主要病因病機(jī)而確立的化痰、祛瘀、解毒治法貫穿乳腺癌治療的始終，創(chuàng)立化瘀丸常加減應(yīng)用于乳腺癌的辨證用藥，收效甚佳?；A(chǔ)研究結(jié)果顯示，化瘀丸含藥血清可提高Lewis細(xì)胞內(nèi)阿霉素藥物濃度，具有化療增效作用[6]。本文從網(wǎng)絡(luò)藥理學(xué)研究出發(fā)，從分子機(jī)制角度，闡述化瘀丸治療乳腺癌的作用機(jī)制，為進(jìn)一步研究提供理論依據(jù)。

通過網(wǎng)絡(luò)藥理學(xué)的預(yù)測及分析，選擇出核心成分，為臨床化瘀丸配伍關(guān)系提供參考依據(jù)。通過預(yù)測，槲皮素可能為化瘀丸治療乳腺癌的有效成分?，F(xiàn)代藥理研究表明，槲皮素在體內(nèi)外能夠有效抑制多種腫瘤細(xì)胞的生長、遷移和轉(zhuǎn)移，以及血管生成，并誘導(dǎo)腫瘤細(xì)胞發(fā)生凋亡、壞死，用于治療肝癌、胃癌、結(jié)腸癌、乳腺癌等[7-8]。同時，槲皮素對放化療具有增敏和緩解不良反應(yīng)的功效[9]。由此可見，槲皮素具有良好的抗癌功效，并且可能成為抗癌增敏劑，從而應(yīng)用于臨床。

本研究篩選化瘀丸的100個活性成分，其中具有潛在治療乳腺癌作用的有7個有效成分，槲皮素位于活性成分的核心位置，有關(guān)槲皮素作用于乳腺癌的相關(guān)機(jī)制分析文獻(xiàn)證實(shí)了本研究的結(jié)果[10-12]。此外，山柰酚、β-谷甾醇也可以在體外抑制乳腺癌細(xì)胞的增殖，誘導(dǎo)乳腺癌細(xì)胞凋亡[13-14]。異鼠李素不僅可以同自噬/有絲分裂抑制劑結(jié)合后用于治療三陰性乳腺癌，而且還能抑制細(xì)胞生長、產(chǎn)生細(xì)胞毒性、增加氧化應(yīng)激、阻滯細(xì)胞周期等[15]。因此，本研究認(rèn)為化瘀丸可以通過多種途徑和方式治療乳腺癌。

通過拓?fù)渚W(wǎng)絡(luò)的構(gòu)建，篩選出排名前5位的候選基因是JUN、FOS、E2F1、E2F2和MYC。原癌基因c-fos作為FOS的同系物，是最早發(fā)現(xiàn)的即早基因之一。它和JUN家族基因以同源或異源二聚體的形式結(jié)合在一起，調(diào)控細(xì)胞的生命活動，其中c-fos和c-Jun組成的異源二聚體是最穩(wěn)定的形式，并在腫瘤發(fā)生、轉(zhuǎn)化、血管生成、侵襲和轉(zhuǎn)移中具有重要作用[16]。E2F轉(zhuǎn)錄因子家族在細(xì)胞的增殖分化及凋亡中起著非常重要的調(diào)節(jié)作用。其中E2F1是E2F家族重要的轉(zhuǎn)錄因子，在細(xì)胞周期的調(diào)節(jié)和凋亡中發(fā)揮重要作用。有相關(guān)文獻(xiàn)報道，E2F1在正常組織中的表達(dá)明顯低于癌組織，在癌組織中的表達(dá)與腫瘤分級、淋巴結(jié)轉(zhuǎn)移及腫瘤TNM分期有關(guān)，可見E2F1在乳腺癌中的表達(dá)與乳腺癌的發(fā)生發(fā)展密切相關(guān)[17]。E2F2蛋白也在乳腺癌組織中高表達(dá)，其表達(dá)水平與乳腺癌的發(fā)生、發(fā)展及生存預(yù)后密切相關(guān)，可以為乳腺癌分子診斷與靶向治療提供參考[18]。MYC癌基因作為近年來分子生物學(xué)的研究熱點(diǎn)之一，已被證實(shí)在人類許多腫瘤的惡性進(jìn)展中起重要作用。MYC在乳腺癌組織中的高表達(dá)與腫瘤直徑、細(xì)胞學(xué)分級、病理分期、激素受體狀態(tài)、HER2狀態(tài)有關(guān)；且MYC陽性患者5年內(nèi)復(fù)發(fā)轉(zhuǎn)移風(fēng)險增高，MYC陽性表達(dá)與乳腺癌預(yù)后相關(guān)[19]。

GO富集分析結(jié)果顯示，化瘀丸的功效與氧化應(yīng)激、轉(zhuǎn)錄因子的激活、趨化因子受體結(jié)合等有關(guān)。轉(zhuǎn)錄因子涉及多種多樣，有研究表明，Krüppel樣轉(zhuǎn)錄因子8（KLF8）在乳腺癌組織及乳腺癌細(xì)胞中均呈高表達(dá)，且與乳腺癌患者的雌激素受體（ER）和孕激素受體（PR）狀態(tài)及預(yù)后生存有關(guān)，可能作為促癌基因參與了乳腺癌的進(jìn)展，或可為乳腺癌預(yù)后判斷及分子靶向治療提供新的切入點(diǎn)[20]。高遷移率族（high-mobility group，HMG）盒轉(zhuǎn)錄因子1（HMG-box transcription factor 1，HBP1）基因的表達(dá)影響乳腺癌的發(fā)生發(fā)展，并與乳腺癌預(yù)后呈正相關(guān)[21]。信號轉(zhuǎn)導(dǎo)和轉(zhuǎn)錄激活因子6（signal transducers and activators of transcription 6，STAT6）的陽性表達(dá)在不同分子分型乳腺癌中無差異，與組織學(xué)分級、臨床分期及腋窩淋巴結(jié)轉(zhuǎn)移相關(guān)；細(xì)胞程序性死亡-配體1（PD-L1）的陽性表達(dá)在三陰型乳腺癌中較高，與組織學(xué)分級和臨床分期相關(guān)；非特殊型浸潤性乳腺癌細(xì)胞STAT6與PD-L1表達(dá)呈正相關(guān)[22]。信號轉(zhuǎn)導(dǎo)和轉(zhuǎn)錄激活因子3（signal transducers and activators of transcription 3，STAT3）是STAT家族成員之一，在細(xì)胞轉(zhuǎn)化、凋亡、分化和生長調(diào)控中具有關(guān)鍵作用[23]。在一組14種乳腺癌細(xì)胞系中，僅在三陰性乳腺癌（TNBC）細(xì)胞中發(fā)現(xiàn)了組成性激活的STAT3。由此可見，STAT3可能是治療TNBC的一個重要靶點(diǎn)[24]。趨化因子及受體對腫瘤細(xì)胞有調(diào)控作用，主要包括調(diào)控腫瘤細(xì)胞生長、腫瘤細(xì)胞代謝、腫瘤細(xì)胞干性、腫瘤細(xì)胞凋亡、腫瘤細(xì)胞侵襲等。同時趨化因子及其受體對腫瘤相關(guān)免疫細(xì)胞有調(diào)節(jié)作用，包括巨噬細(xì)胞、T淋巴細(xì)胞、骨髓來源的免疫抑制細(xì)胞及樹突狀細(xì)胞等[25]。

KEGG富集通路結(jié)果顯示，化瘀丸治療乳腺癌主要涉及病毒感染相關(guān)通路、IL-17信號通路、細(xì)胞衰老信號通路及其他與癌癥相關(guān)的信號通路。多項(xiàng)研究證實(shí)病毒與癌癥的發(fā)生發(fā)展關(guān)系密切，乳腺作為性激素作用的靶器官，人乳頭瘤病毒（HPV）在乳腺癌的發(fā)生發(fā)展中的作用也備受關(guān)注[26]。一項(xiàng)Meta分析顯示，HPV陽性人群乳腺癌的發(fā)病率明顯高于正常人群[27]。IL-17在腫瘤中的臨床意義存在很大爭議，不同研究表明，IL-17在腫瘤發(fā)生發(fā)展中不僅與腫瘤類型有關(guān)，在一定程度上與腫瘤微環(huán)境相關(guān)。IL-17可能通過參與三陰性乳腺癌腫瘤間質(zhì)中炎癥反應(yīng)和微血管形成，從而促進(jìn)腫瘤的生長、侵襲及遠(yuǎn)處轉(zhuǎn)移等惡性生物學(xué)行為的發(fā)生，導(dǎo)致乳腺癌不良的預(yù)后[28]。細(xì)胞衰老在乳腺癌的發(fā)生發(fā)展過程中起重要作用，包括復(fù)制性細(xì)胞衰老、癌基因誘導(dǎo)細(xì)胞衰老、治療誘導(dǎo)的細(xì)胞衰老。腫瘤細(xì)胞衰老相關(guān)分泌表型具有雙面性，既可抑制腫瘤又可促進(jìn)腫瘤進(jìn)展。腫瘤細(xì)胞衰老相關(guān)分泌表型也有雙向調(diào)節(jié)作用，在癌變早期主要發(fā)揮正向調(diào)節(jié)作用，而在腫瘤晚期主要發(fā)揮負(fù)向調(diào)節(jié)作用。不同類型的乳腺癌產(chǎn)生不同程度的衰老相關(guān)分泌表型，重塑腫瘤微環(huán)境，促進(jìn)乳腺癌進(jìn)展[29]。槲皮素可以通過抗骨相關(guān)細(xì)胞衰老作用治療雌激素缺乏骨質(zhì)疏松癥[30]。槲皮素可以下調(diào)IL-6/STAT3信號通路的活性，減少炎性因子水平，阻礙腫瘤細(xì)胞的增殖，促進(jìn)腫瘤細(xì)胞的凋亡，從而保護(hù)結(jié)腸炎性相關(guān)結(jié)腸癌小鼠組織損傷[31]。

本研究通過網(wǎng)絡(luò)藥理學(xué)探討了化瘀丸治療乳腺癌靶點(diǎn)的網(wǎng)絡(luò)關(guān)系，顯示化瘀丸的有效成分與疾病靶點(diǎn)之間的作用機(jī)制，提示化瘀丸治療乳腺癌具有多靶點(diǎn)、多途徑、多通路的特點(diǎn)，為進(jìn)一步實(shí)驗(yàn)研究提供指導(dǎo)。本研究也存在一定的局限性，如有效成分相關(guān)靶點(diǎn)尋找和確認(rèn)方法需要進(jìn)一步優(yōu)化，化瘀丸與IL-17、細(xì)胞衰老、癌癥相關(guān)等信號通路的關(guān)系需要實(shí)驗(yàn)性證據(jù)，本課題組將在后續(xù)的研究中進(jìn)一步進(jìn)行實(shí)驗(yàn)驗(yàn)證。

[1] Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA: A Cancer Journal for Clinicians, 2018,68(6):394-424.

[2] 楊中,徐詠梅,張青,等.化瘀丸聯(lián)合化療治療晚期卵巢癌臨床研究[J].中國中醫(yī)藥信息雜志,2008,15(7):17-19.

[3] GAO L, WANG X D, NIU Y Y, et al. Molecular targets of Chinese herbs: a clinical study of hepatoma based on network pharmacology[J]. Scientific Reports, 2016,6:24944.

[4] 任雪童,狄紫蕊,付冰倩,等.基于網(wǎng)絡(luò)藥理學(xué)聯(lián)合GEO芯片分析百合烏藥湯治療胃癌的分子靶點(diǎn)及作用機(jī)制[J].中醫(yī)藥導(dǎo)報,2021,27(3):125-130.

[5] 赫捷,陳萬青,李霓,等.中國女性乳腺癌篩查與早診早治指南(2021,北京)[J].中華腫瘤雜志,2021,43(4):357-382.

[6] 高磊,郝彩霞,楊國旺,等.化瘀丸含藥血清對Lewis肺癌細(xì)胞內(nèi)阿霉素藥物濃度的影響[J].中醫(yī)雜志,2020,61(9):810-813.

[7] TANG S M, DENG X T, ZHOU J, et al. Pharmacological basis and new insights of quercetin action in respect to its anti-cancer effects[J]. Biomedicine & Pharmacotherapy, 2020,121:109604.

[8] Rauf A, Imran M, Khan I A, et al. Anticancer potential of quercetin: a comprehensive review[J]. Phytotherapy Research, 2018,32(11):2109-2130.

[9] Brito A F, Ribeiro M, Abrantes A M, et al. Quercetin in cancer treatment, alone or in combination with conventional therapeutics[J]. Current Medicinal Chemistry, 2015,22(26):3025-3039.

[10] Elsayed A M, Sherif N M, Hassan N S, et al. Novel quercetin encapsulated chitosan functionalized copper oxide nanoparticles as anti-breast cancer agent via regulating p53 in rat model[J]. International Journal of Biological Macromolecules, 2021,185:134-152.

[11] 蔣多晨,金樂,陳洪曉,等.槲皮素通過靶向GAS5/Notch1信號通路促進(jìn)乳腺癌細(xì)胞凋亡的研究[J].中國藥理學(xué)通報,2021,37(5):637-644.

[12] Schr?der L, Marahrens P, Koch J G, et al. Effects of green tea, matcha tea and their components epigallocatechin gallate and quercetin on MCF-7 and MDA-MB-231 breast carcinoma cells[J]. Oncology Reports, 2019,41(1):387-396.

[13] 李鷹飛,曾世彬,邱華平,等.山柰酚對乳腺癌細(xì)胞增殖及誘導(dǎo)凋亡的作用研究及機(jī)制[J].中國臨床藥理學(xué)雜志,2020,36(22):3679-3682.

[14] 李元鳳.β-谷甾醇-D-葡萄糖苷對乳腺癌細(xì)胞增殖、凋亡和侵襲轉(zhuǎn)移的影響[D].重慶:西南大學(xué),2019.

[15] 閆姍,楊林燕,馬悅,等.異鼠李素在乳腺癌中的研究進(jìn)展[J].青?？萍?2020,27(2):69-72.

[16] 彭洪薇,師銳贊,袁向飛,等.乳腺癌耐藥細(xì)胞中c-fos抗凋亡作用的研究[J].中國癌癥雜志,2014,24(8):581-588.

[17] 王宏坤,鄭繪霞,梁建芳,等.ECT2和E2F1蛋白在乳腺癌中的表達(dá)及臨床意義[J].山西醫(yī)科大學(xué)學(xué)報,2017,48(8):852-855.

[18] 段佳君,張勇,張亞杰.E2F家族蛋白在乳腺癌組織中的表達(dá)特征及其與預(yù)后的相關(guān)性[J].診斷病理學(xué)雜志,2020,27(9):653-657.

[19] 姜丹丹.MYC在乳腺癌組織中的表達(dá)及其對p53缺失的乳腺癌細(xì)胞增殖和凋亡的影響[D].青島:青島大學(xué),2019.

[20] 翟秧秧,師鑫鵬,羅曉勇.Krüppel樣轉(zhuǎn)錄因子8在乳腺癌中的表達(dá)及臨床意義[J].中國普外基礎(chǔ)與臨床雜志,2021,28(10):1320-1327.

[21] 汪軍梅,張毅強(qiáng),黃燕,等.高遷移率族盒轉(zhuǎn)錄因子1與乳腺癌相關(guān)性的生物信息學(xué)分析[J].中國生物制品學(xué)雜志,2021,34(4):437-441.

[22] 石蓉,高琨,楊承綱.信號傳導(dǎo)與轉(zhuǎn)錄活化因子6和程序性死亡配體-1在乳腺癌組織中的表達(dá)及意義[J].中國腫瘤臨床與康復(fù),2021,28(3):263-267.

[23] Aleskandarany M A, Agarwal D, Negm O H, et al. The prognostic significance of STAT3 in invasive breast cancer: analysis of protein and mRNA expressions in large cohorts[J]. Breast Cancer Research and Treatment, 2016,156(1):9-20.

[24] 高碩,王金花.STAT3信號通路在三陰性乳腺癌中的研究進(jìn)展[J].中國細(xì)胞生物學(xué)學(xué)報,2020,42(12):2266-2273.

[25] 王婉玉,呂曉希,胡卓偉,等.趨化因子及其受體在乳腺癌中的研究進(jìn)展[J].藥學(xué)學(xué)報,2021,56(6):1521-1531.

[26] 張煜,劉琨.人乳頭瘤病毒感染與乳腺癌超聲影像學(xué)特征的相關(guān)性研究[J].中國優(yōu)生與遺傳雜志,2020,28(9):1110-1112.

[27] LI N, BI X, ZHANG Y, et al. Human papillomavirus infection and sporadic breast carcinoma risk: a meta-analysis[J]. Breast Cancer Research and Treatment, 2011,126(2):515-520.

[28] 周瓊,董杰,沈夏波,等.白介素-17在三陰性乳腺癌中的表達(dá)及其與新生血管和預(yù)后的相關(guān)性[J].腫瘤學(xué)雜志,2021,27(1):53-57.

[29] 褚玉新,胡欽勇,宋啟斌.衰老相關(guān)分泌表型在乳腺癌進(jìn)展中的作用[J].中國腫瘤,2018,27(5):364-368.

[30] 顧藝婧,傅稼耀,武文婧,等.槲皮素通過抗骨相關(guān)細(xì)胞衰老作用治療雌激素缺乏骨質(zhì)疏松癥的初步研究[J].同濟(jì)大學(xué)學(xué)報(醫(yī)學(xué)版),2019,40(3):274-280.

[31] 蔡劍輝,任約翰,厲金雷,等.槲皮素調(diào)控白細(xì)胞介素-6/信號轉(zhuǎn)導(dǎo)及轉(zhuǎn)錄激活蛋白3信號通路在結(jié)腸炎性相關(guān)結(jié)腸癌小鼠模型的機(jī)制研究[J].中國藥物與臨床,2020,20(6):896-899.

Discussion on the Mechanism and ofPills in the Treatment of Breast Cancer Based on Network Pharmacology Combined with GEO Database and Molecular Docking Technology

LIU Min1, LI Yi-nan2, ZHANG Qing1*

(1. Beijing Hospital of Traditional Chinese Medicine, Capital Medical University, Beijing 100010, China; 2. Beijing University of Chinese Medicine, Beijing 100029, China)

To study the effective components and action mechanism ofPills in the treatment of breast cancer through the combination of GEO database differential gene analysis and molecular docking technology.The TCMSP database was used to retrieve the information of 12 types of Chinese materia medica inPills, and the active components and molecular targets ofPills were obtained.GSE139038 chip data were obtained from the GEO database through searching breast cancer, and the difference genes were analyzed through R software. Targets for the treatment of breast cancer were screened through DrugBank, OMIM, TTD and CTDdatabases.Gene annotation was performed on the target using the UniProt database. The composition-target network map ofPills and breast cancer was constructed by Cytoscape3.2.1 software, and the core target topological network was plotted using plug-ins BisoGenet and CytoNCA. GO enrichment analysis and KEGG enrichment analysis were performed for drugs and disease genes. Finally, the components with high degree of the network were under molecular docking with the core target.From the GEO database, 40 genes with the most significant difference of up-regulated and down-regulated breast cancer were obtained. 7 effective active components ofPills were obtained, including kaempferol, sitosterol, isorhamnetin, hederagenin, quercetin, beta-sitosterol, stigmasterol, which acted on 32 gene targets including JUN, FOS, MYC, E2F1, E2F2, etc. GO function and KEGG pathway enrichment analysis showed that oxidative stress, activation of the transcription factor, chemokinine receptor binding and other aspects were mainly involved in the treatment of breast cancer byPills. The pathways involved in virus infection, cellular senescence signaling pathway, IL-17 signaling pathway, and pathways in hepatocellular carcinoma.The important active components ofPills in the treatment of breast cancer are quercetin, etc., and virus infection, IL-17 and other signaling pathways of cancers are involved.

Pills; breast cancer; network pharmacology; bioinformatics; molecular docking; mechanism

R285；R273.79

A

2095-5707(2022)05-0001-08

10.3969/j.issn.2095-5707.2022.05.001

劉敏,李羿南,張青.基于網(wǎng)絡(luò)藥理學(xué)聯(lián)合GEO數(shù)據(jù)庫及分子對接技術(shù)探討化瘀丸治療乳腺癌的作用機(jī)制[J].中國中醫(yī)藥圖書情報雜志,2022,46(5):1-8.

恩施州科技計劃研究與開發(fā)項(xiàng)目（JCY2019000040）

劉敏，E-mail: liumin624@ccmu.edu.cn

張青，E-mail: zhangqingys@sina.com

（2021-11-24）

（修回日期：2022-01-14；編輯：魏民）