生立嵩，蘇酩，謝亞欣，2，張新軍，彭麗，許祚芝，劉善新*

(1.山東省中醫(yī)藥研究院，山東 濟(jì)南 250014；2.山東中醫(yī)藥大學(xué) 藥學(xué)院，山東 濟(jì)南 250355；3.華潤雙鶴利民藥業(yè)(濟(jì)南)有限公司 質(zhì)量部，山東 濟(jì)南 250200)

半夏厚樸湯始載于漢張仲景《金匱要略》“婦人咽中如有炙臠，半夏厚樸湯主之”，處方為半夏一升、厚樸三兩、茯苓四兩、生姜五兩、蘇葉二兩,有行氣散結(jié)、化痰降逆的功效，近代以來主要治療梅核氣[1]?，F(xiàn)代該方應(yīng)用廣泛，臨床研究用于治療抑郁癥[2]、反流性食管炎[3]、睡眠呼吸暫停綜合征[4]、咳嗽[5]等。前期本課題組結(jié)合臨床應(yīng)用，篩選半夏厚樸湯抗PDE4(磷酸二酯酶4)活性成分，并以液相色譜-質(zhì)譜聯(lián)用分析提取其中成分[6]，但對應(yīng)化合物數(shù)量較少，種類單一，這提示所選靶點(diǎn)及潛在藥效成分存在一定問題。鑒于此，我們以計算機(jī)虛擬篩選與高分辨質(zhì)譜檢測相結(jié)合的方法，應(yīng)用網(wǎng)絡(luò)藥理學(xué)搜集半夏厚樸湯藥材的成分及其作用靶點(diǎn)；采用TCMIO(中藥免疫腫瘤學(xué)數(shù)據(jù)庫)等中藥成分?jǐn)?shù)據(jù)庫和軟件進(jìn)行預(yù)測和分析，預(yù)測半夏厚樸湯主要作用靶點(diǎn)；為研究組方中小分子與蛋白的相互作用機(jī)制，對半夏厚樸湯中已報道化合物[7]進(jìn)行分子對接擬合和親和能打分，篩選親和能值較低的化合物。以高分辨質(zhì)譜法分析半夏厚樸湯提取物的成分及含量，篩選提取物中含量較高且在虛擬篩選中親和能值較低的作為半夏厚樸湯的潛在藥效成分組。

1 儀器與材料

1.1 實(shí)驗儀器

DNA503型千分之一電子天平(南京伯尼塔科學(xué)儀器有限公司)；KDM型控溫電熱套(鄄城華魯電熱儀器有限公司)；KQ-250DA型數(shù)控超聲波清洗機(jī)(昆山市超聲儀器有限公司)；DHG-9146A型電熱恒溫鼓風(fēng)干燥箱(上海精宏實(shí)驗設(shè)備有限公司)；Q-Exactive型超高效液相色譜-四級桿-靜電場軌道阱質(zhì)譜儀(美國賽默飛公司)。

1.2 試藥

茯苓(批號201001)、厚樸(批號200502)、半夏(批號191201)、生姜(批號200702)、蘇葉(批號200501)均購自山東建聯(lián)盛嘉中藥有限公司，由濟(jì)南億民藥業(yè)有限公司冉蓉高級工程師鑒定，均為2020年版《中國藥典》一部收載的正品。

1.3 研究工具

數(shù)據(jù)庫：中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫與分析平臺(TCMSP，https://tcmsp-e.com/)、中藥免疫腫瘤學(xué)數(shù)據(jù)庫(TCMIO，http://tcmio.xielab.net/)、藥物化學(xué)數(shù)據(jù)庫(ChemBL, https://www.ebi.ac.uk/)、RCSB PDB(https://www.rcsb.org/)、STRING(https://www.string-db.org/)；數(shù)據(jù)及化合物處理軟件：Discovery Studio 2016單機(jī)版、AutoDockTools 1.5.6、AutodockVina 1.1.2、Cytoscape 3.6.1；質(zhì)譜分析軟件：Compound Discoverer 3.0、Xcalibur 3.0.63。

2 方法與結(jié)果

2.1 成分?jǐn)?shù)據(jù)庫搜索和建立

