THEORETICAL DENSITY-FUNCTIONAL AND AB-INITIO COMPUTATIONAL STUDY OF THE TAUTOMERIC, VIBRATIONAL AND H-BOND PROPERTIES OF 1,7-DI-CH3-GUANINE- A COMPARISON WITH FT-IR DATA FROM MATRIX-ISOLATION EXPERIMENTS

In this work we report the results of a combined matrix-isolation FT-W and theoretical computational DFT and ab initio study of the tautomer ic and vibrational characteristics of 1,7-dimethyl-guanine (17DMG). Th e HF, MP2 and DFT methods are all in agreement in predicting that the amino-oxo tautomer of 17DMG is the most stable form, and that the imin o-oxo tautomer, which can exist in two different geometrical isomers, is the second most stable form. The energy difference between the amin o-oxo and imino-oxo tautomer is rather targe: 23, 25 and 33 kJ mol(-1) , for the HF, DFT and MP2 methods, respectively. The imino-hydroxy tau tomer is, according to the calculations, by far the least stable form (Delta E > 180 kJ mol(-1)). The FT-IR spectra of 17DMG isolated in an argon matrix are in agreement with the theoretical calculation indicat ing that only the amino-oxo tautomer should be present in the gas phas e in detectable amounts. The experimental IR frequencies agree quite w ell with the theoretical values. When a single scaling factor is appli ed to approximately correct the calculated frequencies for various sys tematic errors in theoretical approach, the mean frequency deviations for the HF/6-31++G*, DFT/6-31G** and DFT/6-31++G** calculated spectra are 17.5 cm(-1), 12.2 cm(-1) and 16.3 cm(-1), respectively. The DFT m ean-frequency-deviations decrease to 8.1 cm(-1) and 10.1 cm-l when var iable scaling factors are applied. The stabilities of the different wa ter complexes of the amino-oxo tautomer of 17DMG are investigated by o ptimizing the complexes at the HF/6-31++G* level of theory and calcul ating the MP2 energies at these optimized geometries. The N-3...HO-H c omplex, which has a closed structure with two H-bonds, is predicted to be the most stable structure. The MP2 calculated energy differences b etween the most stable complex and the C-6=O...HO-K, the H-NH...OH2 an d the H2N...HO-H hetero-dimers are 14.28 kJ mol(-1), 19.54 kJ mol(-1) and 20.73 kJ mol(-1), respectively. (C) 1998 Elsevier Science B.V.