劉福和，陳少軍，倪文娟（浙江醫(yī)藥高等?？茖W(xué)校藥學(xué)院，浙江寧波 315100）

川芎中抗血栓活性成分的計(jì)算機(jī)虛擬篩選研究Δ

劉福和＊，陳少軍，倪文娟（浙江醫(yī)藥高等?？茖W(xué)校藥學(xué)院，浙江寧波 315100）

目的：運(yùn)用計(jì)算機(jī)輔助藥物設(shè)計(jì)尋找川芎抗血栓活性成分。方法：以“血栓癥”為關(guān)鍵詞，從治療靶蛋白數(shù)據(jù)庫中搜索篩選川芎抗血栓靶蛋白；再從蛋白質(zhì)數(shù)據(jù)庫中查詢目標(biāo)靶蛋白并下載其三維結(jié)構(gòu)，采用蛋白質(zhì)預(yù)處理工具確定活性區(qū)域中心坐標(biāo)，應(yīng)用PyRx和Discovery Studio Visualizer軟件對從臺灣中醫(yī)藥資料庫下載的247個(gè)川芎小分子與靶蛋白進(jìn)行分子對接，通過結(jié)合能篩選出活性成分并分析結(jié)合作用力。結(jié)果：篩選出4個(gè)活性成分即新綠原酸、1-H-苯并咪唑-2-胺、3，8-二羥基酰內(nèi)酯、川芎三萜，其分別與凝血酶、抗凝血酶Ⅲ、凝血因子Ⅹa、血栓調(diào)節(jié)蛋白具有較高結(jié)合活性，結(jié)合能分別為－6.1、－4.5、－7.7、－8.6 kJ/mol；分析結(jié)果顯示范德華力、靜電作用力在對接中發(fā)揮著重要作用。結(jié)論：新綠原酸、1-H-苯并咪唑-2-胺、3，8-二羥基酰內(nèi)酯、川芎三萜可能是中藥川芎抗血栓的活性成分。

血栓癥；川芎；計(jì)算機(jī)輔助藥物設(shè)計(jì)；虛擬篩選；靶蛋白；分子對接；活性成分

血栓形成及血栓栓塞性疾病是臨床常見急重癥類疾病，祖國傳統(tǒng)醫(yī)學(xué)將其歸為“血瘀”范疇，治療以益氣活血、祛瘀通絡(luò)為主；而川芎作為常用活血化瘀中藥之一，已被證實(shí)具有良好的抗血栓作用；但其提取物中含有揮發(fā)油、生物堿、酚、有機(jī)酸及苯酞內(nèi)酯等多種化學(xué)物質(zhì)[1-3]，常規(guī)檢測手段難以有效確定活性成分。計(jì)算機(jī)輔助藥物設(shè)計(jì)（Computeraided drug design，CADD）是一種伴隨計(jì)算機(jī)技術(shù)、蛋白質(zhì)組學(xué)以及人類基因組計(jì)劃發(fā)展而形成的新藥開發(fā)技術(shù)，其以計(jì)算機(jī)化學(xué)研究技術(shù)為基礎(chǔ)，通過對藥物與受體生物大分子的關(guān)系進(jìn)行模擬、計(jì)算及預(yù)算，從而設(shè)計(jì)并優(yōu)化合成先導(dǎo)化合物[4]。已有研究顯示，CADD技術(shù)可快速有效地確定中藥有效成分和作用機(jī)制，科學(xué)客觀地對中藥獨(dú)特藥理作用進(jìn)行闡述，已成為中藥新藥設(shè)計(jì)關(guān)鍵技術(shù)之一[5]。在CADD模式下進(jìn)行中藥活性成分篩選，需首先對化合物分子進(jìn)行虛擬篩選，在降低篩選化合物數(shù)量的同時(shí)提高先導(dǎo)化合物發(fā)現(xiàn)效率[6]。川芎中川芎嗪、阿魏酸的抗血栓作用已被廣泛證實(shí)[7]，本次研究采用PyRx軟件的Autodock vina模塊對中藥川芎有抗血栓活性的其他成分進(jìn)行篩選和評價(jià)，擬探討運(yùn)用CADD尋找中藥川芎新的抗血栓活性成分的可行性及其意義，為中藥川芎的現(xiàn)代化研究提供理論基礎(chǔ)和實(shí)踐指導(dǎo)。