從TCMSP、TCMIO等數(shù)據(jù)庫中搜集半夏厚樸湯五味藥材的報道成分，下載所有成分mol2格式文件，采用Autodock自帶的raccoon軟件處理結(jié)構(gòu)，生成pdbqt格式文件，作為對接成分分子數(shù)據(jù)庫。共得到五味藥材所報道化合物402個，主要包括黃酮類、三萜酸類、甾醇類、有機(jī)酸類等。

2.2 網(wǎng)絡(luò)藥理學(xué)研究和靶點(diǎn)收集

本研究采用兩種方法[8-9]搜集半夏厚樸湯的潛在靶點(diǎn)：(1)從藥效成分入手，選擇TCMSP數(shù)據(jù)庫中半夏厚樸湯五味藥材已報道化合物，選擇其中成藥性良好(OB≥30%且DL≥0.18)的化合物對應(yīng)靶點(diǎn)，將這些靶點(diǎn)與Genecards網(wǎng)站中，與咳嗽(關(guān)鍵詞cough)和神志疾病(關(guān)鍵詞mental disease)相關(guān)靶點(diǎn)進(jìn)行Venn交集，重疊靶標(biāo)在STRING數(shù)據(jù)庫中構(gòu)建PPI(protein-protein interaction)網(wǎng)絡(luò)，經(jīng)Cytoscape進(jìn)行拓?fù)浞治龊?，選擇拓?fù)渲递^高的作為半夏厚樸湯候選靶點(diǎn)；(2)從所有報道成分入手，選擇TCMSP數(shù)據(jù)庫中每味藥材所有報道化合物的靶點(diǎn)，將每個化合物所對應(yīng)的靶點(diǎn)進(jìn)行加和統(tǒng)計，選擇相關(guān)度最高的作為半夏厚樸湯潛在靶點(diǎn)。

在TCMSP數(shù)據(jù)庫中，共收集到成藥性良好分子44個，這44個分子對應(yīng)的靶點(diǎn)共有121個，采用Venn 圖分析半夏厚樸湯可能作用的靶點(diǎn)數(shù)量，結(jié)果如圖1所示，藥效成分的預(yù)測靶標(biāo)和咳病(cough)、神志疾病(mental disease)二者的共有靶標(biāo)21個。運(yùn)用STRING 網(wǎng)絡(luò)平臺進(jìn)行重疊靶標(biāo)的PPI網(wǎng)絡(luò)構(gòu)建，同時使用Cytoscape的network analyse功能對PPI網(wǎng)絡(luò)進(jìn)行分析，結(jié)果見表1，度值(degree)排名前2的靶標(biāo)分別是IL6和VEGFA?；诳傮w成分預(yù)測的靶點(diǎn)加和后結(jié)果如表2所示，排名最靠前的是GABR和PTGS兩個靶標(biāo)。綜合上述網(wǎng)絡(luò)藥理學(xué)研究，以這4個靶點(diǎn)作為半夏厚樸湯作用的潛在靶點(diǎn)。

圖1 Venn網(wǎng)絡(luò)圖

表1 基于藥效成分的預(yù)測靶標(biāo)

表2 半夏厚樸湯成分相關(guān)靶蛋白及關(guān)聯(lián)度

2.3 分子對接

對2.2節(jié)中篩選的潛在靶點(diǎn)，參考Autodock軟件[7]的說明進(jìn)行分子處理和對接操作：從RCSB(protein data bank description)網(wǎng)站中搜集相關(guān)靶點(diǎn)蛋白的結(jié)構(gòu)，以蛋白自帶的配體結(jié)構(gòu)和位點(diǎn)作為結(jié)合域，以Autodock軟件對蛋白質(zhì)進(jìn)行去水、加氫、加電荷操作。以Autodock vina 4.1軟件，利用蛋白質(zhì)晶體中自帶的抑制劑位置作為結(jié)合位置，設(shè)置結(jié)合位點(diǎn)和結(jié)合范圍大小，Docking算法采用local search以提高計算速度。利用2.1節(jié)搜集的藥材報道化合物分子進(jìn)行分子對接，將所有結(jié)果的親和能值(affinity，kJ/mol)進(jìn)行排序，選擇親和能值較低的前10%作為潛在有效化合物。

