Department of Chemistry, Memorial University, St. John's, Newfoundland, Canada A1B 3X7.
J Phys Chem B. 2011 Jul 28;115(29):9151-9. doi: 10.1021/jp202098k. Epub 2011 Jun 29.
Oxidation of guanine in DNA yields the nucleobase damage product 8-oxoguanine (8-oxoG), whose further oxidation gives other more stable products. In the present study, the mechanism for the deamination of 8-oxoG with H(2)O, 2H(2)O, H(2)O/OH(-), and 2H(2)O/OH(-) and for protonated 8-oxoG (8-oxoGH(+)) with H(2)O has been investigated using ab initio calculations. All structures were optimized at RHF/6-31G(d), MP2/6-31G(d), and B3LYP with the 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), and 6-31++G(d.p) basis sets. Energies were determined at the G3MP2 level of theory, and solvent calculations were performed using both the polarizable continuum model (PCM) and the solvation model on density (SMD). Intrinsic reaction coordinate calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and, in the final step, an intermediate that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with H(2)O/OH(-) was found to have an overall activation energy of 187, 176, and 156 kJ mol(-1) for the gas phase, PCM, and SMD, respectively, which are ∼50 kJ mol(-1) lower than with H(2)O only. These barriers can be compared to those for the reaction of 8-oxoGH(+) with H(2)O of 248 kJ mol(-1) in the gas phase and 238 kJ mol(-1) in aqueous solution (PCM). The lowest overall activation energies (G3MP2) are for the deamination of 8-oxoG with 2H(2)O/OH(-), 134 kJ mol(-1) in the gas phase and 129 kJ mol(-1) with PCM.
DNA 中鸟嘌呤的氧化产生碱基损伤产物 8-氧代鸟嘌呤(8-oxoG),其进一步氧化产生其他更稳定的产物。在本研究中,使用从头算计算研究了 8-oxoG 与 H(2)O、2H(2)O、H(2)O/OH(-)和 2H(2)O/OH(-)以及质子化 8-oxoG (8-oxoGH(+))与 H(2)O 的脱氨反应机制。所有结构均在 RHF/6-31G(d)、MP2/6-31G(d)和 B3LYP 上进行了优化,采用了 6-31G(d)、6-31+G(d)、6-31G(d,p)、6-31+G(d,p)和 6-31++G(d,p)基组。能量在 G3MP2 理论水平上确定,使用极化连续模型 (PCM) 和密度溶剂模型 (SMD) 进行溶剂计算。进行了内禀反应坐标计算,以表征势能面上的过渡态。还计算了每个反应的热力学性质(ΔE、ΔH 和ΔG)、活化能、焓和 Gibbs 自由能。所有途径都生成初始四面体型中间体,最后一步通过 1,3-质子转移中间体解离生成产物。在 G3MP2 理论水平上,发现 H(2)O/OH(-)脱氨的总活化能分别为气相、PCM 和 SMD 中的 187、176 和 156 kJ mol(-1),比仅用 H(2)O 时低约 50 kJ mol(-1)。这些势垒可以与气相中 8-oxoGH(+)与 H(2)O 的反应的 248 kJ mol(-1)和水溶液中的 238 kJ mol(-1)(PCM)相比较。最低的总活化能(G3MP2)是 8-oxoG 与 2H(2)O/OH(-)的脱氨,气相中为 134 kJ mol(-1),PCM 中为 129 kJ mol(-1)。