1 資料與方法

1.1 資料來源

治療靶蛋白數(shù)據(jù)庫（Therapeutic target database，TTD，新加坡國立大學(xué)開發(fā)，下載網(wǎng)址：http：//bidd.nus. edu.sg/group/ttd/）；蛋白質(zhì)數(shù)據(jù)庫（PDB，美國紐約結(jié)構(gòu)生物信息學(xué)研究實(shí)驗(yàn)室開發(fā)，下載網(wǎng)址：http：//www.rcsb.org/pdb/home/home.do）；蛋白質(zhì)預(yù)處理工具[Protein preparing tool，中國醫(yī)科大學(xué)（臺灣）計(jì)算與系統(tǒng)生物學(xué)實(shí)驗(yàn)室開發(fā)，下載網(wǎng)址：http：//dock.cmu.edu.tw/ligand. php]；臺灣中醫(yī)藥資料庫[Traditional Chinese Medicine Dadabase@Taiwan，中國醫(yī)科大學(xué)（臺灣）計(jì)算與系統(tǒng)生物學(xué)實(shí)驗(yàn)室開發(fā)，下載網(wǎng)址：http：//tcm.cmu.edu.tw/]；PyRx軟件0.8版（美國Scripps研究所開發(fā)）；分子模擬圖形界面軟件（Discovery Studio Visualizer）3.5版（美國Accelrys公司開發(fā)）。

1.2 方法

1.2.1 數(shù)據(jù)挖掘 以“血栓癥”“Thrombosis”作為關(guān)鍵詞，從TTD中搜索與其相關(guān)的蛋白，作為川芎抗血栓的靶蛋白，并且對相應(yīng)的文獻(xiàn)進(jìn)行查詢及篩選，選取試驗(yàn)所需的目標(biāo)靶蛋白[8-11]。

1.2.2 受體結(jié)構(gòu)信息來源及處理 根據(jù)篩選結(jié)果，從PDB中搜索目標(biāo)靶蛋白并下載其三維結(jié)構(gòu)，然后采用Protein preparing tool進(jìn)行預(yù)處理，主要包括刪去自帶的配體小分子、分配相關(guān)電荷及加氫，然后以自帶配體位置為中心，確定活性區(qū)域的中心坐標(biāo)；預(yù)處理完成后蛋白可用于后續(xù)分子對接。以凝血酶為例，預(yù)處理前、后的三維結(jié)構(gòu)圖見圖1。

圖1 凝血酶預(yù)處理前后的三維結(jié)構(gòu)圖Fig 1 Three-dimensional structures of thrombin beforeand after pretreatment

1.2.3 配體小分子化合物庫準(zhǔn)備 從臺灣中醫(yī)藥資料庫中下載本次研究的247個(gè)川芎分子結(jié)構(gòu)，下載的川芎各小分子均符合LipinskiCA提出的“類藥五原則”：（1）分子量小于500；（2）氫鍵給體數(shù)目小于5；（3）氫鍵受體數(shù)目小于10；（4）脂水分配系數(shù)小于5；（5）可旋轉(zhuǎn)鍵的數(shù)量不超過10個(gè)[12]。將各分子mol2結(jié)構(gòu)轉(zhuǎn)換成PyRx軟件能認(rèn)可的pdbqt格式，三維結(jié)構(gòu)均采用MM 2力場進(jìn)行優(yōu)化。應(yīng)用PyRx軟件中的Open babel模塊對川芎各個(gè)小分子進(jìn)行能量最小優(yōu)化。

