Rigid-body molecular dynamics of DNA inside a nucleosome

The role of transcript regions and amino acid choice in nucleosome positioning

Eukaryotic DNA is organized and compacted in a string of nucleosomes, DNA-wrapped protein cylinders. The positions of nucleosomes along DNA are not random but show well-known base pair sequence preferences that result from the sequence-dependent elastic and geometric properties of the DNA double helix. Here, we focus on DNA around transcription start sites, which are known to typically attract nucleosomes in multicellular life forms through their high GC content. We aim to understand how these GC signals, as observed in genome-wide averages, are produced and encoded through different genomic regions (mainly 5' UTRs, coding exons, and introns). Our study uses a bioinformatics approach to decompose the genome-wide GC signal into between-region and within-region signals. We find large differences in GC signal contributions between vertebrates and plants and, remarkably, even between closely related species. Introns contribute most to the GC signal in vertebrates, while in plants the exons dominate. Further, we find signal strengths stronger on DNA than on mRNA, suggesting a biological function of GC signals along the DNA itself, as is the case for nucleosome positioning. Finally, we make the surprising discovery that both the choice of synonymous codons and amino acids contribute to the nucleosome positioning signal.

Predicting the configuration and energy of DNA in a nucleosome by coarse-grain modelling

We present a novel algorithm which uses a coarse-grained model and an energy minimisation procedure to predict the sequence-dependent DNA configuration in a nucleosome together with its energetic cost.

The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study

The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences.

Nucleosomes are biomolecular complexes formed by DNA wrapped around histone proteins. They represent the basic units of Eukaryotic chromosomes, compacting the genome so that it fits into the small nucleus, and regulating important biological processes such as gene expression. Nucleosomes are disassembled during disruptive events such as DNA replication, and re-assembled afterwards to preserve the correct organization of chromatin. However, the molecular details of nucleosome assembly are still not well understood. In particular, experiments found that histones and DNA may associate into a variety of non-canonical complexes, but their precise conformation and role during assembly remain unclear. In this study, we addressed these problems by performing extensive molecular dynamics simulations of nucleosomes undergoing assembly and disassembly. The simulations reveal many insights into the kinetics of assembly, the structure of non-canonical nucleosome intermediates, and the influence of salt concentration and DNA sequence on the assembly process.

Translational nucleosome positioning: A computational study

About three-quarters of eukaryotic DNA is wrapped into nucleosomes; DNA spools with a protein core. The affinity of a given DNA stretch to be incorporated into a nucleosome is known to depend on the base-pair sequence-dependent geometry and elasticity of the DNA double helix. This causes the rotational and translational positioning of nucleosomes. In this study we ask the question whether the latter can be predicted by a simple coarse-grained DNA model with sequence-dependent elasticity, the rigid base-pair model. Whereas this model is known to be rather robust in predicting rotational nucleosome positioning, we show that the translational positioning is a rather subtle effect that is dominated by the guanine-cytosine content dependence of entropy rather than energy. A correct qualitative prediction within the rigid base-pair framework can only be achieved by assuming that DNA elasticity effectively changes on complexation into the nucleosome complex. With that extra assumption we arrive at a model which gives an excellent quantitative agreement to experimental in vitro nucleosome maps, under the additional assumption that nucleosomes equilibrate their positions only locally.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

The base pair-scale diffusion of nucleosomes modulates binding of transcription factors

The structure of promoter chromatin determines the ability of transcription factors (TFs) to bind to DNA and therefore has a profound effect on the expression levels of genes. However, the role of spontaneous nucleosome movements in this process is not fully understood. Here, we developed a single-molecule optical tweezers assay capable of simultaneously characterizing the base pair-scale diffusion of a nucleosome on DNA and the binding of a TF, using the luteinizing hormone β subunit gene (Lhb) promoter and Egr-1 as a model system. Our results demonstrate that nucleosomes undergo confined diffusion, and that the incorporation of the histone variant H2A.Z serves to partially relieve this confinement, inducing a different type of nucleosome repositioning. The increase in diffusion leads to exposure of a TF's binding site and facilitates its association with the DNA, which, in turn, biases the subsequent movement of the nucleosome. Our findings suggest the use of mobile nucleosomes as a general transcriptional regulatory mechanism.

The implication of DNA bending energy for nucleosome positioning and sliding

Nucleosome not only directly affects cellular processes, such as DNA replication, recombination, and transcription, but also severs as a fundamentally important target of epigenetic modifications. Our previous study indicated that the bending property of DNA is important in nucleosome formation, particularly in predicting the dyad positions of nucleosomes on a DNA segment. Here, we investigated the role of bending energy in nucleosome positioning and sliding in depth to decipher sequence-directed mechanism. The results show that bending energy is a good physical index to predict the free energy in the process of nucleosome reconstitution in vitro. Our data also imply that there are at least 20% of the nucleosomes in budding yeast do not adopt canonical positioning, in which underlying sequences wrapped around histones are structurally symmetric. We also revealed distinct patterns of bending energy profile for distinctly organized chromatin structures, such as well-positioned nucleosomes, fuzzy nucleosomes, and linker regions and discussed nucleosome sliding in terms of bending energy. We proposed that the stability of a nucleosome is positively correlated with the strength of the bending anisotropy of DNA segment, and both accessibility and directionality of nucleosome sliding is likely to be modulated by diverse patterns of DNA bending energy profile.

DNA sliding in nucleosomes via twist defect propagation revealed by molecular simulations

While nucleosomes are highly stable structures as fundamental units of chromatin, they also slide along the DNA, either spontaneously or by active remodelers. Here, we investigate the microscopic mechanisms of nucleosome sliding by multiscale molecular simulations, characterizing how the screw-like motion of DNA proceeds via the formation and propagation of twist defects. Firstly, coarse-grained molecular simulations reveal that the sliding dynamics is highly dependent on DNA sequence. Depending on the sequence and the nucleosome super-helical location, we find two distinct types of twist defects: a locally under-twisted DNA region, previously observed in crystal structures, and a locally over-twisted DNA, an unprecedented feature. The stability of the over-twist defect was confirmed via all-atom simulations. Analysis of our trajectories via Markov state modeling highlights how the sequence-dependence of the sliding dynamics is due to the different twist defect energy costs, and in particular how nucleosome regions where defects cannot easily form introduce the kinetic bottlenecks slowing down repositioning. Twist defects can also mediate sliding of nucleosomes made with strong positioning sequences, albeit at a much lower diffusion coefficient, due to a high-energy intermediate state. Finally, we discuss how chromatin remodelers may exploit these spontaneous fluctuations to induce unidirectional sliding of nucleosomes.

The role of DNA sequence in nucleosome breathing

Roughly 3/4 of human genomes are sequestered by nucleosomes, DNA spools with a protein core, dictating a broad range of biological processes, ranging from gene regulation, recombination, and replication, to chromosome condensation. Nucleosomes are dynamical structures and temporarily expose wrapped DNA through spontaneous unspooling from either end, a process called site exposure or nucleosome breathing. Here we ask how this process is influenced by the mechanical properties of the wrapped DNA, which is known to depend on the underlying base pair sequence. Using a coarse-grained nucleosome model we calculate the accessibility profiles for site exposure. We find that the process is very sensitive to sequence effects, so that evolution could potentially tune the accessibility of nucleosomal DNA and would only need a small number of mutations to do so.

Physics behind the mechanical nucleosome positioning code

The positions along DNA molecules of nucleosomes, the most abundant DNA-protein complexes in cells, are influenced by the sequence-dependent DNA mechanics and geometry. This leads to the "nucleosome positioning code", a preference of nucleosomes for certain sequence motives. Here we introduce a simplified model of the nucleosome where a coarse-grained DNA molecule is frozen into an idealized superhelical shape. We calculate the exact sequence preferences of our nucleosome model and find it to reproduce qualitatively all the main features known to influence nucleosome positions. Moreover, using well-controlled approximations to this model allows us to come to a detailed understanding of the physics behind the sequence preferences of nucleosomes.

In silico evidence for sequence-dependent nucleosome sliding

Nucleosomes represent the basic building block of chromatin and provide an important mechanism by which cellular processes are controlled. The locations of nucleosomes across the genome are not random but instead depend on both the underlying DNA sequence and the dynamic action of other proteins within the nucleus. These processes are central to cellular function, and the molecular details of the interplay between DNA sequence and nucleosome dynamics remain poorly understood. In this work, we investigate this interplay in detail by relying on a molecular model, which permits development of a comprehensive picture of the underlying free energy surfaces and the corresponding dynamics of nucleosome repositioning. The mechanism of nucleosome repositioning is shown to be strongly linked to DNA sequence and directly related to the binding energy of a given DNA sequence to the histone core. It is also demonstrated that chromatin remodelers can override DNA-sequence preferences by exerting torque, and the histone H4 tail is then identified as a key component by which DNA-sequence, histone modifications, and chromatin remodelers could in fact be coupled.

Benchmarking and refining probability-based models for nucleosome-DNA interaction

Background

In investigations of nucleosome positioning preferences, a model that assigns an affinity to a given sequence is necessary to make predictions. One important class of models, which treats a nucleosome sequence as a Markov chain, has been applied with success when informed with experimentally measured nucleosomal sequence preferences.

Results

We find that we can also use such models as a fast approximative scheme for computationally expensive biophysical models, vastly increasing their reach. Employing these models in this way also allows us to benchmark them for the first time. Doing so for the approximative in silico models indirectly tells us about the accuracy we can expect of them when applied to real data.

Conclusion

We find that models presented in the literature should perform well, but this performance depends on factors such as the order of the Markov model, the preprocessing of the probability distributions on which the model is based, and the size and quality of the sequence ensemble from which those distributions are calculated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1569-0) contains supplementary material, which is available to authorized users.

