基于網絡藥理學和分子對接探討西黃丸治療宮頸癌的作用機制

鄧鴻飛，楊自力，羅毅

南京中醫藥大學附屬中西醫結合醫院，江蘇 南京 210000

宮頸癌作為威脅女性健康的頭號“殺手”，已成為威脅女性健康的重要因素，全世界每年約有25 萬女性因宮頸癌死亡[1-2]。作為傳統名方的西黃丸已廣泛應用于臨床，研究證明西黃丸輔助治療宮頸癌可以有效減輕西醫治療的不良反應，提高患者的生活質量，提高療效，延長生命[3]。

西黃丸記載于《外科證治全生集》，方由牛黃、麝香、乳香、沒藥4 味中藥組成，具有清熱活血、化痰消徵、解毒散結等功效[4]。根據現代藥理學研究結果，西黃丸具備一系列的藥理效應，這些效應包括抑制惡性腫瘤的增殖和傳播、誘導癌細胞發生凋亡、阻止新血管生成及調節免疫微環境等[3,5-8]。因此，西黃丸常被用于治療多種惡性腫瘤，如乳腺癌和胃癌等。但西黃丸輔助治療宮頸癌的機制尚不明確。

本研究基于網絡藥理學及分子對接的方法[9]來預測西黃丸防治宮頸癌的有效活性成分、作用靶點及參與調控的細胞信號轉導等，闡述西黃丸抗宮頸癌的作用機制，以此為后面的科學研究提供依據。

1 材料與方法

1.1 西黃丸有效成分及目標靶基因的篩選

在TCMSP 數據庫[10]（https://tcmsp-e.com/）中檢索牛黃、乳香、沒藥對應數據庫中的化合物，以口服生物利用度（OB）≥30%，類藥性（DL）≥0.18進行篩選，獲得上述3 味藥的有效成分。同時，在BATMAN-TCM 數據庫[11]（http://bionet.ncpsb.org/batman-tcm/）中，以分數閾值（score cutoff）≥30和P≤0.05 為條件篩選麝香的有效成分及靶點。

1.2 西黃丸有效成分–靶點網絡的構建

利用 UniProt（https://www.uniprot.org/）和BATMAN-TCM 數據庫獲得西黃丸的有效成分作用靶點基因名。通過Excel 構建網絡文件，導入Cytoscape 3.9.1 軟件構活性成分–靶點網絡。

1.3 宮頸癌疾病靶點的篩選

以“cervical cancer”為關鍵詞在GeneCards 人類基因數據庫（https://auth.lifemapsc.com/）和OMIM數據庫（https://www.omim.org/）中檢索宮頸癌的相關基因，合并去重后得到疾病基因靶點。

1.4 西黃丸–宮頸癌共同靶點篩選及蛋白相互作用（PPI）網絡的構建

將西黃丸–宮頸癌共同靶點導入String 12.0 數據庫（https://cn.string-db.org/），物種選擇“人類（homo sapiens）”最低相互作用域值≥0.9，獲得靶點相互作用網絡圖，將該圖導入Cytoscape 3.9.1 軟件構建PPI 網絡。并通過“network analyzer”插件網絡結構進行拓撲分析，根據介數中心性（BC）、緊密度（CC）、節點連接度（degree）篩選得到核心靶點。

1.5 基因本體論（GO）功能富集分析與京都基因與基因組百科全書（KEGG）通路富集分析

將西黃丸–宮頸癌共同靶點以Gene Symbol 的格式導入DAVID 數據庫（https://david.ncifcrf.gov/），進行GO 功能富集分析以及KEGG 基因通路富集分析。通過微生信在線平臺（http://www.bio informatics.com.cn）進行可視化處理，以解釋西黃丸治療宮頸癌的生物功能及相關信號通路過程。

1.6 “西黃丸–靶點–宮頸癌–通路”網絡的構建

通過KEGG 富集分析結果，選取與宮頸癌聯系緊密的通路，通過Cytoscape 3.9.1 軟件構建“西黃丸–靶點–宮頸癌–通路”網絡圖。

1.7 分子對接驗證

