基于網絡藥理學的多西他賽治療絨毛膜癌的作用機制研究

金雪寧, 林 飛, 馮心怡, 李 岑, 崔樹娜

(1. 揚州大學醫學院, 江蘇 揚州, 225000; 2. 揚州大學附屬醫院 婦產科, 江蘇 揚州, 225009)

絨毛膜癌(CC)簡稱絨癌，屬于高度惡性滋養細胞腫瘤，有50%的妊娠CC發生在葡萄胎之后，其余的發生在自然流產和人工流產、正常足月妊娠或早產以及異位妊娠之后[1-2]。目前常用的化療藥物有甲氨蝶呤、環磷酞胺等，但這些化療藥物毒副作用大。多西他賽與紫杉醇都屬于紫杉醇類藥物，但多西他賽抗腫瘤活性更高，抗癌譜更廣，其機制是加強微管蛋白聚合作用和抑制微管解聚作用，導致形成穩定的非功能性微管束，因而破壞腫瘤細胞的有絲分裂[3], 在臨床中被廣泛應用于治療難治性卵巢癌等[4-7]。作者的前期體外實驗[8]證實多西他賽能抑制CC細胞的增殖，誘導細胞凋亡并明顯抑制其侵襲能力。研究[9]報道多西他賽能降低對常規化療藥物不敏感的CC患者血清人絨毛膜促性腺激素(HCG)水平，但多西他賽治療CC的作用機制還缺乏深入研究。網絡藥理學融合了系統生物學、多向藥理學和計算機信息學等多門學科，建立“藥物-成分-基因-靶點-通路-疾病”的多層次網絡關系，探索天然藥物與復雜疾病的相互關系[10-11]。本研究基于前期實驗研究和臨床文獻結果，采用網絡藥理學手段探索多西他賽治療CC的潛在作用機制，現將結果報告如下。

1 材料與方法

1.1 多西他賽、CC的靶點收集

