Sequence-dependent nucleosome sliding in rotation-coupled and uncoupled modes revealed by molecular simulations

Abstract While nucleosome positioning on eukaryotic genome play important roles for genetic regulation, molecular mechanisms of nucleosome positioning and sliding along DNA are not well understood. Here we investigated thermally-activated spontaneous nucleosome sliding mechanisms developing and applying a coarse-grained molecular simulation method that incorporates both long-range electrostatic and short-range hydrogen-bond interactions between histone octamer and DNA. The simulations revealed two distinct sliding modes depending on the nucleosomal DNA sequence. A uniform DNA sequence showed frequent sliding with one base pair step in a rotation-coupled manner, akin to screw-like motions. On the contrary, a strong positioning sequence, the so-called 601 sequence, exhibits rare, abrupt transitions of five and ten base pair steps without rotation. Moreover, we evaluated the importance of hydrogen bond interactions on the sliding mode, finding that strong and weak bonds favor respectively the rotation-coupled and -uncoupled sliding movements. Author summary Nucleosomes are fundamental units of chromatin folding consisting of double-stranded DNA wrapped ∼1.7 times around a histone octamer. By densely populating the eukaryotic genome, nucleosomes enable efficient genome compaction inside the cellular nucleus. However, the portion of DNA occupied by a nucleosome can hardly be accessed by other DNA-binding proteins, obstructing fundamental cellular processes such as DNA replication and transcription. DNA compaction and access by other proteins can simultaneously be achieved via the dynamical repositioning of nucleosomes, which can slide along the DNA sequence. In this study, we developed and used coarse-grained molecular dynamics simulations to reveal the molecular details of nucleosome sliding. We find that the sliding mode is highly dependent on the underlying DNA sequence. Specifically, a sequence with a strong nucleosome positioning signal slides via large jumps by five and ten base pairs, preserving the optimal DNA bending profile. On the other hand, uniform sequences without the positioning signal slide via a screw-like motion of DNA, one base pair at the time. These results show that sequence has a large effect not only on the formation of nucleosomes, but also on the kinetics of repositioning.