茍美玲， 李 玲， 楊文宇， 何宇新

(西華大學(xué)生物工程學(xué)院，四川 成都 610039)

白芷為傘形科多年生草本植物白芷Angelicadahurica(Fisch.ex Hoffm)Benth.et Hook.f或杭白芷A.dahurica(Fisch.ex Hoffm.)Benth.et Hook.f.var.formosana(Boiss)Shan.et Yuan的干燥根。白芷成分復(fù)雜，作用靶點(diǎn)較多，作用方式復(fù)雜，已知化學(xué)成分對(duì)各靶點(diǎn)的作用尚難以用藥效學(xué)實(shí)驗(yàn)一一驗(yàn)證。為全面揭示白芷治療偏頭痛作用機(jī)制，一種理想化的方法是：將代謝預(yù)測(cè)技術(shù)與虛擬篩選方法相結(jié)合，針對(duì)白芷治療偏頭痛所涉及到的作用靶點(diǎn)，用計(jì)算機(jī)軟件逐一模擬體內(nèi)的相互作用，進(jìn)行對(duì)接、評(píng)分，從而初步篩選治療偏頭痛的活性成分，將篩選出的活性化合物再進(jìn)行有針對(duì)性的藥理活性實(shí)驗(yàn)驗(yàn)證，可避免常規(guī)篩選方法的盲目性，能夠大大縮短活性成分篩選時(shí)間，為進(jìn)一步揭示白芷治療偏頭痛作用機(jī)理奠定基礎(chǔ)，并能節(jié)省大量的研究經(jīng)費(fèi)，為天然藥物活性成分篩選建立快捷通道。

本實(shí)驗(yàn)在白芷化學(xué)成分研究的基礎(chǔ)上，以歐前胡素、異歐前胡素、花椒毒素、白當(dāng)歸素等20種化合物[1-6]為配體，借助Molegro Vitual Docker軟件的Docking模塊分別與5-HT、N-甲基-D-天門冬氨酸、多巴胺D3受體、降鈣素基因相關(guān)肽[7-9]、前列腺素、血管緊張素轉(zhuǎn)換酶以及一氧化氮合酶8種治療偏頭痛所涉及的作用靶點(diǎn)(見表1)進(jìn)行分子對(duì)接，根據(jù)這些已知化合物與靶點(diǎn)的對(duì)接親和能大小和對(duì)接位點(diǎn)數(shù)進(jìn)行計(jì)算機(jī)評(píng)分，分析白芷治療偏頭痛的活性成分[10-11]，為進(jìn)一步的動(dòng)物性實(shí)驗(yàn)設(shè)計(jì)提供有效依據(jù)。

表1 20種化學(xué)成分的待對(duì)接靶點(diǎn)

由于硝酸甘油型偏頭痛動(dòng)物模型的病理變化及癥狀表現(xiàn)與人類偏頭痛發(fā)作時(shí)的變化及表現(xiàn)相似性較好，且實(shí)驗(yàn)制作方法簡(jiǎn)便，實(shí)驗(yàn)條件較易控制，在復(fù)制大量的動(dòng)物模型研究使用中能大大縮短時(shí)間，是目前國(guó)內(nèi)研究中運(yùn)用最多的偏頭痛模型，故在本次研究中采用這種模型進(jìn)行研究。

1 材料與方法

1.1 軟件、硬件的準(zhǔn)備

1.1.1 軟件

1.1.1.1 Chemsketch 10.0軟件 繪制白芷中化學(xué)分子結(jié)構(gòu)軟件。

