Multiscale modeling of nucleosome dynamics

Nucleosomes form the fundamental building blocks of chromatin. Subtle modifications of the constituent histone tails mediate chromatin stability and regulate gene expression. For this reason, it is important to understand structural dynamics of nucleosomes at atomic levels. We report a novel multiscale model of the fundamental chromatin unit, a nucleosome, using a simplified model for rapid discrete molecular dynamics simulations and an all-atom model for detailed structural investigation. Using a simplified structural model, we perform equilibrium simulations of a single nucleosome at various temperatures. We further reconstruct all-atom nucleosome structures from simulation trajectories. We find that histone tails bind to nucleosomal DNA via strong salt-bridge interactions over a wide range of temperatures, suggesting a mechanism of chromatin structural organization whereby histone tails regulate inter- and intranucleosomal assemblies via binding with nucleosomal DNA. We identify specific regions of the histone core H2A/H2B-H4/H3-H3/H4-H2B/H2A, termed “cold sites”, which retain a significant fraction of contacts with adjoining residues throughout the simulation, indicating their functional role in nucleosome organization. Cold sites are clustered around H3-H3, H2A-H4 and H4-H2A interhistone interfaces, indicating the necessity of these contacts for nucleosome stability. Essential dynamics analysis of simulation trajectories shows that bending across the H3-H3 is a prominent mode of intranucleosomal dynamics. We postulate that effects of salts on mononucleosomes can be modeled in discrete molecular dynamics by modulating histone-DNA interaction potentials. Local fluctuations in nucleosomal DNA vary significantly along the DNA sequence, suggesting that only a fraction of histone-DNA contacts make strong interactions dominating mononucleosomal dynamics. Our findings suggest that histone tails have a direct functional role in stabilizing higher-order chromatin structure, mediated by salt-bridge interactions with adjacent DNA.     

Rotational and translational positions determine the structural and dynamic impact of a single ribonucleotide incorporated in the nucleosome

Ribonucleotides misincorporated by replicative DNA polymerases are by far the most common DNA lesion. The presence of ribonucleotides in DNA is associated with genome instability, causing replication stress, chromosome fragility, gross chromosomal rearrangements, and other mutagenic events. Furthermore, nucleosome and chromatin assembly as well as nucleosome positioning are affected by the presence of ribonucleotides. Notably, nucleosome formation is significantly reduced by a single ribonucleotide. Single ribonucleotides are primarily removed from DNA by the ribonucleotide excision repair (RER) pathway via the RNase H2 enzyme, which incises the DNA backbone on the 5′-side of the ribonucleotide. While the structural implications of a single ribonucleotide in free duplex DNA have been well studied, how a single ribonucleotide embedded in nucleosomal DNA impacts nucleosome structure and dynamics, and the possible consequent impact on RER, have not been explored. We have carried out 3.5 μs molecular dynamics simulations of a single ribonucleotide incorporated at various translational and rotational positions in a nucleosome core particle. We find that the presence of the 2′−OH group on the ribose impacts the local conformation and dynamics of both the ribonucleotide and nearby DNA nucleotides as well as their interactions with histones; the nature of these disturbances depends on the rotational and translational setting, including whether the ribose faces toward or away from the histones. The ribonucleotide’s preferred C3′-endo pucker is stabilized by interactions with the histones, and furthermore the ribonucleotide can cause dynamic local duplex disturbance involving an abnormal C3′-endo population of the adjacent deoxyribose pucker, minor groove opening, ruptured Watson-Crick pairing, and duplex unwinding that are governed by translation-dependent histone-nucleotide interactions. Possible effects of these disturbances on RER are considered.     

Crystal structures of nucleosome core particles in complex with minor groove DNA-binding ligands

We determined the crystal structures of three nucleosome core particles in complex with site-specific DNA-binding ligands, the pyrrole-imidazole polyamides. While the structure of the histone octamer and its interaction with the DNA remain unaffected by ligand binding, nucleosomal DNA undergoes significant structural changes at the ligand-binding sites and in adjacent regions to accommodate the ligands. Our findings suggest that twist diffusion occurs over long distances through tightly bound nucleosomal DNA. This may be relevant to the mechanism of ATP-dependent and spontaneous nucleosome translocation, and to the effect of bound factors on nucleosome dynamics.     

Histone tail network and modulation in a nucleosome

