陳奇劍 李春燕 肖 先 劉曉龍 薛金濤▲

1.新鄉(xiāng)醫(yī)學(xué)院藥學(xué)院，河南新鄉(xiāng) 453002；2.新鄉(xiāng)醫(yī)學(xué)院三全學(xué)院生物與基礎(chǔ)醫(yī)學(xué)實(shí)驗(yàn)教學(xué)中心，河南新鄉(xiāng) 453003

作為一個(gè)重大公共衛(wèi)生問題，抑郁癥具有高自殺率、高自殘率和高復(fù)發(fā)率的特點(diǎn)。經(jīng)典的抗抑郁藥物起效慢、療效較差，因此尋找快速、有效和安全的藥物用于抑郁癥的治療為當(dāng)前研究的熱點(diǎn)[1-6]。翅果菊（pterocypsela elata）為菊科植物翅果菊Pterocypsela elata（Hemsl.）Shih的根或全草。本研究從翅果菊中共分離得到18個(gè)活性成分，其中萵苣苷B（lactuside B，LB）是一種新型化合物，前期藥理學(xué)研究表明LB和翅果菊具有良好的抗抑郁、改善記憶、抗腦缺血和抗心肌缺血作用[7-10]。鑒于此，翅果菊和LB在抑郁癥的治療中具有廣闊的應(yīng)用前景，但其作用機(jī)制尚不清楚，本研究擬采用網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接方法研究翅果菊發(fā)揮抗抑郁作用的主要活性成分和作用機(jī)制。

1 資料與方法

1.1 軟件

ChemBioOffice軟件（ChemBioDraw和ChemBio3D，Ultra 14.0），Open Babel GUI軟件（Version 2.4.1），Autodock Vina軟件（Version 1.1.2），PyMoL軟件（Version 1.7.2.1）和Autodock Tools軟件（Version 1.5.6）。

1.2 數(shù)據(jù)庫(kù)

DrugBank數(shù) 據(jù) 庫(kù)（Version 5.1.3，https://www.drugbank.ca/），Batman數(shù) 據(jù) 庫(kù) （http://bionet.ncpsb.org/batman-tcm/index.php/Home/Index/index），String數(shù)據(jù)庫(kù) （Version 11.0，http://www.string-db.org/）和The Protein Data Bank（PDB）數(shù) 據(jù) 庫(kù)（http://www.rcsb.org/pdb/home/home.do）。

1.3 網(wǎng)絡(luò)藥理學(xué)

整理課題組從翅果菊中分離得到的化學(xué)分子[7-10]，在ChemBioOffice軟件中分別畫出各化合物化學(xué)結(jié)構(gòu)式，以CDX格式保存。利用Open Babel GUI軟件將CDX格式文件轉(zhuǎn)換為Inchi格式，并上傳至Batman數(shù)據(jù)庫(kù)獲得各成分的作用靶點(diǎn)，選擇與化學(xué)成分?jǐn)M合度高（得分值Z’-score ≥ 23）的靶點(diǎn)進(jìn)行下一步的分析。在DrugBank數(shù)據(jù)庫(kù)中以“Depression”（抑郁癥）為關(guān)鍵詞進(jìn)行檢索，篩選出作用于抑郁癥的活性靶點(diǎn)。將Batman數(shù)據(jù)庫(kù)中獲得的化學(xué)成分作用靶點(diǎn)與DrugBank數(shù)據(jù)庫(kù)中的抗抑郁癥的活性靶點(diǎn)進(jìn)行交叉對(duì)比，篩選翅果菊抗抑郁癥的活性成分和作用靶點(diǎn)。在String數(shù)據(jù)庫(kù)，將物種限定為“Homo sapiens”，對(duì)篩選得到的靶點(diǎn)進(jìn)行基因功能注釋（gene ontology，GO）富集分析和京都基因與基因組百科全書（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路分析，以對(duì)翅果菊的功能分布和通路注釋進(jìn)行研究分析。

1.4 分子對(duì)接

采用ChemBioOffice軟件的ChemBio3D功能將LB和guaianolide Ⅰ（Ⅰ）化學(xué)結(jié)構(gòu)轉(zhuǎn)化成三維格式，并在PDB數(shù)據(jù)庫(kù)檢索并下載目的作用靶點(diǎn)的三維蛋白結(jié)構(gòu)，利用Autodock Tools軟件對(duì)蛋白質(zhì)進(jìn)行常規(guī)處理（加氫、去水、加電荷等），并將蛋白和藥物活性成分轉(zhuǎn)化為pdqpt格式文件保存。采用Autodock Tools軟件進(jìn)行分子對(duì)接參數(shù)的設(shè)置，并進(jìn)一步使用Autodock Vina軟件進(jìn)行分子對(duì)接。最后使用PyMoL軟件處理對(duì)接結(jié)果。