1.2.4 分子對接 川芎活性成分與靶蛋白的匹配情況應(yīng)用PyRx軟件中Autodock vina模塊進(jìn)行分子對接，根據(jù)結(jié)合能的大小進(jìn)行打分。結(jié)合位點(diǎn)坐標(biāo)參照各蛋白自帶配體的坐標(biāo)，運(yùn)用Grid box（盒子）將蛋白質(zhì)嵌入一個(gè)三維網(wǎng)格中，Grid box大小設(shè)定為9埃米（?，1?＝10－10m）×9?×9?；另血栓調(diào)節(jié)蛋白（1DQB）晶體結(jié)構(gòu)因不自帶配體，結(jié)合位點(diǎn)及Grid box大小采用默認(rèn)值，可以涵蓋各個(gè)活性位點(diǎn)。除具體坐標(biāo)外的其他參數(shù)均采用默認(rèn)值。

1.2.5 數(shù)據(jù)分析與圖像處理 采用Discovery Studio Visualizer對分子對接結(jié)果進(jìn)行分析并制作三維和二維圖像。

2 結(jié)果

2.1 數(shù)據(jù)挖掘結(jié)果

在TTD中，共檢索到30余種靶蛋白，涉及止血系統(tǒng)和纖溶系統(tǒng)等多個(gè)環(huán)節(jié)。選擇與血栓癥密切關(guān)聯(lián)的4個(gè)靶蛋白，即凝血酶、抗凝血酶Ⅲ、凝血因子Ⅹa、血栓調(diào)節(jié)蛋白，而這些靶蛋白也是川芎抗血栓作用研究的熱點(diǎn)。

2.2 分子對接結(jié)果

下載247個(gè)川芎小分子配體與靶蛋白進(jìn)行對接，結(jié)果顯示新綠原酸（Neochlorogenic acid）、1-H-苯并咪唑-2-胺（1-H-benzim idazole-2-amine）、3，8-二羥基酰內(nèi)酯（3，8-dihydrodiligustilide）、川芎三萜（Chuanxiongterpene）分別與凝血酶、抗凝血酶Ⅲ、凝血因子Ⅹa、血栓調(diào)節(jié)蛋白結(jié)合具有較低結(jié)合能，可作為靶蛋白抑制劑先導(dǎo)化合物加以研究。靶蛋白的靶標(biāo)名稱及篩選結(jié)合位點(diǎn)坐標(biāo)見表1，川芎活性成分與靶蛋白對接的結(jié)合能見表2。

表1 靶蛋白的靶標(biāo)名稱及篩選結(jié)合位點(diǎn)坐標(biāo)Tab 1 Target nam e of target protein and coordinates of screening binding site

表2 川芎活性成分與靶蛋白對接的結(jié)合能（kJ/mol）Tab 2 Binding energy of active ingredientsof chuanxiong rhizoma docking w ith target protein（k J/mol）

2.3 對接結(jié)果分析

2.3.1 新綠原酸與凝血酶的相互作用 川芎的247個(gè)小分子配體與凝血酶進(jìn)行分子對接后，新綠原酸與凝血酶的結(jié)合能最低，為－6.1 kJ/mol；其對接效果見圖2（圖C中“HIS”等為氨基酸縮寫，下同）。

圖2 新綠原酸與凝血酶分子對接效果圖Fig 2 Docking renderings of chlorogenic acid w ith throm binm olecules

從圖2C可以看出，范德華力在此對接中起重要作用；同時(shí)GLY233與新綠原酸形成了氫鍵。

2.3.2 1-H-苯并咪唑-2-胺與抗凝血酶Ⅲ的相互作用 川芎的247個(gè)小分子配體與抗凝血酶Ⅲ進(jìn)行分子對接后，1-H-苯并咪唑-2-胺與抗凝血酶Ⅲ結(jié)合能最低，為－4.5 kJ/mol；其對接效果見圖3。

圖3 1-H-苯并咪唑-2-胺與抗凝血酶Ⅲ分子對接效果圖Fig 3 Docking renderings of 1-H-benzim idazole-2-am inew ith antithrombinⅢ

從圖3C可以看出，靜電作用力在此對接中起著重要作用；同時(shí)，LYS193、NAG861與1-H-苯并咪唑-2-胺形成了鹽橋，這是離子π-π相互作用引起的。

