劉靖麗,郭 勇

(1.陜西中醫(yī)藥大學(xué) 藥學(xué)院,陜西 咸陽(yáng) 712046; 2.四川大學(xué) 化學(xué)學(xué)院,四川 成都 610064)

青蒿素及其衍生物抗瘧活性的密度泛函理論研究

劉靖麗1*,郭 勇2

(1.陜西中醫(yī)藥大學(xué) 藥學(xué)院,陜西 咸陽(yáng) 712046; 2.四川大學(xué) 化學(xué)學(xué)院,四川 成都 610064)

為了深入研究青蒿素及其衍生物的抗瘧活性和分子結(jié)構(gòu)之間的關(guān)系，運(yùn)用密度泛函理論B3LYP方法，以6-31G*為基組對(duì)青蒿素及其衍生物二氫青蒿素、蒿甲醚和青蒿琥酯進(jìn)行了優(yōu)化計(jì)算.從分子的幾何構(gòu)型、NBO電荷及前線軌道能等方面分析了青蒿素及其衍生物的抗瘧活性與結(jié)構(gòu)之間的關(guān)系.青蒿素及其衍生物結(jié)構(gòu)中的過(guò)氧橋鍵是其抗瘧作用的活性位，O17和O20帶負(fù)電荷越多、ΔE(LUMO-HOMO)越低、EHOMO能級(jí)越高，分子的抗瘧活性越強(qiáng).結(jié)果表明，4種化合物的抗瘧活性順序?yàn)椋呵噍锼?二氫青蒿素<蒿甲醚<青蒿琥酯，與臨床實(shí)驗(yàn)結(jié)果相吻合.

青蒿素；抗瘧活性；密度泛函理論

青蒿為菊科艾屬一年生草本植物黃花蒿，性寒，味苦、辛，歸肝、膽經(jīng)，具有清熱解暑、除蒸、截瘧之功效,青蒿自古以來(lái)就被用于發(fā)熱和瘧疾的治療[1].上世紀(jì)70年代，屠呦呦等從青蒿中分離出一種具有抗瘧活性的倍半萜內(nèi)酯化合物，命名為青蒿素(artemisinin).通過(guò)化學(xué)反應(yīng)、光譜數(shù)據(jù)和X射線單晶衍射方法證明青蒿素結(jié)構(gòu)中含有獨(dú)特的過(guò)氧橋鍵[2-3].后來(lái)人們對(duì)青蒿素化學(xué)結(jié)構(gòu)進(jìn)行了改造，以青蒿素母體為基礎(chǔ)，在青蒿素的碳氧雙鍵上加上不同基團(tuán)，合成了二氫青蒿素(dihydroartemisinin)、蒿甲醚(artemether)和青蒿琥酯(artesunate)等衍生物[4-6].

實(shí)驗(yàn)表明[4]，二氫青蒿素的活性比青蒿素高1倍，但不穩(wěn)定，易分解.脂溶性的蒿甲醚也具有抗瘧活性，其活性高于青蒿素，現(xiàn)已應(yīng)用于臨床.青蒿琥酯是目前唯一有效的青蒿素水溶性衍生物，是經(jīng)我國(guó)衛(wèi)生部批準(zhǔn)上市用于抗瘧的一類(lèi)新藥[6].近年來(lái)的相關(guān)研究表明[7-9]，青蒿素除具有抗瘧疾作用外，在心血管疾病、抗腫瘤、免疫調(diào)節(jié)等多方面具重要的藥理作用和潛在的應(yīng)用價(jià)值.但是有關(guān)其結(jié)構(gòu)和活性關(guān)系的理論研究甚少，本文中我們利用密度泛函(DFT)方法對(duì)上述4個(gè)目標(biāo)化合物進(jìn)行研究，試圖從微觀層次上闡明其結(jié)構(gòu)與活性之間的關(guān)系，為設(shè)計(jì)和合成抗瘧活性更強(qiáng)的青蒿素類(lèi)藥物分子提供有用的理論信息.