2 結(jié)果

2.1 活性成分整理

本研究前期研究中從翅果菊中提取分離得到18個(gè) 化 學(xué) 成 分[7-10]，分 別 是LB，11β，13-二氫萵苣內(nèi)酯乙酸酯，11β,13-二氫萵苣內(nèi)酯，GI，guaianolide Ⅱ，guaianolide Ⅲ，（24R）-5α-豆甾-7，22（E）-二烯-3α-醇，齊墩果酸，β-香樹脂），β-谷甾醇，胡蘿卜苷，3，3′，4-三甲氧基鞣花酸，硬脂酸，3-oxo-friedlane，stigmasterol，aralia cerebroside，ursa-12-ene-11-one-3-oloctocosate和hexacosanoic acid，其中LB從翅果菊中分離得到一種新型化合物。

2.2 抗抑郁活性成分和作用靶點(diǎn)的篩選

通過統(tǒng)計(jì)分析Batman數(shù)據(jù)庫(kù)篩選得到的活性靶點(diǎn)，翅果菊18個(gè)活性成分共得到409個(gè)作用靶點(diǎn)。在DrugBank數(shù)據(jù)庫(kù)中共檢索得到292個(gè)抗抑郁作用靶點(diǎn)。交叉對(duì)比化學(xué)成分的作用靶點(diǎn)與抗抑郁作用靶點(diǎn)，篩選得到1作用靶點(diǎn)γ-氨基丁酸受體亞單位β-3（γ-aminobutyric acid receptor subunit beta-3，GABRB3，ID：P28472），翅果菊中共有2個(gè)藥效活性成分（LB和GI）作用于該靶點(diǎn)。

2.3 GO富集分析和KEGG通路分析

在上述對(duì)翅果菊活性成分和作用靶點(diǎn)篩選的基礎(chǔ)上，采用GO富集分析和KEGG代謝通路分析進(jìn)一步進(jìn)行網(wǎng)絡(luò)藥理學(xué)研究。GO富集分析主要包括生物過程、分子功能和細(xì)胞組分3個(gè)部分，并以P.Adust分析某一指標(biāo)和功能對(duì)富集的重要性，P.Adust值<0.05表示相應(yīng)指標(biāo)具有顯著性影響，P.Adust值越小，表示越具有顯著性影響，統(tǒng)計(jì)分析中一般取其對(duì)數(shù)的負(fù)值：-log10（P. Adust），故其值越大，表示越具有顯著性。

圖1A，GO富集分析共發(fā)現(xiàn)有顯著性影響的生物過程17個(gè)，其中顯著性意義較大的有γ-氨基丁酸信號(hào)通路（P.Adust=8.11E-26）、氯離子跨膜轉(zhuǎn)運(yùn)（P.Adust=8.05E-21）和G蛋白耦聯(lián)受體信號(hào)通路（P.Adust=1.05E-12）等。圖1B，共發(fā)現(xiàn)有顯著性影響的分子功能7個(gè)，其中顯著性意義較大的有GABA-A受體活性（P.Adust = 4.55E-27）、細(xì)胞外配體門控離子通道活性（P.Adust=2.67E-22）和氯離子通道活性（P.Adust=3.64E-22）等。圖1C，共發(fā)現(xiàn)有顯著性影響的細(xì)胞組分5個(gè)，其中顯著性意義較大的有氯離子通道復(fù)合體（P.Adust=1.08E-23）、GABA-A受 體復(fù) 合物（P.Adust =1.08E-23）和 突觸后膜（P.Adust=1.69E-17）等。圖2，KEGG通路分析結(jié)果表明，翅果菊活性成分作用于抑郁癥的KEGG通路中有顯著性影響的共有7個(gè)，其中顯著性意義較大的有尼古丁成癮（P.Adust=3.47E-25）、GABA能 突 觸（P.Adust=2.36E-22）和 嗎 啡 成 癮（P.Adust=2.36E-22）等。

圖1 作用靶點(diǎn)GO富集分析結(jié)果

