A Course Grain Model of Histones and DNA in the Nucleosome

Abstract

There have been significant advances and interest in developing coarse grained models of DNA based on elastic rod theories. These models have focused largely on describing the physical-mechanical behavior of the DNA itself, without protein. Here we utilize coarse-grain models for both DNA and protein molecules to investigate the interaction of DNA with protein-complexes. We have focused our attention on the level of detail required in the coarse graining to achieve large scale deformations of the DNA. For this purpose we have chosen the nucleosome, a histone-DNA complex as our example system because (a) there is significant deformation of DNA in the bound state as compared to the free state (i.e., 1.7 turns of a superhelix), (b) the DNA in the complex is on the order of one persistence length, and (c) there are no sequence specific contacts between the protein and the DNA. These three features perhaps make this system a right candidate to be described by the physical-mechanical behavior of DNA.In our model of the nucleosome the DNA is represented, numerically, as a continuous elastic rod with a continuous charge distribution interacting with a multipolar representation of the electrostatic potential for the nucleosome's histone core. A contact (non-penetration) constraint is also included to account for the histone excluded volume. Using this model we investigate the physical basis of superhelical formation and the dependency of nucleosome stability on DNA sequence.