核酸堿基互變異構體的結構、穩(wěn)定性及其物理化學性質

曹國進 鄭衛(wèi)軍

(1清華大學化學系，有機光電子與分子工程教育部重點實驗室，北京100084;

2中國科學院化學研究所，分子反應動力學國家重點實驗室，北京分子科學國家實驗室，北京100190)

1 引言

核酸堿基是DNA及RNA分子的重要組成部分，它們在基因遺傳信息的傳遞方面起著主導作用.1核酸堿基在DNA及RNA分子中主要以標準異構體形式存在，但是當它們與金屬離子結合或成為激發(fā)態(tài)時，核酸堿基的標準異構體的分子內會發(fā)生質子轉移，從而導致核酸堿基的異構化.2這些異構化后的核酸堿基稀有異構體通常能夠導致堿基對的錯配對，進而形成不規(guī)則的DNA分子結構，最終很可能導致DNA或RNA分子雙螺旋發(fā)生自發(fā)突變.3例如稀有亞胺型異構體在金屬陽離子存在的條件下變得更加穩(wěn)定，這些異構體會進一步誘導有機體的突變.4，5因此，對核酸堿基的互變異構體進行系統的研究，有助于人們深入認識DNA和RNA分子的結構和性質.國際上有很多研究小組已經通過實驗和理論對核酸堿基互變異構體的結構、相對能量及其性質進行了研究.因為在DNA或RNA分子中，核酸堿基的互變異構體會受到各種環(huán)境的干擾，所以，現在很多研究小組采用獨立于DNA或RNA的分子骨架來研究核酸堿基的互變異構體.在實驗方面，人們利用X射線晶體衍射實驗技術測定了核酸堿基標準異構體的結構，并且還通過光譜技術包括紫外光譜、紅外和拉曼光譜以及微波光譜等與理論計算相結合的方法確定了部分稀有異構體的相對含量和穩(wěn)定性順序.6理論計算方面，半經驗、微擾理論以及密度泛函理論等方法經常被采用.很早之前人們就采用半經驗方法對核酸堿基互變異構體進行了研究，但其結果與實驗值偏差較大，后來人們較少采用這種方法.到目前為止密度泛函方法是研究核酸堿基最普遍采用的計算方法.van Mourik等7的研究表明密度泛函方法在研究尿嘧啶不同異構體的結構和相對穩(wěn)定性方面值得信賴.此外，基組的選取同樣會對堿基互變異構體的穩(wěn)定性順序產生一定的影響.盡管前人進行了很多研究，但是，以前的實驗和理論研究一般針對于某一特定堿基的一部分互變異構體，還有一些互變異構體沒有考慮到，研究仍不全面，對于核酸堿基的互變異構體的結構、相對穩(wěn)定性缺乏系統的歸納總結，不同研究小組之間的研究結果也存在互相矛盾的地方.因此，本文對文獻中涉及核酸堿基互變異構體的實驗和理論研究進行了綜述.在對前人研究進行歸納總結的基礎上，我們還利用密度泛函方法對核酸堿基的所有可能的互變異構體進行了系統的理論計算，并根據它們的相對能量對所有可能的互變異構體進行了排序.同時，我們將計算結果與文獻中已有的實驗值進行了比較，計算結果與實驗結果比較一致，因此，我們計算所得到的核酸堿基的互變異構體的結構、相對能量及排序是可信的，可以為人們研究核酸堿基、DNA、RNA分子提供重要參考.此外，因為核酸堿基的電子親和能、電離能、質子親和能等物理化學性質可以為生物、化學、材料等方面的研究提供重要的基礎性信息，因此我們還對它們的相關研究進行了綜述.

2 核酸堿基的互變異構體

2.1 尿嘧啶

尿嘧啶是構成RNA分子的一種重要堿基，在DNA分子中與它對應的堿基是胸腺嘧啶.國際上有多個研究小組對尿嘧啶的互變異構體開展了實驗研究.1982年，Beak和White8通過量熱實驗發(fā)現尿嘧啶的羰基-烯醇(U3)異構體比它的二羰基(U1)異構體的能量高(79.5±25.1)kJ·mol?1.1986年，Ferenczy等9通過電子衍射實驗測定了尿嘧啶的標準異構體的精確結構.1986年和1988年，Fujii和合作者10，11發(fā)現尿嘧啶的U2異構體要比U1異構體能量高40.2 kJ·mol?1.他們采用激光誘導熒光技術研究了超聲分子束噴射冷卻條件下的尿嘧啶，11發(fā)現尿嘧啶的U1和U3異構體在氣相中是共同存在的.同樣在1988年，Brown等12研究了尿嘧啶的微波譜，發(fā)現在氣相中U1異構體是最穩(wěn)定的異構體，并沒有發(fā)現其他互變異構體.此外，尿嘧啶的光電子能譜實驗、13紫外光譜實驗、14，15紅外和拉曼光譜實驗16?18也證實了U1異構體在氣相和凝聚相中為最穩(wěn)定的結構.

國際上有不少研究小組對尿嘧啶互變異構體的結構和相對能量進行了理論計算.191970年，Bodor等20利用半經驗方法中的Hückel分子軌道近似在沒有優(yōu)化結構只計算單點能的情況下得出尿嘧啶異構體的穩(wěn)定性順序為U1＞U2＞U11＞U3＞U4.1979年，Czermiński等21利用半經驗MINDO/3方法對尿嘧啶異構體結構進行了優(yōu)化并計算了能量，得到穩(wěn)定性順序為 U1＞U3＞U11＞U2＞U4＞U10.1984年，Scanlan等22利用Hartree-Fock方法和3-21G基組計算得到尿嘧啶異構體的穩(wěn)定性順序是U1＞U2＞U3＞U4＞U10＞U11.1987年，Saunders等23利用半經驗算法如AM1、MINDO/3、INDO等預測了U1異構體為最穩(wěn)定的異構體.1991年，Katritzky等24利用AM1方法對尿嘧啶的互變異構體做了勢能面計算，發(fā)現U1異構體為勢能曲線能量最低點，第二個次穩(wěn)定異構體是能量較標準異構體高39 kJ·mol?1的U2異構體.1992年，Leszczynski25利用較為近似的計算方法和較小的STO-3G基組以及高精度的MP4方法和較大的6-31G(d，p)基組對尿嘧啶互變異構體進行了計算，發(fā)現U1異構體均為最低能量異構體，而后面互變異構體的相對能量受方法和基組的影響比較大.1993年，Boughton和Pulay26運用MP4(SDQ)/6-31G**研究了尿嘧啶六個異構體的結構和相對能量，得到穩(wěn)定性順序是 U1＞U2＞U3＞U4＞U10＞U11.1994年，Estrin等27運用了局域和梯度校正的密度泛函(LDA+ZPE，NL1+ZPE，NL2+ZPE)計算方法得到相同的結果.1999年，Tian等28應用密度泛函理論及不同交換相關函數(SL，BLYP，B3LYP以及BHandH)研究了氣相中尿嘧啶的八個異構體，得到了它們的穩(wěn)定性順序為 U1＞U2＞U3＞U4＞U5＞U9＞U11＞U12.雖然之前這些研究小組利用半經驗方法或從頭算對尿嘧啶的互變異構體進行了理論計算，但是他們的計算僅涉及了最多8個尿嘧啶異構體.22，25，29?31實際上，2個氫原子還能夠結合到尿嘧啶的N1、N3、O8和O10位點，另外，O―H鍵的方向不同會產生不同的順反異構體，這樣得到尿嘧啶的異構體的總數為13個，包括了1個二羰基異構體、4個二羥基異構體、4個2-羥基-4-羰基異構體、4個2-羰基-4-羥基平面互變異構體.Kryachko等32，33用B3LYP/6-31+G(d，p)對全部13個尿嘧啶的異構體結構和勢能曲線進行了計算，發(fā)現了之前沒計算的5個互變異構體，這5個異構體的相對能量在前面Tian等27計算得到的8個異構體的相對能量之間，此外他們還計算得到了其中一些異構體之間的勢能曲線.