圖2 作用靶點(diǎn)的KEGG通路分析結(jié)果

2.4 分子對(duì)接

分子對(duì)接方法廣泛應(yīng)用于藥物活性分子與靶蛋白的對(duì)接分析[11-16]。表1在PDB數(shù)據(jù)庫(kù)以GA B R B 3為關(guān)鍵詞檢索相關(guān)蛋白受體，共得到11個(gè)蛋白受體，將LB、GI分別與11個(gè)蛋白受體文件進(jìn)行分子對(duì)接。一般認(rèn)為分子與蛋白的結(jié)合自由能（Affinity）越小，兩者結(jié)合越緊密牢固，結(jié)合能<-5.0 kcal/mol時(shí)有較好的結(jié)合活性，結(jié)合能<-7.0 kcal/mol時(shí)提示有強(qiáng)烈的結(jié)合活性[11-12]。表1 LB與除7A5V外的其余10個(gè)蛋白受體均具有強(qiáng)烈的結(jié)合活性（Affinity<-7.0 kcal/mol），其 中 與6I53的 結(jié) 合活性最強(qiáng)（Affinity=-9.2 kcal/mol），其分子對(duì)接結(jié)果如圖3A，LB與6I53的結(jié)合位點(diǎn)在其蛋白內(nèi)部的空腔中（圖中數(shù)字為L(zhǎng)B與6I53中相近的氨基酸位點(diǎn)的分子間距離，單位：?）。G I與11個(gè)蛋白受體均具有強(qiáng)烈的結(jié)合活性（Affinity<-7.0 kcal/mol），其 中 與6I53和6HUG的結(jié)合活性最強(qiáng)，結(jié)合自由能均為9.5 kcal/mol，其對(duì)接結(jié)果見圖3B和圖3C，GI與6I53和6HUG的結(jié)合位點(diǎn)均在其蛋白外側(cè)的空腔中。

表1 活性成分和作用靶點(diǎn)對(duì)接結(jié)果

圖3 活性成分與靶點(diǎn)對(duì)接結(jié)果

3 討論

網(wǎng)絡(luò)藥理學(xué)結(jié)果表明，在翅果菊18個(gè)活性成分中篩選得到2個(gè)活性成分（LB和GI）及其作用靶點(diǎn)GABRB3，并采用分子對(duì)接技術(shù)對(duì)LB和GI與GABRB3進(jìn)行了對(duì)接驗(yàn)證，結(jié)果表明，LB和GI均與GABRB3的11個(gè)蛋白受體有良好的結(jié)合能力。GABRB3是γ-氨基丁酸A（γ-aminobutyric acid A，GABA）受體中功能顯著的一種，而GABA是中樞神經(jīng)系統(tǒng)中最重要的抑制性神經(jīng)遞質(zhì)之一，研究表明GABRB3與多種神經(jīng)類疾病有關(guān)，其異常表達(dá)可引起包括抑郁癥等多種精神疾病[17-18]。辛軍彩[17]發(fā)現(xiàn)GABRB3基因是抑郁癥的一種遺傳相關(guān)基因。張嬙[18]研究發(fā)現(xiàn)GABRB3等基因的異常表達(dá)可能是導(dǎo)致抑郁癥的中樞機(jī)制之一，糾正該基因的異常表達(dá)則對(duì)抑郁癥具有良好的干預(yù)作用。

GO富集和KEGG通路分析表明，翅果菊及其活性成分的作用機(jī)制與γ-氨基丁酸信號(hào)通路、GABA-A受體活性、GABA-A受體復(fù)合體和GABA能突觸等聯(lián)系密切，可見GABA的水平和功能對(duì)翅果菊及其活性成分發(fā)揮抗抑郁作用具有重要意義。劉向六等[19]對(duì)抑郁癥大鼠強(qiáng)迫游泳模型中研究發(fā)現(xiàn)，氯胺酮的抗抑郁作用與大鼠前額皮層 GABA下調(diào)有關(guān)。周月等[20]采用GABA和茶氨酸聯(lián)合干預(yù)抑郁癥大鼠，給藥后對(duì)抑郁癥有一定改善作用。

本研究從分子網(wǎng)絡(luò)層面揭示LB在抑郁癥的治療中具有良好的應(yīng)用前景，并分析了作用靶點(diǎn)和信號(hào)通路，為其作用機(jī)制的深入研究和藥理學(xué)研究應(yīng)用提供了研究依據(jù)和線索。