金属离子(Na~+、K~+、Ca~(2+)、Mg~(2+)、Zn~(2+))与鸟嘌呤异构体配合物的稳定性

在B3LYP/6-311++G**水平上用极化连续介质模型(PCM)系统研究了金属离子(M+/2+=Na+,K+,Ca2+,Mg2+,Zn2+)和十三种鸟嘌呤异构体形成的配合物GnxM+/2+(n为鸟嘌呤异构体的编号,x表示M+/2+与鸟嘌呤异构体的结合位点)在气(g)液(a)两相中的稳定性顺序.着重探讨了液相中配合物的稳定性差异,并且从溶质-溶剂效应、结合能、形变能及异构体的相对能量等几个方面分析了造成稳定顺序发生变化的原因.报道了溶液中这五种金属离子与鸟嘌呤异构体结合形成的六种基态配合物:aG1N2,N3Na+,aG1N2,N3K+,aG1O6,N7Ca2+,aG1N2,N3Mg2+(aG1O6,N7Mg2+),aG2N3,N9Zn2+.可以看出,除了在Zn2+配合物中鸟嘌呤异构体为G2外,构成其余四种金属离子配合物的鸟嘌呤异构体都是G1,但结合位点不同.同时对气相中各类配合物稳定性也进行了系统的排序,并报道了几种较稳定的配合物,如:gG3N1,O6K+,gG5N1,O6K+,gG3N1,O6Ca2+/Mg2+,gG4O6,N7Ca2+/Mg2+.     

溶液中Zn2+与腺嘌呤异构体间相互作用的理论研究

用B3LYP/6-311++G**方法和PCM及Onsager模型研究了Zn2+与腺嘌呤异构体在溶液中的11种较为稳定配合物. 结果显示, 这些配合物在溶液中的稳定性顺序与气相中明显不同, 其结合位点表现出如下的规律性, 在亚氨基类配合物中, Zn2+与腺嘌呤的N7、N6位结合比与N1、N6位结合形成的配合物更稳定; 氨基类配合物中, Zn2+以"双齿"形式与腺嘌呤异构体上的氮结合时的优先顺序为(N3和N9)>(N7和N6)>(N1和N6). 研究表明, 不论气相还是溶液相, 孤立的腺嘌呤分子内的质子转移较困难, 结合Zn2+后也不能明显降低关键步骤的活化能; 结合Cu2+却能明显地降低气相中关键步骤的活化能, 但溶剂效应却不利于Cu2+引发腺嘌呤分子内的质子转移.     

槲皮素与腺嘌呤氢键作用的最佳位点

本文优化得到了16个由槲皮素与腺嘌呤形成的氢键复合物的稳定结构,并计算了它们的结合能.研究发现,在气相和水相中,槲皮素均通过qu1位点与腺嘌呤作用形成稳定的氢键复合物.比较了腺嘌呤与槲皮素形成的氢键复合物、腺嘌呤与胸腺嘧啶形成的Watson-Crick碱基对的相对稳定性.在气相条件下Watson-Crick碱基对更稳定,在水相条件下腺嘌呤与槲皮素形成的氢键复合物更稳定,说明水相条件下腺嘌呤与槲皮素之间的相互作用强于与胸腺嘧啶之间的相互作用.基于标准反应Gibbs自由能变的计算结果估算了水相条件下腺嘌呤与槲皮素形成的氢键复合物和Watson-Crick碱基对的相对平衡浓度.     

碱(土)金属离子与2-(3′-羟基-2′-吡啶基)苯并噁唑阳离子-π复合物及其分子内质子转移过程的理论研究

采用密度泛函B3LYP方法,在6-311++G(d,p)基组水平上对碱(土)金属离子(Li+,Na+,K+,Be2+,Mg2+和Ca2+)与2-(3’-羟基-2’-吡啶基)苯并噁唑(HPyBO)的36种阳离子-π复合物的初始构型进行了几何全优化,并计算了其相互作用能.结果表明,碱(土)金属离子与HPyBO复合物有较强的阳离子-π相互作用,部分复合物甚至达到了化学键的强度.相对能量的变化表明碱(土)金属离子的作用能改变HPyBO分子内质子转移过程的能垒,甚至可以导致优势构型反转.当考虑水的溶剂效应后,各质子转移异构体的相对能量及质子转移的能垒均有一定程度的改变.另外,应用分子中的原子(AIM)方法对复合物分子内氢键的键临界点性质进行了分析.     

2,6-二巯基嘌呤互变异构体热力学稳定性的密度泛函理论研究