我們利用B3LYP/6-31++G(d，p)方法對尿嘧啶所有13個異構體幾何結構和相對能量進行了計算，在計算過程中考慮了零點能的校正，結果列于圖1中.B3LYP雜化泛函方法是一種較為經濟的方法，該方法得到的結果與Kryachko等32，33的計算結果非常接近，最優(yōu)結構參數也與實驗值較為一致.

2.2 胸腺嘧啶

胸腺嘧啶可看作尿嘧啶C5配位氫被甲基所取代的產物.有一些研究小組已經對胸腺嘧啶互變異構體進行了實驗研究.Ozeki等34早在1969年就測定了胸腺嘧啶標準異構體(T1)的晶體結構.1989年，Brown等35研究了胸腺嘧啶的微波譜，發(fā)現在氣相中胸腺嘧啶的二羰基標準異構體(T1)是最主要的異構體，這與Ozeki得到的晶體結構較為一致.1999年，Morsy等36研究了胸腺嘧啶在室溫下的熒光光譜，發(fā)現在水溶液中胸腺嘧啶標準異構體和稀有異構體的相對穩(wěn)定性強烈依賴于溶液pH值的大小.

國際上一些理論研究小組已經對胸腺嘧啶的異構體進行了計算，但是這些計算只考慮了部分異構體.Scanlan和Hillier22利用從頭算中的3-21G基組計算了其中6個異構體的相對能量，發(fā)現2，4-二羰基標準異構體(T1)是最低能量異構體，其幾何結構參數也與實驗值較為一致.實際上，胸腺嘧啶的互變異構體總共有13個，包括1個二羰基異構體、4個二羥基異構體、4個2-羥基-4-羰基異構體、4個2-羰基-4-羥基異構體.37，38雖然一些稀有異構體具有較高的相對能量，在一般條件下難以穩(wěn)定存在，但是在熱平衡條件改變情況下也可能會以較低濃度存在.Ha和Gunthard38運用HF/6-31G(d，p)和MP2/HF/6-31G(d，p)對所有13個異構體進行了結構優(yōu)化和相對能量計算，并用HF/6-31G(d，p)預測了從反式轉變?yōu)轫樖綍r的能量變化.例如2-羥基-反式-4-羰基-N1H異構體(T7)到 2-羥基-順式-4-羰基-N1H(T13)異構體轉變時的能量變化約為27.2(±3.3)kJ·mol?1.

我們利用B3LYP/6-31++G(d，p)方法對胸腺嘧啶所有13個異構體的幾何結構和相對能量進行了計算，在計算過程中同樣考慮了零點能的校正，結果列于圖2中.此能量順序與Ha等38的計算結果非常接近，最優(yōu)幾何結構與上述實驗測量值也較為一致.

2.3 胞嘧啶

圖1 尿嘧啶的互變異構體以及通過B3LYP/6-31++G(d，p)計算得到的相對能量Fig.1 All the tautomers of uracil and their relative energies calculated by B3LYP/6-31++G(d，p)

對于胞嘧啶的互變異構體，人們在實驗方面做了大量的工作.早在20世紀六七十年代時，Ueda39和Dreyfus40等發(fā)現當胞嘧啶以固相存在時以及在水溶液中，氨基-羰基-N1H異構體(C1)都是其最主要的異構體.1978年，Yu等41通過胞嘧啶的光電子能譜實驗證實在氣相中C1異構體是最主要的異構體.1988年和1989年，Szczesniak42和Nowak43等通過胞嘧啶的基質隔離紅外光譜實驗證實，在220°C高溫蒸發(fā)條件下，胞嘧啶氨基-羥基-順式異構體(C2)和異構體C1同時存在，其中C2異構體以較高濃度存在，原因可能是在基質隔離樣品中紫外光照射影響了互變異構體的相對濃度.同樣在1989年，Brown等44利用特殊噴嘴對胞嘧啶進行噴射冷卻后，研究了它的分子束微波光譜，結果發(fā)現胞嘧啶除了存在上述兩個異構體外，還包括少量的亞氨-反式-羰基異構體(C4).胞嘧啶這三個異構體的相對穩(wěn)定順序是異構體C4＜C2＜C1.2002年，Nir等45報道了胞嘧啶的共振增強多光子電離(REMPI)光譜，發(fā)現胞嘧啶在激光脫附和噴射冷卻條件下存在C1、C2或C3異構體.2005年，Tomic等46認為激光濺射條件下未能發(fā)現亞氨-羰基異構體的原因可能是因為樣品由固體粉末壓制得到，而固相中胞嘧啶僅以異構體C1存在.此外，激光濺射過程是一個高度非平衡過程.