2.3.3 3，8-二羥基酰內(nèi)酯與凝血因子Ⅹa的相互作用 川芎的247個(gè)小分子配體與凝血因子Ⅹa進(jìn)行分子對接后，3，8-二羥基酰內(nèi)酯與凝血因子Ⅹa的結(jié)合能最低，為－7.7 kJ/mol；其對接效果見圖4。

圖4 3，8-二羥基酰內(nèi)酯與凝血因子Ⅹa對接效果圖Fig 4 Docking renderings of 3，8-dihydrodiligustilidew ith coagulation factorⅩa

從圖4C可以看出，范德華力在此對接中起著重要作用。

2.3.4 川芎三萜與血栓調(diào)節(jié)蛋白的相互作用 川芎的247個(gè)小分子配體與血栓調(diào)節(jié)蛋白進(jìn)行分子對接后，川芎三萜與血栓調(diào)節(jié)蛋白的結(jié)合能最低，為－8.6 kJ/mol；其對接效果見圖5。

圖5 川芎三萜與血栓調(diào)節(jié)蛋白的對接效果圖Fig 5 Docking renderings of chuanxiongterpene w ith thrombomodulin

從圖5C可以看出，范德華力與靜電作用力在此對接中起著重要作用。

3 討論

傳統(tǒng)藥物研究方法為中醫(yī)藥現(xiàn)代化研究奠定了良好基礎(chǔ)，但由于中藥自身具有多組分的特點(diǎn)，涉及有效藥物成分、配伍、人體內(nèi)吸收及代謝等問題存在，常規(guī)試驗(yàn)研究技術(shù)用于中藥活性成分檢測具有較大盲目性[13]。首先，中藥中含有多類微量成分故無法通過傳統(tǒng)分離篩選方法確定其作用機(jī)制及藥理作用；其次，中藥篩選目前多采用基于疾病特征的動物實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行，難以從分子機(jī)制上闡明中藥作用機(jī)制；最后，目前尚無可描述中藥傳統(tǒng)功效及科學(xué)內(nèi)涵的符合中醫(yī)藥理論評價(jià)的模型[14]。

臺灣中醫(yī)藥資料庫是全世界目前最大的中醫(yī)藥小分子化合物數(shù)據(jù)庫之一，數(shù)據(jù)庫中包含了453種中藥的小分子化合物，共含有超過20 000種的小分子化學(xué)成分；數(shù)據(jù)庫提供了中藥小分子活性成分虛擬篩選、分子對接及分子動力學(xué)模擬等應(yīng)用工具，可有效解決中西藥分歧[15]，故本試驗(yàn)采用此庫作為依據(jù)。

相較于傳統(tǒng)藥物研發(fā)方法，CADD技術(shù)能夠顯著提高藥物設(shè)計(jì)速度，減少藥物創(chuàng)新工作所需人力和物力。而虛擬篩選技術(shù)應(yīng)用可在生物活性篩選前以計(jì)算機(jī)對化合物分子進(jìn)行提前篩選，主要目的為降低實(shí)際篩選化合物數(shù)目，提高先導(dǎo)化合物發(fā)現(xiàn)效率；其中分子對接技術(shù)應(yīng)用最為普遍，即從小分子數(shù)據(jù)庫中檢索能與之匹配的候選化合物，通過計(jì)算機(jī)對接技術(shù)評價(jià)已知藥物分子與受體大分子間是否可相互結(jié)合，并對兩者間相互作用親和力及結(jié)合模式進(jìn)行評價(jià)[16]。分子對接由于從整體上考慮配體與受體的結(jié)合效果，因而能更好地避免其他方法的弊端，例如局部作用較好而整體結(jié)合欠佳的狀況[17]。本次研究中采用的PyRx軟件是應(yīng)用較為廣泛的對接軟件之一，其能夠自動模擬配體在受體活性位點(diǎn)的作用情況，記錄最優(yōu)相互作用方式，并具有配體三維數(shù)據(jù)庫搜索功能。