The structural unit of eukaryotic chromatin is a nucleosome, comprising two histone H2A/H2B heterodimers and one histone (H3/H4)₂ tetramer, wrapped around by ∼146-bp core DNA and linker DNA. Flexible histone tails sticking out from the core undergo posttranslational modifications that are responsible for various epigenetic functions. Recently, the functional dynamics of histone tails and their modulation within the nucleosome and nucleosomal complexes have been investigated by integrating NMR, molecular dynamics simulations, and cryo-electron microscopy approaches. In particular, recent NMR studies have revealed correlations in the structures of histone N-terminal tails between H2A and H2B, as well as between H3 and H4 depending on linker DNA, suggesting that histone tail networks exist even within the nucleosome.     

A Mechanism for Histone Chaperoning Activity of Nucleoplasmin: Thermodynamic and Structural Models

Nucleoplasmin (NP), a histone chaperone, acts as a reservoir for histones H2A-H2B in Xenopus laevis eggs and can displace sperm nuclear basic proteins and linker histones from the chromatin fiber of sperm and quiescent somatic nuclei. NP has been proposed to mediate the dynamic exchange of histones during the expression of certain genes and assists the assembly of nucleosomes by modulating the interaction between histones and DNA. Here, solution structural models of full-length NP and NP complexes with the functionally distinct nucleosomal core and linker histones are presented for the first time, providing a picture of the physical interactions between the nucleosomal and linker histones with NP core and tail domains. Small-angle X-ray scattering and isothermal titration calorimetry reveal that NP pentamer can accommodate five histones, either H2A-H2B dimers or H5, and that NP core and tail domains are intimately involved in the association with histones. The analysis of the binding events, employing a site-specific cooperative model, reveals a negative cooperativity-based regulatory mechanism for the linker histone/nucleosomal histone exchange. The two histone types bind with drastically different intrinsic affinity, and the strongest affinity is observed for the NP variant that mimicks the hyperphosphorylated active protein. The different “affinity windows” for H5 and H2A-H2B might allow NP to fulfill its histone chaperone role, simultaneously acting as a reservoir for the core histones and a chromatin decondensing factor. Our data are compatible with the previously proposed model where NP facilitates nucleosome assembly by removing the linker histones and depositing H2A-H2B dimers onto DNA.     

Characterization of Dnmt1 Binding and DNA Methylation on Nucleosomes and Nucleosomal Arrays

The packaging of DNA into nucleosomes and the organisation into higher order structures of chromatin limits the access of sequence specific DNA binding factors to DNA. In cells, DNA methylation is preferentially occuring in the linker region of nucleosomes, suggesting a structural impact of chromatin on DNA methylation. These observations raise the question whether DNA methyltransferases are capable to recognize the nucleosomal substrates and to modify the packaged DNA. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the maintenance DNA methyltransferase Dnmt1. Our binding studies show that Dnmt1 has a DNA length sensing activity, binding cooperatively to DNA, and requiring a minimal DNA length of 20 bp. Dnmt1 needs linker DNA to bind to nucleosomes and most efficiently recognizes nucleosomes with symmetric DNA linkers. Footprinting experiments reveal that Dnmt1 binds to both DNA linkers exiting the nucleosome core. The binding pattern correlates with the efficient methylation of DNA linkers. However, the enzyme lacks the ability to methylate nucleosomal CpG sites on mononucleosomes and nucleosomal arrays, unless chromatin remodeling enzymes create a dynamic chromatin state. In addition, our results show that Dnmt1 functionally interacts with specific chromatin remodeling enzymes to enable complete methylation of hemi-methylated DNA in chromatin.     

Breaths,Twists, and Turns of Atomistic Nucleosomes

Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145–147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.     

DNA sequence-dependent contributions of core histone tails to nucleosome stability: differential effects of acetylation and proteolytic tail removal

