活血化瘀中藥抗同種異體移植排斥潛在活性成分與機(jī)制的網(wǎng)絡(luò)藥理學(xué)研究

曾玉群 邵麗娜 何強(qiáng) 戴振華

器官移植手術(shù)后的排斥反應(yīng)是困擾器官移植發(fā)展的主要難題之一。免疫抑制劑的使用是器官移植界的重大突破，然而長期服用免疫抑制劑會導(dǎo)致多種毒副反應(yīng)，例如導(dǎo)致感染、誘發(fā)腫瘤等[1]。因此，探索移植排斥的免疫機(jī)制，尋找新的藥物作用靶點(diǎn)是抑制移植排斥和誘導(dǎo)移植物長期存活或免疫耐受的重要策略。

近幾十年我國器官移植飛速發(fā)展，針對抗移植排斥的藥物研究，祖國醫(yī)學(xué)在動物與臨床試驗(yàn)中取得了較大進(jìn)展。多項(xiàng)研究顯示中醫(yī)藥在抗移植免疫排斥以及減少免疫抑制劑不良反應(yīng)方面具有良好的療效[2]。目前中醫(yī)理論認(rèn)為器官移植后排斥反應(yīng)會使機(jī)體處于氣虛、血瘀、痰濕互結(jié)狀態(tài)，相關(guān)中藥被證明能抑制動物器官移植排斥反應(yīng)，尤其是多個活血化瘀中藥的抗排斥作用[3]。但中藥成分復(fù)雜，靶點(diǎn)眾多，機(jī)制不明，影響了進(jìn)一步開發(fā)及其臨床使用。

網(wǎng)絡(luò)藥理學(xué)是2007年英國藥理學(xué)家Hopkins提出的概念，其基于系統(tǒng)生物學(xué)將藥物分子、作用靶點(diǎn)及疾病3者的網(wǎng)絡(luò)關(guān)系相連，從而研究藥物治療疾病的作用機(jī)制[4]。網(wǎng)絡(luò)藥理學(xué)為我國中醫(yī)藥的研發(fā)提供了新思路和新策略。本研究將利用網(wǎng)絡(luò)藥理學(xué)去挖掘活血化瘀中藥作用于同種異體移植排斥的活性成分及其機(jī)制，為中醫(yī)藥抗器官移植排斥提供理論依據(jù)。

1 資料和方法

1.1 數(shù)據(jù)來源 活血化瘀中藥來源于《中藥學(xué)》課本第十二章活血化瘀藥[5]，共包括川芎、延胡索、郁金、姜黃、乳香、沒藥、丹參、紅花、桃仁、益母草、澤蘭、牛膝、雞血藤、王不留行、月季花、凌霄花、鳳仙花、土鱉蟲、自然銅、蘇木、骨碎補(bǔ)、血竭、兒茶、劉寄奴、莪術(shù)、三棱、水蛭、虻蟲、斑蝥、穿山甲32味中藥。中藥活性成分來自BATMAN-TCM在線數(shù)據(jù)庫[6]，同種異體移植排斥靶基因來自 PALM-IST[7]及 GeneCards數(shù)據(jù)庫[8]。

1.2 中藥活性化合物及靶標(biāo)的篩選 采用BATMANTCM在線數(shù)據(jù)庫[6]（http://bionet.ncpsb.org/batman-tcm/）檢索獲取中藥活性成分及相關(guān)靶基因。BATMAN-TCM在線數(shù)據(jù)庫是第一個專門為中藥分子機(jī)制設(shè)計的在線生物信息學(xué)分析工具，目前被中醫(yī)藥領(lǐng)域廣泛使用，其基于藥物與藥物、蛋白與蛋白相似性算法來預(yù)測每一種中藥及方劑的潛在靶點(diǎn)。檢索參數(shù)選擇中藥各化合物成分得分（Score Cut off）≥48 分，P≤0.05。