理論方面，雖然目前有很多關于胞嘧啶的互變異構體的相對穩(wěn)定性的計算，22，27，43，47?58但是之前的計算只考慮了部分的異構體.由于這些異構體能量差別非常小，采用不同的理論方法往往會得到不同的結果，22因此，這些異構體在氣相中的穩(wěn)定性順序還存在爭議.1988年，Kwiatkowski等58在MP2/3-21G優(yōu)化基礎上采用MP2/6-31G(d)計算了異構體C1、C2和C4的單點能.它們的穩(wěn)定性順序和Brown等44的實驗結果一致.1996年，Colominas等48采用不同方法對前6個異構體在優(yōu)化結構的基礎上進行了單點能的計算，使用如下6種方法，A:HF/6-31G(d)//HF/6-31G(d);B:HF/6-311++G(d，p)//HF/6-31G(d);C:MP2/6-311++G(d，p)//HF/6-31G(d);D:MP4/6-311++G(d，p)//MP2/6-31G(d);E:B3LYP/6-31G(d)//MP2/6-31G(d);F:B3LYP/6-311++G(d，p)//MP2/6-31G(d)，得到的結果不完全相同.前4個計算得到的最低能量異構體為C2，后兩個得到的最優(yōu)異構體為標準異構體C1.1997年，Fogarasi47利用MP2/DZP或MBPT(2)/DZP方法對胞嘧啶異構體的結構進行優(yōu)化的基礎上，采用HF、MP2和CCSD方法以及不同基組計算了胞嘧啶的5個最低異構體的單點能.其中在MBPT(2)/DZP方法優(yōu)化結構的基礎上應用CCSD/DZP方法計算這些異構體的單點能，得到它們的相對 穩(wěn)定性順 序為 C5＜C1＜C3＜C4＜C2.1998 年，Kobayashi53用CCSD/cc-pvtz(-f)計算了胞嘧啶的互變異構體的相對穩(wěn)定性，發(fā)現穩(wěn)定性順序為C4＜C5＜C1＜C3＜C2.2000年，Sambrano等57利用MP2/6-31++G(d，p)方法計算了胞嘧啶的互變異構體的相對穩(wěn)定性，發(fā)現C2異構體為最穩(wěn)定結構，但是利用B3LYP/6-31++G(d，p)方法得到的結果顯示C1異構體最為穩(wěn)定.

圖2 胸腺嘧啶的互變異構體以及通過B3LYP/6-31++G(d，p)計算得到的相對能量Fig.2 All the tautomers of thymine and their relative energies calculated by B3LYP/6-31++G(d，p)

我們同樣利用B3LYP/6-31++G(d，p)方法對胞嘧啶所有14個異構體的幾何結構和相對能量進行了計算，在計算過程中也考慮了零點能的校正，結果列于圖3中.此結果與Sambrano等57的計算結果一致，最優(yōu)結構也與實驗得到的結果一致.

2.4 腺嘌呤

人們通過大量的實驗認識到腺嘌呤的互變異構體隨著實驗條件不同而有所不同.1975年，Chenon等59在二甲基亞砜(DMSO)溶液中利用腺嘌呤13C核磁共振譜對其互變異構體進行了研究，結果顯示腺嘌呤的氨基-N9H標準異構體(A1)是最為穩(wěn)定的異構體，其含量約為85%.同年，Dreyfus等60利用溫度躍變弛豫測量法(T-jump relaxation measurements)證實了在20°C水溶液中，A3異構體濃度為A1異構體的28%.1983年，Gonnella等61利用15N核磁共振譜發(fā)現A1和氨基-N7H(以下簡稱A3)異構體在水溶液中的含量分別約為80%和20%，而在DMSO溶液中A3異構體含量降到了13.5%.此外他們估算了這兩個異構體能量差約為4.2 kJ·mol?1，這只有對腺嘌呤異構體進行高精度的量子化學計算才能區(qū)分.之后Nowak等62，63在腺嘌呤的N9或N7位上進行15N同位素取代，再測定其基質隔離紅外光譜，發(fā)現A1標準異構體在低溫基質中是含量最多的異構體.2001年，Laxer等64利用15N核磁共振譜發(fā)現腺嘌呤在DMSO溶液中除了以A1和A3異構體形式存在，還少量以氨基-N3H異構體(A2)形式存在.2002年，Plützer和Kleinermanns65利用共振雙光子電離譜(R2PI)以及紅外-紫外(IR-UV)光譜證實了通過加熱蒸發(fā)方法得到的氣相腺嘌呤中存在著A1和A3異構體.2003年，Cohen等66利用腺嘌呤水溶液中的超快瞬態(tài)吸收光譜研究發(fā)現在室溫條件下腺嘌呤存在有A1和A3兩種異構體，其中A3異構體的含量為22%±4%.同時在DMSO溶液中腺嘌呤13C-1H核磁共振自旋-自旋耦合常數的測定表明A3異構體在這兩種異構體中所占比例為19%.直到2009年，Vogt等67才利用氣相電子衍射實驗測定了氣相中A1標準異構體的精確結構.

圖3 胞嘧啶的互變異構體以及通過B3LYP/6-31++G(d，p)計算得到的相對能量Fig.3 All the tautomers of cytosine and their relative energies calculated by B3LYP/6-31++G(d，p)

文獻中有很多關于腺嘌呤異構體相對能量及其結構參數的計算，采用的理論計算方法包括HF/6-31G、68B3LYP/ccpVDZ、69BP86/QZ4P、70RI-MP2/TZVPP、71MP2/aug-ccpVDZ、72CASSCF(4/4)/6-31G、68CCSD/aug-cc-pVDZ72以及 MCSCF(12，10)/DZP73等.BP86/QZ4P、70MP2/6-311++G(d，p)、74RI-MP2/TZVPP71以及MCQDPT2/MCSCF1473等計算方法表明A1異構體為最低能量異構體，A3和A2異構體能量比A1異構體能量要高29.3?33.5 kJ·mol?1，兩者之間存在一個高度為263.8 kJ·mol?1的能壘.1990年，Karplus和Wiorkiewicz-Kuczera75利用HF/4-21G對A1和A3異構體的結構進行了優(yōu)化并模擬了它們的振動譜圖，并通過與實際譜圖對比，發(fā)現多晶腺嘌呤的紅外光譜中僅有A1異構體的吸收峰，而在腺嘌呤的低溫基質隔離紅外光譜中存在A1和A3兩種異構體的吸收峰.2004年，Hanus等71使用了RIMP2/TZVPP計算了不同腺嘌呤異構體的相對能量，并得到了A3異構體和A2異構體在DMSO溶液和水溶液中都可能存在的結論.2007年，Kim等76進一步研究了A1、A3和A2腺嘌呤異構體在水溶液中互變異構化的機理.2010年，Aidas等77利用DFT/MM方法模擬了腺嘌呤的15N和13C核磁共振譜，發(fā)現A1和A3異構體在溶液中的含量分別約為83%和17%.

我們利用B3LYP/6-31++G(d，p)方法對腺嘌呤所有14個異構體幾何結構和相對能量進行了計算，在計算過程中也考慮了零點能的校正，結果列于圖4中.此計算得到的最低能量結構為腺嘌呤A1異構體，其幾何參數與腺嘌呤N9H異構體氣相電子衍射實驗測定的實驗值67較為一致.

