The Role of Aromaticity, Hybridization, Electrostatics, and Covalency in Resonance-Assisted Hydrogen Bonds of Adenine-Thymine (AT) Base Pairs and Their Mimics

Hydrogen bonds play a crucial role in many biochemical processes and in supramolecular chemistry. In this study, we show quantum chemically that neither aromaticity nor other forms of assistance are responsible for the enhanced stability of the hydrogen bonds in adenine-thymine (AT) DNA base pairs. This follows from extensive bonding analyses of AT and smaller analogs thereof, based on dispersion-corrected density functional theory (DFT). Removing the aromatic rings of either A or T has no effect on the Watson-Crick bond strength. Only when the smaller mimics become saturated, that is, when the hydrogen-bond acceptor and donor groups go from sp(2) to sp(3), does the stability of the resulting model complexes suddenly drop. Bonding analyses based on quantitative Kohn-Sham molecular orbital theory and corresponding energy decomposition analyses (EDA) show that the stronger hydrogen bonds in the unsaturated model complexes and in AT stem from stronger electrostatic interactions as well as enhanced donor-acceptor interactions in the sigma-electron system, with the covalency being responsible for shortening the hydrogen bonds in these dimers.