Jones D B, Wang F, Winkler D A, Brunger M J
School of Chemistry, Physics and Earth Sciences, Flinders University, GPO Box 2100, Adelaide, SA 5001, Australia.
Biophys Chem. 2006 May 20;121(2):105-20. doi: 10.1016/j.bpc.2005.12.006. Epub 2006 Feb 7.
Electronic structural signatures of the guanine-7H and guanine-9H tautomers have been investigated on an orbital by orbital basis using dual space analysis. A combination of density functional theory (B3LYP/TZVP), the statistical average of model orbital potentials (SAOP/TZ2P) method and outer valence Green's function theory (OVGF/TZVP) has been used to generate optimal tautomer geometries and accurate ionization energy spectra for the guanine tautomer pair. The present work found that the non-planar form for both of the guanine keto pair possesses lower energies than their corresponding planar counterparts, and that the canonical form of the guanine-7H tautomer has slightly lower total energy than guanine-9H. This latter result is in agreement with previous experimental and theoretical findings. In the planar guanine pair the geometric parameters and anisotropic molecular properties are compared, focusing on changes caused by the mobile proton transfer. It is demonstrated that the mobile proton only causes limited disturbance to isotropic properties, such as geometry and the energetics, of the guanine keto tautomer pair. The exception to this general statement is for related local changes such as the N((7))-C((8)) and C((8))-N((9)) bond length resonance between the single and double bonds, reflecting the nitrogen atom being bonded with the mobile proton in the tautomers. The mobile proton distorts the electron distribution of the tautomers, which leads to significant changes in the molecular anisotropic properties. The dipole moment of guanine-7H is altered by about a factor of three, from 2.23 to 7.05 D (guanine-9H), and the molecular electrostatic potentials also reflect significant electron charge distortion. The outer valence orbital momentum distributions, which were obtained using the plane wave impulse approximation (PWIA), have demonstrated quantitatively that the outer valence orbitals of the tautomer pair can be divided into three groups. That is orbitals 1a''-7a'' and 18a', which do not have visible alternations in the tautomeric process (which consist of either pi orbitals or are close to the inner valence shell); a second group comprising orbitals 19a'-22a', 25a', 26a', 28a', 29a' and 31a', which show small perturbations as a result of the mobile hydrogen locations; and group three, orbitals 23a', 24a', 27a', 30a' and 32a', which demonstrate significant changes due to the mobile proton transfer and are therefore considered as signature orbitals of the G-7H/G-9H keto tautomeric process.
利用双空间分析,在逐个轨道的基础上研究了鸟嘌呤-7H和鸟嘌呤-9H互变异构体的电子结构特征。结合密度泛函理论(B3LYP/TZVP)、模型轨道势的统计平均值(SAOP/TZ2P)方法和外层价电子格林函数理论(OVGF/TZVP),生成了鸟嘌呤互变异构体对的最佳互变异构体几何结构和精确的电离能谱。目前的研究发现,鸟嘌呤酮式对的两种非平面形式的能量均低于其相应的平面形式,并且鸟嘌呤-7H互变异构体的标准形式的总能量略低于鸟嘌呤-9H。后一结果与先前的实验和理论发现一致。在平面鸟嘌呤对中,比较了几何参数和各向异性分子性质,重点关注由移动质子转移引起的变化。结果表明,移动质子仅对鸟嘌呤酮式互变异构体对的各向同性性质(如几何结构和能量学)造成有限的干扰。这一普遍结论的例外情况是相关的局部变化,如单键和双键之间的N((7))-C((8))和C((8))-N((9))键长共振,这反映了互变异构体中与移动质子相连的氮原子。移动质子使互变异构体的电子分布发生畸变,从而导致分子各向异性性质发生显著变化。鸟嘌呤-7H的偶极矩改变了约三倍,从2.23 D变为7.05 D(鸟嘌呤-9H),分子静电势也反映出显著的电荷畸变。利用平面波脉冲近似(PWIA)获得的外层价轨道动量分布定量地表明,互变异构体对的外层价轨道可分为三组。即轨道1a''-7a''和18a',在互变异构过程中没有明显变化(由π轨道组成或接近内层价壳层);第二组包括轨道19a'-22a'、25a'、26a'、28a'、29a'和31a',由于移动氢的位置而表现出小的扰动;第三组是轨道23a'、24a'、27a'、30a'和32a',由于移动质子转移而表现出显著变化,因此被视为G-7H/G-9H酮式互变异构过程的特征轨道。