Modulation of nucleosome stability in chromatin plays an important role in eukaryotic gene expression. The core histone N-terminal tail domains are believed to modulate the stability of wrapping nucleosomal DNA and the stability of the chromatin filament. We analyzed the contribution of the tail domains to the stability of nucleosomes containing selected DNA sequences that are intrinsically straight, curved, flexible, or inflexible. We find that the presence of the histone tail domains stabilizes nucleosomes containing DNA sequences that are intrinsically straight or curved. However, the tails do not significantly contribute to the free energy of nucleosome formation with flexible DNA. Interestingly, hyperacetylation of the core histone tail domains does not recapitulate the effect of tail removal by limited proteolysis with regard to nucleosome stability. We find that acetylation of the tails has the same minor effect on nucleosome stability for all the selected DNA sequences. A comparison of histone partitioning between long donor chromatin, acceptor DNA, and free histones in solution shows that the core histone tails mediate internucleosomal interactions within an H1-depleted chromatin fiber amounting to an average free energy of about 1 kcal/mol. Thus, such interactions would be significant with regard to the free energies of sequence-dependent nucleosome positioning. Last, we analyzed the contribution of the H2A/H2B dimers to nucleosome stability. We find that the intact nucleosome is stabilized by 900 cal/mol by the presence of the dimers regardless of sequence. The biological implications of these observations are discussed.     

The use of Quantitative Agarose Gel Electrophoresis for rapid analysis of the integrity of protein–DNA complexes

Recent biochemical studies evaluated the affinity of histones to DNA in the context of nucleosome core particle (NCP). These have indicated a concentration-dependence for nucleosome stability. However, when studying chromatin the preferred templates are nucleosome arrays (NA) and not the NCP. Biochemical methods are poorly suited for structural analysis of chromatin. To overcome that technical hindrance, and investigate the effect of concentration on stability of the histone–DNA interactions, we have applied the multigel Quantitative Agarose Gel Electrophoresis (QAGE) method to in vitro-assembled nucleosomal arrays. The results demonstrated the method to be extremely valuable for the evaluation of the effect of low concentration on NA. However, QAGE is a fairly time-demanding and complex method. To maximize the efficiency of use of this technology, we devised a protocol that allowed for multiple sets of templates to be analyzed simultaneously. Briefly, samples can be loaded at regular intervals and analyzed individually for their molecular composition. The technique presented in this study describes the calibration steps and proof of concept necessary to validate the use of multiple loading of multigel to evaluate the composition of nucleosomal arrays as a function of concentration.     

Strategies for crystallizing a chromatin protein in complex with the nucleosome core particle

The molecular details of how chromatin factors and enzymes interact with the nucleosome are critical to understanding fundamental genetic processes including cell division and gene regulation. A structural understanding of such processes has been hindered by the difficulty in producing diffraction-quality crystals of chromatin proteins in complex with the nucleosome. We describe here the steps used to grow crystals of the 300-kDa RCC1 chromatin factor/nucleosome core particle complex that diffract to 2.9-Å resolution. These steps include both pre- and postcrystallization strategies potentially useful to other complexes. We screened multiple variant RCC1/nucleosome core particle complexes assembled using different RCC1 homologs and deletion variants, and nucleosomes containing nucleosomal DNA with different sequences and lengths, as well as histone deletion variants. We found that using RCC1 from different species produced different crystal forms of the RCC1/nucleosome complex consistent with key crystal packing interactions mediated by RCC1. Optimization of postcrystallization soaks to dehydrate the crystals dramatically improved the diffraction quality of the RCC1/nucleosome crystal from 5.0- to 2.9-Å resolution.     

Laser Raman spectra of calf thymus chromatin and its constituents.

Extensive Raman measurements have been made on calf thymus chromatin, core chromatin, the (H3,H4)/DNA complex, and isolated DNA. The results indicate that the alpha-helical content of the nucleosomal histones gradually increases as they form the heterocomplexes that lead to the formation of the octameric nucleosome core. The secondary structure of the latter is not modified as it binds to DNA. The spectra indicate that the DNA essentially retains its B conformation in nucleosomes, although slight changes probably occur in the ribose-phosphate backbone. No specific interactions between the nucleosomal histones and DNA can be established from the spectra, but histone H1 possibly interacts selectively with the thymine bases.     

Studies on chromatin

A new model for the fine structure of the chromatin subunit (or 'nucleosome') is proposed. The model is based on previous experimental findings [1-14] and on two new suggestions, namely: (1) Eight histones form a toroidal-shaped histone coe of nucleosome and are arranged in the following ciruclar sequence: (see article). (2) DNA is 'kinked' around a toroidal-shaped histone core in a 'solenoid-like' mode, each kink occurring every 10 base pairs along DNA. The electron microscopic evidence for a toroidal shape of the nucleosome is described in the preceding paper [13]. The possibility of the existence of kinks in the DNA double helix was considered recently by Crick and Klug [14]. The proposed model of the nucleosome, being more detailed than earlier models permits us to explain in direct structural terms the yet unordered set of data bearing on the pattern of histone-histone interactions in chromatin, the results of a mild deoxyribonuclease digestion of DNA within the nucleosomal particle and also the quantitative data on the unwinding of the DNA duplex upon formation of the nucleosome.     