1 計(jì)算方法

運(yùn)用DFT(B3LYP)方法，以6-31G*為基組對(duì)4個(gè)目標(biāo)化合物青蒿素(1)、二氫青蒿素(2)、蒿甲醚(3)和青蒿琥酯(4)的幾何構(gòu)型進(jìn)行全優(yōu)化，得到4個(gè)目標(biāo)化合物的最優(yōu)幾何結(jié)構(gòu)，并在此基礎(chǔ)上進(jìn)行頻率計(jì)算.結(jié)果顯示，優(yōu)化得到的結(jié)構(gòu)均無(wú)虛頻，說(shuō)明我們得到的是能量極小結(jié)構(gòu).最后進(jìn)行NBO(Natural Bond Orbital)電荷和分子前線軌道分析，所有的計(jì)算均由Gaussian 03程序完成[10].

2 結(jié)果與討論

2.1 分子的幾何平衡構(gòu)型

采用B3LYP/6-31G*方法對(duì)4個(gè)目標(biāo)化合物的結(jié)構(gòu)進(jìn)行了幾何全優(yōu)化，獲得了各個(gè)化合物的幾何平衡構(gòu)型，如圖1所示.由此結(jié)構(gòu)計(jì)算所得到的頻率中沒(méi)有虛頻，說(shuō)明理論計(jì)算幾何結(jié)構(gòu)為穩(wěn)定結(jié)構(gòu)，表1列出了計(jì)算所得的部分幾何結(jié)構(gòu)參數(shù)，并與青蒿素的晶體結(jié)構(gòu)數(shù)據(jù)[3]進(jìn)行了對(duì)比.

圖1 青蒿素及其衍生物的結(jié)構(gòu)和原子編號(hào)Fig.1 The geometry and atom numbering of artemisinin and its derivatives

SpeciesExp[3]1234SpeciesExp[3]1234C4-O30.1390.1380.1430.1430.142O20-C5-C6106.2106.4105.6105.6105.5C2-O30.1440.1440.1430.1430.143O20-C5-C8104.4104.7103.3103.4103.4C2-O170.1400.1410.1420.1420.142O20-C5-C4111.0110.1112.2112.0112.2O17-O200.1480.1460.1460.1460.146C2-O17-O20-C547.446.143.243.244.0C5-O200.1460.1460.1460.1460.146C(5-4-14-12)-24.534.935.836.1C14-O150.1200.1210.1400.1390.143C(1-7-6-5)--35.6-39.7-39.4-38.9C14-O130.1340.1360.1440.1440.142C(5-8-10-11)--0.051.11.11.2O13-C40.1450.1440.1400.1400.141C(6-11-10-9)--53.0-53.8-53.7-53.9C10-C2-O17112.3111.8111.4111.4111.4C(6-5-8-9)-52.7951.351.351.4C2-O17-O20107.5108.8108.8108.9108.9O13-C14-C8-C5--31.3-52.0-51.8-53.1O17-O20-C5111.1111.9113.2113.1112.9C4-O13-C14-C12-14.465.064.966.7

從表1中可以看出，采用B3LYP/6-31G*方法計(jì)算所得4個(gè)目標(biāo)化合物的結(jié)構(gòu)參數(shù)與其母體化合物青蒿素的實(shí)驗(yàn)值[3]非常接近.由分子價(jià)鍵理論可知，鍵長(zhǎng)越長(zhǎng)，鍵能越小，鍵越容易斷裂，反應(yīng)活性越高.可見(jiàn)分子的鍵強(qiáng)度與鍵長(zhǎng)有密切關(guān)系，所以我們可以通過(guò)判斷鍵長(zhǎng)來(lái)推斷反應(yīng)活性.正常的過(guò)氧鍵長(zhǎng)為0.13 nm左右[11].4個(gè)目標(biāo)化合物的過(guò)氧鍵長(zhǎng)均大于0.13 nm，說(shuō)明4個(gè)化合物中的過(guò)氧鍵比較容易斷裂，從一定程度上反映出青蒿素及其衍生物的抗瘧活性與過(guò)氧鍵有密切的關(guān)系.另外從二面角的數(shù)值可以看出，包括過(guò)氧鍵的七元環(huán)存在不同程度的扭曲，這種扭曲顯然是由環(huán)與環(huán)之間不同程度的排斥所導(dǎo)致的.環(huán)己烷是一個(gè)椅式的穩(wěn)定構(gòu)象，六元環(huán)內(nèi)酯是一個(gè)扭曲的椅式構(gòu)象，理論計(jì)算結(jié)果與其母體青蒿素的晶體結(jié)構(gòu)相吻合[3].