通過PDB 數據庫（https://www.rcsb.org/）下載核心靶點排名前10 位蛋白質結構，利用PyMol、Autodock-vina 軟件處理大分子蛋白結構，后保存為“pdbqt”格式文件。通過PubChem 數據庫（https://pubchem.ncbi.nlm.nih.gov/）得到西黃丸核心有效成分小分子的3D 結構，并通過Autodock Tools 軟件處理后保存為“pdbqt”格式文件。最后利用Autodock Vina 1.5.6 軟件將大分子和小分子進行對接，利用PyMol 軟件處理對接結果，建立對接相互作用圖。

2 結果

2.1 西黃丸活性成分及其靶點

檢索到西黃丸有效活性成分共70 個，對應的靶點共428 個，經篩選后，牛黃、麝香、乳香、沒藥有效成分分別為5、25、6、34 個，靶點分別有7、280、13、195 個，見表1。

2.2 西黃丸活性成分–靶點網絡圖

采用Cytoscape 3.9.1 軟件構建西黃丸活性成分–靶點網絡，共有502 個節點，1 221 條邊。degree值顯示活性成分中槲皮素（quercetin，degree＝147）、表雄酮（3beta-hydroxy-5alpha-androstan-17-one，degree＝64）、雄甾酮（androsterone，degree＝64）、3α-羥基-5α-17-酮（3alpha-hydroxy-5alpha-androstan-17-one，degree＝64）、5-順-環十五烷-L-酮（5-ciscyclopentadecen-1-one，degree＝63）等在西黃丸中具有主要作用，見圖1、表2。

圖1 西黃丸“藥物–活性成分–作用靶點”網絡圖Fig.1 Xihuang Pills “drug -active ingredient -target” network diagram

表2 西黃丸主要活性成分Table 2 Main active ingredient of Xihuang Pills

2.3 共同靶點及PPI 網絡

在GeneCards 和OMIM 數據庫中分別找到相關靶點1 062、495 個，合并去重后得到宮頸癌疾病靶點1 449 個，將西黃丸靶點與宮頸癌疾病靶點進行篩選得到交集基因靶點116 個；通過String 數據庫，最終得到114 個節點和1 121 條邊。通過Cytoscape 3.9.1 軟件生成PPI 網絡，該網絡共包含114 個節點和1 121 條邊。顏色越深、節點越大，degree 值越高，提示這些蛋白越重要。用插件Centiscape 2.2 對PPI 網絡進行拓撲分析，以CC 大于中位數0.556 064 306、BC 大于中位數0.003 372 877、degree 大于2 倍中位數值33，篩選得到腫瘤蛋白P53（TP53）、表皮生長因子受體（EGFR）、蛋白激酶B1（Akt1）、白細胞介素（IL）-6、原癌基因（MYC）、轉錄因子AP-1（JUN）、半胱氨酸蛋白酶3（CASP3）、IL-1B、腫瘤壞死因子（TNF）、雌激素受體1（ESR1）、B 淋巴細胞2（Bcl-2）、低氧誘導因子-1A（HIF-1A）、基質金屬蛋白酶9（MMP9）、V-Rel 網狀內皮增生病毒癌基因同源物A（RELA）等22 個核心靶點，見圖2～4。根據degree 排名前10 位的基因見表3。這些靶點在西黃丸治療宮頸癌中發揮主要作用。

圖2 西黃丸–宮頸癌PPI 網絡Fig.2 Xihuang Pills– cervical cancer PPI network

圖3 基于Cytoscape 構建的PPI 網絡圖Fig.3 PPI network diagram based on Cytoscape

圖4 西黃丸–宮頸癌PPI 關鍵靶點Fig.4 Key target PPI of Xihuang Pills– cervical cancer

表3 交集基因靶點信息Table 3 Intersection gene target information

2.4 GO 生物過程和KEGG 通路富集分析

交集基因通過DAVID 數據庫進行富集分析，獲得1 018 個GO 條目，包括生物過程（BP）805個，主要包括細胞凋亡過程正負調控、RNA 聚合酶啟動子的轉錄調控、對刺激反應的調控、細胞增殖調控、基因表達的調控等方面；細胞組成（CC）70個，涉及轉錄因子復合體、細胞質、線粒體、細胞核、染色體等；分子功能（MF）143 個，涉及蛋白質結合、泛素蛋白連接酶結合、轉錄因子結合、蛋白質異構化活性等、蛋白激酶活性等。根據P＜0.05，分別取排名前10 位的的條目進行可視化處理，見圖5。