2.5 鳥嘌呤

與胞嘧啶等核酸堿基類似，鳥嘌呤互變異構體的相對穩(wěn)定性受分子所處環(huán)境影響較大.到目前為止，國際上有很多關于鳥嘌呤互變異構體的實驗研究.早在1980年，LeBreton和合作者78研究了氣相中通過加熱方法得到的鳥嘌呤及其甲烷取代物光電子能譜，發(fā)現在氣相中羰基-N7H異構體(G1)是鳥嘌呤最穩(wěn)定的異構體.Nir等79，80在1999年和2001年分別報道了通過分子束噴射冷卻的鳥嘌呤的共振增強多光子電離(REMPI)光譜，紫外-紫外(UV-UV)和紅外-紫外雙共振譜(IR-UV double resonance spectroscopy)，發(fā)現了羰基-N9H(G2)、羰基-N7H(G1)、烯醇-N9H(G3)3個異構體.2001年，Piuzzi等81通過鳥嘌呤異構體的紫外光譜實驗，同樣證實在激光脫附及分子束噴射冷卻條件下存在著G1、G2和烯醇-N9H-反式(G4)三種鳥嘌呤異構體.在與Nir等人相似的研究中，Mons等82證實分子束噴射冷卻條件下鳥嘌呤存在G1、G2、G3以及烯醇-N7H(G5)四種異構體.這兩種研究的主要差別是在譜圖中34755 cm?1峰開始僅被指認為G3異構體，但后來作者重新指認為G5和G3異構體.2006年，Choi和Miller83基于隔離在氦納米液滴中的鳥嘌呤紅外激光光譜實驗，結合MP2方法以及6-311++G(d，p)和aug-cc-pVDZ基組，發(fā)現有G1、G2、G3和G4鳥嘌呤異構體存在.受Choi和Miller啟發(fā)，Mons等84又對之前氣相中得到的鳥嘌呤共振增強雙光子電離(R2PI)光譜重新指認為 G1、G5、羰基-N7H-亞氨-順式(G8)及其反式(G6)異構體.令人驚奇的是，盡管亞氨基異構體的能量遠遠高于標準異構體的能量，但還是被發(fā)現存在于激光脫附實驗中，很可能是由于鳥嘌呤激發(fā)態(tài)中出現了互變異構體依賴的超快弛豫現象.2009年，Zhou等85通過真空紫外(VUV)光電離技術研究了鳥嘌呤的互變異構體.他們發(fā)現熱蒸發(fā)和激光脫附兩種技術制備得到的異構體有很大的不同，前一種制備方法得到了G1、G2、G3、G4和G5鳥嘌呤異構體，后一種制備方法得到了G6和G8鳥嘌呤異構體.在理論計算方面，有一些研究小組已經對鳥嘌呤互變異構體進行了研究.1990年，Sabio等86利用半經驗AM1及Hartree-Fock方法和3-21G、STO-3G基組對鳥嘌呤的36個互變異構體進行了計算，結果顯示G2是最主要的異構體.1994年，Hobza和Sponer50運用MP2/6-31G(d，p)等方法計算得到鳥嘌呤非平面型結構比平面型結構能量約低4.2 kJ·mol?1.1999年，Ha等54利用HF/6-31G(d，p)和MP2/6-31G(d，p)方法對鳥嘌呤的36個互變異構體進行了計算，結果顯示G2和G3是兩個最穩(wěn)定的異構體.2005年，Gorb等87利用MP2和CCSD(T)方法及不同的基組對鳥嘌呤異構體進行計算，結果顯示鳥嘌呤四個最低能量異構體G1、G2、G3和G4相對能量差在4.2 kJ·mol?1以內，第五個較低能量異構體比最低能量異構體能量高12.6 kJ·mol?1以上.2006年，Shukla和Leszczynski88利用HF、B3LYP、MP2方法和6-311G(d，p)、6-311++G(d，p)基組對鳥嘌呤的28個互變異構體進行了計算，結果顯示G1是氣相中最為穩(wěn)定的異構體，G2比G1能量高約1.7 kJ·mol?1.同樣在2006年，Liang等89利用密度泛函理論中的B3LYP/6-311+G(d，p)基組對鳥嘌呤的36個互變異構體的相對能量以及部分堿基異構體相互轉化的116個過渡態(tài)進行了計算，結果與Shukla等88的較為接近，同樣顯示G1是氣相中最為穩(wěn)定的異構體，G2比G1能量高約1.7 kJ·mol?1.

圖4 腺嘌呤的互變異構體以及通過B3LYP/6-31++G(d，p)計算得到的相對能量Fig.4 All the tautomers of adenine and their relative energies calculated by B3LYP/6-31++G(d，p)

圖5 鳥嘌呤的互變異構體以及通過B3LYP/6-31++G(d，p)計算得到的相對能量Fig.5 All the tautomers of guanine and their relative energies calculated by B3LYP/6-31++G(d，p)

我們利用B3LYP/6-31++G(d，p)方法對鳥嘌呤所有36個異構體幾何結構和相對能量進行了計算，在計算過程中也考慮了零點能的校正，結果列于圖5中.結果表明G1異構體在氣相中是最低能量結構，與上述實驗結果一致.

3 電子親和能

由于核酸堿基是極性分子，并且它們的偶極矩比較大，例如尿嘧啶(～1.5×10?29C·m)，胸腺嘧啶(～1.5×10?29C·m)，胞嘧啶(～2.3×10?29C·m)，30因此電子結合到核酸堿基分子上形成的核酸堿基陰離子主要通過靜電－偶極相互作用.它們脫附電子后對應的中性分子結構與陰離子相比幾乎沒有改變.因此實驗測得的電子親和能一般是偶極束縛(dipolebound)而非價電子束縛(valence-bound)的電子親和能(EA).901996年，Desfrancois等91利用堿基與里德堡態(tài)的氙原子碰撞的實驗方法測得尿嘧啶、胸腺嘧啶和腺嘌呤的EA值分別為(54±35)、(68±20)、(12±5)meV.同年，Bowen小組92測定了尿嘧啶和胸腺嘧啶的光電子能譜，得到它們的EA值分別為(93±7)、(69±7)meV.1998年，Schiedt等90測定了經過分子束噴射冷卻的尿嘧啶、胸腺嘧啶以及胞嘧啶C1和C2異構體的光脫附-光電子能譜，結果顯示它們的EA值分別為(86±8)、(62±8)、(230±8)、(85±8)meV.理論計算方面，Dolgounitcheva等93計算得到尿嘧啶負離子的絕熱電子脫附能(ADE)和垂直脫附能(VDE)都為 57 meV.Wesolowski等94運用 B3LYP/TZ2P++方法對所有DNA及RNA堿基的電子親和能進行了計算，得到U、T、C、A和G的絕熱電子親和能為分別0.19、0.16、?0.02、?0.17和0.07 eV.Li等95采用B3LYP密度泛函方法結合多種基組進行計算得到DNA及RNA堿基的電子親和能順序為U≈T＞C＞A＞G.

