Protonation of thymine, cytosine, adenine, and guanine DNA nucleic acid bases: Theoretical investigation into the framework of density functional theory

Gradient-corrected density functional computations with triple-zeta-type basis sets were performed to determine the preferred protonation site and the absolute gas-phase proton affinities of the most stable tautomer of the DNA bases thymine (T), cytosine (C), adenine (A), and guanine (G). Charge distribution, bond orders, and molecular electrostatic potentials were considered to rationalize the obtained results. The vibrational frequencies and the contribution of the zero-point energies were also computed. Significant geometrical changes in bond lengths and angles near the protonation sites were found. At 298 K, proton affinities values of 208.8 (T), 229.1 (C), 225.8 (A), and 230.3 (G) kcal/mol were obtained in agreement with experimental results.