圖5 GO 功能富集分析Fig.5 GO functional enrichment analysis histogram

通過KEGG 富集分析共獲得164 條通路，根據P＜0.01，取排名前40 位的通路進行可視化處理，見圖6?？v軸表示富集途徑名稱，橫軸表示基因比例。氣泡顏色深淺與P值呈正比，顏色越深P值越??；氣泡大小與基因數量呈正比，氣泡越大，數量越多。結果顯示涉及的通路主要有胰腺癌（pancreatic cancer）、膀胱癌（bladder cancer）、前列腺癌（prostate cancer）、結直腸癌（colorectal cancer）、HIF-1、磷脂酰肌醇3-激酶（PI3K）/Akt、TP53 等信號通路。

圖6 KEGG 通路富集分析Fig.6 KEGG pathway enrichment analysis bubble diagram

選取前50 位富集通路中的IL-17、TNF、PI3K/Akt、p53、HIF-1 信號通路，構建“西黃丸–靶點–宮頸癌–通路”網絡，見圖7。

圖7 西黃丸–靶點–宮頸癌–通路互作網絡Fig.7 Xihuang Pills– target– cervical cancer– pathway interaction network

2.5 分子對接驗證

對篩選得到的10 個核心成分及10 個核心靶點進行分子對接。并用結合能結果作圖，分子對接結合能熱圖見圖8（在BATMAN-TCM 數據庫里顯示雄甾酮與3α-羥基-5α-17-酮顯示為麝香的2 種有效成分，在PubChem 數據庫顯示其為同一化學結構）。最低結合能≤?1.2 kcal/mol（1 cal＝4.4 J），即表明藥效分子與蛋白對接效果良好，結合能越低，分子與蛋白結合能力越好。其中結合自由能排名靠前的為表雄酮與ESR1（?10.45 kcal/mol）、雄甾酮（3α-羥基-5α-17-酮）與ESR1（?9.97 kcal/mol）、4-雄烯二酮與ESR1（?9.95 kcal/mol）、雄甾酮（3α-羥基-5α-17-酮）與EGFR（?9.55 kcal/mol）、雌二醇與ESR1（?9.13 kcal/mol）、4-雄烯二酮與 Akt1（?8.82 kcal/mol）。通過Pymol 軟件對上述的結果進行處理，見圖9。

圖8 分子對接結合能熱圖Fig.8 Molecular docking binding energy heat map

圖9 分子對接模式圖Fig.9 Molecular docking model diagram

3 討論

通過活性成分–靶點網絡圖得知西黃丸關鍵活性成分有表雄酮、槲皮素、雄甾酮、3α-羥基-5α-17-酮、5-順-環十五烷-L-酮等?，F代研究表明，槲皮素是一種具有高抗癌活性的天然類黃酮，特別是對于與人乳頭瘤病毒（HPV）相關的宮頸癌等癌癥具有治療作用[12]。槲皮素可通過誘導ESR 通路促進宮頸癌HELA 細胞的凋亡，達到抗腫瘤目的[13]。研究表明，槲皮素具有作為化學預防/抗腫瘤劑的內在潛力，并通過逆轉上皮–間充質轉化（EMT）信號傳導來減緩腫瘤進展[14]。惡性腫瘤中的重組人泛素結合酶E2-S（UBE2S）通過EMT 信號傳導有助于細胞運動，并能被槲皮素逆轉，且UBE2S 可能有助于宮頸癌中HIF-1A 信號傳導。結果表明，槲皮素通過降低UBE2S表達對宮頸癌的轉移抑制作用[15]。這些研究不僅為抑制宮頸癌的轉移提供了新的見解，也提示槲皮素在治療宮頸癌中具有重要價值。表雄酮可以在機體內與脫氫表雄酮進行轉化，已經有研究推測脫氫表雄酮可能通過抑制腫瘤細胞的增殖和遷移來預防或治療宮頸癌[16-17]?，F在研究已經表明，雄甾酮具有很好的細胞毒活性，具有潛在的抗癌作用。因此，西黃丸能夠通過相關通路調控腫瘤細胞的增殖轉移，誘導凋亡。