上述理論和實驗上的核酸堿基電子親和能(EA)數據90?96總結于表1中.

4 電離能

核酸堿基的電離是DNA或RNA分子受到輻射損傷的初始過程，因此研究它們的電離能在生物方面有著重要的意義.97Hush和Cheung98很早就通過光電子能譜實驗測得了核酸堿基的第一垂直電離能，之后Orlov等99通過氣相中的光致電離質譜實驗測得了核酸堿基的第一絕熱電離能，他們的結果均顯示腺嘌呤和鳥嘌呤是兩個最容易被氧化的核酸堿基.Sevilla等100采用了MP2/6-31+G(d)//6-31G*方法得到的垂直電離能和絕熱電離能與實驗結果較為一致.Hutter和Clark101采用B3LYP/6-31G(d)方法對DNA堿基電離能進行了計算，它們計算的結果均比實驗值低0.20?0.34 eV.Le?o等102利用AM1半經驗方法計算了胸腺嘧啶的6個異構體、胞嘧啶的6個異構體、腺嘌呤的8個異構體以及鳥嘌呤的15個異構體的電離能.他們的結果顯示同一堿基互變異構體的電離能差別約在0.7 eV以內.

表1 DNA及RNA堿基的電子親和能(EAs)Table 1 Electron affinities(EAs)of the DNAand RNA nucleobases

表2 DNA及RNA堿基的電離能(IPs)Table 2 Ionization potentials(IPs)of the DNAand RNA nucleobases

