楊夢(mèng)霞,毛 昀,周 磊,朱世杰,王寧軍,蘆殿榮
(1.北京中醫(yī)藥大學(xué),北京 100029 2.中國(guó)中醫(yī)科學(xué)院望京醫(yī)院,北京 100102)
骨轉(zhuǎn)移癌是乳腺、肺和前列腺癌等惡性腫瘤常見(jiàn)并發(fā)癥,可引起疼痛、脊髓壓迫、病理性骨折等一系列臨床表現(xiàn),給晚期腫瘤患者的生存質(zhì)量帶來(lái)了嚴(yán)重不良影響[1]。目前多采用手術(shù)、放療、雙膦酸鹽類(lèi)藥物治療等延緩骨轉(zhuǎn)移,但總體療效欠佳,且不良反應(yīng)明顯[2]。近年來(lái)中醫(yī)藥以其副作用少、療效佳等優(yōu)勢(shì)在骨轉(zhuǎn)移癌治療中取得了較好的臨床效果,在臨床上受到了廣泛關(guān)注[3]。骨轉(zhuǎn)移癌屬于中醫(yī)學(xué)骨痹、骨蝕、骨瘤、脛陰疽等范疇,其病機(jī)多為正虛邪實(shí)。其中正虛多為腎虛、陽(yáng)虛,邪實(shí)多為寒、痰、血瘀和癌毒等,故在治療上應(yīng)重視扶正,并聯(lián)合溫陽(yáng)活血、散寒化痰,攻毒抑癌等方法[4]?;诖?,我院自擬溫陽(yáng)活血湯治療陰寒凝滯型骨轉(zhuǎn)移癌,取得了較好的臨床療效,但尚未對(duì)其作用機(jī)制進(jìn)行研究。近年來(lái)網(wǎng)絡(luò)藥理學(xué)因其融合了系統(tǒng)藥理、生物學(xué)及計(jì)算機(jī)信息學(xué)等內(nèi)容的獨(dú)特優(yōu)勢(shì),發(fā)展迅速,可通過(guò)多層次相互作用的中藥、化合物、靶點(diǎn)、疾病基因和信號(hào)通路網(wǎng)絡(luò)示意圖的構(gòu)建來(lái)探討中醫(yī)藥作用機(jī)理,是較為全新的一種方法。本研究擬運(yùn)用網(wǎng)絡(luò)藥理學(xué)方法探討溫陽(yáng)活血湯治療陰寒凝滯型骨轉(zhuǎn)移癌的分子機(jī)制,以期能夠在基礎(chǔ)研究中為后來(lái)學(xué)者提供新思路,在臨床上為醫(yī)生用藥提供一定參考價(jià)值。
1.1溫陽(yáng)活血湯有效成分的篩選:在數(shù)據(jù)庫(kù)TCMSP中,以口服生物利用度OB、類(lèi)藥性指數(shù)DL分別不小于30%和0.18為篩選條件,收集溫陽(yáng)活血湯8味中藥化合物成分;此外,以TCMID、TCM Database@Taiwan為補(bǔ)充數(shù)據(jù)庫(kù),并聯(lián)合文獻(xiàn)數(shù)據(jù)庫(kù)如知網(wǎng)、萬(wàn)方等檢索溫陽(yáng)活血湯相關(guān)成分。
1.2藥物有效成分作用靶點(diǎn)的預(yù)測(cè):使用數(shù)據(jù)庫(kù)TCMSP查找主要有效成分相關(guān)靶蛋白,并結(jié)合Uniprot數(shù)據(jù)庫(kù)查詢(xún)其對(duì)應(yīng)的Gene Symbol;在PubChem中獲取化合物有效成分的SMILES結(jié)構(gòu),并在數(shù)據(jù)庫(kù)Swiss Target Prediction中將上述獲得的結(jié)構(gòu)輸入,且以“Homo sapiens”為篩選條件進(jìn)行相關(guān)靶點(diǎn)的查詢(xún),保留匹配度高的靶點(diǎn),之后將兩個(gè)數(shù)據(jù)庫(kù)中獲得的靶點(diǎn)進(jìn)行去重,并在MCODE插件中進(jìn)行核心基因的篩選。
1.3骨轉(zhuǎn)移癌潛在治療靶點(diǎn)篩選:本研究將“Gene”和“Homo sapiens”設(shè)置為檢索條件,在疾病基因數(shù)據(jù)庫(kù)DisGeNet和Drugbank中尋找骨轉(zhuǎn)移癌相關(guān)基因,最后去重。
1.4溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的主要潛在靶點(diǎn):將藥物相關(guān)作用靶點(diǎn)和疾病基因在Venny2.1中進(jìn)行映射,篩選出溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的靶點(diǎn)有105個(gè)。之后在STRING11.0數(shù)據(jù)庫(kù)中,輸入上述靶點(diǎn)并將“Homo sapiens”作為選擇物種,可獲得PPI互作網(wǎng)絡(luò),并在Cytoscape3.6.0中進(jìn)行可視化。
1.5靶基因的功能分析:在DAVID數(shù)據(jù)庫(kù)中對(duì)溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的潛在靶點(diǎn)做GO和KEGG富集分析,其中GO分析包含了生物過(guò)程BP、細(xì)胞成分CC和分子功能MF。
2.1溫陽(yáng)活血湯主要活性成分:在數(shù)據(jù)庫(kù)TCMSP、TCMID和TCM Database@Taiwan中檢索到溫陽(yáng)活血湯所有化學(xué)成分共計(jì)707個(gè)。以O(shè)B、DL分別不小于30%和0.18為篩選閾值,并結(jié)合文獻(xiàn)檢索,合計(jì)有效成分共109個(gè),見(jiàn)表1。
表1 溫陽(yáng)活血湯的主要有效成分
leonticineSexangularetin13-methylpalmatrubinebeta-sitosterolN-methyllaurotetaninesitosterolnorglaucing黨參(21個(gè))poriferasta-7,22E-dien-3beta-olPontevedrinePerlolyrinePseudocoptisineDiop24240-05-9ZINC03978781saulatineStigmasterolstylopine7-Methoxy-2-methylisoflavoneTetrahydrocorysamineSpinasteroltetrahydroprotopapaverineChrysanthemaxanthinST057701FrutinoneA2,3,9,10-tetramethoxy-13-methyl-5,6-dihydroisoquinolino[2,1-b]isoquinolin-8-oneluteolinStigmasterolTaraxerolPalmatinestigmast-7-enolFumarine3-beta-HydroxymethyllenetanshiquinoneIsocorypalminemethylicosa-11,14-dienoatebicuculline5alpha-Stigmastan-3,6-dioneC093677-(beta-Xylosyl)cephalomannine_qtquercetinDaturilin山慈菇(3個(gè))beta-sitosterolglyciteinStigmasterolSpinosideA2-methoxy-9,10-dihydrophenanthrene-4,5-diol(8S,9S,10R,13R,14S,17R)-17-[(E,2R,5S)-5-ethyl-6-methylhept-3-en-2-yl]-10,13-dimethyl-1,2,4,7,8,9,11,12,14,15,16,17-dodecahydrocyclo-penta[a]phenanthren-3-one補(bǔ)骨脂(3個(gè))Daucosterol11-Hydroxyrankinidineisobavachin炙黃芪(1個(gè))beta-sitosterolStigmasterol
2.2溫陽(yáng)活血湯潛在治療靶點(diǎn):在數(shù)據(jù)庫(kù)TCMSP和Swiss Target Prediction中檢索溫陽(yáng)活血湯109個(gè)化合物的相關(guān)靶點(diǎn),去重后保留高匹配度靶點(diǎn),最終得到相關(guān)靶點(diǎn)374個(gè)。在數(shù)據(jù)庫(kù)STRING11.0中獲得374個(gè)相關(guān)靶點(diǎn)PPI網(wǎng)絡(luò),并在Cytoscape3.6.0中對(duì)此網(wǎng)絡(luò)作可視化分析,發(fā)現(xiàn)此網(wǎng)絡(luò)中包含節(jié)點(diǎn)369個(gè)和邊數(shù)7623條。在MCODE插件中進(jìn)行基因篩選,得到158個(gè)核心基因,位于A(yíng)、B、C三個(gè)子簇,見(jiàn)圖1。三個(gè)子簇包含的節(jié)點(diǎn)數(shù)分別為66、59和33;組成邊數(shù)分別為1686、426和182。
圖1 溫陽(yáng)活血湯潛在治療靶點(diǎn)核心基因
2.3溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的主要潛在靶點(diǎn):以“Bone metastasis”為檢索詞,在疾病基因數(shù)據(jù)庫(kù)DisGeNet和Drugbank中最終得到773個(gè)靶基因,與藥物374個(gè)相關(guān)靶點(diǎn)進(jìn)行映射,見(jiàn)圖2;最后獲得105個(gè)溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的靶點(diǎn),見(jiàn)圖3A。節(jié)點(diǎn)連接度Degree數(shù)值越大,說(shuō)明有越多條邊與該節(jié)點(diǎn)相連,即該靶點(diǎn)越重要。據(jù)此,篩選出Degree值位于前25的關(guān)鍵靶點(diǎn),見(jiàn)圖3。
圖2 溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的潛在治療靶點(diǎn)和骨轉(zhuǎn)移癌疾病靶點(diǎn)的韋恩圖
圖3 溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的105個(gè)靶基因(A)及核心基因(B)
2.4溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌靶基因的GO功能分析:在數(shù)據(jù)庫(kù)DAVID中對(duì)105個(gè)潛在治療靶點(diǎn)進(jìn)行GO功能富集分析,共得到143條結(jié)果,其中BP、CC和MF分別涉及88項(xiàng)、11項(xiàng)和44項(xiàng)。以錯(cuò)誤發(fā)現(xiàn)率FDR<0.05為篩選條件,獲得100個(gè)條目,其中BP63項(xiàng)、CC5項(xiàng)、MF32項(xiàng),將BP、MF的前10項(xiàng)和CC的前5項(xiàng)繪制成氣泡圖來(lái)進(jìn)行富集分析,見(jiàn)圖4。在BP方面,包括脂多糖介導(dǎo)的信號(hào)通路、5-羥色胺受體信號(hào)通路和細(xì)胞對(duì)血管內(nèi)皮生長(zhǎng)因子刺激的反應(yīng)等;在CC中,涉及樹(shù)突、內(nèi)溶酶體腔、血小板α顆粒腔、受體復(fù)合物、蛋白質(zhì)細(xì)胞外基質(zhì)等;在MF方面,涉及G蛋白偶聯(lián)的5-羥色胺受體活性、絲氨酸型內(nèi)肽酶活性、蛋白酪氨酸激酶活性和蛋白激酶活性等。
圖4 溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌靶基因的GO功能富集分析
2.5溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌靶基因KEGG通路分析:通過(guò)對(duì)105個(gè)潛在治療靶點(diǎn)做KEGG通路分析發(fā)現(xiàn),與溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌作用機(jī)理相關(guān)的通路有92條,以FDR<0.05為富集篩選條件,得到84條符合該條件的通路,取排名前25的通路在氣泡圖中進(jìn)行富集功能分析,見(jiàn)圖5。相關(guān)通路主要有HIF-1信號(hào)通路、TNF信號(hào)通路、VEGF信號(hào)通路、ErbB信號(hào)通路、Toll樣受體信號(hào)通路、Rap1信號(hào)通路等。
圖5 溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌靶基因的KEGG通路富集分析
現(xiàn)代醫(yī)學(xué)認(rèn)為骨轉(zhuǎn)移癌發(fā)生機(jī)制與骨微環(huán)境和腫瘤細(xì)胞的相互作用有關(guān),在轉(zhuǎn)移過(guò)程中,成骨細(xì)胞、破骨細(xì)胞兩者之間動(dòng)態(tài)失衡所產(chǎn)生的環(huán)境,更加有利于腫瘤細(xì)胞的轉(zhuǎn)移[5]。中醫(yī)學(xué)認(rèn)為罹患腫瘤的患者多處于邪正并存的狀態(tài),故在治則上應(yīng)當(dāng)扶正與驅(qū)邪并重,以改善機(jī)體狀態(tài)[6]。自擬溫陽(yáng)活血湯中以干姜9g、肉桂10g和補(bǔ)骨脂20g溫陽(yáng);以紅花10g和延胡索20g活血;以黨參30g、炙黃芪30g和山慈菇20g扶正以驅(qū)邪,全方八位中藥共奏溫陽(yáng)活血之功,但其作用機(jī)制尚未明確。
通過(guò)篩選溫陽(yáng)活血湯的主要有效成分并構(gòu)建多層次相互作用的成分與靶點(diǎn)相互作用網(wǎng)絡(luò)示意圖,發(fā)現(xiàn)槲皮素、豆甾醇、木犀草素、β-谷甾醇、山奈酚、黃芩素和β-胡蘿卜素等為其核心化合物,且大多數(shù)能抑癌、抗炎和調(diào)節(jié)免疫系統(tǒng)。其中槲皮素作為關(guān)鍵化合物,一方面可以通過(guò)抑制RANK-L、RANK等蛋白以及IL-8,TNF-α等炎癥介質(zhì)的表達(dá),增加OPG的水平,進(jìn)一步降低PAR2/TRPV1通路中的相關(guān)分子(如PAR2,TRPV1等)和TNF-α和IL-1β的水平,最終使破骨細(xì)胞活化受抑、骨破壞程度得以下降[7]。豆甾醇和β-谷甾醇均可抑制脂多糖誘導(dǎo)的巨噬細(xì)胞產(chǎn)生和分泌多種炎癥介質(zhì),如活性氧,一氧化氮和TNF-α,主要與NF-κB信號(hào)通路相關(guān):其中TNF-α可通過(guò)對(duì)NF-κB激酶IKKβ的抑制功能來(lái)增加NF-κB,而NF-κB3抑制劑的泛素化可減少NF-κB[8]。木犀草素可通過(guò)抑制PCNA蛋白和MMP-9蛋白的表達(dá)及PI3K/AKT通路活化來(lái)阻礙宮頸癌HeLa細(xì)胞的增殖和遷移[9]。
通過(guò)Cytoscape3.6.0軟件構(gòu)建溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的蛋白互作網(wǎng)絡(luò)示意圖,發(fā)現(xiàn)VEGFA、AKT1、IL6、TP53、ALB、CXCL8、CASP3、TNF、EGFR、EGF和MAPK3等為其關(guān)鍵靶點(diǎn)。VEGFA具有調(diào)節(jié)血管生成和促進(jìn)癌細(xì)胞增殖的作用,其通路受到HIF和JAK/STAT通路的調(diào)節(jié),且VEGFA信號(hào)級(jí)聯(lián)與MAPK通路有關(guān)[10]。促骨因子IL-6和TNF-α都可以刺激成骨細(xì)胞產(chǎn)生RANKL,可與RANK即破骨細(xì)胞前體表面的受體結(jié)合,刺激下游信號(hào)分子,包括MAPK和PI3K/AKT,并可對(duì)破骨細(xì)胞前體向功能性破骨細(xì)胞的轉(zhuǎn)化起促進(jìn)作用。CXCL8,也稱(chēng)為IL-8,通過(guò)與CXCR-1和CXCR2作用起調(diào)節(jié)炎癥反應(yīng)的作用。此外,有研究發(fā)現(xiàn)在肝癌組織中CXCL8上調(diào),可增加ERK、p-ERK和凋亡抑制基因Survivin的表達(dá),降低CASP3和BAX的表達(dá),從而使癌細(xì)胞增殖和遷移能力增加[11]。
根據(jù)GO和KEGG富集分析結(jié)果可以發(fā)現(xiàn),溫陽(yáng)活血湯治療骨轉(zhuǎn)移癌的作用機(jī)制主要包含抑癌,調(diào)節(jié)骨微環(huán)境和腫瘤細(xì)胞等方面。當(dāng)骨微環(huán)境缺氧時(shí),HIF-1可啟動(dòng)缺氧基因的轉(zhuǎn)錄,這些基因?qū)ρ苌?、腫瘤細(xì)胞凋亡和生長(zhǎng)因子/細(xì)胞因子活性至關(guān)重要。此外,HIF-1還可與生長(zhǎng)因子/細(xì)胞因子(如VEGF,EGF等)相互作用,這些因子在癌細(xì)胞上的表達(dá)可進(jìn)一步促進(jìn)癌細(xì)胞的增殖和轉(zhuǎn)移[12]。有研究表明,在缺血性壞死骨組織中TNF-α表達(dá)顯著上調(diào),可激活p38MAPK/NF-κB信號(hào)通路。而抑制p38mapk/NF-κB通路可減弱TNF-α誘導(dǎo)的自噬作用,并增強(qiáng)TNF-α誘導(dǎo)的成骨細(xì)胞凋亡[13]。此外,也有研究發(fā)現(xiàn),長(zhǎng)非編碼RNA的SPRY4內(nèi)含子轉(zhuǎn)錄1可通過(guò)TNF信號(hào)通路促進(jìn)肝細(xì)胞癌的進(jìn)展和轉(zhuǎn)移[14]。Toll樣受體(TLR)可調(diào)節(jié)機(jī)體免疫功能和細(xì)胞代謝,其中TLR4是TLR家族中最具代表性的成員,包含兩條信號(hào)通路:一是MyD88依賴(lài)型通路,可通過(guò)催化TRAF6泛素化啟動(dòng)NF-κB和MAPK表達(dá),促進(jìn)炎性細(xì)胞因子釋放;二是MyD88非依賴(lài)型通路,該通路可通過(guò)活化NF-κB以促進(jìn)炎性因子釋放,也可激活干擾素誘導(dǎo)基因IP-10與調(diào)節(jié)活化正常T細(xì)胞表達(dá)與分泌的趨化因子[15]。
綜上所述,基于網(wǎng)絡(luò)藥理學(xué)研究方法,可發(fā)現(xiàn)溫陽(yáng)活血湯有效化合物為槲皮素、豆甾醇、木犀草素等,可通過(guò)HIF-1、TNF、VEGF等信號(hào)通路調(diào)控VEGFA、IL6、TNF等靶點(diǎn),進(jìn)而發(fā)揮抑癌、抗骨轉(zhuǎn)移和調(diào)節(jié)免疫系統(tǒng)的作用,對(duì)臨床上治療骨轉(zhuǎn)移癌用藥具有一定參考意義,但仍需后期基礎(chǔ)實(shí)驗(yàn)加以驗(yàn)證。