Computer simulation of the 30-nanometer chromatin fiber

A new Monte Carlo model for the structure of chromatin is presented here. Based on our previous work on superhelical DNA and polynucleosomes, it reintegrates aspects of the "solenoid" and the "zig-zag" models. The DNA is modeled as a flexible elastic polymer chain, consisting of segments connected by elastic bending, torsional, and stretching springs. The electrostatic interaction between the DNA segments is described by the Debye-Huckel approximation. Nucleosome core particles are represented by oblate ellipsoids; their interaction potential has been parameterized by a comparison with data from liquid crystals of nucleosome solutions. DNA and chromatosomes are linked either at the surface of the chromatosome or through a rigid nucleosome stem. Equilibrium ensembles of 100-nucleosome chains at physiological ionic strength were generated by a Metropolis-Monte Carlo algorithm. For a DNA linked at the nucleosome stem and a nucleosome repeat of 200 bp, the simulated fiber diameter of 32 nm and the mass density of 6.1 nucleosomes per 11 nm fiber length are in excellent agreement with experimental values from the literature. The experimental value of the inclination of DNA and nucleosomes to the fiber axis could also be reproduced. Whereas the linker DNA connects chromatosomes on opposite sides of the fiber, the overall packing of the nucleosomes leads to a helical aspect of the structure. The persistence length of the simulated fibers is 265 nm. For more random fibers where the tilt angles between two nucleosomes are chosen according to a Gaussian distribution along the fiber, the persistence length decreases to 30 nm with increasing width of the distribution, whereas the other observable parameters such as the mass density remain unchanged. Polynucleosomes with repeat lengths of 212 bp also form fibers with the expected experimental properties. Systems with larger repeat length form fibers, but the mass density is significantly lower than the measured value. The theoretical characteristics of a fiber with a repeat length of 192 bp where DNA and nucleosomes are connected at the core particle are in agreement with the experimental values. Systems without a stem and a repeat length of 217 bp do not form fibers.