2.2 NBO電荷布局?jǐn)?shù)

在藥理學(xué)中，當(dāng)藥物分子與受體相互作用時(shí)，必須通過(guò)兩者之間電子的相互作用來(lái)進(jìn)行.因此，藥物分子的電荷分布對(duì)藥物活性有重要的影響.分析藥物分子內(nèi)電荷分布可以揭示其與受體的作用位點(diǎn).表2列出了4個(gè)目標(biāo)化合物中起主要作用的C和O原子的NBO電荷.

表2 C原子和O原子的NBO電荷Table 2 NBO charge on the carbon atoms and oxygen atoms

從表2中可以看出，4個(gè)目標(biāo)化合物的電荷分布規(guī)律極為相似，除了青蒿琥酯.負(fù)電荷集中分布在5個(gè)O原子上，而正電荷則主要分布在與O原子相連的C上，由此分析O原子主要與受體的正電性區(qū)域相作用，而C14、C2、C4及C5則與受體的負(fù)電性區(qū)域相作用.對(duì)青蒿素分子結(jié)構(gòu)與其抗瘧活性的構(gòu)效關(guān)系研究表明過(guò)氧鍵是分子中的關(guān)鍵活性部位[12]，有活性的青蒿素類(lèi)似物，如果把過(guò)氧基-O-O-變成-O-、-O-C-O-都將失去活性[13].隨后進(jìn)一步研究證實(shí)，青蒿素與瘧原蟲(chóng)體內(nèi)的血紅素作用時(shí)，中心原子Fe2+的空軌道可結(jié)合富電子的原子和基團(tuán)[14].SHUKLA等[15]利用分子模擬方法將青蒿素與血紅素進(jìn)行對(duì)接研究，結(jié)果提示青蒿素的分子內(nèi)過(guò)氧橋結(jié)構(gòu)緊靠于血紅素的中心鐵原子上.因此，青蒿素在發(fā)揮抗瘧活性時(shí)，主要以供電子為主，過(guò)氧橋鍵是此類(lèi)分子的主要供電子區(qū)，過(guò)氧橋鍵的給電子能力直接影響化合物的抗瘧活性.這一結(jié)論與程永浩等[16]通過(guò)半經(jīng)驗(yàn)的MNDO方法研究青蒿素的電子結(jié)構(gòu)與生物活性關(guān)系的計(jì)算結(jié)果一致.由此可見(jiàn)，過(guò)氧原子O17和O20所帶的負(fù)電荷越多，青蒿素與血紅素的相互作用越強(qiáng)，分子的抗瘧活性越強(qiáng).根據(jù)表2中的電荷分布推測(cè)4個(gè)目標(biāo)化合物的抗瘧活性順序可能為：青蒿素<二氫青蒿素≤蒿甲醚<青蒿琥酯.電荷分布分析不夠充分，還需要考慮分子的前線軌道.

2.3 前線分子軌道及軌道能