1.3 關(guān)鍵活性成分靶基因集的構(gòu)建 利用PALM-IST[7]（http://www.hpppi.iicb.res.in/ctm）及GeneCards[8]（https://www.genecards.org/）在線數(shù)據(jù)庫，以“同種異體移植排斥（allograft rejection）”作為檢索關(guān)鍵詞，獲取相關(guān)靶基因作為疾病病理靶標(biāo)。同時與中藥活性成分所預(yù)測的靶基因進(jìn)行匹配，采用Cytoscape 3.3.0軟件[9]構(gòu)建活血化瘀中藥關(guān)鍵活性化合物的靶基因集。最后使用STRING數(shù)據(jù)庫[10]（http://string-db.org/），選擇“homo sapiens（人類）”物種進(jìn)行靶基因間相互作用分析。PALM-IST數(shù)據(jù)庫是基于Pubmed的文獻(xiàn)資源進(jìn)行基因/蛋白名稱、藥物、疾病等生物實(shí)體及其相互關(guān)系提取的檢索工具，其可根據(jù)文獻(xiàn)的摘要進(jìn)行數(shù)據(jù)檢索挖掘。GeneCards數(shù)據(jù)庫是目前使用最多的人類基因靶點(diǎn)集成數(shù)據(jù)庫，其自動集成了150多個網(wǎng)絡(luò)來源的基因數(shù)據(jù)，包括了基因組學(xué)、蛋白質(zhì)組學(xué)、轉(zhuǎn)錄組學(xué)以及疾病等信息，其界面友好，方便用戶使用。STRING數(shù)據(jù)庫是目前使用最廣泛的進(jìn)行已知和預(yù)測蛋白質(zhì)-蛋白質(zhì)相互之間作用關(guān)系的數(shù)據(jù)庫。

1.4 靶基因基因本體（gene ontology，GO）和基因組百科全書（kyoto encyclopedia of genes and genomes，KEGG）富集分析利用DAVID[11]（https://david.ncifcrf.gov/）及STRING數(shù)據(jù)庫將獲取的所有關(guān)鍵活性成分靶標(biāo)進(jìn)行生物信息學(xué)在線分析，包括GO分析和通路分析：對靶標(biāo)基因進(jìn)行GO注釋，并計算每條GO項(xiàng)目的P值，以P＜0.05為標(biāo)準(zhǔn)篩選出顯著性富集的GO項(xiàng)目，獲得GO相關(guān)生物學(xué)過程、分子功能、細(xì)胞成分等；對靶標(biāo)進(jìn)行KEGG通路注釋，得到其所參與的通路及每條通路的P值，以P＜0.05為標(biāo)準(zhǔn)篩選出顯著性通路，最終得到與同種異體移植排斥的相關(guān)通路。DAVID數(shù)據(jù)庫是目前使用最廣泛的可對大量基因背后的生物學(xué)含義進(jìn)行全面功能注釋的在線工具。

1.5 中藥-活性成分-靶標(biāo)-移植排斥的網(wǎng)絡(luò)構(gòu)建及分析 對活血化瘀中藥活性成分與同種異體移植排斥靶標(biāo)之間的復(fù)雜關(guān)系進(jìn)行網(wǎng)絡(luò)分析。利用蛋白質(zhì)-蛋白質(zhì)相互作用信息，將得到的活血化瘀中藥關(guān)鍵化學(xué)活性成分、關(guān)鍵活性成分的靶標(biāo)、重要的同種異體移植排斥靶標(biāo)及相互作用的蛋白連接在一起，以構(gòu)建中藥-活性成分-靶標(biāo)-移植排斥網(wǎng)絡(luò)。使用Cytoscape 3.3.0對網(wǎng)絡(luò)進(jìn)行可視化和分析，并使用Cytoscape軟件的NetworkAnalyzer APP計算每個節(jié)點(diǎn)的拓?fù)涮卣鳌⒕哂休^高“連接度（degree）”、“介數(shù)（betweenness centrality）”的網(wǎng)絡(luò)節(jié)點(diǎn)識別為同種異體移植排斥的候選活血化瘀中藥的核心靶標(biāo)[9]。

2 結(jié)果

2.1 活血化瘀中藥活性成分的篩選 32味活血化瘀中藥經(jīng)BATMAN-TCM檢索后總共獲得922個相關(guān)活性化合物，得到503個獨(dú)立的具有化學(xué)結(jié)構(gòu)的化學(xué)活性成分，其中最主要的9種活性成分具體見表1，可知活血化瘀類中藥具有一些共同的化學(xué)活性成分，且姜黃、郁金、川芎3種中藥是活血化瘀中藥含有主要活性成分的代表中藥。

表1 活血化瘀中藥主要含有的活性成分

2.2 活血化瘀中藥關(guān)鍵活性成分靶標(biāo)的篩選 經(jīng)BATMAN-TCM檢索后共獲取相關(guān)化學(xué)活性成分預(yù)測的相關(guān)靶標(biāo)5 224個，獨(dú)立靶標(biāo)共554個。經(jīng)PALM-IST工具檢索到與同種異體移植排斥相關(guān)的靶基因共2 397個（共10 960篇摘要），經(jīng)GeneCards檢索到相關(guān)的靶基因共1 679個，將得到的靶標(biāo)去掉重復(fù)基因最后得到同種異體移植排斥靶基因共704個。經(jīng)Cytoscape 3.3.0軟件將704個同種異體移植排斥靶基因與554個活血化瘀中藥活性成分靶基因匹配后共篩選出67個關(guān)鍵靶基因，選擇較高的連接度中位數(shù)的3倍（連接度≥12）和介數(shù)獲得關(guān)鍵17個靶基因，見表2。進(jìn)一步使用STRING工具獲取17個關(guān)鍵靶標(biāo)之間的相互作用關(guān)系（選擇最低互相關(guān)系得分0.4），共獲取44個互相關(guān)系，見圖1。