Multiplexing Genetic and Nucleosome Positioning Codes: A Computational Approach

Eukaryotic DNA is strongly bent inside fundamental packaging units: the nucleosomes. It is known that their positions are strongly influenced by the mechanical properties of the underlying DNA sequence. Here we discuss the possibility that these mechanical properties and the concomitant nucleosome positions are not just a side product of the given DNA sequence, e.g. that of the genes, but that a mechanical evolution of DNA molecules might have taken place. We first demonstrate the possibility of multiplexing classical and mechanical genetic information using a computational nucleosome model. In a second step we give evidence for genome-wide multiplexing in Saccharomyces cerevisiae and Schizosacharomyces pombe. This suggests that the exact positions of nucleosomes play crucial roles in chromatin function.

Nucleosome dynamics: Sequence matters

About three quarter of all eukaryotic DNA is wrapped around protein cylinders, forming nucleosomes. Even though the histone proteins that make up the core of nucleosomes are highly conserved in evolution, nucleosomes can be very different from each other due to posttranslational modifications of the histones. Another crucial factor in making nucleosomes unique has so far been underappreciated: the sequence of their DNA. This review provides an overview of the experimental and theoretical progress that increasingly points to the importance of the nucleosomal base pair sequence. Specifically, we discuss the role of the underlying base pair sequence in nucleosome positioning, sliding, breathing, force-induced unwrapping, dissociation and partial assembly and also how the sequence can influence higher-order structures. A new view emerges: the physical properties of nucleosomes, especially their dynamical properties, are determined to a large extent by the mechanical properties of their DNA, which in turn depends on DNA sequence.

