In silicio stretching of chromatin

We present Monte-Carlo (MC) simulations of the stretching of a single 30 nm chromatin fiber. The model approximates the DNA by a flexible polymer chain with Debye-H\"uckel electrostatics and uses a two-angle zig-zag model for the geometry of the linker DNA connecting the nucleosomes. The latter are represented by flat disks interacting via an attractive Gay-Berne potential. Our results show that the stiffness of the chromatin fiber strongly depends on the linker DNA length. Furthermore, changing the twisting angle between nucleosomes from 90 deg to 130 deg increases the stiffness significantly. An increase in the opening angle from 22 deg to 34 deg leads to softer fibers for small linker lengths. We observe that fibers containing a linker histone at each nucleosome are stiffer compared to those without the linker histone. The simulated persistence lengths and elastic moduli agree with experimental data. Finally, we show that the chromatin fiber does not behave as an isotropic elastic rod, but its rigidity depends on the direction of deformation: chromatin is much more resistant to stretching than to bending.