孫 穎 任愛民 閔春剛 鄒陸一 任雪峰

(吉林大學理論化學研究所,理論化學計算國家重點實驗室,長春 130023)

海螢熒光素衍生物發(fā)光反應關鍵步驟的理論研究

孫 穎 任愛民*閔春剛 鄒陸一 任雪峰

(吉林大學理論化學研究所,理論化學計算國家重點實驗室,長春 130023)

用不同取代基對化學發(fā)光物質(zhì)6-芳基-2-甲基咪唑[1,2-α]吡嗪-3-(7H)酮環(huán)(海螢發(fā)光的類似物)的芳基位進行取代,形成系列海螢熒光素類似物MIPa-MIPd;并采用B3LYP/6-311+G(d,p)方法,通過電子抽取能(EEP)的計算和自然電荷布居分析(NPA),研究了在氣相、二甲亞砜(DMSO)和二甘醇二甲醚(DG)中海螢發(fā)光的類似物從陰離子變化到自由基過程中取代基的作用.結果表明:在這個過程中,吲哚在DG中作為取代基(MIPb)時的EEP最小,電荷變化最大,說明這種取代基有利于反應的進行.

密度泛函理論; 海螢熒光素; 化學發(fā)光; 電子抽取能; 自然電荷布居分析

海螢由于高的發(fā)光效率而成為研究熱點.海螢把熒光素和熒光素酶噴射到海水中產(chǎn)生藍色熒光(發(fā)光效率0.3),海螢熒光素的晶體結構在1957年由Shimomura等[1]測出.由于生物發(fā)光的靈敏度很高,所以海熒熒光素的研究對生物分子的快速敏感檢測和對功能活性分子以及細胞的活動規(guī)律進行原位跟蹤觀察發(fā)揮著非常重要的作用[2].海熒熒光素在生物醫(yī)學研究[3]、臨床診斷、食品[4]以及環(huán)境檢測[5-6]等方面的應用也越來越廣泛,相關產(chǎn)品也逐步實現(xiàn)產(chǎn)業(yè)化和規(guī)?；痆7].并且海螢熒光素及其熒光素酶都很穩(wěn)定,可保存多年,這可用于某些特殊情況下的照明[8].但是海螢熒光素高效率的發(fā)光問題至今還沒有清楚的理論描述.海螢生物發(fā)光通常的反應機理是海螢在體內(nèi)熒光素酶和氧氣等的作用下,發(fā)光分子(通常稱為熒光素)被氧化成激發(fā)態(tài),在退激的過程中發(fā)出可見光,所以生物發(fā)光本質(zhì)也是化學發(fā)光,不同之處在于,在生物發(fā)光中激發(fā)態(tài)的形成往往伴隨酶的作用[7,9].1966年研究者發(fā)現(xiàn),一個咪唑并吡嗪酮環(huán)的存在是發(fā)光反應的核心,另外海螢熒光素的同系物也被合成出來,并用來研究海螢生物發(fā)光的機理[9-13].Takahashi等[14]對海螢生物發(fā)光的研究表明,咪唑并吡嗪酮環(huán)的陰離子(MIP-)和與氧氣反應生成咪唑并吡嗪酮環(huán)的自由基(MIP·)的過程(如示意圖1所示)對整個發(fā)光反應起到了關鍵的作用,并且決定了反應的速率[14-15].本文以海螢熒光素活性中心分子結構咪唑并吡嗪酮環(huán)為核心,在C6位上用不同的芳基取代(見圖1),用密度泛函理論(DFT)[16-17],B3LYP方法[18-20]結合6-311+G(d,p)基組從理論上研究這些海螢熒光素衍生物的熒光活性中心分子的電子結構,從而從分子水平上理解C6位上取代基的變化對該發(fā)光反應決速步驟的影響.

1 理論方法