相關核酸堿基的電離能的實驗值和理論值(文獻97中采用P3//MBPT(2)/6-311G**方法，文獻100中采用ROMP2/6-31+G(d)//ROHF/6-31G*方法，文獻101中采用B3LYP/D95*//UHF/6-31G*方法)總結于表2中.

5 質子親和能

核酸堿基的質子化反應在穩(wěn)定DNA三螺旋結構、酶催化反應以及有機體誘導突變等過程中起著重要的作用，測定它們的質子親和能可以為生物、化學等方面的研究提供重要的基礎性信息.1979年，Meot-Ner103就在脈沖式高壓質譜儀中利用氣相中的質子轉移動力學測定了胸腺嘧啶、胞嘧啶和腺嘌呤的質子親和能.Lias104在1984年對DNA堿基的質子親和能進行了較為全面的總結，該數據主要是通過測定氣相條件下質子轉移反應的平衡常數，然后再利用公式計算獲得.Greco等105利用質量分離離子動能(MIKE)譜圖重新測定了DNA堿基的質子親和能，得到了鳥嘌呤質子親和能的準確結果.理論計算方面，Colominas，48Podolyan，106Russo107以及Chandra108等分別采用MP4/6-311++G(d，p)//MP2/6-31G(d)、MP2/6-311++G(d，p)//MP2/6-31+G(d，p)、VWN/TZVP以及B3LYP/6-31+G(d，p)等方法計算得到了DNA堿基的質子親和能，他們的結果均與實驗值較為接近.相關的核酸堿基質子親和能的實驗值和理論值總結于表3中.

表3 DNA及RNA堿基的質子親和能(PAs)Table 3 Proton affinities(PAs)of the DNAand RNA nucleobases

6 總結和展望

我們對有關核酸堿基互變異構體的實驗和理論研究的進展進行了綜述.理論研究表明尿嘧啶、胸腺嘧啶均有13個互變異構體，包括1個二羰基異構體，4個二羥基異構體、4個2-羥基-4-羰基異構體、4個2-羰基-4-羥基互變異構體，而實驗研究中僅證實了2，4-二羰基、2-羥基-順式-4-羰基-N3H以及2-羰基-4-羥基、順式-N1H三種異構體形式存在.文獻中只對胞嘧啶部分互變異構體進行了計算和總結，而實際上胞嘧啶共有14個互變異構體.實驗研究已證實胞嘧啶存在著氨基-羰基-N1H(C1)、氨基-羥基-順式(C2)和亞氨-反式-羰基(C4)三種異構體.腺嘌呤共有14個互變異構體，現有實驗已證實存在著腺嘌呤N9H(A1)、氨基-N3H(A2)和氨基-N7H(A3)三種異構體.理論研究顯示鳥嘌呤共有36個互變異構體，而實驗上先后已證實鳥嘌呤存在著羰基-N7H(G1)、羰基-N9H(G2)、烯醇-N9H(G3)、烯醇-N9H-反式 (G4)、羰基-N7H-亞氨-順式(G8)及其反式(G6)異構體.所有核酸堿基互變異構體的穩(wěn)定性順序在理論計算中會強烈依賴于所用的方法和基組，在實驗中會隨著所處環(huán)境的變化而不同.總的說來，無論是理論還是實驗中，除了鳥嘌呤，其他核酸堿基最穩(wěn)定的異構體均為其標準異構體.而對于鳥嘌呤，標準異構體(G2)和羰基-N7H異構體(G1)在能量上差別非常小，幾乎是簡并的.在對前人研究進行歸納總結的基礎上，我們還利用密度泛函方法對核酸堿基所有可能的互變異構體進行了系統的理論計算，并根據它們的相對能量對所有可能的互變異構體進行了排序.同時，我們將計算結果與文獻中已有的實驗值進行了比較，計算結果與實驗值比較一致，因此，我們計算所得到的核酸堿基的互變異構體的結構、相對能量的排序是可信的，可以為人們研究核酸堿基、DNA、RNA分子提供重要參考.我們還對文獻中有關核酸堿基的電子親和能、電離能、質子親和能的研究進行了綜述.

近年來有關核酸堿基互變異構體的研究主要集中于核酸堿基互變異構體與金屬、109多聚氫化物110的相互作用以及它們的異構化、111，112質子化和去質子化、113激發(fā)態(tài)動力學114等方面.到目前為止引起核酸堿基異構化的原因以及利用不規(guī)則的DNA分子結構對納米材料進行自組裝的機理等方面還不是特別清楚.隨著實驗技術和理論計算方法的提高，有關核酸堿基的互變異構體及其所構成的DNA及RNA分子在生物醫(yī)學、材料、化學等領域的應用可能會成為人們關注的熱點.

(1) Blackburn，G.M.Nucleic Acids in Chemistry and Biology;Royal Society of Chemistry:Cambridge，2006;p 13.

(2) Wojnarowska，Z.;Paluch，M.;Wlodarczyk，P.;Dulski，M.;Wrzalik，R.;Roland，C.M.J.Phys.Chem.Lett.2012，3，2288.doi:10.1021/jz300541t

(3) Lee，H.;Popodi，E.;Tang，H.X.;Foster，P.L.Proc.Natl.Acad.Sci.U.S.A.2012，109，E2774.

(4) Lippert，B.;Gupta，D.Dalton Transactions 2009，4619.

(5) Ai，H.Q.;Chen，J.P.;Zhang，C.J.Phys.Chem.B 2012，116，13624.doi:10.1021/jp308937k

(6) Neidle，S.Principles of Nucleic Acid Structure;Elsevier:London，2008;p 1.

(7) van Mourik，T.;Danilov，V.I.;Dailidonis，V.V.;Kurita，N.;Wakabayashi，H.;Tsukamoto，T.Theor.Chem.Acc.2010，125，233.doi:10.1007/s00214-009-0630-0

(8) Beak，P.;White，J.M.J.Am.Chem.Soc.1982，104，7073.doi:10.1021/ja00389a032

(9) Ferenczy，G.;Harsányi，L.;Rozsondai，B.;Hargittai，I.J.Mol.Struct.1986，140，71.doi:10.1016/0022-2860(86)80148-X

(10) Fujii，M.;Tamura，T.;Mikami，N.;Ito，M.Chem.Phys.Lett.1986，126，583.doi:10.1016/S0009-2614(86)80178-6

(11) Tsuchiya，Y.;Tamura，T.;Fujii，M.;Ito，M.J.Phys.Chem.1988，92，1760.doi:10.1021/j100318a013

(12) Brown，R.D.;Godfrey，P.D.;McNaughton，D.;Pierlot，A.P.J.Am.Chem.Soc.1988，110，2329.doi:10.1021/ja00215a069

(13) Dougherty，D.;Wittel，K.;Meeks，J.;McGlynn，S.P.J.Am.Chem.Soc.1976，98，3815.doi:10.1021/ja00429a013

(14) Brady，B.B.;Peteanu，L.A.;Levy，D.H.Chem.Phys.Lett.1988，147，538.doi:10.1016/0009-2614(88)80264-1

(15) Kubota，M.;Kobayashi，T.J.Electron.Spectrosc.Relat.Phenom.1996，82，61.doi:10.1016/S0368-2048(96)03047-2

(16) Colarusso，P.;Zhang，K.;Guo，B.;Bernath，P.F.Chem.Phys.Lett.1997，269，39.doi:10.1016/S0009-2614(97)00245-5

(17) Szczesniak，M.;Nowak，M.J.;Rostkowska，H.;Szczepaniak，K.;Person，W.B.;Shugar，D.J.Am.Chem.Soc.1983，105，5969.doi:10.1021/ja00357a002

(18) Chin，S.;Scott，I.;Szczepani，K.;Person，W.B.J.Am.Chem.Soc.1984，106，3415.doi:10.1021/ja00324a006

(19) Basch，H.;Garmer，D.R.;Jasien，P.G.;Krauss，M.;Stevens，W.J.Chem.Phys.Lett.1989，163，514.doi:10.1016/0009-2614(89)85179-6

(20) Bodor，N.;Dewar，M.J.S.;Harget，A.J.J.Am.Chem.Soc.1970，92，2929.doi:10.1021/ja00713a001

(21) Czermiński，R.;Lesyng，B.;Pohorille，A.Int.J.Quantum Chem.1979，16，605.

(22) Scanlan，M.J.;Hillier，I.H.J.Am.Chem.Soc.1984，106，3737.doi:10.1021/ja00325a005

(23) Saunders，M.;Webb，G.;Tute，M.J.Mol.Struct.1987，158，69.doi:10.1016/0022-2860(87)80005-4

(24) Katritzky，A.R.;Karelson，M.J.Am.Chem.Soc.1991，113，1561.doi:10.1021/ja00005a017

(25) Leszczynski，J.J.Phys.Chem.1992，96，1649.doi:10.1021/j100183a029

(26) Boughton，J.W.;Pulay，P.Int.J.Quantum Chem.1993，47，49.

(27) Estrin，D.A.;Paglieri，L.;Corongiu，G.J.Phys.Chem.1994，98，5653.doi:10.1021/j100073a014

(28) Tian，S.X.;Zhang，C.F.;Zhang，Z.J.;Chen，X.J.;Xu，K.Z.Chem.Phys.1999，242，217.doi:10.1016/S0301-0104(99)00009-9

(29) Hobza，P.;Sponer，J.Chem.Rev.1999，99，3247.doi:10.1021/cr9800255

(30) Jasien，P.G.;Fitzgerald，G.J.Chem.Phys.1990，93，2554.doi:10.1063/1.458894

(31) Les，A.;Adamowicz，L.J.Phys.Chem.1990，94，7021.doi:10.1021/j100381a020

(32) Kryachko，E.S.;Nguyen，M.T.;Zeegers-Huyskens，T.J.Phys.Chem.A 2001，105，1288.doi:10.1021/jp001031j

(33) Kryachko，E.S.;Nguyen，M.T.;Zeegers-Huyskens，T.J.Phys.Chem.A 2001，105，1934.doi:10.1021/jp0019411

(34) Ozeki，K.;Sakabe，N.;Tanaka，J.Acta Crystallogr.，Sect.B:Struct.Sci.1969，25，1038.doi:10.1107/S0567740869003505

(35) Brown，R.D.;Godfrey，P.D.;McNaughton，D.;Pierlot，A.P.J.Chem.Soc.，Chem.Commun.1989，37.

(36) Morsy，M.;Al-Somali，A.;Suwaiyan，A.J.Phys.Chem.B 1999，103，11205.doi:10.1021/jp990858e

(37) Rejnek，J.;Hanus，M.;Kabelac，M.;Ryjacek，F.;Hobza，P.Phys.Chem.Chem.Phys.2005，7，2006.doi:10.1039/b501499a

(38) Ha，T.K.;Gunthard，H.H.J.Am.Chem.Soc.1993，115，11939.doi:10.1021/ja00078a036

(39) Ueda，T.;Fox，J.J.J.Am.Chem.Soc.1963，85，4024.doi:10.1021/ja00907a026

(40) Dreyfus，M.;Bensaude，O.;Dodin，G.;Dubois，J.E.J.Am.Chem.Soc.1976，98，6338.doi:10.1021/ja00436a045

(41) Yu，C.;Peng，S.;Akiyama，I.;Lin，J.;LeBreton，P.R.J.Am.Chem.Soc.1978，100，2303.doi:10.1021/ja00476a006

(42) Szczesniak，M.;Szczepaniak，K.;Kwiatkowski，J.;KuBulat，K.;Person，W.J.Am.Chem.Soc.1988，110，8319.doi:10.1021/ja00233a006

(43) Nowak，M.J.;Lapinski，L.;Fulara，J.Spectrochim.Acta，Part A:Mol.Spectrosc.1989，45，229.doi:10.1016/0584-8539(89)80129-1

(44) Brown，R.D.;Godfrey，P.D.;McNaughton，D.;Pierlot，A.P.J.Am.Chem.Soc.1989，111，2308.doi:10.1021/ja00188a058

(45) Nir，E.;Müller，M.;Grace，L.I.;de Vries，M.S.Chem.Phys.Lett.2002，355，59.doi:10.1016/S0009-2614(02)00180-X

(46) Tomic，K.;Tatchen，J.;Marian，C.M.J.Phys.Chem.A 2005，109，8410.doi:10.1021/jp051510o

(47) Fogarasi，G.J.Mol.Struct.1997，413，271.

(48) Colominas，C.;Luque，F.J.;Orozco，M.J.Am.Chem.Soc.1996，118，6811.doi:10.1021/ja954293l

(49) Kwiatkowski，J.S.;Leszczyński，J.J.Phys.Chem.1996，100，941.doi:10.1021/jp9514640

(50) Sponer，J.;Hobza，P.J.Phys.Chem.1994，98，3161.doi:10.1021/j100063a019

(51) Gould，I.R.;Burton，N.A.;Hall，R.J.;Hillier，I.H.J.Mol.Struct:Theochem 1995，331，147.doi:10.1016/0166-1280(94)03887-Q

(52) Les，A.;Adamowicz，L.;Bartlett，R.J.J.Phys.Chem.1989，93，4001.doi:10.1021/j100347a028

(53) Kobayashi，R.J.Phys.Chem.A 1998，102，10813.doi:10.1021/jp9829546

(54) Ha，T.K.;Keller，H.J.;Gunde，R.;Gunthard，H.H.J.Phys.Chem.A 1999，103，6612.doi:10.1021/jp984564p

(55) Fogarasi，G.J.Phys.Chem.A 2002，106，1381.doi:10.1021/jp013067x

(56) van Mourik，T.;Benoit，D.M.;Price，S.L.;Clary，D.C.Phys.Chem.Chem.Phys.2000，2，1281.doi:10.1039/a909183a

(57) Sambrano，J.R.;de Souza，A.R.;Queralt，J.J.;Andrés，J.Chem.Phys.Lett.2000，317，437.doi:10.1016/S0009-2614(99)01394-9

(58) Kwiatkowski，J.S.;Bartlett，R.J.;Person，W.B.J.Am.Chem.Soc.1988，110，2353.doi:10.1021/ja00216a001

(59) Chenon，M.T.;Pugmire，R.J.;Grant，D.M.;Panzica，R.P.;Townsend，L.B.J.Am.Chem.Soc.1975，97，4636.doi:10.1021/ja00849a028

(60) Dreyfus，M.;Dodin，G.;Bensaude，O.;Dubois，J.E.J.Am.Chem.Soc.1975，97，2369.doi:10.1021/ja00842a011

(61) Gonnella，N.C.;Nakanishi，H.;Holtwick，J.B.;Horowitz，D.S.;Kanamori，K.;Leonard，N.J.;Roberts，J.D.J.Am.Chem.Soc.1983，105，2050.doi:10.1021/ja00345a063

(62) Nowak，M.J.;Lapinski，L.;Kwiatkowski，J.S.;Leszczynski，J.J.Phys.Chem.1996，100，3527.doi:10.1021/jp9530008

(63) Nowak，M.J.;Rostkowska，H.;Lapinski，L.;Kwiatkowski，J.S.;Leszczynski，J.J.Phys.Chem.1994，98，2813.doi:10.1021/j100062a015

(64) Laxer，A.;Major，D.T.;Gottlieb，H.E.;Fischer，B.J.Org.Chem.2001，66，5463.doi:10.1021/jo010344n

(65) Plützer，C.;Kleinermanns，K.Phys.Chem.Chem.Phys.2002，4，4877.doi:10.1039/b204595h

(66) Cohen，B.;Hare，P.M.;Kohler，B.J.Am.Chem.Soc.2003，125，13594.doi:10.1021/ja035628z

(67) Vogt，N.;Dorofeeva，O.V.;Sipachev，V.A.;Rykov，A.N.J.Phys.Chem.A 2009，113，13816.doi:10.1021/jp905755u

(68) Broo，A.J.Phys.Chem.A 1998，102，526.doi:10.1021/jp9713625

(69) Mennucci，B.;Toniolo，A.;Tomasi，J.J.Phys.Chem.A 2001，105，4749.doi:10.1021/jp0045843

(70) Guerra，C.F.;Bickelhaupt，F.M.;Saha，S.;Wang，F.J.Phys.Chem.A 2006，110，4012.doi:10.1021/jp057275r

(71) Hanus，M.;Kabelac，M.;Rejnek，J.;Ryjacek，F.;Hobza，P.J.Phys.Chem.B 2004，108，2087.

(72) Roca-Sanjuán，D.;Merchán，M.;Serrano-Andrés，L.;Rubio，M.J.Chem.Phys.2008，129，095104.doi:10.1063/1.2958286

(73) Salter，L.M.;Chaban，G.M.J.Phys.Chem.A 2002，106，4251.doi:10.1021/jp014620d

(74) Choi，M.Y.;Dong，F.;Han，S.W.;Miller，R.E.J.Phys.Chem.A 2008，112，7185.doi:10.1021/jp8012688

(75)Wiorkiewicz-Kuczera，J.;Karplus，M.J.Am.Chem.Soc.1990，112，5324.doi:10.1021/ja00169a045

(76) Kim，H.S.;Ahn，D.S.;Chung，S.Y.;Kim，S.K.;Lee，S.J.Phys.Chem.A 2007，111，8007.doi:10.1021/jp074229d

(77)Aidas，K.;Mikkelsen，K.V.;Kongsted，J.Phys.Chem.Chem.Phys.2010，12，761.doi:10.1039/b915604f

(78) Lin，J.;Yu，C.;Peng，S.;Akiyama，I.;Li，K.;Lee，L.K.;LeBreton，P.R.J.Phys.Chem.1980，84，1006.doi:10.1021/j100446a015

(79) Nir，E.;Grace，L.;Brauer，B.;de Vries，M.S.J.Am.Chem.Soc.1999，121，4896.doi:10.1021/ja984088g

(80) Nir，E.;Janzen，C.;Imhof，P.;Kleinermanns，K.;de Vries，M.S.J.Chem.Phys.2001，115，4604.doi:10.1063/1.1391443

(81) Piuzzi，F.;Mons，M.;Dimicoli，I.;Tardivel，B.;Zhao，Q.Chem.Phys.2001，270，205.doi:10.1016/S0301-0104(01)00393-7

(82) Mons，M.;Dimicoli，I.;Piuzzi，F.;Tardivel，B.;Elhanine，M.J.Phys.Chem.A 2002，106，5088.doi:10.1021/jp0139742

(83) Choi，M.Y.;Miller，R.E.J.Am.Chem.Soc.2006，128，7320.doi:10.1021/ja060741l

(84) Mons，M.;Piuzzi，F.;Dimicoli，I.;Gorb，L.;Leszczynski，J.J.Phys.Chem.A 2006，110，10921.doi:10.1021/jp063738x

(85) Zhou，J.;Kostko，O.;Nicolas，C.;Tang，X.;Belau，L.;de Vries，M.S.;Ahmed，M.J.Phys.Chem.A 2009，113，4829.doi:10.1021/jp811107x

(86)Sabio，M.;Topiol，S.;Lumma，W.C.J.Phys.Chem.1990，94，1366.doi:10.1021/j100367a032

(87) Gorb，L.;Kaczmarek，A.;Gorb，A.;Sadlej，A.J.;Leszczynski，J.J.Phys.Chem.B 2005，109，13770.doi:10.1021/jp050394m

(88) Shukla，M.K.;Leszczynski，J.Chem.Phys.Lett.2006，429，261.doi:10.1016/j.cplett.2006.08.037

(89) Liang，W.;Li，H.R.;Hu，X.B.;Han，S.J.Chem.Phys.2006，328，93.doi:10.1016/j.chemphys.2006.06.025

(90) Schiedt，J.;Weinkauf，R.;Neumark，D.M.;Schlag，E.W.Chem.Phys.1998，239，511.doi:10.1016/S0301-0104(98)00361-9

(91) Desfrancois，C.;Abdoul Carime，H.;Schermann，J.J.Chem.Phys.1996，104，7792.doi:10.1063/1.471484

(92) Hendricks，J.;Lyapustina，S.;De Clercq，H.;Snodgrass，J.;Bowen，K.J.Chem.Phys.1996，104，7788.doi:10.1063/1.471482

(93) Dolgounitcheva，O.;Zakrzewski，V.G.;Ortiz，J.V.Chem.Phys.Lett.1999，307，220.doi:10.1016/S0009-2614(99)00492-3

(94) Wesolowski，S.S.;Leininger，M.L.;Pentchev，P.N.;Schaefer，H.F.J.Am.Chem.Soc.2001，123，4023.doi:10.1021/ja003814o

(95) Li，X.;Cai，Z.;Sevilla，M.D.J.Phys.Chem.A 2002，106，1596.doi:10.1021/jp013337b

(96) Desfrancois，C.;Periquet，V.;Bouteiller，Y.;Schermann，J.P.J.Phys.Chem.A 1998，102，1274.doi:10.1021/jp9728417

(97) Dolgounitcheva，O.;Zakrzewski，V.G.;Ortiz，J.V.J.Phys.Chem.A 2002，106，8411.doi:10.1021/jp020080o

(98) Hush，N.S.;Cheung，A.S.Chem.Phys.Lett.1975，34，11.doi:10.1016/0009-2614(75)80190-4

(99) Orlov，V.;Smirnov，A.;Varshavsky，Y.M.Tetrahedron Lett.1976，17，4377.doi:10.1016/0040-4039(76)80120-7

(100) Sevilla，M.D.;Besler，B.;Colson，A.O.J.Phys.Chem.1995，99，1060.doi:10.1021/j100003a032

(101) Hutter，M.;Clark，T.J.Am.Chem.Soc.1996，118，7574.doi:10.1021/ja953370+

(102) Le?o，M.B.C.;Longo，R.L.;Pav?o，A.C.J.Mol.Struct.-Theochem 1999，490，145.doi:10.1016/S0166-1280(99)00095-0

(103) Meot-Ner，M.J.Am.Chem.Soc.1979，101，2396.doi:10.1021/ja00503a027

(104) Lias，S.J.Phys.Chem.Ref.Data 1984，13，695.doi:10.1063/1.555719

(105) Greco，F.;Liguori，A.;Sindona，G.;Uccella，N.J.Am.Chem.Soc.1990，112，9092.doi:10.1021/ja00181a009

(106) Podolyan，Y.;Gorb，L.;Leszczynski，J.J.Phys.Chem.A 2000，104，7346.doi:10.1021/jp000740u

(107) Russo，N.;Toscano，M.;Grand，A.;Jolibois，F.J.Comput.Chem.1998，19，989.

(108) Chandra，A.K.;Nguyen，M.T.;Uchimaru，T.;Zeegers-Huyskens，T.J.Phys.Chem.A 1999，103，8853.doi:10.1021/jp990647+

(109) Cao，G.J.;Xu，H.G.;Li，R.Z.;Zheng，W.J.Chem.Phys.2012，136，014305.doi:10.1063/1.3671945

(110) Esteruelas，M.A.;Garcia-Raboso，J.;Olivan，M.Inorg.Chem.2012，51，9522.doi:10.1021/ic3013238

(111) Javan，M.J.;Tehrani，Z.A.;Fattahi，A.;Jamshidi，Z.Struct.Chem.2012，23，1843.doi:10.1007/s11224-012-9993-2

(112)Lin，Y.X.;Wang，H.Y.;Gao，S.M.;Wu，Y.X.;Li，R.H.Acta Phys.-Chim.Sin.2013，29，1233.[林月霞，王紅艷，高思敏，吳穎曦，李汝虎.物理化學學報，2013，29，1233.]doi:10.3866/PKU.WHXB201304022

(113) Palamarchuk，G.V.;Shishkin，O.V.;Gorb，L.;Leszczynski，J.J.Phys.Chem.B 2013，117，2841.doi:10.1021/jp311363c

(114) Villabona-Monsalve，J.P.;Noria，R.;Matsika，S.;Peon，J.J.Am.Chem.Soc.2012，134，7820.doi:10.1021/ja300546x