利用密度泛函(DFT)B3LYP/6-311G(d,p)方法,水相计算采用自洽反应场(SCRF)中的Onsager模型,对气相和水相中可能存在的13种2,6-二巯基嘌呤互变异构体进行了全优化,并计算了各异构体的热力学参数、偶极矩及原子净电荷。计算结果表明,不论是气相还是水相,二硫酮DTP(1,3,7)是最稳定的异构体。溶剂化效应使各异构体的稳定性均增强,偶极矩大者其稳定性显著增大。溶剂化吉布斯自由能与异构体在两相中偶极矩之差存在相关性。二硫酮DTP(1,3,7)在水相中与致癌物BPDE进行亲核取代反应时,二硫酮DTP(1,3,7)中的S10原子优先进攻亲电试剂BPDE.     

6-巯基嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G水平下，对8种气相和水相中可能存在的6-巯基嘌呤异构体进行了几何构型的全自由度优化，并计算出它们的总能量、焓、熵、吉布斯自由能．Onsager反应场溶剂模型用于水相的计算．计算结果表明，6-巯基嘌呤在气相和水相中主要以硫酮形式存在．在气相中，硫酮．N(7)(H)要比硫酮-N(9)(H)更稳定，而在水相中，则硫酮-N(9)(H)要比硫酮-N(7)(H)更稳定．计算结果同已有实验结果一致．6-巯基嘌呤的熵效应小，对互变异构平衡几乎没有显著的影响，而焓变对互变异构产生了主要的影响．较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响，溶剂化吉布斯自由能与异构体的气相偶极矩存在相关性．     

5-氟尿嘧啶和5-氯尿嘧啶及其互变异构体的理论计算研究

采用HF/3-21G方法, 对6种气相和水相中可能存在的5-氟尿嘧啶(和5-氯尿嘧啶)互变异构体进行了构象分析.采用B3LYP/6-311＋G**方法对处于优势构象时的各互变异构体进行了几何全优化, 并计算出它们的总能量、焓、熵、吉布斯自由能. Onsager反应场溶剂模型用于水相的计算. 计算结果表明, 5-氟尿嘧啶和5-氯尿嘧啶在气相中和水相中主要以双酮形式存在. 5-氟尿嘧啶和5-氯尿嘧啶的熵效应小, 对互变异构平衡没有显著的影响, 而焓变对互变异构产生了主要的影响. 讨论了水溶剂化作用对异构体的能量、电荷分布和偶极矩的影响. 溶剂化自由能与异构体的气相偶极矩存在相关性. 另外, 详细地将5-氟尿嘧啶和5-氯尿嘧啶与尿嘧啶进行了对比, 获得三者最稳定异构体间电子结构异同的有用信息.     

6-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平下,对14种气相和水相中可能存在的6-硫代黄嘌呤异构体进行了几何构型的全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能。Onsager反应场溶剂模型用于水相的计算.计算结果表明,6-硫代黄嘌呤在气相中和水相中主要以硫酮的形式存在.在气相和水相中,硫酮-N7(H)均比硫酮-N9(H)更稳定.计算结果同已有实验结果一致.6-硫代黄嘌呤异构化的熵效应小,对互变异构平衡几乎没有显著的影响,而焓变对互变异构产生了主要的影响.较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响.     

密度泛函研究不同价态金属离子与谷胱甘肽相互作用

采用密度泛函DFT/B3LYP方法,研究了在气相和生物环境内稳定存在的2种构型的还原型谷胱甘肽（GSH）与不同价态金属铬离子(Cr2+,Cr3+,Cr6+)相互作用。 金属离子的电荷越高、半径越小,与GSH结合能越大,使GSH的变形程度也越大。 金属Cr6+在气相和液相条件与GSH作用均促使了GSH的骨架断裂,末端羧基发生脱羧。 Cr3+和Cr2+与气相中性和液相两性离子结构的GSH分子相互作用均形成了9种稳定的复合物,与气相计算结果相比,考虑溶剂化效应之后,金属离子与GSH两性离子作用的结合能要比与在气相条件下中性的GSH相互作用能大大降低。     

2-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平上,对14种气相和水相中可能存在的2-硫代黄嘌呤互变异构体进行了几何构型全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能.Onsager反应场溶剂模型用于水相的计算.计算结果表明,2-硫代黄嘌呤在气相和水相中主要以酮式结构形式存在,与已有实验结果一致.在气相和水相中,酮式结构—N(7)(H)均比酮式结构—N(9)(H)更稳定.2-硫代黄嘌呤互变异构的熵效应小,对互变异构平衡没有显著的影响,而焓变对互变异构却产生了主要的影响.水溶剂化自由能与异构体的气相偶极矩存在相关性.另外,较详细地考察了2-硫代黄嘌呤与6-硫代黄嘌呤的相对稳定性.     

Bond energies and attachments sites of sodium and potassium cations to DNA and RNA nucleic acid bases in the gas phase

