Dissecting DNA-histone interactions in the nucleosome by molecular dynamics simulations of DNA unwrapping

The nucleosome complex of DNA wrapped around a histone protein octamer organizes the genome of eukaryotes and regulates the access of protein factors to the DNA. We performed molecular dynamics simulations of the nucleosome in explicit water to study the dynamics of its histone-DNA interactions. A high-resolution histone-DNA interaction map was derived that revealed a five-nucleotide periodicity, in which the two DNA strands of the double helix made alternating contacts. On the 100-ns timescale, the histone tails mostly maintained their initial positions relative to the DNA, and the spontaneous unwrapping of DNA was limited to 1–2 basepairs. In steered molecular dynamics simulations, external forces were applied to the linker DNA to investigate the unwrapping pathway of the nucleosomal DNA. In comparison with a nucleosome without the unstructured N-terminal histone tails, the following findings were obtained: 1), Two main barriers during unwrapping were identified at DNA position ±70 and ±45 basepairs relative to the central DNA basepair at the dyad axis. 2), DNA interactions of the histone H3 N-terminus and the histone H2A C-terminus opposed the initiation of unwrapping. 3), The N-terminal tails of H2A, H2B, and H4 counteracted the unwrapping process at later stages and were essential determinants of nucleosome dynamics. Our detailed analysis of DNA-histone interactions revealed molecular mechanisms for modulating access to nucleosomal DNA via conformational rearrangements of its structure.     

Dynamics of Forced Nucleosome Unraveling and Role of Nonuniform Histone-DNA Interactions

A coarse-grained model of the nucleosome is introduced to investigate the dynamics of force-induced unwrapping of DNA from histone octamers. In this model, the DNA is treated as a charged, discrete worm-like chain, and the octamer is treated as a rigid cylinder carrying a positively charged superhelical groove that accommodates 1.7 turns of DNA. The groove charges are parameterized to reproduce the nonuniform histone/DNA interaction free energy profile and the loading rate-dependent unwrapping forces, both obtained from single-molecule experiments. Brownian dynamics simulations of the model under constant loading conditions reveal that nucleosome unraveling occurs in three distinct stages. At small extensions, the flanking DNA exhibits rapid unwrapping-rewrapping (breathing) dynamics and the octamer flips ～180° and moves toward the pulling axis. At intermediate extensions, the outer turn of DNA unwraps gradually and the octamer swivels about the taut linkers and flips a further ～90° to orient its superhelical axis almost parallel to the pulling axis. At large extensions, a portion of the inner turn unwraps abruptly with a notable rip in the force-extension plot and a >90° flip of the octamer. The remaining inner turn unwraps reversibly to leave a small portion of DNA attached to the octamer despite extended pulling. Our simulations further reveal that the nonuniform histone/DNA interactions in canonical nucleosomes serve to: stabilize the inner turn against unraveling while enhancing the breathing dynamics of the nucleosome and prevent dissociation of the octamer from the DNA while facilitating its mobility along the DNA. Thus, the modulation of the histone/DNA interactions could constitute one possible mechanism for regulating the accessibility of the nucleosome-wound DNA sequences.     

The effect of epigenetic modifications on the secondary structures and possible binding positions of the N-terminal tail of histone H3 in the nucleosome: a computational study

The roles of histone tails as substrates for reversible chemical modifications and dynamic cognate surfaces for the binding of regulatory proteins are well established. Despite these crucial roles, experimentally derived knowledge of the structure and possible binding sites of histone tails in chromatin is limited. In this study, we utilized molecular dynamics of isolated histone H3 N-terminal peptides to investigate its structure as a function of post-translational modifications that are known to be associated with defined chromatin states. We observed a structural preference for α-helices in isoforms associated with an inactive chromatin state, while isoforms associated with active chromatin states lacked α-helical content. The physicochemical effect of the post-translational modifications was highlighted by the interaction of arginine side-chains with the phosphorylated serine residues in the inactive isoform. We also showed that the isoforms exhibit different tail lengths, and, using molecular docking of the first 15 N-terminal residues of an H3 isoform, identified potential binding sites between the superhelical gyres on the octamer surface, close to the site of DNA entry/exit in the nucleosome. We discuss the possible functional role of the binding of the H3 tail within the nucleosome on both nucleosome and chromatin structure and stability.     

Sequence-dependent nucleosome structural and dynamic polymorphism. Potential involvement of histone H2B N-terminal tail proximal domain