密度泛函理論研究化學發(fā)光現(xiàn)象的可靠性已得到廣泛證實[21-22].本文運用Gaussian 03程序包[23],采用密度泛函理論DFT中的B3LYP泛函方法對類海螢熒光素發(fā)光過程中關鍵步驟(示意圖1中的反應式(1))涉及的陰離子、自由基的結構進行了優(yōu)化,考慮所研究分子6-芳基-2-甲基咪唑[1,2-α]吡嗪-3-(7H)酮環(huán)上7位H原子的重要性和相應自由基的活潑性,計算中采用了6-311+G(d,p)基組.文中對6-芳基-2-甲基咪唑[1,2-α]吡嗪-3-(7H)酮環(huán)(海螢發(fā)光的相似物)芳基位不同取代基分子(MIP)從陰離子到自由基的變化過程(示意圖1中的反應式(1))的電子抽取能(EEP)和自然電荷布居進行了計算研究.在這里電子抽取能(EEP)定義為陰離子失去一個電子形成自由基所需的垂直電離能.運用自然鍵軌道(NBO)分析程序得到所計算分子中的原子自然電荷布居,分析各種分子內(nèi)、分子間電荷變化.溶劑作用選用極化連續(xù)模型(PCM)[24]來模擬.溶劑二甲亞楓(DMSO)和二甘醇二甲醚(DG)的相對介電常數(shù)分別采用46.7和6.25.

2 結果與討論

本文研究的化學發(fā)光反應過程(見示意圖1)即MIP在堿性溶劑作用下失去 H+離子形成 MIP-, MIP-在氧氣等的作用下,發(fā)生電子轉移,即O2抽取了一個電子形成O-·20,使陰離子MIP-變成自由基MIP·,它們繼續(xù)作用產(chǎn)生了發(fā)光分子(DT),發(fā)光分子被氧化成激發(fā)態(tài)(1AAP*),1AAP*在退激的過程中發(fā)出可見光.其中MIP-→MIP·過程為整個發(fā)光反應的決速步驟[14].對這個決速步驟有影響的因素除了骨架分子結構咪唑并吡嗪酮環(huán)外,取代基的影響是值得研究的.已有的報道結果[12-13]顯示:在C6位上給電子基團取代咪唑并吡嗪酮環(huán)有較好的熒光效率,所以我們選擇了三種給電子基團(如圖1所示的a-c)和一種吸電子基團(d)對咪唑并吡嗪中心環(huán)的C6位進行取代,形成四種反應物MIPa,MIPb,MIPc, MIPd.

2.1 分子幾何優(yōu)化

示意圖 1給出了MIP→MIP-→MIP·變化過程[14-15,25-27]的示意圖,圖1給出了C6位上不同芳基取代的MIP衍生物的具體分子結構,帶有不同取代基的MIP分子陰離子和自由基優(yōu)化結構如圖2所示.

表1 MIPa-MIPd陰離子和自由基的C—C6、C═O鍵長及其差值Table 1 C—C6,C═O bond lengths and bond length differences of the anions and free radicals for compounds MIPa-MIPd

把不同取代基取代咪唑并吡嗪酮化物的中性自由基分子和陰離子作對比,結果表明:在氣相中,給受電子基團與MIP間基本上保持在同一平面(二面角在179.15°-179.90°之間,而吸電子基團為173.85°),取代基與MIP連接處的C—C6鍵長有明顯變化.有研究表明[14],當咪唑并吡嗪酮環(huán)在溶劑DG中時,其發(fā)光效率接近于海螢生物體本身,所以我們在前面優(yōu)化的幾何構型下,采用PCM模型在DG溶劑中再次優(yōu)化其構型.部分優(yōu)化結果如表1所示,結果表明:在MIP-變?yōu)镸IP·時吸電子取代基團與咪唑并吡嗪酮環(huán)的角度發(fā)生的變化相對給電子基團偏大,因而吸電子基團取代時共軛減弱.陰離子與自由基中的C—C6鍵長差值的關系是MIPb=MIPc≈MIPa＞MIPd,說明取代基為給電子性質(zhì)時,從陰離子到自由基時C—C6鍵長變化比較大并且都縮短了(氣相中縮短了0.0009-0.0010 nm,DG中縮短了0.0012-0.0013 nm),表明此時的自由基相對陰離子來說,取代基與咪唑并吡嗪酮環(huán)間共軛加強,而由吸電子基取代時C—C6鍵長也縮短了,但是幅度比較小(氣相中縮短了0.0005 nm,DG中縮短了0.0006 nm).這樣無論在氣相還是溶劑中都呈現(xiàn)一個規(guī)律:具有給電子基的陰離子a,b,c相對吸電子基的陰離子d鍵長變化大,給電子基團與中心結構的π共軛相互作用加強,并且在溶劑中C—C6鍵長變化幅度都大于在氣相中,說明溶劑增強了給電子基團與中心結構片段之間的相互作用.