Gas-phase metal affinities of DNA and RNA bases for the Na(+) and K(+) ions were determined at density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of nucleic acid bases and the alkali ions on the different binding sites, were considered. Structural features of the sodium and potassium complexes were found to be similar except in some uracil and thymine compounds in which the tendency of potassium ion toward monocoordination appeared evident. B3LYP bond energies for both metal ions were in agreement with the available experimental results in the cases of uracil and thymine for which the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base. For adenine, although the interaction of the ions with the most stable free tautomer generated the least stable molecular complex, the best agreement with experiment was found in just this case. For the remaining cytosine and guanine bases, our calculations indicated that the metal ion affinity value closest to experiment should be determined taking into account the role played by the different tautomers of the free bases with similar energy and all the possible complexes obtained by them.     

Protonated adenine: Tautomers,solvated clusters,and dissociation mechanisms

Ab initio and density functional theory calculations at the B3-MP2 and CCSD(T)/6-311 + G(3df,2p) levels of theory are reported that address the protonation of adenine in the gas phase, water clusters, and bulk aqueous solution. The calculations point to N-1-protonated adenine (1+) as the thermodynamically most stable cationic tautomer in the gas phase, water clusters, and bulk solution. This strongly indicates that electrospray ionization of adenine solutions produces tautomer 1+ with a specificity calculated as 97-90% in the 298-473 K temperature range. The mechanisms for elimination of hydrogen atoms and ammonia from 1+ have also been studied computationally. Ion 1+ is calculated to undergo fast migrations of protons among positions N-1, C-2, N-3, N-10, N-7, and C-8 that result in an exchange of five hydrogens before loss of a hydrogen atom forming adenine cation radical at 415 kJ mol(-1) dissociation threshold energy. The elimination of ammonia is found to be substantially endothermic requiring 376-380 kJ mol(-1) at the dissociation threshold and depending on the dissociation pathway. The overall dissociation is slowed by the involvement of ion-molecule complexes along the dissociation pathways. The competing isomerization of 1+ proceeds by a sequence of ring opening, internal rotations, imine flipping, ring closures, and proton migrations to effectively exchange the N-1 and N-10 atoms in 1+, so that either can be eliminated as ammonia. This mechanism explains the previous N-1/N-10 exchange upon collision-induced dissociation of protonated adenine.     

Ca2+‐complex stability of GAPAGPLIVPY peptide in gas and aqueous phase,investigated by affinity capillary electrophoresis and molecular dynamics simulations and compared to mass spectrometric results

Strong, sequence‐specific gas‐phase bindings between proline‐rich peptides and alkaline earth metal ions in nanoESI‐MS experiments were reported by Lehmann et al. (Rapid Commun. Mass Spectrom. 2006, 20, 2404–2410), however its relevance for physiological‐like aqueous phase is uncertain. Therefore, the complexes should also be studied in aqueous solution and the relevance of the MS method for binding studies be evaluated. A mobility shift ACE method was used for determining the binding between the small peptide GAPAGPLIVPY and various metal ions in aqueous solution. The findings were compared to the MS results and further explained using computational methods. While the MS data showed a strong alkaline earth ion binding, the ACE results showed nonsignificant binding. The proposed vacuum state complex also decomposed during a molecular dynamic simulation in aqueous solution. This study shows that the formed stable peptide–metal ion adducts in the gas phase by ESI‐MS does not imply the existence of analogous adducts in the aqueous phase. Comparing peptide–metal ion interaction under the gaseous MS and aqueous ACE conditions showed huge difference in binding behavior.     

复合物Ca+-RNA碱基的结构与稳定性

Gas-phase metal ion affinities and optimized structures of RNA nucleic acid bases for the Ca+ were determined at a density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of RNA nucleic acid and Ca+ on the different binding sites, were considered. For Cytosine, the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base tautomers. As to thymine, the bond energy of the ion with the most stable tautomer of the free nucleic acid base is the weakest among the three tautomer’s complexes, and that of the ion with least stable tautomer of the free nucleic acid base is the strongest . Uracil is similar to thymine. The two kinds of relation, bond energy and total energy for the complex, are in disagreement, as the metal affinities of RNA bases for the Ca+ depend on binding sites, and total energy of complex (Ca+-RNA base) relies on all atoms and their relative positions in the complex.     

水溶液中氨基酸侧链对G∶C碱基对间氢键影响的理论研究