1.1.1.2 Corina在線轉(zhuǎn)換程序(http://www.molecular-networks.com) 對(duì)白芷中化學(xué)分子結(jié)構(gòu)進(jìn)行三維優(yōu)化。

1.1.1.3 Molegro Virtual Docker 2007 Trial(Molegro ApS，MVD 2007，Denmark， Aarhus，http://www.molegro. com) 具有高精度對(duì)接程序虛擬對(duì)接軟件。

1.1.2 設(shè)備 戴爾Dimension E530臺(tái)式計(jì)算機(jī)(intel core2 E4500處理器，1.80 GHz，2.00 Gb內(nèi)存， Microsoft Windows 7 Professional操作系統(tǒng))。

1.1.3 藥物 白當(dāng)歸素、東莨菪苷、硬脂酸、水合氧化前胡素、歐前胡素、白芷內(nèi)酯，蒸餾水(含0.5%吐溫-80)配成所需質(zhì)量濃度的混懸液備用。

1.1.4 試劑 硝酸甘油注射液，河南潤(rùn)弘制藥股份有限公司；麥角胺咖啡因片，哈爾濱泰華藥業(yè)股份有限公司。

1.1.5 動(dòng)物 SD大鼠，雌雄兼用，體質(zhì)量200 g左右，由成都中醫(yī)藥大學(xué)實(shí)驗(yàn)動(dòng)物中心提供，動(dòng)物合格證號(hào)scxl(川2013-05)。

1.2 活性靶點(diǎn)結(jié)構(gòu)和參考配體信息來源及處理 從蛋白質(zhì)數(shù)據(jù)庫(http://www.rcsb.org/pdb/home/home.do)下載治療偏頭痛作用涉及的以下靶點(diǎn)三維結(jié)構(gòu)模型：5-HT、N-甲基-D-天門冬氨酸受體、多巴胺D3受體、降鈣素基因相關(guān)肽、前列腺素、血管緊張轉(zhuǎn)換酶以及一氧化氮合酶。將以上8個(gè)靶點(diǎn)的單晶衍射圖中的乙酸分子、水分子、金屬離子用Molegro Virtual Docker軟件除去，只保留這8個(gè)蛋白分子與原始配體，最后得到可用于對(duì)接的5-HT、N-甲基-D-天門冬氨酸受體、多巴胺D3受體、降鈣素基因相關(guān)肽、前列腺素、血管緊張素轉(zhuǎn)換酶以及一氧化氮合酶的三維立體結(jié)構(gòu)。

1.3 配體化合物的三維結(jié)構(gòu)優(yōu)化 將白芷中已知的20種化合物作為對(duì)接所需的配體小分子化合物庫，首先利用ChemSketch 10.0軟件的繪圖功能，繪制出化合物的平面結(jié)構(gòu)圖(見表2)，然后從CORINA—Fast Generation of High-Quality 3D Molecular Models (http://www.molecular-networks.com/products/corina)中產(chǎn)生三維結(jié)構(gòu)并進(jìn)行優(yōu)化。

1.4 配體與活性靶點(diǎn)的分子對(duì)接 首先啟動(dòng)Molegro Virtual Docker 2007軟件，移入治療偏頭痛靶點(diǎn)三維結(jié)構(gòu)，同時(shí)檢測(cè)作用靶點(diǎn)活性結(jié)構(gòu)，然后移入已進(jìn)行三維結(jié)構(gòu)優(yōu)化的20種化合物分子。設(shè)置Molegro Virtual Docker 2007軟件中各種對(duì)接運(yùn)行參數(shù)，接著運(yùn)行對(duì)接模塊、進(jìn)行虛擬對(duì)接運(yùn)算，保存運(yùn)算結(jié)果。最后分析運(yùn)算結(jié)果和各種參數(shù)的相互關(guān)系，預(yù)測(cè)對(duì)接結(jié)果中各配體分子結(jié)合方式的分子對(duì)接評(píng)分，根據(jù)分子對(duì)接評(píng)分絕對(duì)值進(jìn)行排序，與參考配體評(píng)分進(jìn)行對(duì)比，從而篩選出有潛在活性的成分[12]。

表2 20種化學(xué)成分的二維結(jié)構(gòu)

1.5 虛擬篩選確定的活性成分對(duì)硝酸甘油型偏頭痛動(dòng)物模型的影響[13]取SD大鼠90只，雌雄各半，按體質(zhì)量隨機(jī)分層分為9組，每組10只，分別為空白組、模型組、陽性組、白當(dāng)歸素組、東莨菪苷組、硬脂酸組、水合氧化前胡素組、歐前胡素組、白芷內(nèi)酯組。除空白組外其余8組動(dòng)物進(jìn)行硝酸甘油型偏頭痛模型造模(于大鼠皮下注射硝酸甘油10 mg/kg，以出現(xiàn)前肢頻繁搔頭、雙耳發(fā)紅等預(yù)示造模成功)，同時(shí)記錄下各組大鼠耳紅出現(xiàn)時(shí)間。造模15 min后，陽性組予麥角胺咖啡因溶液灌胃，實(shí)驗(yàn)組則予上述虛擬篩選所得活性成分混懸液，空白組和模型組予等容積實(shí)驗(yàn)時(shí)用蒸餾水(含0.5%吐溫-80助溶)灌胃。給藥后4 h內(nèi)，觀察并對(duì)各組大鼠耳紅消失時(shí)間進(jìn)行記錄，同時(shí)以每0.5 h為一個(gè)時(shí)間段，觀察各組大鼠在每0.5 h爬籠次數(shù)和內(nèi)撓頭次數(shù)，以大鼠前肢抓籠蓋為爬籠行為癥狀的記錄標(biāo)志，以大鼠后肢撓頭為撓頭行為記錄標(biāo)志。

1.6 統(tǒng)計(jì)學(xué)方法 采用PEMS 3.1統(tǒng)計(jì)軟件包進(jìn)行。兩兩比較選用t檢驗(yàn)，多組比較選用方差分析，差齊性選用 LSD 或 Tukey 等多重比較，方差不齊選用Tamhane＇s T2多重比較。

2 結(jié)果與分析

2.1 白芷治療偏頭痛活性成分的虛擬篩選結(jié)果與分析

2.1.1 針對(duì)一氧化氮合酶虛擬篩選的結(jié)果 參考配體L-NAME的評(píng)分為-79.174，白芷的化學(xué)成分中與一氧化氮合酶對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于L-NAME評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表3。

表3 與一氧化氮合酶靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

2.1.2 針對(duì)血管緊張素轉(zhuǎn)換酶的虛擬篩選結(jié)果 參考配體血管緊張素轉(zhuǎn)換酶抑制劑的評(píng)分為-131.852，白芷的化學(xué)成分中與血管緊張素轉(zhuǎn)換酶對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值排序前6的分子見表4。

表4 與血管緊張素轉(zhuǎn)換酶靶點(diǎn)分子對(duì)接評(píng)分排序前6的化合物

2.1.3 針對(duì)腎上腺素-β虛擬篩選的結(jié)果 參考配體普萘洛爾的評(píng)分為-104.136，白芷的化學(xué)成分中與腎上腺素-β對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于普萘洛爾評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表5。

表5 與腎上腺素-β靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

2.1.4 針對(duì)前列腺素的虛擬篩選結(jié)果 參考配體吲哚美辛的評(píng)分為-110.224，白芷的化學(xué)成分中與前列腺素對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值排序前6的分子見表6。

表6 與前列腺素靶點(diǎn)分子對(duì)接評(píng)分排序前6的化合物

2.1.5 針對(duì)降鈣素基因相關(guān)肽虛擬篩選的結(jié)果 參考配體吳茱萸次堿的評(píng)分為-108.399，白芷的化學(xué)成分中與降鈣素基因相關(guān)肽對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于吳茱萸次堿評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表7。

2.1.6 針對(duì)多巴胺D3受體虛擬篩選的結(jié)果 參考配體BP-897的評(píng)分為-130.229，白芷的化學(xué)成分中與多巴胺D3受體對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于MK-801評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表8。

表7 與降鈣素基因相關(guān)肽靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

表8 與多巴胺D3受體靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

2.1.7 針對(duì)N-甲基-D天門冬氨酸受體虛擬篩選的結(jié)果 參考配體MK-801的評(píng)分為-92.204，白芷的化學(xué)成分中與N-甲基-D天門冬氨酸受體對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于MK-801評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表9。

2.1.8 針對(duì)5-HT受體虛擬篩選的結(jié)果 參考配體阿米替林的評(píng)分為-92.204，白芷的化學(xué)成分中與5-HT受體對(duì)接時(shí)分子對(duì)接評(píng)分絕對(duì)值大于阿米替林評(píng)分90%的分子且MolDockScore評(píng)分排序前10的分子見表10。

表9 與N-甲基-D天門冬氨酸受體靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

表10 與5-HT受體靶點(diǎn)分子對(duì)接評(píng)分排序前10的化合物

2.1.9 分析 以上虛擬篩選結(jié)果說明：(1)針對(duì)上述8個(gè)作用靶點(diǎn)篩出的理論活性成分可能在一定程度上反映了白芷治療偏頭痛的物質(zhì)基礎(chǔ)，值得進(jìn)一步進(jìn)行藥理實(shí)驗(yàn)驗(yàn)證；(2)白芷治療偏頭痛可能主要針對(duì)多巴胺受體D3、降鈣素基因相關(guān)肽、前列腺素、血管緊張素轉(zhuǎn)換酶4個(gè)靶點(diǎn)，這需要進(jìn)一步擴(kuò)大靶點(diǎn)篩選范圍進(jìn)行確證；(3)部分化學(xué)成分可能具有多靶點(diǎn)活性，見表11。

表11 具有多靶點(diǎn)活性的化學(xué)成分

表11結(jié)果顯示：通過白芷中20種化學(xué)成分針對(duì)8個(gè)治療偏頭痛相關(guān)靶點(diǎn)的虛擬篩選，并與參考配體評(píng)分絕對(duì)值進(jìn)行對(duì)比，發(fā)現(xiàn)白芷中的白當(dāng)歸素、東莨菪苷、硬脂酸、水合氧化前胡素、歐前胡素、白芷內(nèi)酯6種成分對(duì)多個(gè)治療偏頭痛相關(guān)靶點(diǎn)有抑制活性。正是這種可能存在的單分子的多靶點(diǎn)作用，為中藥針對(duì)單靶點(diǎn)藥效不強(qiáng)提供了可能的解決方法，使得白芷治療偏頭痛作用表現(xiàn)出整體治療、調(diào)制并舉的特點(diǎn)。現(xiàn)目前中醫(yī)藥研究領(lǐng)域的難點(diǎn)之一就有用當(dāng)代醫(yī)學(xué)語言解釋白芷這種整體治療的特色與原理，因此這一結(jié)果為進(jìn)一步研究白芷整體治療偏頭痛提供了有利依據(jù)。

2.2 藥效學(xué)實(shí)驗(yàn)結(jié)果與分析

2.2.1 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠耳紅癥狀的影響結(jié)果，見表12。

表12結(jié)果顯示：虛擬篩選所得活性成分白當(dāng)歸素、東莨菪苷、水合氧化前胡素、歐前胡素可明顯縮短耳紅消失時(shí)間，與模型組比較有顯著差異(P<0.05)，硬脂酸、白芷內(nèi)酯也存在相同的作用趨勢(shì)，但無統(tǒng)計(jì)學(xué)意義。

2.2.2 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠撓頭次數(shù)的影響結(jié)果，見表13。

表12 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠耳紅癥狀的影響

表13 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠撓頭次數(shù)的影響

表13結(jié)果顯示：虛擬篩選所得活性成分白當(dāng)歸素、東莨菪苷、水合氧化前胡素、歐前胡素可明顯減少偏頭痛模型大鼠撓頭次數(shù)，與模型組比較有顯著差異(P<0.05)，硬脂酸、白芷內(nèi)酯也存在相同的作用趨勢(shì)，但無統(tǒng)計(jì)學(xué)意義。

2.2.3 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠爬籠次數(shù)的影響結(jié)果，見表14。

表14 虛擬篩選確定的活性成分對(duì)硝酸甘油致偏頭痛大鼠爬籠次數(shù)的影響

表14結(jié)果顯示：虛擬篩選所得活性成分白當(dāng)歸素、東莨菪苷、硬脂酸、水合氧化前胡素、歐前胡素、白芷內(nèi)酯對(duì)大鼠爬籠次數(shù)無明顯影響。

2.2.4 分析 上述硝酸甘油型偏頭痛動(dòng)物模型的研究結(jié)果顯示：虛擬篩選所得活性成分白當(dāng)歸素、東莨菪苷、水合氧化前胡素、歐前胡素對(duì)硝酸甘油型偏頭痛模型有明顯影響，硬脂酸、白芷內(nèi)酯對(duì)硝酸甘油型偏頭痛模型的影響不明顯。

3 結(jié)論

3.1 本實(shí)驗(yàn)所用的MVD分子對(duì)接軟件的基本原理是利用各種參數(shù)主要包括電性參數(shù)、立體參數(shù)、疏水參數(shù)等對(duì)配體分子與靶點(diǎn)分子的相互作用進(jìn)行模擬，從而給出白芷化學(xué)成分與所選取的8個(gè)治療偏頭痛相關(guān)靶點(diǎn)相互作用的評(píng)分，根據(jù)與參考配體評(píng)分絕對(duì)值比較預(yù)測(cè)配體的活性。這些所用的的參數(shù)都是建立在大量實(shí)驗(yàn)數(shù)據(jù)基礎(chǔ)之上的，并且決定著對(duì)接軟件精度。

3.2 虛擬篩選能為隨后的藥效學(xué)驗(yàn)證實(shí)驗(yàn)提供有用的線索和基礎(chǔ)，將虛擬篩選與藥效學(xué)實(shí)驗(yàn)研究結(jié)合將能更快速和全面地明確白芷治療偏頭痛有效成分和多種藥效指標(biāo)間的復(fù)雜對(duì)應(yīng)關(guān)系，從而為進(jìn)一步揭示白芷治療偏頭痛作用機(jī)制奠定基礎(chǔ)。

3.3 基于虛擬篩選的白芷藥效物質(zhì)基礎(chǔ)研究是在中醫(yī)理論指導(dǎo)下的研究，現(xiàn)目前整體治療偏頭痛作用還不能被白芷中的任何一種化學(xué)成分全面反映和體現(xiàn)，這也是與大多數(shù)學(xué)者發(fā)現(xiàn)中藥活性成分越純，但藥理作用越差的重要原因所在，因此在虛擬篩選中除了考慮配體與受體的相互作用外，還要考慮分子的化學(xué)穩(wěn)定性、類藥性、毒性和合成的難易性以及體內(nèi)轉(zhuǎn)運(yùn)和代謝分布等藥代動(dòng)力學(xué)方面的性質(zhì)。

3.4 硝酸甘油型偏頭痛動(dòng)物模型的研究結(jié)果顯示：白芷中白當(dāng)歸素、東莨菪苷、水合氧化前胡素、歐前胡素對(duì)硝酸甘油型偏頭痛的影響較明顯，通過藥效學(xué)實(shí)驗(yàn)對(duì)虛擬篩選出的白芷治療偏頭痛活性成分的有效性進(jìn)行了初步驗(yàn)證。

參考文獻(xiàn)：

[1] 黃 欣,蘇樂群,邵 偉,等.白芷香豆素類化合物的成分分析及CO2超臨界萃取工藝的研究[J].中國(guó)藥房,2005,16(11):824-826.

[2] 賈 英,孫振蛟,包文芳.正交設(shè)計(jì)法優(yōu)選白芷的提取工藝[J].沈陽藥科大學(xué)學(xué)報(bào),2005,22(3):217-219.

[3] 鄧捷圓,高廣慧,趙春杰,等.HPLC法同時(shí)測(cè)定不同產(chǎn)地白芷中2種香豆素的含量[J].沈陽藥科大學(xué)學(xué)報(bào),2004,21(5):354-357.

[4] 王夢(mèng)月,賈敏如,馬逾英,等.不同入藥部分及不同加工方法對(duì)白芷香豆素類成分含量的影響[J].中藥材,2004,27(11):826-828.

[5] 楊 松,賈 英,畢開順.RP-HPLC法測(cè)定桂枝中4種活性成分的含量[J].藥物分析雜志,2004,24(2):143-146.

[6] 趙 越,段秀俊,郭旭輝,等.對(duì)紅花、白芷兩種中藥TLC鑒別方法的改進(jìn)[J].中醫(yī)藥研究,2000,16(2):54-55.

[7] 連亞軍,騰軍放,王左生,等.偏頭痛患者血漿一氧化氮和降鈣素基因相關(guān)肽含量測(cè)定及發(fā)病機(jī)制探討[J].中華神經(jīng)外科雜志,1999，32(2):94-95.

[8] 鐘廣偉,李 煒,鄧干初,等.針刺對(duì)偏頭痛大鼠腦內(nèi)降鈣素基因相關(guān)肽基因表達(dá)的影響[J].中國(guó)臨床康復(fù),2002,6(24):3680-3681.

[9] 任 旭,牛爭(zhēng)平.降鈣素基因相關(guān)肽在偏頭痛發(fā)病機(jī)制中的研究進(jìn)展[J].中西醫(yī)結(jié)合心腦血管病雜志,2004,2(4):229-231.

[10] 吳 敏,陳美娟.偏頭痛發(fā)病機(jī)制的研究進(jìn)展[J].四川生理科學(xué)雜志,2007,29(4):173-176.

[11] 魯啟洪,徐金枝.偏頭痛發(fā)病機(jī)制研究進(jìn)展[J].國(guó)際中醫(yī)中藥雜志,2008,30(2):125-132.

[12] 楊文宇,萬德光,楊鑫嵎.虛擬篩選輔助揭示中藥藥效物質(zhì)基礎(chǔ)的思路與初步實(shí)踐[J].中草藥,2011,42(9):1655-1672.

[13] 秦旭華,李祖?zhèn)?朱果夫,等.川白芷總香豆素治療偏頭痛的部分藥效學(xué)實(shí)驗(yàn)研究[J].四川中醫(yī),2008,6(15):4-6.