Unravelling Protein–DNA Interactions at Molecular Level: A DFT and NCI Study

Histone–DNA interactions were probed computationally at a molecular level, by characterizing the bimolecular clusters constituted by selected amino acid derivatives with polar (asparagine and glutamine), nonpolar (alanine, valine, and isoleucine), and charged (arginine) side chains and methylated pyrimidinic (1-methylcytosine and 1-methylthymine) and puric (9-methyladenine and 9-methylguanine) DNA bases. The computational approach combined different methodologies: a molecular mechanics (MMFFs forced field) conformational search and structural and vibrational density-functional calculations (M06-2X with double and triple-ζ Pople’s basis sets). To dissect the interactions, intermolecular forces were analyzed with the Non-Covalent Interactions (NCI) analysis. The results for the 24 different clusters studied show a noticeable correlation between the calculated binding energies and the propensities for protein–DNA base interactions found in the literature. Such correlation holded even for the interaction of the selected amino acid derivatives with Watson and Crick pairs. Therefore, the balance between hydrogen bonds and van der Waals interactions (specially stacking) in the control of the final shape of the investigated amino acid–DNA base pairs seems to be well reproduced in dispersion-corrected DFT molecular models, reinforcing the idea that the specificity between the amino acids and the DNA bases play an important role in the regulation of DNA.