Journal of Physical Chemistry A, Vol.108, No.17, 3817-3828, 2004
Bader's and reactivity descriptors' analysis of DNA base pairs
The characterization of hydrogen bonds (H bonds) present in various canonical and noncanonical base pairs of DNA in terms of the electron density (rho(r(c))). Laplacian of electron density (del(2) rho(r(c))), and integrated atomic properties of hydrogen atoms involved in hydrogen bonding has been made. The existence of hydrogen bond critical points (HBCP) and the values of rho(r(c)) and del(2)rho(r(c)) at HBCP indicate the presence of hydrogen bonding and the closed-shell kind of interactions governing DNA base pairing. It has been observed that the calculated total rho(r(c)) at HBCPs of the base pairs varies linearly with their interaction energy. The total del(2)rho(r(c)) also exhibits a similar trend with the interaction energy. Integrated atomic properties such as charge, energy, volume, the first moment of the H-bonded hydrogen atoms in the DNA bases and in the H-bonded base pairs, and the change (Delta) in the corresponding properties arising upon pairing have also been calculated and compared. The calculated global chemical hardness of various DNA bases and DNA base pairs has been used to establish a linear relation with interaction energy. The charge transfer involved in the formation of DNA base pairs has also been qualitatively determined using Parr's formula and also from AIM-derived charges.