2.2 電子抽取能的計算

由Takahashi等[14]的研究可知,咪唑并吡嗪酮環(huán)的陰離子(MIP-)與氧氣反應生成咪唑并吡嗪酮環(huán)的自由基(MIP·)的過程(如示意圖1所示)對整個發(fā)光反應起到了關鍵的作用,并且決定了反應的速率[15].電子抽取能EEP可評估一個分子失去電子的能力,EEP越小說明失去電子時需要的能量越小,即越容易失去電子.對于海螢衍生物熒光過程、性質(zhì)和熒光效率來說,除了發(fā)光體本身分子結構起決定作用外,外界環(huán)境如溶劑的影響也不容忽視.已有實驗報道[27-28]表明在非質(zhì)子溶劑中海螢類似物的熒光效率較高,所以,我們選擇了在DMSO和DG非質(zhì)子溶劑中及氣相中對MIPa-MIPd分子進行深入的理論研究.首先我們對反應物 MIPa--MIPd-的HOMO-LUMO能隙做了分析,結果列于表2.結果表明,它們的能隙在氣相和DMSO中的順序為MIPa-＜MIPd-＜MIPc-＜MIPb-,在溶劑 DG中變?yōu)镸IPa-＜MIPb-＜MIPc-＜MIPd-,顯然MIPb-在DG中的動力學穩(wěn)定性降低,說明在類海螢環(huán)境DG中MIPb-在整個發(fā)光反應的決速步驟中優(yōu)先于其他取代物MIPc--MIPd-.

其次分別計算了MIPa-MIPd分子的電子抽取能(表3),通過表3可知:無論是在氣相中還是兩種溶劑(DMSO和DG)中,給電子基取代物MIPa-MIPc的EEP總小于吸電子基取代物MIPd的,EEP越小,反應越容易進行.表3結果說明給電子基團作為取代基時更有助于反應快速進行.在DMSO中的EEP總體小于在DG溶劑中的EEP,說明海螢衍生物的化學發(fā)光反應決速步驟(1)在DMSO中比在DG中更快速,這點與實驗報道[25]相一致.但在DG中,對不同取代基化合物的EEP比較發(fā)現(xiàn),吲哚(b)取代物MIPb-最容易失去一個電子轉化為自由基MIPb·,即EEP最小,反應較其他取代物更易進行,從而從理論上表明在類海螢環(huán)境DG下整個發(fā)光反應的決速步驟中吲哚取代時優(yōu)先于其他取代MIP,并且b作為MIP分子的取代基時更接近海螢的原型[29-31],因此這是海螢在DG中發(fā)光效率比較高[27-28]的主要原因之一.

表2 反應物MIPa--MIPd-在不同介質(zhì)中的HOMOLUMO能隙(ΔEH-L)Table 2 HOMO-LUMO energy gaps(ΔEH-L)for compounds MIPa--MIPd-in different media

表3 MIPa-MIPd分子在不同介質(zhì)中的電子抽取能(EEP)Table 3 Electron extraction potentials(EEP)for compounds MIPa-MIPd in different media

為了驗證在不同溶劑中,對于反應過程MIP-→MIP·,咪唑并吡嗪酮環(huán)C6位上不同取代基Ar對決速步驟(1)的貢獻作用,我們還比較了陰離子和自由基在各自最優(yōu)構型下的能量差,結果列于表4.從表4中可以看出:取代基為a,b,c時的咪唑并吡嗪酮環(huán)的ΔE值一般都小于d,這說明給體取代時MIP-到MIP·這個過程中更容易進行.在DMSO中a,b,c取代物的ΔE值總體小于在DG中的ΔE,并且在DG中,吲哚取代(b)時能量差ΔE最小.這些結果與上面EEP計算結果一致,說明MIP自由基的弛豫能在非質(zhì)子化溶劑中較小,可忽略,用EEP值可以判斷不同取代基對這個反應的影響.