Elastic model for dinucleosome structure and energy

The equilibrium structure of a dinucleosome is studied using an elastic model that takes into account the force and torque balance conditions. Using the proper boundary conditions, it is found that the conformational energy of the problem does not depend on the length of the linker DNA. In addition it is shown that the two histone octamers are almost perpendicular to each other, and the linker DNA in short lengths is almost straight. These findings could shed some light on the role of DNA elasticity in the chromatin structure.

Why Do Nucleosomes Unwrap Asymmetrically?

Nucleosomes, DNA spools with a protein core, engage about three-quarters of eukaryotic DNA and play a critical role in chromosomal processes, ranging from gene regulation, recombination, and replication to chromosome condensation. For more than a decade, micromanipulation experiments where nucleosomes are put under tension, as well as the theoretical interpretations of these experiments, have deepened our understanding of the stability and dynamics of nucleosomes. Here we give a theoretical explanation for a surprising new experimental finding: nucleosomes wrapped onto the 601 positioning sequence (the sequence used in most laboratories) respond highly asymmetrically to external forces by always unwrapping from the same end. Using a computational nucleosome model, we show that this asymmetry can be explained by differences in the DNA mechanics of two very short stretches on the wrapped DNA portion. Our finding suggests that the physical properties of nucleosomes, here the response to forces, can be tuned locally by the choice of the underlying base-pair sequence. This leads to a new view of nucleosomes: a physically highly varied set of DNA-protein complexes whose properties can be tuned on evolutionary time scales to their specific function in the genomic context.