分子對接主要參數(shù)如表3所示，每個靶點(diǎn)結(jié)合能最低的前10個化合物結(jié)果如表4所示，4個靶點(diǎn)中對接較好的化合物為黃酮類和黃酮苷類，其次為三萜酸類和有機(jī)酸酯類。其中有3個黃酮類化合物在3個靶點(diǎn)中對接表現(xiàn)較好，這可能是PAINS[10]類化合物，也可能是藥物潛在有效分子。其中打分最優(yōu)化合物與靶點(diǎn)對接的二維圖如圖2所示。

表3 分子對接參數(shù)

表4 半夏厚樸湯與候選靶點(diǎn)分子對接結(jié)果

注：灰色線代表小分子結(jié)構(gòu)圖；粉色、紫色、深粉色虛線代表疏水作用力；深綠色虛線代表氫鍵作用力；淺綠色虛線代表范德華力；深紅色代表unfavorable donor-donor 作用。

2.4 飲片提取和樣品制備

表5 樣品提取方法

2.5 高分辨質(zhì)譜及CD分析

色譜柱Thermo Scientific Hypersil(100 mm×2.1mm)；流動相為0.5%甲酸溶液(A)-甲醇(B)，梯度洗脫(0～27 min，流動相A 10%～90%)，條件如下所示：流速0.3 mL/min；質(zhì)譜條件掃描范圍80～1 200m/z；源內(nèi)誘導(dǎo)碰撞解離電壓(CID)0.0 eV；分辨率70 000；流速0.3 mL/min；柱溫30 ℃；正/負(fù)離子檢測，進(jìn)樣量5 μL。

參考文獻(xiàn)[11]方法使用Compound Discoverer軟件分析質(zhì)譜數(shù)據(jù)，采用chemspider、m/zcloud等質(zhì)譜數(shù)據(jù)庫匹配化合物。使用Xcalibur軟件對檢測數(shù)據(jù)進(jìn)行處理。選擇m/zcloud best match 打分≥70.0分的化合物，認(rèn)為有較高的匹配度。

將高分辨質(zhì)譜數(shù)據(jù)經(jīng)Compound Discover 3.0軟件分析后共收集高匹配度化合物176個，經(jīng)比對，發(fā)現(xiàn)有19個化合物在半夏厚樸湯報道成分?jǐn)?shù)據(jù)庫中，這些成分在既往研究中已證實(shí)存在于半夏厚樸湯成分藥材中，又能在提取物中被檢測到，因此可認(rèn)為是半夏厚樸湯的潛在有效成分組，具體結(jié)果見表6，其中有3個在分子對接中表現(xiàn)較好(分子對接打分前10%)：野黃芩苷，芹菜素，熊果酸。含上述化合物的高分辨質(zhì)譜圖結(jié)果如圖3所示。

表6 半夏厚樸湯潛在有效成分組的親和能

注：①—野黃芩苷；②—芹菜素；③—熊果酸。

以19個對應(yīng)化合物的親和能和離子峰面積賦權(quán)作為指標(biāo)，一般認(rèn)為，親和能和結(jié)合能力呈對數(shù)關(guān)系，所以賦權(quán)公式按下式計算：

式中，S為以潛在有效成分指標(biāo)的分子對接親合能賦權(quán)值，i為化合物編號，A為離子峰面積，Af為親和能。計算結(jié)果見表7。

表7 潛在有效成分的親和能賦權(quán)值

3 討論和結(jié)論