根據(jù)前線分子軌道理論，最高的電子占據(jù)軌道(HOMO)和最低的電子未占軌道(LUMO)在化學(xué)反應(yīng)中起著主要作用，前線軌道往往與分子的生物活性有著直接的聯(lián)系.因此，研究前線軌道的性質(zhì)可以為確定活性部分以及探討作用機(jī)制等提供重要信息.HOMO具有優(yōu)先提供電子的重要作用，HOMO能級(jí)越高，該物質(zhì)越易失去電子參與反應(yīng)；LUMO具有優(yōu)先接受電子的重要作用，LUMO能級(jí)越低，該物質(zhì)越易得到電子參與反應(yīng)[17-19].圖2給出了4個(gè)目標(biāo)化合物的HOMO軌道和LUMO軌道圖，由圖2可以看出4個(gè)化合物的HOMO成份主要集中在過(guò)氧橋鍵O17和O20原子上，另外C5、C2、C4、O15、O13原子也有貢獻(xiàn).進(jìn)一步說(shuō)明，青蒿素及其衍生物二氫青蒿素、蒿甲醚和青蒿琥酯在發(fā)揮抗瘧活性時(shí)，均以供電子為主，過(guò)氧橋鍵是決定其抗瘧活性的關(guān)鍵幾何結(jié)構(gòu).

表3分別列出了4個(gè)目標(biāo)化合物的前線軌道能級(jí)數(shù)據(jù)，表中ΔE(LUMO-HOMO)為最高占據(jù)軌道能量EHOMO與最低未占據(jù)軌道能量ELUMO的能級(jí)差.ΔE(LUMO-HOMO)的大小反映了電子從最高占據(jù)軌道向最低空軌道躍遷需要的能量大小，能級(jí)差越小，電子越容易被激發(fā)，與電子得失相關(guān)的反應(yīng)活性越強(qiáng).對(duì)于4個(gè)目標(biāo)化合物來(lái)說(shuō)，ΔE(LUMO-HOMO)越小，過(guò)氧橋鍵的氧原子越容易給出電子，抗瘧活性越強(qiáng).由表3ΔE的大小可以得出4個(gè)化合物的抗瘧活性順序?yàn)椋呵噍锼?二氫青蒿素<蒿甲醚<青蒿琥酯，與前面電荷分布的分析結(jié)果一致.同時(shí)EHOMO能級(jí)越高，活性越高，由此得出的活性順序?yàn)椋呵噍锼?二氫青蒿素<蒿甲醚<青蒿琥酯，與前面得到的活性順序符合.無(wú)論是ΔE還是EHOMO，都表明4個(gè)目標(biāo)化合物中，青蒿琥酯的活性是最強(qiáng)的，而二氫青蒿素與蒿甲醚的活性均高于青蒿素的抗瘧活性.

圖2 4個(gè)化合物的HOMO(a)與LUMO(b)軌道圖Fig.2 HOMO (a) and LUMO (b) orbitals of four compounds

表3 4個(gè)化合物的最高占據(jù)軌道能和最低空軌道能以及兩者之差(單位:eV)Table 3 Energies of HOMO and LUMO and the energy gap between HOMO and LUMO of four compounds (in eV)

3 結(jié)論

采用密度泛函理論B3LYP方法對(duì)青蒿素及其衍生物(二氫青蒿素、蒿甲醚和青蒿琥酯)分子進(jìn)行了優(yōu)化計(jì)算，得到了4個(gè)化合物的理論平衡幾何結(jié)構(gòu)，與其母體青蒿素的晶體結(jié)構(gòu)吻合的較好.青蒿素及其衍生物的抗瘧活性與其結(jié)構(gòu)中的過(guò)氧鍵密切相關(guān)，從分子的幾何結(jié)構(gòu)、NBO電荷和前線軌道等方面分析了4個(gè)化合物的抗瘧活性與結(jié)構(gòu)的關(guān)系，結(jié)果表明，過(guò)氧橋鍵是其發(fā)揮抗瘧活性的關(guān)鍵幾何結(jié)構(gòu).從電子結(jié)構(gòu)的角度對(duì)青蒿素及其衍生物抗瘧活性的差異作出了合理的解釋?zhuān)瑸榍噍锼丶捌溲苌锏慕Y(jié)構(gòu)改造提供了理論指導(dǎo).