本研究運(yùn)用PyRx軟件的Autodock vina模塊對中藥川芎中除川芎嗪和阿魏酸外的其他抗血栓活性成分進(jìn)行篩選，從247種川芎成分中篩選出活性分子新綠原酸、1-H-苯并咪唑-2-胺、3，8-二羥基酰內(nèi)酯、川芎三萜分別與凝血酶、抗凝血酶Ⅲ、凝血因子Ⅹa、血栓調(diào)節(jié)蛋白結(jié)合活性較高；然后運(yùn)用Discovery Studio Visualizer軟件對分子對接最匹配的蛋白復(fù)合物相互作用力進(jìn)行分析，其結(jié)合能分別為－6.1、－4.5、－7.7、－8.6 kJ/mol，證實(shí)范德華力和靜電作用力在各自的對接中起著重要作用。已有研究顯示，新綠原酸可以減少剪切應(yīng)激引起的血小板聚集；凝血酶GLY233對酶水解專一性具有決定作用，并可識別和定位底物[18]。從圖2C可以看出，新綠原酸與凝血酶核心氨基酸均有相互作用。以上均說明川芎活性成分與靶蛋白的核心氨基酸存在相互的作用力。而二維圖可以幫助了解川芎活性成分與受體對接的作用力和位置結(jié)合等情況，有助于闡述抗血栓作用機(jī)制。

綜上所述，活性分子新綠原酸、1-H-苯并咪唑-2-胺、3，8-二羥基酰內(nèi)酯、川芎三萜可能是中藥川芎抗血栓的活性成分。本文建立的方法彌補(bǔ)了傳統(tǒng)藥理學(xué)實(shí)驗(yàn)方法的不足，提高了先導(dǎo)化合物的發(fā)現(xiàn)效率，在中藥藥理學(xué)研究中具有可行性。

[1] 穆桂榮.川芎嗪注射液的臨床應(yīng)用進(jìn)展[J].中國基層醫(yī)藥，2011，18（9）：1283-1284.

[2] 李偉霞，唐于平，郭建明，等.比較評價(jià)當(dāng)歸川芎配伍對急性血瘀大鼠血液流變學(xué)及凝血功能的影響[J].中國中西醫(yī)結(jié)合雜志，2012，32（6）：806-811.

[3] 陳斌，賈曉斌，黃洋.通脈復(fù)方微丸多元釋藥系統(tǒng)中川芎酚酸速釋單元的制備工藝研究[J].中國藥房，2012，23（39）：3676-3678.

[4] 吳瓊，康宏，王歡歡，等.AT1受體的中藥活性成分篩選模型及其作用機(jī)理研究[J].化學(xué)學(xué)報(bào)，2012，70（6）：796-802.

[5] 宋芳嬌，曾克武，王學(xué)美，等.中藥復(fù)方有效成分組相關(guān)研究方法的研究進(jìn)展[J].環(huán)球中醫(yī)藥，2012，5（12）：951-955.

[6] Liang XL，Liao ZG，Zhu JY，etal.The absorption characterization effects and mechanism of Radix Angelicae dahuricae extracts on baicalin in Radix Scutellariae using in vivo and in vitro absorption models[J].J Ethnopharmacol，2012，139（1）：52-57.

[7] 徐筱杰.計(jì)算機(jī)輔助藥物分子設(shè)計(jì)[M].北京：化學(xué)工業(yè)出版社，2004：3-10.

[8] 成春英，孫勇，文志斌，等.川芎嗪對凝血酶誘導(dǎo)血管內(nèi)皮細(xì)胞組織因子表達(dá)的影響[J].南方醫(yī)科大學(xué)學(xué)報(bào)，2009，29（8）：1743-1747.

[9] 丁紅香，王明山，王良興.低O2高CO2肺動脈高壓及川芎嗪保護(hù)的大鼠血漿抗凝血酶-Ⅲ變化及意義[J].心肺血管病雜志，2000，19（2）：138-140.

[10] 成春英，孫勇，文志斌，等.川芎嗪抑制凝血酶誘導(dǎo)血管內(nèi)皮細(xì)胞組織因子表達(dá)的機(jī)制[J].中國動脈硬化雜志，2010，18（3）：184-188.