本文采用整合藥理學(xué)方法，結(jié)合半夏厚樸湯傳統(tǒng)和現(xiàn)代臨床應(yīng)用分析其潛在靶點(diǎn)，這些靶點(diǎn)涉及神經(jīng)系統(tǒng)、免疫系統(tǒng)，采用虛擬篩選和分子對接技術(shù)，篩選半夏厚樸湯潛在有效成分，在篩選中表現(xiàn)良好的化合物主要包括黃酮類、三萜酸類等，其中熊果酸[12-13]、芹菜素[14]、野黃芩苷[15]等均有文獻(xiàn)報道參與IL6、PTGS等靶點(diǎn)和信號通路的作用。采用高分辨質(zhì)譜方法和化合物匹配軟件，分析了不同提取方法下半夏厚樸湯的成分，其中有121個匹配打分超過70.0分，有19個與半夏厚樸湯藥材報道成分吻合，有3個虛擬篩選打分靠前(前10%)。結(jié)合中藥傳統(tǒng)應(yīng)用方式，認(rèn)為這些化合物可能共同形成半夏厚樸湯潛在藥效成分組。

在本研究中，采用多種方法開展網(wǎng)絡(luò)藥理學(xué)靶點(diǎn)富集。最終搜集的4個靶點(diǎn)中，有3個為常見的炎癥、免疫相關(guān)信號通路，而GABR則為中醫(yī)神志病常見靶點(diǎn)，這與半夏厚樸湯用于抑郁癥的臨床應(yīng)用契合較好。說明本研究中網(wǎng)絡(luò)藥理學(xué)分析方法可靠，結(jié)果可信，在一定程度上彌補(bǔ)了近年來腫瘤免疫學(xué)研究的不足，為后續(xù)的小分子虛擬篩選提供了較為可靠的靶點(diǎn)結(jié)論。

中藥成分復(fù)雜，藥理作用廣泛，在特定疾病的治療中具備獨(dú)特優(yōu)勢。中藥成分發(fā)揮藥理作用，往往需要達(dá)到特定的劑量，本文研究過程中，對所涉及潛在有效化合物定量研究未能開展，發(fā)揮藥效的劑量因素，還要進(jìn)一步查閱文獻(xiàn)。本文以不同溶劑提取半夏厚樸湯，采用高分辨質(zhì)譜分析，可比較其中特定成分的大小差異，為后續(xù)的提取工藝研究奠定基礎(chǔ)。

目前的中藥系統(tǒng)藥理學(xué)和網(wǎng)絡(luò)藥理學(xué)研究中，對中藥的生物信息學(xué)研究較多，而對發(fā)揮藥理作用的物質(zhì)基礎(chǔ)部分，即成分的研究重視程度還不夠，在篩選得到潛在有效靶點(diǎn)和相關(guān)化合物后的分子對接、靶點(diǎn)的結(jié)合域選擇、親和能等主要指標(biāo)相對比較籠統(tǒng)，也往往不對靶蛋白自帶的抑制劑小分子進(jìn)行對照研究以保證分子對接的可靠性。在本研究中，為克服對接分子數(shù)量少帶來的數(shù)據(jù)偏差，引入虛擬篩選功能，對文獻(xiàn)報道的全部四百多個化合物均進(jìn)行對接運(yùn)算，這樣既能兼顧靶點(diǎn)的全面性，又能有效地衡量化合物針對特定靶點(diǎn)的有效性。相較而言，采用本研究方法更能篩選得到有較高活性的潛在有效化合物。

本研究確定了半夏厚樸湯發(fā)揮藥效的潛在活性成分(組)，但這些成分是否能有效代表半夏厚樸湯的絕大多數(shù)藥理活性還不得而知。后續(xù)的藥理實(shí)驗和分子生物學(xué)分析有待開展，希望通過后續(xù)藥效學(xué)驗證，對本文中所發(fā)現(xiàn)的潛在活性成分組進(jìn)行深入研究，驗證成分組是否能真正代替半夏厚樸湯整體藥效，以及這些成分的具體藥理活性。