[1] 屠呦呦.青蒿及青蒿素類(lèi)藥物[M].北京:化學(xué)工業(yè)出版社,2009:5.

TU Y Y.Artemisinin and artemisinin-based drugs [M].Beijing:Chemical Industry Press,2009:5.

[2] 青蒿結(jié)構(gòu)研究協(xié)作組.一種新型的倍半萜內(nèi)酯-青蒿素[J].科學(xué)通報(bào),1977,22(3):123.

Research Group of Artemisinin Structure.A new kind of half terpene lactones-artemisinin [J].Chinese Science Bulletin,1977,22(3):123.

[3] 中國(guó)科學(xué)院生物物理研究所青蒿素協(xié)作組.青蒿素的晶體結(jié)構(gòu)及其絕對(duì)構(gòu)型[J].中國(guó)科學(xué),1979,11:1114-1128.

Institute of Biophysics,Chinese Academy of Sciences,Artemisinin Group.The crystal structure of artemisinin and its absolute configuration [J].Science China,1979,11:1114-1128.

[4] 李英,虞佩琳,陳一心,等.青蒿素衍生物的合成[J].科學(xué)通報(bào),1979,24(14):667-669.

LI Y,YU P L,CHEN Y X,et al.Synthesis of some derivatives of artemisinin [J].Chinese Science Bulletin,1979,24(14):667-669.

[5] 梁潔,李英.青蒿素芳香醚類(lèi)衍生物的合成[J].中國(guó)藥物化學(xué)雜志,1996,19:22-25.

LIANG J,LI Y.Synthesis of the arylether derivatives of artemisinin [J].Chinese Journal of Medicinal Chemistry,1996,19:22-25.

[6] 劉旭.青蒿素衍生物的研究[J].藥學(xué)通報(bào),1980,15(4):183.

LIU X.Research of artemisinin derivatives [J].Chinese Pharmaceutical Bulletin,1980,15(4):183.

[7] 郭燕,王俊,陳正堂.青蒿素類(lèi)藥物的藥理作用新進(jìn)展[J].中國(guó)臨床藥理學(xué)與治療學(xué),2006,11(6):615-620.

GUO Y,WANG J,CHEN Z T.Recent advancement in pharmacological effects of artemisinin and its derivatives [J].Chinese Journal of Clinical Pharmacology and Therapeutics,2006,11(6):615-620.

[8] 李琛琛,尹昆,閆歌.抗瘧藥青蒿素及其衍生物相關(guān)藥理作用研究進(jìn)展[J].中國(guó)病原生物雜志,2016,11(2):185-187.

LI C C,YIN K,YAN G.Advances in research on the pharmacological action of the antimalarial artemisinin and its derivatives [J].Journal of Pathogen Biology,2016,11(2):185-187.

[9] 張楠.青蒿素類(lèi)藥物的主要研究進(jìn)展[J].中國(guó)藥物評(píng)價(jià),2013,30(1):13-16.

ZHANG N.Research progress of the artemisinin drugs [J].Chinese Journal of Drug Evaluation,2013,30(1):13-16.

[10] FRISCH M J,TRUCKS G W,SCHLEGEL H B,et al.Gaussian 03 [CP].Revision A.01,Wallingford CT:Gaussian,Inc.,2004.

[11] HNYK D,MACHACEK J,ARGUELLO G A,et al.Stucture and conformational properties of bis(trifluoromethyl) peroxydicarbonate,CF3OC(O)O-OC(O)OCF3[J].The Journal of Physical Chemistry A,2003,107(6):847-851.

[12] 陳揚(yáng),朱世民,陳洪淵.青蒿素類(lèi)抗瘧藥研究進(jìn)展[J].藥學(xué)學(xué)報(bào),1998,33(3):234-239.

CHEN Y,ZHU S M,CHEN H Y.The research progress of artemisinin-based antimalarial drugs [J].Acta Pharmaceutica Sinica,1998,33(3):234-239.