Partial Unwrapping and Histone Tail Dynamics in Nucleosome Revealed by Coarse-Grained Molecular Simulations

Nucleosomes, basic units of chromatin, are known to show spontaneous DNA unwrapping dynamics that are crucial for transcriptional activation, but its structural details are yet to be elucidated. Here, employing a coarse-grained molecular model that captures residue-level structural details up to histone tails, we simulated equilibrium fluctuations and forced unwrapping of single nucleosomes at various conditions. The equilibrium simulations showed spontaneous unwrapping from outer DNA and subsequent rewrapping dynamics, which are in good agreement with experiments. We found several distinct partially unwrapped states of nucleosomes, as well as reversible transitions among these states. At a low salt concentration, histone tails tend to sit in the concave cleft between the histone octamer and DNA, tightening the nucleosome. At a higher salt concentration, the tails tend to bound to the outer side of DNA or be expanded outwards, which led to higher degree of unwrapping. Of the four types of histone tails, H3 and H2B tail dynamics are markedly correlated with partial unwrapping of DNA, and, moreover, their contributions were distinct. Acetylation in histone tails was simply mimicked by changing their charges, which enhanced the unwrapping, especially markedly for H3 and H2B tails.

Nucleosomes, folding units of chromatin, wrap DNA about 1.75 turns and provide bottlenecks for transcription. Recent experiments showed that nucleosomes are not rigid but dynamic, showing spontaneous and partial unwrapping which is thus important for transcriptional activation. Experimentally, however, one cannot directly watch DNA unwrapping at high resolution. On the other hand, molecular dynamics simulations have high spatio-temporal resolution and thus can be powerful and complementary to experiments. Here, we put forward coarse-grained modeling of protein-DNA interactions at residue-level resolution, which is rather generic and thus can be applied to any protein-DNA complexes. By this method, we could reveal spontaneous and salt-concentration dependent partial unwrapping of DNA from nucleosomes. In addition to consistency with single molecule experiments, the simulation showed multiple and distinct intermediate states of unwrapping. Interestingly, partial unwrapping of DNA is correlated with certain parts of histone tail dynamics. Deleting positive charges in histone tails that mimics histone acetylation facilitated partial unwrapping, most significantly for H3 and H2B.

Computational study of remodeling in a nucleosomal array

Chromatin remodeling complexes utilize the energy of ATP hydrolysis to change the packing state of chromatin, e.g. by catalysing the sliding of nucleosomes along DNA. Here we present simple models to describe experimental data of changes in DNA accessibility along a synthetic, repetitive array of nucleosomes during remodeling by the ACF enzyme or its isolated ATPase subunit, ISWI. We find substantial qualitative differences between the remodeling activities of ISWI and ACF. To understand better the observed behavior for the ACF remodeler, we study more microscopic models of nucleosomal arrays.

Inferring coarse-grain histone-DNA interaction potentials from high-resolution structures of the nucleosome