2.3 前線軌道及自然電荷布居分析

首先對四種取代基參與前線軌道的情況加以分析,所研究分子在氣相和兩種溶劑中的前線軌道用Gaussview軟件繪出并示于圖3中.

如圖3所示,從LUMO軌道看,各種取代基的陰離子,其LUMO軌道在有溶劑和無溶劑參與時顯示出特別明顯的差別,LUMO在氣相中基本由取代基組成,但在DMSO和DG中則完全由核結構MIP組成.而自由基分子的LUMO組成受溶劑影響不大,有無溶劑存在都完全由核結構MIP組成.由HOMO軌道看,氣相中各種取代基參與軌道構成較小,只有在溶劑協(xié)助下,取代基才參與最高占據(jù)軌道的組成中.a,b,c三種給電子取代基在自由基分子的HOMO中所占成分遠遠大于吸電子取代基d,顯然吸電子的苯唑呋喃基在MIP-→MIP·過程的活化作用不如其他三個給電子取代基.對于三個給電子基的活化自由基作用從軌道分布圖上無法區(qū)分,需要進行電荷布居定量分析進行比較.

表4 MIPa-MIPd陰離子和自由基最優(yōu)構型下的能量差Table 4 Energy differences between the optimized anions and free radicals for compounds MIPa-MIPd

利用自然鍵軌道(NBO)程序計算了MIPa、MIPb、MIPc分子中的原子布居數(shù),采用自然電荷布居分析(NPA)法分析每一個原子的電荷,從而得到分子片段的電荷分布.分別對海螢核心結構的陰離子和自由基做了電荷布居分析,目的是研究MIP-→MIP·過程電荷轉移的情況.

根據(jù)前面的計算,我們知道在DG中,在MIP-到MIP·這個過程中當吲哚(b)作為咪唑并吡嗪酮環(huán)C6位上的取代基時的EEP最小,其電子轉移能力最強,在整個發(fā)光反應的決速步驟中起到了關鍵作用.在DG中咪唑并吡嗪酮環(huán)上C6位上的取代基是怎樣影響中心環(huán)的電荷分布呢?我們對其進行了NBO計算分析,其結果示于圖4.

其中的片段電荷是與取代芳基連接后海螢熒光素母體部分的電荷量.由圖4可知,MIPa--MIPc-在DG中一個負電荷主要分布在中心環(huán)的三個N原子、酮基氧和甲基C上.離取代基最遠端的N1、酮基O上聚集的電荷在b取代(MIPb-)時最多.b取代時,甲基C原子及酮基O上收集負電荷最多,失去一個負電荷變?yōu)镸IP·時甲基C上獲得負電荷(0.01895e)比a,c(0.01815e,0.01726e)取代時多,使得MIP·的2位C原子形成相對于MIP-更為缺電子的狀態(tài),利于O-·20的進攻.

圖4也表現(xiàn)了MIP陰離子的靜態(tài)誘導效應較強,導致酮基氧、甲基碳、N7位上和自由基C2上的電荷分布較多;MIP自由基的共軛效應增強,主要表現(xiàn)為C═O鍵長、取代基與MIP連接處的C—C6鍵長都變短.這樣,從MIP陰離子到自由基過程中,酮基氧、N7和C2位上的電荷數(shù)變化相對其他原子要明顯,尤其是酮基氧變化更加明顯(見圖4).在氣相中從MIP陰離子到MIP自由基的過程中,MIPa-MIPc上C═O鍵長沒有變化,但是在DG溶劑中時就發(fā)現(xiàn)當吲哚作為MIP取代基團時C═O鍵長相對其他取代基團變化大很多,而且b取代相對于a, c取代時整個MIP片段的電荷轉移數(shù)最大,即b取代基對母體的影響最大,說明取代基為b時對這步反應過程起到了積極的作用.

