Matsuda and coworkers demonstrated that imidazopyridopyrimidine nucelobases (N N , O O , N O , tO O , and O N ) can mimic Watson–Crick nucleobase in forming H-bonds in DNA double helix. In the present study, we address the question about the strengths of the H-bonds in imidazopyridopyrimidine base pairs compared to those in Watson–Crick ones by focusing particularly on the nature of these interactions. Optimized structures of imidazopyridopyrimidine, imidazopyridopyrimidine–Watson–Crick, and Watson–Crick base pairs are obtained at the DFTB3LYP/6-311++G (d,p). The nature and strength of the intramolecular H-bonds in these base pairs have been investigated based on natural bond orbital (NBO method) to consider the effect of charge transfer, “atoms-in-molecules” (AIM) topological parameters, and decomposition of the interaction energies using the energy decomposition analysis (EDA). These investigations imply that N N –O O and N O –O N can form base pairs with four H-bonds (most stable than those of Watson–Crick base pairs) when they incorporated into DNA double helix. Furthermore, it can be deduced that O N and N N nucleobases form energetically more favorable pairs with adenine and guanine than the normal Watson–Crick counter parts. These results can be helpful for the stabilization and regulation of a variety of new base-pairing motif of DNA structures.