[11] 熊石龍，文志斌，何曉凡，等.川芎嗪抑制腫瘤壞死因子α誘導(dǎo)人臍靜脈血管內(nèi)皮細(xì)胞組織因子表達(dá)的機(jī)制研究

[J].廣東醫(yī)學(xué)，2008，29（4）：545-548.

[12] Lu W J，Xu C，Pei Z，et al.The tamoxifenmetabolite norendoxifen is a potentand selective inhibitor of aromatase（CYP19）and a potential lead compound for novel therapeutic agents[J].Breast Cancer Res Treat，2012，133（1）：99-109.

[13] Chikan NA，Bhavaniprasad V，Anbarasu K，et al.From natural products to drugs for epimutation computer-aided drug design[J].Appl Biochem Biotechnol，2013，170（1）：164-175.

[14] 劉磊，馬英，王潤玲，等.基于神經(jīng)氨酸酶結(jié)構(gòu)的流感病毒抑制劑的計(jì)算機(jī)輔助藥物設(shè)計(jì)[J].天津醫(yī)科大學(xué)學(xué)報(bào)，2013，19（3）：192-195.

[15] Tsai TY，Chang KW，Chen CY.iScreen：world’s first cloud-computing web server for virtual screening and de novo drug design based on TCM database@Taiwan[J].J ComputAided MolDes，2011，25（6）：525-531.

[16] Singh R，Balupuri A，Sobhia ME，et al.Development of 3D-pharmacophoremodel followed by successive virtual screening，molecular docking and ADME studies for the design of potent CCR2 antagonists for inflammation-driven diseases[J].Molecular Simulation，2013，39（1）：49-58.

[17] Gomes FV，A lves FH，Guimar?es FS，et al.Cannabidiol administration into the bed nucleus of the stria terminalis alters cardiovascular responses induced by acute restraint stress through 5-HT1A receptor[J].Eur Neuropsychopharmacol，2013，23（9）：1096-1104.

[18] Lao YM，Jiang JG，Yan L，et al.Application ofmetabonom ic analytical techniques in themodernization and toxi-

Study on the Computer Virtual Screening of Antithrombotic Active Ingredients in Chuanxiong Rhizoma

LIU Fuhe，CHEN Shaojun，NIWenjuan（School of Pharmacy，Zhejiang Pharmaceutical College，Zhejiang Ningbo 315100，China）

OBJECTIVE：To find the active ingredient of on antithrombotic chuanxiong rhizoma using computer aided drug design.METHODS：Using“thrombosis”as keyword，thrombosis related proteinswere searched and screened in therapeutic target database；target proteins’three-dimensional structure were downloaded in protein database，then the protein preparing tool were used to determine the coordinates of the active area center.PyRx software and Discovery Studio Visualizer were used to match the 247 smallmolecules of chuanxiong rhizoma w ith target protein that downloaded from Taiwan traditional Chinesemedicine database.The active molecules were screened and binding force was analyzed.RESULTS：Active molecules of neochlorogenic acid，1-H-benzimidazole-2-amine，3，8-dihydrodiligustilide，chuanxiongterpene were selected by blinding energy，and there were high binding activity among these activemolecules，thrombin，antithrombinⅢ，coagulation factorⅩa and thrombomodulin，and the binding energy were－6.1，－4.5，－7.7，－8.6 kJ/mol.Analysis results showed van Edward force and electrostatic interactions played an important role in their respective docking.CONCLUSIONS：Neochlorogenic acid，1-H-benzim idazole-2-am ine，3，8-dihydrodiligusti-lide，chuanxiongterpenemay be the antithrombotic activity ingredients of Chuanxiong rhizoma.

Thrombogenesis；Chuanxiong rhizoma；Computer aided drug design；Virtual screening；Target protein；Molecular docking；Active ingredients

R9-39

A

1001-0408（2017）16-2182-05 DOI 10.6039/j.issn.1001-0408.2017.16.06

浙江省自然科學(xué)基金資助項(xiàng)目（No.LY15H280009）

＊講師，碩士。研究方向：中藥藥理。電話：0574-88223170。E-mail：89994226@qq.com