Feasibility study of hydrogen-bonded nucleic acid base pairs in gas and water phases — A theoretical study

To investigate the base pair binding probabilities for nucleic acid bases, numerous models were studied for contacts between adenine, thymine, guanine, cytosine, and uracil using density functional theory (DFT) in combination with the 6–311G* basis set. We obtained an assessment for the energy given by our calculations in gas and aqueous phases, which showed that it should be incorporated into hydrogen bonding and propeller rotational energies. The 42 complexes of base pairs (5 regular and 37 irregular base pairs) were proposed and their hydrogen-bonding (H-bonding) properties were verified. The hydrogen bonds in some irregular base pairs, including CC, UU, and TT (series 1), were stronger than in regular GC and AT base pairs. Also, the strength of the hydrogen bonds in the proposed base pairs, including CU, GG, GU, and TU (series 2), were similar to regular base pairs from an energetic point of view. The propeller rotations revealed a higher rotational barrier energy (6–7.5 kcal/mol; 1 cal = 4.184 J) for irregular base pairs (series 1 and 2) than regular GC and AT ones (1–3 kcal/mol). Nevertheless, the trend in these affinities of the complex contact probabilities and their biological properties were confirmed by our calculations.