Ladik János, Bende Attila, Bogár Ferenc
Chair for Theoretical Chemistry and Laboratory of the National Foundation for Cancer Research, Friedrich-Alexander-University-Erlangen-Nürnberg, Egerlandstr 3, Erlangen, Germany.
J Chem Phys. 2008 Mar 14;128(10):105101. doi: 10.1063/1.2832860.
Using the ab initio Hartree-Fock crystal orbital method in its linear combination of atomic orbital form, the energy band structure of the four homo-DNA-base stacks and those of poly(adenilic acid), polythymidine, and polycytidine were calculated both in the absence and presence of their surrounding water molecules. For these computations Clementi's double zeta basis set was applied. To facilitate the interpretation of the results, the calculations were supplemented by the calculations of the six narrow bands above the conduction band of poly(guanilic acid) with water. Further, the sugar-phosphate chain as well as the water structures around poly(adenilic acid) and polythymidine, respectively, were computed. Three important features have emerged from these calculations. (1) The nonbase-type or water-type bands in the fundamental gap are all close to the corresponding conduction bands. (2) The very broad conduction band (1.70 eV) of the guanine stack is split off to seven narrow bands in the case of poly(guanilic acid) (both without and with water) showing that in the energy range of the originally guanine-stack-type conduction band, states belonging to the sugar, to PO(4)(-), to Na(+), and to water mix with the guanine-type states. (3) It is apparent that at the homopolynucleotides with water in three cases the valence bands are very similar (polycytidine, because it has a very narrow valence band, does not fall into this category). We have supplemented these calculations by the computation of correlation effects on the band structures of the base stacks by solving the inverse Dyson equation in its diagonal approximation taken for the self-energy the MP2 many body perturbation theory expression. In all cases the too large fundamental gap decreased by 2-3 eV. In most cases the widths of the valence and conduction bands, respectively, decreased (but not in all cases). This unusual behavior is most probably due to the rather large complexity of the systems. From all this emerges the following picture for the charge transport in DNA: There is a possibility in short segments of the DNA helix of a Bloch-type conduction of holes through the nucleotide base stacks of DNA combined with hopping (and in a lesser degree with tunneling). The motivation of this large scale computation was that recently in Zurich (ETH) they have performed high resolution x-ray diffraction experiments on the structure of the nucleosomes. The 8 nucleohistones in them are wrapped around by a DNA superhelix of 147 base pairs in the DNA B form. The most recent investigations have shown that between the DNA superhelix (mostly from its PO(4) (-) groups) there is a charge transfer to the positively charged side chains (first of all arginines and lysines) of the histones at 120 sites of the superhelix. This would cause a hole conduction in DNA and an electronic one in the proteins.
采用从头算哈特里 - 福克晶体轨道方法的原子轨道线性组合形式,计算了四种同型DNA碱基堆积以及聚腺苷酸、聚胸苷和聚胞苷在有无周围水分子存在时的能带结构。对于这些计算,应用了克莱门蒂的双ζ基组。为便于结果解释,通过计算聚鸟苷酸与水体系导带上方的六个窄带对计算进行了补充。此外,分别计算了聚腺苷酸和聚胸苷周围的糖 - 磷酸链以及水结构。这些计算得出了三个重要特征。(1)基本能隙中的非碱基型或水型能带都靠近相应的导带。(2)鸟嘌呤堆积的非常宽的导带(1.70 eV)在聚鸟苷酸(无论有无水)的情况下分裂为七个窄带,这表明在原本鸟嘌呤堆积型导带的能量范围内,属于糖、磷酸根(PO(4)(-))、钠离子(Na(+))和水的态与鸟嘌呤型态混合。(3)显然,在三种有水的同聚核苷酸情况下,价带非常相似(聚胞苷因为其价带非常窄,不属于此类别)。我们通过求解以MP2多体微扰理论表达式为自能的对角近似下的逆戴森方程,计算了碱基堆积能带结构的相关效应,对这些计算进行了补充。在所有情况下,过大的基本能隙减小了2 - 3 eV。在大多数情况下,价带和导带的宽度分别减小了(但并非在所有情况下)。这种不寻常的行为很可能是由于系统相当大的复杂性。由此得出了关于DNA中电荷传输的如下图景:在DNA螺旋的短片段中,有可能通过DNA的核苷酸碱基堆积进行布洛赫型空穴传导,并伴有跳跃(以及在较小程度上伴有隧穿)。进行这种大规模计算的动机是最近在苏黎世(瑞士联邦理工学院)他们对核小体的结构进行了高分辨率x射线衍射实验。其中的8个核组蛋白被B型DNA的147个碱基对的超螺旋缠绕。最新研究表明,在DNA超螺旋(主要来自其磷酸根(PO(4) (-))基团)与超螺旋120个位点处组蛋白带正电的侧链(首先是精氨酸和赖氨酸)之间存在电荷转移。这将导致DNA中的空穴传导和蛋白质中的电子传导。