在Genecards數據庫(https://www.genecards.org/)中分別以“Docetaxol” “choriocarcinoma”為關鍵詞檢索靶點，并將其導入Excel表格中進行標準化處理，進行統計分析。

1.2 關鍵靶點篩選

將多西他賽所對應的靶點基因和CC相關的靶點基因進行映射，獲得兩者的共同基因，即多西他賽治療CC的作用靶點。

1.3 蛋白質互作(PPI)網絡構建

將篩選出的基因輸入String 11.0數據庫，將蛋白種類設置為人類(homo sapiens), “minimum required interaction score”參數設置為“highest confidence(0.900)”, 其余參數保持初始默認設置，導出PPI數據文件，建立蛋白與蛋白之間的網絡。將結果導入Cytoscape軟件，對結果進行可視化。

1.4 Hub基因篩選

通過Cytoscape軟件的Cytohubba插件中的MCC算法，計算并篩選排名前10位的多西他賽治療CC的Hub基因。

1.5 基因本體論(GO)富集和京都基因和基因組百科全書(KEGG)通路富集

將共同基因輸入DAVID 6.7數據庫(http://david.ncifcrf.gov/)進行GO富集分析和KEGG通路分析，研究多西他賽治療CC的關鍵靶點的生物功能以及參與的信號通路，以P<0.05和錯誤發現率(FDR)<0.05作為篩選條件，篩選具有顯著差異的生物過程和靶點通路。

2 結 果

2.1 靶點蛋白PPI網絡的構建

在Genecards數據庫分別以“Docetaxol”和“choriocarcinoma”為關鍵詞檢索靶點，分別獲得953和1 137個靶點，將1 137個疾病靶點與953個藥物靶點輸入Venny2.1 軟件繪制韋恩圖，篩選出294個共同靶點，見圖1。將共同靶點在String 11.0 平臺分析的數據導入Cytoscape3.6.1軟件，構建PPI網絡模型。以Combined_score>0.9為篩選條件繪制PPI網絡圖，其中節點為靶點蛋白，邊為各蛋白間互作關系，節點連接邊數越多代表其在網絡中作用越重要，見圖2。通過Cytoscape 軟件插件Cytohubba計算排名前10位的Hub基因，其中TP53、AKT1、CASP8、STAT3等是多西他賽治療CC的關鍵基因。見圖3。

圖1 藥物和疾病靶點韋恩圖

圖2 交集靶點相互作用網絡關系

圖3 通過MCC算法計算排名前10位的多西他賽治療CC的hub基因

2.2 GO功能和KEGG富集分析

將294個預測靶點導入DAVID數據庫中進行GO富集分析和KEGG通路分析。GO富集分析中，參與的生物過程有409個，主要包括調控細胞增殖(細胞增殖調節)、細胞凋亡(細胞凋亡負調控、細胞死亡的調控、程序性細胞死亡的調控)、信號傳導(信號傳導的正調控)、調節分子代謝(大分子代謝及合成過程的正調控、有機物質的反應等)。KEGG通路富集分析確定了28條相關信號通路，主要涉及癌癥通路、細胞凋亡通路、Focal Adhesion通路和p53、ErbB、MAPK等信號通路。見圖4、圖5。

圖4 多西他賽潛在靶點的富集分析(生物過程)

原點由小到大代表分值逐漸升高,靶點富集數目多,通路顯著性逐漸增強; 顏色由藍變紅代表P值逐漸減小。圖5 多西他賽潛在靶點的KEGG通路注釋圖

3 討 論

CC屬于惡性滋養細胞癌，是由滋養層滋養細胞增殖引起的惡性轉化，導致這種轉化的分子機制尚未確定。CC可繼發于葡萄胎妊娠，也可繼發于非葡萄胎妊娠，病情發展迅速，病死率高，主要通過血行播散，轉移早且易擴散，肺是CC最易轉移的部位，其次為陰道、盆腔等，轉移部位皆有局部出血灶[1, 12]。

多西他賽也稱多西紫杉醇，是紫杉醇類化學藥物的衍生物，臨床應用廣泛，常被用于治療難治性卵巢癌等，但其應用于CC方面的研究較少[4-6]。作者前期實驗研究發現，多西他賽可以抑制CC細胞增殖的能力，并能誘導CC細胞凋亡，其作用濃度與時間呈正相關，并能明顯抑制CC細胞體外侵襲的能力[8]。對高度化療耐藥的CC細胞，每周使用單劑多西他賽可快速降低HCG水平(超過95%), 且對人體毒性最小，但其作用機制有待進一步深入研究。

本研究結果發現，多西他賽可靶向P53、AKT1、STAT3、VEGFA等294個基因，參與細胞凋亡與增殖的調控、分子合成與代謝的調控等主要生物過程。抑癌基因TP53編碼p53腫瘤蛋白，其能調節細胞周期和避免細胞癌變發生。任何導致TP53失活的突變都會引起其抑癌功能的喪失而導致癌癥發生[13]。研究[14]發現女性TP53突變攜帶者罹患CC的風險更高，且易產生化療耐藥，因此靶向調控p53對于多西他賽治療CC具有重要的臨床價值。蛋白激酶B1(AKT1)是AKT的一個亞型，是生存信號通路PI3K/AKT的關鍵分子，與腫瘤細胞的生長、增殖、促進細胞侵襲和轉移以及促進血管形成密切相關。研究[15]發現，在CC細胞系中，逆轉AKT1表達可能通過抑制PI3K/AKT信號通路及靶向負調控VEGFA表達，最終抑制CC細胞增殖、遷移和侵襲。信號轉導與轉錄激活因子3(STAT3)是一類具有信號轉導及轉錄調控雙重功能的蛋白質家族，激活后與DNA分子結合發揮作用[16]。多項研究[17-18]證實STAT3是參與惡性細胞轉化的信號轉導介質, IL-6/STAT3通路與CC細胞增殖密切相關，阻斷 STAT3信號通路能夠抑制CC細胞的增殖，并誘導CC細胞凋亡。

研究[19-20]表明, p53信號通路與腫瘤發展及預后密切相關，其可以激活多條信號通路，阻礙細胞周期，誘導細胞凋亡，促進細胞衰老以及DNA修復過程，最終抑制腫瘤生長。研究[21]發現p53和p21CIP1是p53信號通路中重要的靶基因，在CC JEG-3細胞中，細胞周期蛋白依懶性激酶抑制劑p21CIPI表達被激活后，使G2細胞周期蛋白依賴性激酶CDK1表達顯著降低，同時p53絲氨酸磷酸化水平增加, p53信號通路激活后使細胞周期阻滯在G2/M期，進而激活DNA損傷反應通路，抑制細胞增殖。

本研究借助網絡藥理學方法對多西他賽和CC進行靶標挖掘和通路分析，得到294個共同靶向基因和CC多條信號通路，包括TP53、AKT1、STAT3、VEGAF等基因以及細胞凋亡通路、p53通路等信號通路，證明了多西他賽可能通過調節多個靶基因、調控蛋白表達和信號傳導通路而發揮對CC的治療作用，初步揭示了多西他賽治療CC的作用機制，為后續臨床研究提供了新思路。