Comprehensive analysis unveils altered binding kinetics of 5-/6-methylCytosine/adenine modifications in R2R3-DNA system

Abstract

Recent studies have shown that DNA methylation is an important epigenetic marker. Two prominent forms are methylation of the C5 position of cytosine and methylation of the C6 position of adenine. Given the vital significance of DNA methylation, investigating the mechanisms that influence protein binding remains a compelling pursuit. This study used molecular dynamics simulations to investigate the binding patterns of R2R3 protein and four differentially methylated DNAs. The alanine scanning combined with interaction entropy method was used to identify key residues that respond to different methylation patterns. The order of protein binding ability to DNA is as follows: unmethylated DNA > A11 methylation (5′-A6mAC-3′) (6m2A system) > A10 methylation (5′-6mAAC-3′) (6m1A system) > both A10 and A11 methylation (5′-6mA6mAC-3′) (6mAA system) > C12 methylation (5′-AA5mC-3′) (5mC system). All methylation systems lead to the sixth α helix (H6) (residues D105 to L116) moving away from the binding interface, and in the 5mC and 6m1A systems, the third α helix (H3) (residues G54 to L65) exhibits a similar trend. When the positively charged amino acids in H3 and H6 move away from the binding interface, their electrostatic and van der Waals interactions with the negatively charged DNA are weakened. Structural changes induced by methylation contributed to the destabilization of the hydrogen bond network near the original binding site, except for the 6m2A system. Moreover, there is a positive correlation between the number of methylated sites and the probability of distorting the DNA structure. Our study explores how different methylation patterns affect binding and structural adaptability, and have implications for drug discovery and understanding diseases related to abnormal methylation.