The histone-DNA interaction in the nucleosome is a fundamental mechanism of genomic compaction and regulation, which remains largely unknown despite increasing structural knowledge of the complex. In this paper, we propose a framework for the extraction of a nanoscale histone-DNA force-field from a collection of high-resolution structures, which may be adapted to a larger class of protein-DNA complexes. We applied the procedure to a large crystallographic database extended by snapshots from molecular dynamics simulations. The comparison of the structural models first shows that, at histone-DNA contact sites, the DNA base-pairs are shifted outwards locally, consistent with locally repulsive forces exerted by the histones. The second step shows that the various force profiles of the structures under analysis derive locally from a unique, sequence-independent, quadratic repulsive force-field, while the sequence preferences are entirely due to internal DNA mechanics. We have thus obtained the first knowledge-derived nanoscale interaction potential for histone-DNA in the nucleosome. The conformations obtained by relaxation of nucleosomal DNA with high-affinity sequences in this potential accurately reproduce the experimental values of binding preferences. Finally we address the more generic binding mechanisms relevant to the 80% genomic sequences incorporated in nucleosomes, by computing the conformation of nucleosomal DNA with sequence-averaged properties. This conformation differs from those found in crystals, and the analysis suggests that repulsive histone forces are related to local stretch tension in nucleosomal DNA, mostly between adjacent contact points. This tension could play a role in the stability of the complex.

The influence of DNA shape fluctuations on fluorescence resonance energy transfer efficiency measurements in nucleosomes

Fluorescence resonance energy transfer (FRET) measurements allow one to observe site exposure in nucleosomes, i.e. the transient unwrapping of a part of the wrapped DNA from the histone octamer. In such experiments one can typically distinguish between a closed state and an open state but in principle one might hope to detect several states, each corresponding to a certain number of open binding sites. Here we show that even in an ideal FRET setup it would be hard to detect unwrapping states with intermediate levels of FRET efficiencies. As the unwrapped DNA molecule, modelled here as a wormlike chain, has a finite stiffness, shape fluctuations smear out FRET signals completely from such intermediate states.

A simple model for DNA bridging proteins and bacterial or human genomes: bridging-induced attraction and genome compaction

We present computer simulations of the phase behaviour of an ensemble of proteins interacting with a polymer, mimicking non-specific binding to a piece of bacterial DNA or eukaryotic chromatin. The proteins can simultaneously bind to the polymer in two or more places to create protein bridges. Despite the lack of any explicit interaction between the proteins or between DNA segments, our simulations confirm previous results showing that when the protein-polymer interaction is sufficiently strong, the proteins come together to form clusters. Furthermore, a sufficiently large concentration of bridging proteins leads to the compaction of the swollen polymer into a globular phase. Here we characterise both the formation of protein clusters and the polymer collapse as a function of protein concentration, protein-polymer affinity and fibre flexibility.

DNA damage may drive nucleosomal reorganization to facilitate damage detection

One issue in genome maintenance is how DNA repair proteins find lesions at rates that seem to exceed diffusion-limited search rates. We propose a phenomenon where DNA damage induces nucleosomal rearrangements which move lesions to potential rendezvous points in the chromatin structure. These rendezvous points are the dyad and the linker DNA between histones, positions in the chromatin which are more likely to be accessible by repair proteins engaged in a random search. The feasibility of this mechanism is tested by considering the statistical mechanics of DNA containing a single lesion wrapped onto the nucleosome. We consider lesions which make the DNA either more flexible or more rigid by modeling the lesion as either a decrease or an increase in the bending energy. We include this energy in a partition function model of nucleosome breathing. Our results indicate that the steady state for a breathing nucleosome will most likely position the lesion at the dyad or in the linker, depending on the energy of the lesion. A role for DNA binding proteins and chromatin remodelers is suggested based on their ability to alter the mechanical properties of the DNA and DNA-histone binding, respectively. We speculate that these positions around the nucleosome potentially serve as rendezvous points where DNA lesions may be encountered by repair proteins which may be sterically hindered from searching the rest of the nucleosomal DNA. The strength of the repositioning is strongly dependent on the structural details of the DNA lesion and the wrapping and breathing of the nucleosome. A more sophisticated evaluation of this proposed mechanism will require detailed information about breathing dynamics, the structure of partially wrapped nucleosomes, and the structural properties of damaged DNA.

Confinement dynamics of a semiflexible chain inside nano-spheres

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.