MIP-到MIP·過程中咪唑并吡嗪酮環(huán)在DG溶劑中的電荷分布差值分別為0.88872e(a),0.90799e (b)和0.90321e(c).在溶劑DG中進行MIP-→MIP·的反應時,中心片段咪唑并吡嗪酮環(huán)(母體)的電荷差絕對值越大說明取代基給電子能力越強.因此,b和a,c相比是最好的給電子體.這個結果與電子抽取能的計算結論相一致.在MIP-→MIP·過程中,給電子基團吲哚取代基與咪唑并吡嗪酮環(huán)相連時單電子轉移效率最大,促使反應MIP-→MIP·快速進行.

3 結 論

從理論上研究了海螢衍生物即以咪唑并吡嗪酮環(huán)為骨架,在C6位上用不同的芳基進行取代形成的一系列衍生物.對其發(fā)光反應的關鍵步驟中涉及的陰離子和自由基做了幾何結構、鍵長變化、HOMOLUMO能隙分析、電子抽取能的計算和自然電荷布居分析.通過陰離子和自由基中的C—C6以及C═O鍵長分析表明,給電子基團與中心結構的π共軛相互作用都強于吸電子基團,并且溶劑效應增強了給電子基團與中心結構片段之間的相互作用.在DG中MIPb-的HOMO-LUMO能隙值變化比較大,吲哚作為咪唑并吡嗪酮環(huán)芳基(C6)位取代基時電子抽取能最小、陰離子和自由基之間的電荷差最大,表明在這種溶劑中其給電子能力增強、穩(wěn)定性降低,更容易實現(xiàn)MIP-→MIP·這個過程,這對提高海熒類似物熒光效率起到了積極的作用.我們將在接下來的工作里研究海螢發(fā)光機理的全過程.總之,本研究從理論計算上驗證了在DG這種溶劑中,咪唑并吡嗪酮環(huán)C6位的取代基為吲哚時類似海螢的原型,有較高的電荷轉移效率和較低的HOMO-LUMO能隙和電子抽取能,增加海螢類似物發(fā)光反應過程的速率,有助于產(chǎn)生較高的發(fā)光效率.

1 Shimomura,O.;Goto,T.;Hirata,Y.Chem.Soc.Jpn.,1957,30: 929

2 Yang,Y.B.;Ji,S.M.;Zhou,F.K.;Zhao,J.Z.Biosensors and Bioelectronics,2009,24:3442

3 Wu,C.;Kawasaki,K.Anal.Chem.,2007,79:1634

4 Tan,D.J.;Gu,R.X.Food and Machinery,1994,3:25 [譚道經(jīng),顧瑞霞.食品與機械,1994,3:25]

5 Hu,Z.;Yin,J.G.;Jiang,N.;Chen,Q.J.Chin.J.Vet.Sci.,2008, 28:732 [胡 哲,尹繼剛,姜 寧,陳啟軍.中國獸醫(yī)學報, 2008,28:732]

6 Lin,H.Envi.Sci.,1981,6:67 [林 海.環(huán)境科學,1981,6:67]

7 Yang,X.L.Continuing Medical Education,2007,21:56 [楊曉林.繼續(xù)醫(yī)學教育,2007,21:56]

8 Zhang,J.C.;Wang,D.Y.Modern photochemistry.Beijing: Beijing Industry Press,2006:319-320 [張建成,王奪元.現(xiàn)代光化學.北京:北京工業(yè)出版社,2006:319-320]

9 Katsunori,T.Bioorg.Chem.,2007,35:82

10 Shimomura,O.;Johnson,F.H.;Masugi,T.Nature,1969,164: 1299

11 Shimomura,O.;Johnson,F.H.Methods Enzymol.,1978,57:331

12 Kishi,Y.;Goto,T.;Hirata,Y.;Shimomura,O.;Johnson,F.H. Tetrahedron Lett.,1966,7:3427

13 Kishi,Y.;Goto,T.;Inoue,S.;Sugiura,S.;Kishimoto,H. Tetrahedron Lett.,1966,7:3445

14 Takahashi,Y.;Kondo,H.;Maki,S.;Niwa,H.;Hiroshi,I.;Takashi, H.Tetrahedron Lett.,2006,47:6057

15 Kondo,H.;Igarashi,T.;Maki,S.;Niwa,H.;Ikeda,H.;Hirano,T. Tetrahedron Lett.,2005,46:7701