2.3 生物GO富集及通路分析 經(jīng)DAVID工具對17個預(yù)測關(guān)鍵靶基因進(jìn)行GO富集后，共獲得52條生物學(xué)含義富集結(jié)果，其中包括生物學(xué)過程41項(xiàng)，分子功能8項(xiàng)，細(xì)胞成分3項(xiàng)，主要的GO生物學(xué)過程富集結(jié)果見表3，可見多條相關(guān)生物學(xué)過程均與免疫反應(yīng)、炎癥反應(yīng)、細(xì)胞凋亡等過程相關(guān)。進(jìn)一步進(jìn)行KEGG的通路分析，得到主要的16條通路見圖2，可見核因子κB（NF-κB）信號通路、IL-17信號通路、TNF信號通路和過氧化物酶體增殖物激活受體（PPAR）信號通路主要參與活血化瘀中藥活性成分作用于同種異體移植排斥的相關(guān)藥理機(jī)制。

2.4 中藥-活性成分-靶標(biāo)-移植排斥網(wǎng)絡(luò)分析 根據(jù)篩選出來的活血化瘀中藥關(guān)鍵26種活性成分的17個核心靶標(biāo)及其網(wǎng)絡(luò)拓?fù)鋮?shù)，利用Cytoscape 3.3.0軟件構(gòu)建中藥-活性成分-靶標(biāo)-同種異體移植排斥網(wǎng)絡(luò)圖（圖3）。進(jìn)一步利用藥理網(wǎng)絡(luò)分析關(guān)鍵活血化瘀活性成分作用于同種異體移植排斥的相關(guān)藥理機(jī)制可能是活血化瘀中藥的活性成分主要是萜類化合物包括如莪術(shù)二酮（Curdione）、香芹酮（Carvone）、丹參新酮（Miltirone）、百里酚（Thymol）、姜黃酮（Turmerone）等通過 NF-κB 信號通路、IL-17信號通路、TNF信號通路以及PPAR信號通路作用于 PTGS2、LTA、ATM、TNF、APOA2、PPARG、IL-13、IL-1B靶標(biāo)從而改善同種異體移植排斥反應(yīng)，見圖4。

表2 關(guān)鍵活性成分17個核心靶基因

圖1 關(guān)鍵17個靶基因間相互作用關(guān)系

表3 關(guān)鍵靶基因主要生物學(xué)過程

圖2 關(guān)鍵靶標(biāo)的潛在信號通路

圖3 核心中藥-活性成分-靶標(biāo)-移植排斥網(wǎng)絡(luò)藥理圖（關(guān)鍵17個靶標(biāo)）

圖4 活血化瘀中藥活性成分作用于同種異體移植排斥的藥理機(jī)制

3 討論

眾所周知，免疫排斥反應(yīng)是限制同種異體器官移植的主要原因，細(xì)胞免疫是免疫排斥反應(yīng)的主要啟動機(jī)制。細(xì)胞免疫中T淋巴細(xì)胞在移植免疫排斥中起著關(guān)鍵作用。T細(xì)胞中的輔助性T（T help，Th）細(xì)胞包含不同的亞群[12-13]，主導(dǎo)細(xì)胞免疫，主要分泌IL、干擾素（Interferon，IFN）等細(xì)胞因子，其中Th1主要分泌IFN-γ等，Th2主要分泌 IL-4、IL-5等，Th17主要分泌 IL-17A、IL-22以及IL-13等[14]。Th1細(xì)胞是參與移植免疫應(yīng)答的主要免疫細(xì)胞，而Th2細(xì)胞對移植排斥反應(yīng)有一定的調(diào)節(jié)作用，Th17細(xì)胞可能在一定程度上也參與了移植排斥反應(yīng)，IL-17在移植排斥中大量分泌也促使移植排斥反應(yīng)[15-16]。因此靶向IL-17可抑制同種異體移植排斥反應(yīng)。而由抗體介導(dǎo)的體液免疫也參與移植排斥反應(yīng)，其中包括慢性排斥系統(tǒng)炎癥反應(yīng)，例如NF-κB通路也參與其中等[17-18]。另外，目前的相關(guān)研究顯示TNF受體及相關(guān)分子可控制淋巴器官的發(fā)育，影響T細(xì)胞和抗原遞呈細(xì)胞的活化和存活，并改變細(xì)胞因子和趨化因子的產(chǎn)生[19]。且越來越多的數(shù)據(jù)表明腫瘤壞死因子受體（TNFR）/TNF超家族成員對移植免疫中發(fā)揮作用的CD8+T細(xì)胞的功能尤為重要，可見其在移植排斥中具有重要作用[20]。因此靶向TNFR/TNF超家族分子的策略似乎很有可能用于控制或預(yù)防移植器官和組織的排斥反應(yīng)。可見本研究挖掘的IL-17、TNF以及NF-κB信號通路及相關(guān)靶標(biāo)如IL-13、IL-1B、TNF等與移植排斥機(jī)制緊密相關(guān)，作用于其信號通路及靶標(biāo)將成為抗同種異體移植排斥藥物開發(fā)的重要策略。