[13] BUTLER A R,WU Y L.Artemisinin (qinghaosu):a new type of antimalarial drug [J].Chemical Society Reviews,1992,21(2):85-90.

[14] GOLENSER J,WAKNINE J H,KRUGLIAK M,et al.Current perspectives on the mechanism of action of artemisinins [J].International Journal for Parasitology,2006,36(14):1427-1441.

[15] SHUKLA K L,GUND T M,MESHNICK S R.Molecular modeling studies of the artemisinin (qinghaosu)-hemin interaction:Docking between the antimalarial agent and its putative receptor [J].Journal of Molecular Graphics,1995,13(4):215-222.

[16] 程永浩,李曉艷,鄭世均.青蒿素及其衍生物的量子化學(xué)及構(gòu)效關(guān)系研究[J].計(jì)算機(jī)與應(yīng)用化學(xué),2000,17(21):21-22.

CHENG Y H,LI X Y,ZHENG S J.Quantum chemistry and structure-activity relationship studies of artemisinin and its derivatives [J].Computers and Applied Chemistry,2000,17(21):21-22.

[17] 張?chǎng)?楊英杰,呂慶章.4種甘草黃酮類(lèi)化合物抗氧化活性的密度泛函理論研究[J].計(jì)算機(jī)與應(yīng)用化學(xué),2012,29(6):656-660.

ZHANG X,YANG Y J,Lü Q Z.Density functional theory calculations on antioxidation activity of four flavones from radix glycyrrhizae [J].Computers and Applied Chemistry,2012,29(6):656-660.

[18] 張?chǎng)?楊英杰,呂慶章.黃芪異黃酮類(lèi)化合物抗氧化活性的密度泛函理論研究[J].化學(xué)研究與應(yīng)用,2012,24(11):1662-1669.

ZHANG X,YANG Y J,Lü Q Z.Density functional theory calculations on antioxidation activity of the isoflavone compounds from astragalus [J].Chemical Research and Application,2012,24(11):1662-1669.

[19] 陳瑩,徐抗震,宋紀(jì)蓉,等.酚酸抗氧化活性的理論計(jì)算[J].食品科學(xué),2011,32(9):36-39.

CHEN Y,XU K Z,SONG J R,et al.Theoretical study on antioxidant activity of phenolic acids [J].Food Science,2011,32(9):36-39.

[責(zé)任編輯：吳文鵬]

Density functional theory studies on antimalarial activity of artemisinin and its derivatives

LIU Jingli1*,GUO Yong2

(1.CollegeofPharmacy,ShaanxiUnivserityofChineseMedicine,Xianyang712046,Shaanxi,China; 2.CollegeofChemistry,SichuanUniversity,Chengdu,610064,Sichuan,China)

In order to study the relationship between the structure and antimalarial activity of artemisinin and its derivatives,the density functional theory (DFT) B3LYP method with 6-31G*as a basis set was used to optimize the structure of artemisinin,dihydroartemisinin,artemether and artesunate.The properties of the artemisinin and its derivatives have been discussed in detail based on their molecular structures,NBO charges and the energies of the frontier molecular orbitals.The results showed that the endoperoxy bridge is likely the most active site and essential for antimalarial activity of artemisinin and its derivatives.The more negative charge of O17 and O20,the lower enengy of ΔE(LUMO-HOMO)and the higher energy ofEHOMOresulted in the stronger antimalarial activity of the compounds.It was revealed that the antimalarial activity of the compounds goes along with the order:artemisinin

artemisinin; antimalarial activity; density functional theory (DFT)

2016-09-08.

陜西中醫(yī)藥大學(xué)重點(diǎn)培育項(xiàng)目(2015PY11).

劉靖麗(1984-),女,講師,研究方向?yàn)榱孔踊瘜W(xué).*

,E-mail:sculiujl@126.com.

O641

A

1008-1011(2017)01-0087-05