利用量子化学方法研究了气相和水溶液下,氨基酸侧链与鸟嘌呤和胞嘧啶间的氢键作用.应用B3LYP/6-31+G(d,p)方法优化复合物几何结构,使用MP2/aug-cc-p VDZ方法进行复合物能量、自然键轨道(NBO)电荷和二阶稳定化能的计算.结果表明,水溶液可使氨基酸侧链与碱基或碱基对之间氢键键能显著减小;带电复合物气相和水溶液氢键键能之差范围为50.63~146.48 k J/mol,中性为0.17~24.94 k J/mol;电荷的转移量与氢键键能成正比,电荷转移量越多,复合物越稳定;二阶稳定化能与氢键键长成反比,与电荷转移量成正比,且气相与水溶液氢键二阶稳定化能之比约为两相的电荷转移量之比.水溶液对该类体系中氢键作用具有明显影响.     

Adenine radicals in the gas phase: an experimental and computational study of hydrogen atom adducts to adenine

The elusive hydrogen atom adduct to the N-1 position in adenine, which is thought to be the initial intermediate of chemical damage, was specifically generated in the gas phase and characterized by neutralization-reionization mass spectrometry. The N-1 adduct, 1,2-dihydroaden-2-yl radical (1), was generated by femtosecond electron transfer to N-1-protonated adenine that was selectively produced by electrospray ionization of adenine in aqueous-methanol solution. Radical 1 is an intrinsically stable species in the gas phase that undergoes specific loss of the N-1-hydrogen atom to form adenine, but does not isomerize to the more stable C-2 adduct, 1,2-dihydroaden-1-yl radical (5). Radicals 1 that are formed in the fifth and higher electronically excited states of DeltaE > or = 2.5 eV can also undergo ring-cleavage dissociations resulting in expulsion of HCN. The relative stabilities, dissociation, and transition state energies for several hydrogen atom adducts to adenine have been established computationally at highly correlated levels of theory. Transition state theory calculations of 298 K rate constants in the gas phase, including quantum tunnel corrections, indicate the branching ratios for H-atom additions to C-8, C-2, N-3, N-1, and N-7 positions in adenine as 0.68, 0.20, 0.08, 0.03, and 0.01, respectively. The relative free energies of adenine radicals in aqueous solution point to the C-8 adduct as the most stable tautomer, which is predicted to be the predominating (>99.9%) product at thermal equilibrium in solution at 298 K.     

Na+, Mg2+, and Zn2+ binding to all tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine: a correlated ab initio quantum chemical study

Interactions of adenine, cytosine, guanine, and thymine with Na(+), Mg(2+), and Zn(2+) cations were studied using an approximate resolution of identity correlated second-order MP2 (RI-MP2) method with the TZVPP ([5s3p2d1f/3s2p1d]) basis set. All existing tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine were considered. Cations bind mostly in a bidentate manner, and stabilization energies of these complexes are larger than those in the case when cations bind in a unidentate manner. The cation...Y (Y equal to N or O) distances for divalent metals are shorter than those for Na(+) and for Zn(2+) are mostly shorter than the Mg(2+)...Y distance. The intermolecular distances between the cation and the base for complexes containing adenine and cytosine are systematically shorter than those for complexes containing guanine and thymine. Only for cytosine the canonical keto/amino tautomer structure with ions represents the global minimum. For guanine, the metalated canonical form is again the most stable, but its stabilization energy is within less than 5% of the stabilization energies of the two other rare tautomers, which indicates that the canonical form and these two rare tautomers could coexist. The canonical structures of adenine and thymine in the presence of ions are considerably less stable (by more than 10%) than the complexes of the rare tautomers. It can be concluded that the interaction of Na(+), Mg(2+), and Zn(2+) cations with cytosine in the gas phase will not induce the change of the canonical form to the rare tautomeric form. In the case of isolated guanine, the equilibrium of the canonical form with rare tautomers can be found. For isolated adenine and thymine the presence of rare tautomers is highly probable.     

Influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions. Threshold collision-induced dissociation and theoretical studies

The influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions is investigated both experimentally and theoretically. Bond dissociation energies of alkali metal ion-halouracil complexes, M+(XU), are determined using threshold collision-induced dissociation techniques in a guided ion beam mass spectrometer, where M+ = Li+, Na+, and K+ and XU = 5-fluorouracil, 5-chlorouracil, 6-chlorouracil, 5-bromouracil, and 5-iodouracil. The structures and theoretical bond dissociation energies of these complexes are determined from ab initio calculations. Theoretical calculations are also performed to examine the influence of halogenation on the acidities, proton affinities, and Watson-Crick base pairing energies. Halogenation of uracil is found to produce a decrease in the proton affinity, an increase in the alkali metal ion binding affinities, an increase in the acidity, and stabilization of the A::U base pair. In addition, alkali metal ion binding is expected to lead to an increase in the stability of nucleic acids by reducing the charge on the nucleic acid in a zwitterion effect as well as through additional noncovalent interactions between the alkali metal ion and the nucleobases.