Premium
Density functional study of guanine and uracil quartets and of guanine quartet/metal ion complexes
Author(s) -
Meyer Michael,
Steinke Thomas,
Brandl Maria,
Sühnel Jürgen
Publication year - 2000
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/1096-987x(20010115)22:1<109::aid-jcc11>3.0.co;2-5
Subject(s) - guanine , chemistry , hydrogen bond , ion , crystallography , stacking , uracil , density functional theory , metal , base (topology) , interaction energy , metal ions in aqueous solution , stereochemistry , computational chemistry , atomic physics , molecule , physics , dna , nucleotide , mathematical analysis , biochemistry , mathematics , organic chemistry , gene
The structures and interaction energies of guanine and uracil quartets have been determined by B3LYP hybrid density‐functional calculations. The total interaction energy Δ E T of the C 4h ‐symmetric guanine quartet consisting of Hoogsteen‐type base pairs with two hydrogen bonds between two neighbor bases is −66.07 kcal/mol at the highest level. The uracil quartet with C6 H6O4 interactions between the individual bases has only a small interaction energy of −20.92 kcal mol −1 , and the interaction energy of −24.63 kcal/mol for the alternative structure with N3H3O4 hydrogen bonds is only slightly more negative. Cooperative effects contribute between 10 and 25% to all interaction energies. Complexes of metal ions with G‐quartets can be classified into different structure types. The one with Ca 2+ in the central cavity adopts a C 4h ‐symmetric structure with coplanar bases, whereas the energies of the planar and nonplanar Na + complexes are almost identical. The small ions Li + , Be 2+ , Cu + , and Zn 2+ prefer a nonplanar S 4 ‐symmetric structure. The lack of coplanarity prevents probably a stacking of these base quartets. The central cavity is too small for K + ions and, therefore, this ion favors in contrast to all other investigated ions a C 4 ‐symmetric complex, which is 4.73 kcal/mol more stable than the C 4h ‐symmetric one. The distance 1.665 Å between K + and the root‐mean‐square plane of the guanine bases is approximately half of the distance between two stacked G‐quartets. The total interaction energy of alkaline earth ion complexes exceeds those with alkali ions. Within both groups of ions the interaction energy decreases with an increasing row position in the periodic table. The B3LYP and BLYP methods lead to similar structures and energies. Both methods are suitable for hydrogen‐bonded biological systems. Compared with the before‐mentioned methods, the HCTH functional leads to longer hydrogen bonds and different relative energies for two U‐quartets. Finally, we calculated also structures and relative energies with the MMFF94 forcefield. Contrary to all DFT methods, MMFF94 predicts bifurcated CHO contacts in the uracil quartet. In the G‐quartet, the MMFF94 hydrogen bond distances N2H22N7 are shorter than the DFT distances, whereas the N1H1O6 distances are longer. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 109–124, 2001