In silico studies to explore the mutagenic ability of 5-halo/oxy/li-oxy-uracil bases with guanine of DNA base pairs.

DNA nucleobases are reactive in nature and undergo modifications by deamination, oxidation, alkylation, or hydrolysis processes. Many such modified bases are susceptible to mutagenesis when formed in cellular DNA. The mutagenesis can occur by mispairing with DNA nucleobases by a DNA polymerase during replication. We have performed a study of mispairing of DNA bases with unnatural bases computationally. 5-Halo uracils have been studied as mispairs in mutagenesis; however, the reports on their different forms are scarce in the literature. The stability of mispairs with keto form, enol form, and ionized form of 5-halo-uracil has been computed with the M06-2X/6-31+G** level of theory. The enol form of 5-halo-uracil showed remarkable stability toward DNA mispair compared to the corresponding keto and ionized forms. (F)U-G mispair showed the highest stability in the series and (Halo)(U(enol/ionized)-G mispair interactions energies are more stable than the natural G-C basepair of DNA. To enhance the stability of DNA mispairs, we have introduced the hydroxyl group in the place of halogen atoms, which provides additional hydrogen-bonding interactions in the system while forming the 5-membered ring. The study has been further extended with lithiated 5-hydroxymethyl-uracil to stabilize the DNA mispair. (CH2OLi)U(ionized)-G mispair has shown the highest stability (ΔG = -32.4 kcal/mol) with multi O-Li interactions. AIM (atoms in molecules) and EDA (energy decomposition analysis) analysis has been performed to examine the nature of noncovalent interactions in such mispairs. EDA analysis has shown that electrostatic energy mainly contributes toward the interaction energy of mispairs. The higher stability achieved in these studied mispairs can play a pivotal role in the mutagenesis and can help to attain the mutation for many desired biological processes.