我國臨床器官移植已開展五十余年，中醫(yī)藥及其活性化合物被用于抗器官移植排斥及減少免疫抑制劑的毒副反應(yīng)，且取得一定的療效。常用的中藥如雷公藤中分離出來的雷公藤多甙包括萜類化合物雷公藤甲素以及冬蟲夏草中提取的FTY-720均證明可有效抑制臨床患者移植排斥反應(yīng)，可協(xié)同免疫抑制劑發(fā)揮抑制器官移植排斥且同時降低免疫抑制劑的毒副反應(yīng)，已成為新的免疫抑制劑[21-23]。另外多個主要活血化瘀中藥活性成分也都在抗器官移植排斥中發(fā)揮作用，如蘇木提取物蘇木素A可延長大鼠心臟移植物的存活，從而減弱大鼠心臟同種異體移植的急性排斥反應(yīng)[24]；川芎嗪可聯(lián)合小劑量環(huán)孢素（cyclosporine，CsA）抑制小腸移植急性排斥反應(yīng)更佳[25]；姜黃素具有免疫抑制作用，在BN-WKY模型大鼠中單用姜黃素14d可將其平均生存時間由9.1d延長至20.5～24.5d，聯(lián)合姜黃素和CsA可延長至28.5～35.6d，比單用CsA效果好[26]；筆者也證實(shí)了相關(guān)活血化瘀中藥活性成分山奈酚（Kaempferol）可通過提高叉頭狀螺旋轉(zhuǎn)錄因子3表達(dá)而增加調(diào)節(jié)性T細(xì)胞（Treg）數(shù)量以協(xié)助小劑量CsA更好發(fā)揮免疫抑制作用，從而最大限度降低CsA的不良反應(yīng)[27]；Lu等[28]也發(fā)現(xiàn)丹參素（Tanshinol）可顯著延遲小鼠同種異體移植心臟排斥反應(yīng)，并與雷帕霉素（Rapamycin）發(fā)揮協(xié)同作用延長小鼠同種異體心臟移植物的存活時間。因此，本研究挖掘出了作用于同種異體排斥反應(yīng)的主要活血化瘀中藥包括川芎、丹參、姜黃、莪術(shù)、郁金等，且預(yù)測了其可能發(fā)揮作用的主要活性化合物，這將有助于進(jìn)一步挖掘出對器官移植排斥有效的免疫抑制新藥，以及可改善免疫抑制劑毒副反應(yīng)的輔助新藥。

為了揭示對疾病和涉及藥物的活性物質(zhì)的相關(guān)作用，本研究進(jìn)行了網(wǎng)絡(luò)藥理學(xué)分析，以研究活血化瘀中藥作用于同種異體移植排斥的相關(guān)藥理機(jī)制。通過網(wǎng)絡(luò)藥理學(xué)方法，筆者從分子或系統(tǒng)水平構(gòu)建活性化合物與蛋白質(zhì)/基因靶標(biāo)之間相互作用的生物網(wǎng)絡(luò)，預(yù)測了與同種異體移植排斥相關(guān)的活血化瘀中藥活性成分和潛在靶標(biāo)。預(yù)測結(jié)果表明活血化瘀中藥26種活性成分與17個與同種異體移植排斥相關(guān)的靶標(biāo)相互作用。結(jié)果提示一些活性成分可以同時調(diào)節(jié)相同的靶標(biāo)，這很好地解釋了中藥的協(xié)同機(jī)制。值得一提的是，計算模型只能給出提示與預(yù)測，仍然需要通過實(shí)驗(yàn)進(jìn)行驗(yàn)證。對中藥藥理機(jī)制的理解是其現(xiàn)代化的必要條件，網(wǎng)絡(luò)藥理學(xué)的綜合和整體思想與中藥整體觀的本質(zhì)是一致的。因此，網(wǎng)絡(luò)藥理學(xué)可以被認(rèn)為是預(yù)測和發(fā)現(xiàn)治療特定疾病的新藥研發(fā)的有效方法，這將為揭示中醫(yī)藥的作用機(jī)制提供新的策略。