The Effect of Water in Four Adenine-Thymine and Three Guanine-Cytosine Pairs: Combining Quantum and Statistical Mechanics

The molecular interactions between the nucleic acid bases and water molecules are important in organism. Despite Adenine-Thymine Hoogsteen base pair and Guanine-Cytosine Watson-Crick base pair have been demonstrated to be most stable in a gas phase, the effect of water on the stability of these base pairs remains elusive. Here we report the structural and thermodynamic characteristics on possible Adenine-Thymine and Guanine-Cytosine base pairs in a gas phase as well as in an aqueous phase by using quantum mechanical method and statistical mechanical calculations. First, we optimized the direct base-pair interaction energies of four Adenine-Thymine base pairs (Hoogsteen base pair, reverse Hoogsteen base pair, Watson-Crick base pair, and reverse Watson-Crick base pair) and three Guanine-Cytosine base pairs (GC1 base pair, GC2 base pair, and Watson Crick base pair) in a gas phase at the $B3LYP/6-31+G^{**}$ level. Then, the effect of solvent was quantified by the electronic reorganization energy and the solvation free energy by statistical mechanical calculations. Thereby, we discuss the effect of water on the stability of Adenine-Thymine and Guanine-Cytosine base pairs, and argue why Adenine-Thymine Watson-Crick base pair and Guanine-Cytosine Watson-Crick base pair are most stable in an aqueous environment.