16 Hohenberg,P.;Kohn,W.Phys.Rev.B,1964,136:864

17 Kohn,W.;Sham,L.J.Phys.Rev.A,1965,140:1133

18 Becke,A.D.Phys.Rev.A,1988,38:3098

19 Becke,A.D.J.Chem.Phys.,1993,98:5648

20 Lee,C.;Yang,W.;Parr,R.G.Phys.Rev.B,1988,37:785

21 Navizet,I.;Liu,Y.J.;Ferre,N.;Xiao,H.Y.;Fang,W.H.;Lindh, R.J.Am.Chem.Soc.,2010,132:706

22 Min,C.G.;Ren,A.M.;Guo,J.F.;Li,Z.W.;Zou,L.Y.;Goddard, J.D.;Feng,J.K.ChemPhysChem,2010,11:251

23 Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03. Revision B.01.Pittsburgh,PA:Gaussian Inc.,2003

24 Mennucci,B.;Tomasi,J.J.Chem.Phys.,1997,106:5151

25 Hirano,T.;Takahashi,Y.;Kondo,H.;Maki,S.;Kojima,S.; Hiroshi,I.;Haruki,N.Photochem.Photobiol.Sci.,2008,7:197

26 Okajima,T.;Tokuda,K.;Ohsaka,T.Bioelectrochem.Bioenerg., 1996,41:205

27 Goto,T.Pure Appl.Chem.,1968,17:421

28 Johnso,F.H.;Satchel,H.D.;Taylor,E.C.;Shimomura,O. Chemiluminescence and fluorescence of cypridina luciferin and of some new indole compounds in dimethylsufoxide//Bioluminescence in Progress.Princeton:Princeton University Press,1966:67

29 Goto,T.;Inoue,S.;Sugiura,S.;Nishikawa,K.;Isobe,M.;Abe,Y. Tetrahedron Lett.,1968,37:4035

30 Toya,Y.;Kayano,T.;Sato,K.;Goto,Y.Chem.Soc.Jpn.,1992, 65:2475

31 Nakamura,H.;Aizawa,M.;Takeuchi,D.;Murai,A.;Shimomura, O.Tetrahedron Lett.,2000,41:2185

Theoretical Investigation of the Key Reaction for the Chemiluminescence of Cypridina Luciferin Analogues

SUN Ying REN Ai-Min*MIN Chun-Gang ZOU Lu-Yi REN Xue-Feng
(State Key Laboratory of Theoretical and Computational Chemistry,Institute of Theoretical Chemistry,Jilin University, Changchun 130023,P.R.China)

A series of 6-aryl-2-methylimidazo[1,2-α]pyrazin-3-(7H)-ketone based derivatives MIPa-MIPd(cypridina luciferin analogues)were investigated using the B3LYP/6-311+G(d,p)method,and we modeled a 6-site replacement by different substituents.The effect of different substitutions on cypridina luciferin analogues in their transition from anions to free radicals was studied by electron extraction potentials(EEP)and natural charge population analysis(NPA) in gas phase,dimethyl sulfoxide(DMSO),and diglyme(DG).The calculated results indicated that MIPb(3-indolyl as the substituent)had a lower EEP and a larger natural charge population change occurred at the MIP moiety in DG when transitioning from the anion to its free radical than in the other solvents.We also showed that 3-indolyl,as a substituent,accelerated the luminescent reaction of MIP.

Density functional theory; Cypridina luciferin; Chemiluminescence; Electron extraction potential; Natural charge population analysis

O641

Received:June 7,2010;Revised:July 17,2010;Published on Web:September 3,2010.

*Corresponding author.Email:aimin_ren@yahoo.com;Tel:+86-431-88498567.

The project was supported by the National Natural Science Foundation of China(20673045,20973078),Start Fund for Returned Overseas People, China((2008)890),and Major State Basic Research Development Program of China(973)(2002CBN613406).

國家自然科學基金(20673045,20973078),留學回國人員啟動基金(外交司留(2008)890號)和國家重點基礎研究發(fā)展計劃項目(973)

(2002CBN613406)資助

?Editorial office of Acta Physico-Chimica Sinica