The use of Quantitative Agarose Gel Electrophoresis for rapid analysis of the integrity of protein-DNA complexes

Recent biochemical studies evaluated the affinity of histones to DNA in the context of nucleosome core particle (NCP). These have indicated a concentration-dependence for nucleosome stability. However, when studying chromatin the preferred templates are nucleosome arrays (NA) and not the NCP. Biochemical methods are poorly suited for structural analysis of chromatin. To overcome that technical hindrance, and investigate the effect of concentration on stability of the histone-DNA interactions, we have applied the multigel Quantitative Agarose Gel Electrophoresis (QAGE) method to in vitro-assembled nucleosomal arrays. The results demonstrated the method to be extremely valuable for the evaluation of the effect of low concentration on NA. However, QAGE is a fairly time-demanding and complex method. To maximize the efficiency of use of this technology, we devised a protocol that allowed for multiple sets of templates to be analyzed simultaneously. Briefly, samples can be loaded at regular intervals and analyzed individually for their molecular composition. The technique presented in this study describes the calibration steps and proof of concept necessary to validate the use of multiple loading of multigel to evaluate the composition of nucleosomal arrays as a function of concentration.     

Salt-dependent intra- and internucleosomal interactions of the H3 tail domain in a model oligonucleosomal array

The core histone tail domains are known to be key regulators of chromatin structure and function. The tails are required for condensation of nucleosome arrays into secondary and tertiary chromatin structures, yet little is known regarding tail structures or sites of tail interactions in chromatin. We have developed a system to test the hypothesis that the tails participate in internucleosomal interactions during salt-dependent chromatin condensation, and here we used it to examine interactions of the H3 tail domain. We found that the H3 tail participates primarily in intranucleosome interactions when the nucleosome array exists in an extended "beads-on-a-string" conformation and that tail interactions reorganize to engage in primarily internucleosome interactions as the array successively undergoes salt-dependent folding and oligomerization. These results indicated that the location and interactions of the H3 tail domain are dependent upon the degree of condensation of the nucleosomal array, suggesting a mechanism by which alterations in tail interactions may elaborate different structural and functional states of chromatin.     

Spontaneous access to DNA target sites in folded chromatin fibers

DNA wrapped in nucleosomes is sterically occluded from many protein complexes that must act on it; how such complexes gain access to nucleosomal DNA is not known. In vitro studies on isolated nucleosomes show that they undergo spontaneous partial unwrapping conformational transitions, which make the wrapped nucleosomal DNA transiently accessible. Thus, site exposure might provide a general mechanism allowing access of protein complexes to nucleosomal DNA. However, existing quantitative analyses of site exposure focused on single nucleosomes, while the presence of neighbor nucleosomes and concomitant chromatin folding might significantly influence site exposure. In this work, we carried out quantitative studies on the accessibility of nucleosomal DNA in homogeneous nucleosome arrays. Two striking findings emerged. Organization into chromatin fibers changes the accessibility of nucleosomal DNA only modestly, from ∼ 3-fold decreases to ∼ 8-fold increases in accessibility. This means that nucleosome arrays are intrinsically dynamic and accessible even when they are visibly condensed. In contrast, chromatin folding decreases the accessibility of linker DNA by as much as ∼ 50-fold. Thus, nucleosome positioning dramatically influences the accessibility of target sites located inside nucleosomes, while chromatin folding dramatically regulates access to target sites in linker DNA.     

Principles of nucleosome recognition by chromatin factors and enzymes

The recent torrent of structures of chromatin complexes determined by cryoelectron microscopy provides an opportunity to discern general principles for how chromatin factors and enzymes interact with their nucleosome substrate. We find that many chromatin proteins use a strikingly similar arginine anchor and variant arginine interactions to bind to the nucleosome acidic patch. We also observe that many chromatin proteins target the H3 and H2B histone fold α1-loop1 elbows and the H2B C-terminal helix on the nucleosomal histone face. These interactions with the histones can be complemented with interactions with and distortions of nucleosomal DNA.