Sequence-dependent Kink-and-Slide Deformations of Nucleosomal DNA Facilitated by Histone Arginines Bound in the Minor Groove

Abstract

In addition to bending and twisting deformabilities, the lateral displacements of the DNA axis (Kink-and-Slide) play an important role in DNA wrapping around the histone core (M. Y. Tolstorukov, A. V. Colasanti, D. M. McCandlish, W. K. Olson, V. B. Zhurkin, J. Mol. Biol. 371, 725-738 (2007)). Here, we show that these Kink-and-Slide deformations are likely to be stabilized by the arginine residues of histones interacting with the minor groove of DNA. The arginines are positioned asymmetrically in the minor groove, being closer to one strand. The asymmetric arginine-DNA interactions facilitate lateral displacement of base pairs across the DNA grooves, thus leading to a stepwise accumulation of the superhelical pitch of nucleosomal DNA. To understand the sequence dependence of such Kink-and-Slide deformations, we performed all-atom calculations of DNA hexamers with the YR and RY steps in the center. We found that when the unrestrained DNA deformations are allowed, the YR steps tend to bend into the major groove, and RY steps bend into the minor groove. However, when the nucleosomal Kink-and-Slide deformation is considered, the YR steps prove to be more favorable for bending into the minor groove. Overall, the Kink-and-Slide deformation energy of DNA increases in the order TA < CA < CG < GC < AC < AT. We propose a simple stereochemical model accounting for this sequence dependence. Our results agree with experimental data indicating that the TA step most frequently occurs in the minor-groove kink positions in the most stable nucleosomes. Our computations demonstrate that the Kink-and-Slide distortion is accompanied by the BI to BII transition. This fact, together with irregularities in the two-dimensional (Roll, Slide) energy contour maps, suggest that the Kink-and-Slide deformations represent a nonharmonic behavior of the duplex. This explains the difference between the two estimates of the DNA deformation energy in nucleosome - the earlier one made using knowledge-based elastic energy functions, and the current one based on all-atom calculations. Our findings are useful for refining the score functions for the prediction of nucleosome positioning. In addition, the reverse bending behavior of the YR and RY steps revealed under the Kink-and-Slide constraint is important for understanding the molecular mechanisms of binding transcription factors (such as p53) to DNA exposed on the surface of nucleosome.