通過西黃丸–宮頸癌靶點網絡圖推測西黃丸治療宮頸癌的主要靶點包括TP53、EGFR、Akt1、IL-6、MYC 基因、JUN、CASP3、TNF 等。TP53 基因的體細胞突變是人類癌癥中最常見的改變之一，TP53 突變通過影響DNA 損傷修復途徑，減緩細胞周期，造成染色質重塑和細胞凋亡，從而影響腫瘤的發展過程，是腫瘤治療潛在的預后和預測標志物，也是藥物干預的靶標之一[18-19]。宮頸癌的發生與TP53 基因密切相關[20]，EGFR 水平的升高能夠促進實體瘤生長[21]。EGFR/Akt 信號傳導的過度激活使NANOG 腫瘤細胞對細胞毒性T 淋巴細胞（CTL）殺傷具有抗性。NANOG/LC3B-p/EGFR 軸存在于各種類型的人類癌癥中，并與宮頸癌患者的總生存率呈負相關[22]。研究表明，通過抑制EGFR 基因的表達水平，能夠抑制宮頸癌細胞的細胞活力、增殖、遷移和癌癥干細胞樣表型，從而達到治療宮頸癌的復發和轉移的目的[23]。研究表明，circ-Akt1在宮頸癌樣品和細胞系中的高表達，體內檢測證實circ-Akt1 會促進宮頸癌腫瘤生長[24]。研究表明，IL-6 的表達與腫瘤大小、婦科腫瘤分期標準組織學分級和腫瘤分化顯著相關。IL-6 表達陽性與預后不良顯著相關，其被認為是宮頸癌預后和潛在治療靶點的有價值的生物標志物[25]。有學者提出MYC 基因的過表達促進了侵襲性宮頸癌的轉移[26]。有研究表明，TNF-α 啟動子區域參與宮頸癌的發展[27]。

通過GO 富集分析可知，西黃丸治療宮頸癌的生物功能主要涉及基因轉錄調控、信號調節、RNA轉錄調控及對藥物刺激的反應等；參與了蛋白質結合和酶結合等分子功能。KEGG 富集分析顯示，西黃丸治療宮頸癌涉及了多條信號通路，包括IL17、TNF、PI3K/Akt、p53、HIF-1 等信號通路。IL-17 被認為是炎癥的重要因素。試驗研究證明，IL-17 通過激活Janus 激酶2（JAK2）/信號轉導與轉錄激活因子3（STAT3）和PI3K/Akt/核因子-κB（NF-κB）信號傳導促進宮頸癌的進展[28]。TNF 信號通路調控著細胞的生長和侵襲轉移，對宮頸癌疾病的發展有著重要的作用[29]。PI3K/Akt 信號通路調控著人宮頸癌細胞系（CaSki、HELA）的細胞活力和遷移，從而影響宮頸癌的復發和轉移[30]。研究表明，T 細胞免疫球蛋白黏蛋白-1（TIM-1）過表達通過調節PI3K/Akt/p53 和PI3K/Akt/mTOR 信號通路促進細胞遷移和侵襲，抑制宮頸癌細胞凋亡，可能是該病的候選診斷生物標志物[31]。研究表明，HIF 信號通路可能在宮頸癌中起關鍵作用[32]，參與宮頸癌細胞的增殖和調亡[33-35]。

本研究采用網絡藥理學結合分子對接的方法，篩選分析西黃丸和宮頸癌的相關靶點，并對藥物的活性成分和重要靶點進行分子對接，對接結果均良好，由此可以推斷西黃丸發揮主要作用的活性成分有槲皮素、表雄酮、雄甾酮等；潛在治療靶點是TP53、EGFR、Akt1、IL-6、MYC、TNF 等；相關調控通路有IL17、TNF、PI3K/Akt、p53、HIF-1、癌癥等信號通路。由此得知西黃丸治療宮頸癌的作用靶點是多途徑、多層面的，為西黃丸的后續實驗研究提供了更多的理論依據。

利益沖突所有作者均聲明不存在利益沖突