Relaxation of nucleosomes on an homologous series (pBR) of ca 350-370 bp DNA minicircles originating from plasmid pBR322 was recently used as a tool to study their structure and dynamics. These nucleosomes thermally fluctuated between three distinct DNA conformations within a histone N-terminal tail-modulated equilibrium: one conformation was canonical, with 1.75 turn wrapping and negatively crossed entering and exiting DNAs; another was also "closed", but with these DNAs positively crossed; and the third was "open", with a lower than 1.5 turn wrapping and uncrossed DNAs. In this work, a new minicircle series (5S) of similar size was used, which contained the 5S nucleosome positioning sequence. Results showed that DNA in pBR nucleosomes was untwisted by approximately 0.2 turn relative to 5S nucleosomes, which DNase I footprinting confirmed in revealing a approximately 1 bp untwisting at each of the two dyad-distal sites where H2B N-terminal tails pass between the two gyres. In contrast, both nucleosomes showed untwistings at the dyad-proximal sites, i.e. on the other gyre, which were also observed in the high-resolution crystal structure. 5S nucleosomes also differ with respect to their dynamics: they hardly accessed the positively crossed conformation, but had an easier access to the negatively crossed conformation. Simulation showed that such reverse effects on the conformational free energies could be simply achieved by slightly altering the trajectories of entering and exiting DNAs. We propose that this is accomplished by H2B tail untwisting at the distal sites through action at a distance ( approximately 20 bp) on H3-tail interactions with the small groove at the nucleosome entry-exit. These results may help to gain a first glimpse into the two perhaps most intriguing features of the high-resolution structure: the alignment of the grooves on the two gyres and the passage of H2B and H3 N-terminal tails between them.     

Structure of the nucleosome core particle at 8 A resolution

The x-ray crystallographic structure of the nucleosome core particle has been determined using 8 A resolution diffraction data. The particle has a mean diameter of 106 A and a maximum thickness of 65 A in the superhelical axis direction. The longest chord through the histone core measures 85 A and is in a non-axial direction. The 1.87 turn superhelix consists of B-DNA with about 78 base pairs or 7.6 helical repeats per superhelical turn. The mean DNA helical repeat contains 10.2 +/- 0.05 base pairs and spans 35 A, slightly more than standard B-DNA. The superhelix varies several Angstroms in radius and pitch, and has three distinct domains of curvature (with radii of curvature of 60, 45 and 51 A). These regions are separated by localized sharper bends +/- 10 and +/- 40 base pairs from the center of the particle, resulting in an overall radius of curvature about 43 A. Compression of superhelical DNA grooves on the inner surface and expansion on the outer surface can be seen throughout the DNA electron density. This density has been fit with a double helical ribbon model providing groove width estimates of 12 +/- 1 A inside vs. 19 +/- 1 A outside for the major groove, and 8 +/- 1 A inside vs. 13 +/- 1 A outside for the minor groove. The histone core is primarily contained within the bounds defined by the superhelical DNA, contacting the DNA where the phosphate backbone faces in toward the core. Possible extensions of density between the gyres have been located, but these are below the significance level of the electron density map. In cross-section, a tripartite organization of the histone octamer is apparent, with the tetramer occupying the central region and the dimers at the extremes. Several extensions of histone density are present which form contacts between nucleosomes in the crystal, perhaps representing flexible or "tail" histone regions. The radius of gyration of the histone portion of the electron density is calculated to be 30.4 A (in reasonable agreement with solution scattering values), and the histone core volume in the map is 93% of its theoretical volume.     

DNA methylation cues in nucleosome geometry,stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3–4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.     

Nucleosome structural dynamics and super spiral linking number paradox

The minireview presents recent results obtained in the experiments with DNA minicircles containing reconstituted nucleosomes. This system allows one to register and to characterize conformational dynamics of nucleosomes and of subnucleosomal particles containing histone tetramer (H3-H4)2. In particular, it has revealed an important role of the histone N-terminal tails in this dynamics. Solution of the linking number paradox and relevance of the results obtained to chromatin structural dynamics are discussed.     

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.     

Nucleosome Structural Dynamics and Linking Number Paradox

The minireview presents recent results obtained in the experiments with DNA minicircles containing reconstituted nucleosomes. This system allows one to register and to characterize the conformational dynamics of nucleosomes and of subnucleosomal particles containing histone tetramer (H3-H4)₂. In particular, it has revealed an important role of the histone N-terminal tails in this dynamics. Solution of the linking number paradox and the relevance of the results obtained to chromatin structural dynamics are discussed.     

Nucleosome composition regulates the histone H3 tail conformational ensemble and accessibility

Hexasomes and tetrasomes are intermediates in nucleosome assembly and disassembly. Their formation is promoted by histone chaperones, ATP-dependent remodelers, and RNA polymerase II. In addition, hexasomes are maintained in transcribed genes and could be an important regulatory factor. While nucleosome composition has been shown to affect the structure and accessibility of DNA, its influence on histone tails is largely unknown. Here, we investigate the conformational dynamics of the H3 tail in the hexasome and tetrasome. Using a combination of NMR spectroscopy, MD simulations, and trypsin proteolysis, we find that the conformational ensemble of the H3 tail is regulated by nucleosome composition. As has been found for the nucleosome, the H3 tails bind robustly to DNA within the hexasome and tetrasome, but upon loss of the H2A/H2B dimer, we determined that the adjacent H3 tail has an altered conformational ensemble, increase in dynamics, and increase in accessibility. Similar to observations of DNA dynamics, this is seen to be asymmetric in the hexasome. Our results indicate that nucleosome composition has the potential to regulate chromatin signaling and ultimately help shape the chromatin landscape.     

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.     

