Ca, Cd, Zn, and their ions interacting with Cytosine: a theoretical study.

Metal atoms play a major role in the chemical behavior of biological systems. In this work, known issues of the metal-base interactions, such as the stabilization of different tautomers of cytosine that could be incompatible with the DNA double helix, are researched using DFT methods. Ca-, Zn-, and Cd-cytosine in neutral and ionic forms were studied at the B3LYP/LANL2DZ level. Several neutral and ionic isomers were found within an interval of 10 kcal/mol of relative stability, with the most stable isomer in each group being a compound derived from the canonical isomer of cytosine, except for the dications where two iso-energetic isomers were found. Interatomic lengths from each metal atom to the nearest atoms in cytosine's ring were larger than 2 A, discouraging the possibility of a covalent interaction, as supported by additional evidence from molecular orbitals. The interaction between metal and cytosine, electrostatic in nature, is reinforced with the increase of the metal's nuclear charge. Additionally, the ionization energies of the metal-cytosine compounds exhibit a significant reduction (below 6 eV) compared with that of plain cytosine (8.7 eV), posing an interesting possibility with respect to the experimental determination of the photoelectron spectra of these compounds. Analyses of the energetics of the global reactions to form cationic species show that metal cations bind more strongly to neutral cytosine than to neutral metals. Metal dications form the most stable compounds with neutral cytosine, and the stabilities of these systems decrease as (Zn-cyt)2+ > (Cd-cyt)2+ > (Ca-cyt)2+. Aromaticities computed via the HOMA indexes also support the observation regarding the greater affinity of cytosine for metal cations.