DNA methylation cues in nucleosome geometry, stability, and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported in both prokaryotes and eukaryotes, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here we investigated the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations (five microsecond long trajectories, 60 microseconds in total), we generated extensive atomic level conformational full nucleosome ensembles. Our results revealed that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, shifting the population equilibrium of sugar-phosphate backbone geometry. These conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microsecond time scale.