Structural Studies of Overlapping Dinucleosomes in Solution

An overlapping dinucleosome (OLDN) is a structure composed of one hexasome and one octasome and appears to be formed through nucleosome collision promoted by nucleosome remodeling factor(s). In this study, the solution structure of the OLDN was investigated through the integration of small-angle x-ray and neutron scattering (SAXS and SANS, respectively), computer modeling, and molecular dynamics simulations. Starting from the crystal structure, we generated a conformational ensemble based on normal mode analysis and searched for the conformations that reproduced the SAXS and SANS scattering curves well. We found that inclusion of histone tails, which are not observed in the crystal structure, greatly improved model quality. The obtained structural models suggest that OLDNs adopt a variety of conformations stabilized by histone tails situated at the interface between the hexasome and octasome, simultaneously binding to both the hexasomal and octasomal DNA. In addition, our models define a possible direction for the conformational changes or dynamics, which may provide important information that furthers our understanding of the role of chromatin dynamics in gene regulation.     

30 nm chromatin fibre decompaction requires both H4-K16 acetylation and linker histone eviction

The mechanism by which chromatin is decondensed to permit access to DNA is largely unknown. Here, using a model nucleosome array reconstituted from recombinant histone octamers, we have defined the relative contribution of the individual histone octamer N-terminal tails as well as the effect of a targeted histone tail acetylation on the compaction state of the 30 nm chromatin fiber. This study goes beyond previous studies as it is based on a nucleosome array that is very long (61 nucleosomes) and contains a stoichiometric concentration of bound linker histone, which is essential for the formation of the 30 nm chromatin fiber. We find that compaction is regulated in two steps: Introduction of H4 acetylated to 30% on K16 inhibits compaction to a greater degree than deletion of the H4 N-terminal tail. Further decompaction is achieved by removal of the linker histone.     

Linker histone interaction shows divalent character with both supercoiled and linear DNA

The interaction of linker histone H1 with both linear and superhelical double-stranded DNA has been investigated at low ionic strengths. Gel mobility retardation experiments demonstrate strikingly different behavior for the two forms of DNA. First, the experiments strongly suggest that linker histone binds to superhelical DNA in a negatively cooperative mode. In contrast, binding of linker histone to linear DNA under the conditions employed here shows no cooperativity. Second, binding of linker histone to linear DNA results in aggregation of histone-DNA complexes, even at very low levels of input histone H1. Because H1 has been shown to interact as a monomer, this aggregation is evidence of the divalent character of the linker histone, for without H1's ability to bind to two duplex strands of DNA, aggregation could not occur. Although aggregation can be made to occur with superhelical DNA, it can do so only at near-saturation levels of input histone H1. Finally, in direct competition, linker histone binds to superhelical DNA to the complete exclusion of linear DNA, indicating that the linker histone's function is related to the crossover structures that differentiate superhelical DNA from linear DNA. We develop a model that explains the observed behavior of binding of linker histone to superhelical DNA that is consistent with both the divalent character of the linker histone and the negative cooperativity by which linker histone and superhelical DNA interact.     

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl₂, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5–12) carrying three positive charges on lysine⁺ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.     

Conformational and thermodynamic properties of supercoiled DNA.

We used Monte Carlo simulations to investigate the conformational and thermodynamic properties of DNA molecules with physiological levels of supercoiling. Three parameters determine the properties of DNA in this model: Kuhn statistical length, torsional rigidity and effective double-helix diameter. The chains in the simulation resemble strongly those observed by electron microscopy and have the conformation of an interwound superhelix whose axis is often branched. We compared the geometry of simulated chains with that determined experimentally by electron microscopy and by topological methods. We found a very close agreement between the Monte Carlo and experimental values for writhe, superhelix axis length and the number of superhelical turns. The computed number of superhelix branches was found to be dependent on superhelix density, DNA chain length and double-helix diameter. We investigated the thermodynamics of supercoiling and found that at low superhelix density the entropic contribution to superhelix free energy is negligible, whereas at high superhelix density, the entropic and enthalpic contributions are nearly equal. We calculated the effect of supercoiling on the spatial distribution of DNA segments. The probability that a pair of DNA sites separated along the chain contour by at least 50 nm are juxtaposed is about two orders of magnitude greater in supercoiled DNA than in relaxed DNA. This increase in the effective local concentration of DNA is not strongly dependent on the contour separation between the sites. We discuss the implications of this enhancement of site juxtaposition by supercoiling in the context of protein-DNA interactions involving multiple DNA-binding sites.     

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.