Nucleosome dynamics: Sequence matters

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

Sequence Dependence in Nucleosome Dynamics

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, gene expression, to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse grained force field for energy estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the Self-Organized Polymer model with Side Chains for proteins and the Three Interaction Site model for DNA in a modular fashion, without re-calibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , the crystallographic B-factors for DNA-protein complexes, with varying sizes and topologies, are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts the scattering profiles that are in quantitative agreement with SAXS experiments as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which does not alter the balance of electrostatic interactions, but affects chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.

Specific DNA binding of archaeal histones HMfA and HMfB

In archaea, histones play a role in genome compaction and are involved in transcription regulation. Whereas archaeal histones bind DNA without sequence specificity, they bind preferentially to DNA containing repeats of alternating A/T and G/C motifs. These motifs are also present on the artificial sequence "Clone20," a high-affinity model sequence for binding of the histones from Methanothermus fervidus. Here, we investigate the binding of HMfA and HMfB to Clone20 DNA. We show that specific binding at low protein concentrations (<30 nM) yields a modest level of DNA compaction, attributed to tetrameric nucleosome formation, whereas nonspecific binding strongly compacts DNA. We also demonstrate that histones impaired in hypernucleosome formation are still able to recognize the Clone20 sequence. Histone tetramers indeed exhibit a higher binding affinity for Clone20 than nonspecific DNA. Our results indicate that a high-affinity DNA sequence does not act as a nucleation site, but is bound by a tetramer which we propose is geometrically different from the hypernucleosome. Such a mode of histone binding might permit sequence-driven modulation of hypernucleosome size. These findings might be extrapolated to histone variants that do not form hypernucleosomes. Versatile binding modes of histones could provide a platform for functional interplay between genome compaction and transcription.

Multiplexing mechanical and translational cues on genes

The genetic code gives precise instructions on how to translate codons into amino acids. Due to the degeneracy of the genetic code-18 out of 20 amino acids are encoded for by more than one codon-more information can be stored in a basepair sequence. Indeed, various types of additional information have been discussed in the literature, e.g., the positioning of nucleosomes along eukaryotic genomes and the modulation of the translating efficiency in ribosomes to influence cotranslational protein folding. The purpose of this study is to show that it is indeed possible to carry more than one additional layer of information on top of a gene. In particular, we show how much translation efficiency and nucleosome positioning can be adjusted simultaneously without changing the encoded protein. We achieve this by mapping genes on weighted graphs that contain all synonymous genes, and then finding shortest paths through these graphs. This enables us, for example, to readjust the disrupted translational efficiency profile after a gene has been introduced from one organism (e.g., human) into another (e.g., yeast) without greatly changing the nucleosome landscape intrinsically encoded by the DNA molecule.

OCT4 interprets and enhances nucleosome flexibility

Pioneer transcription factors are proteins that induce cellular identity transitions by binding to inaccessible regions of DNA in nuclear chromatin. They contribute to chromatin opening and recruit other factors to regulatory DNA elements. The structural features and dynamics modulating their interaction with nucleosomes are still unresolved. From a combination of experiments and molecular simulations, we reveal here how the pioneer factor and master regulator of pluripotency, Oct4, interprets and enhances nucleosome structural flexibility. The magnitude of Oct4’s impact on nucleosome dynamics depends on the binding site position and the mobility of the unstructured tails of nucleosomal histone proteins. Oct4 uses both its DNA binding domains to propagate and stabilize open nucleosome conformations, one for specific sequence recognition and the other for nonspecific interactions with nearby regions of DNA. Our findings provide a structural basis for the versatility of transcription factors in engaging with nucleosomes and have implications for understanding how pioneer factors induce chromatin dynamics.

Reconstituted TAD-size chromatin fibers feature heterogeneous nucleosome clusters

Large topologically associated domains (TADs) contain irregularly spaced nucleosome clutches, and interactions between such clutches are thought to aid the compaction of these domains. Here, we reconstituted TAD-sized chromatin fibers containing hundreds of nucleosomes on native source human and lambda-phage DNA and compared their mechanical properties at the single-molecule level with shorter ‘601’ arrays with various nucleosome repeat lengths. Fluorescent imaging showed increased compaction upon saturation of the DNA with histones and increasing magnesium concentration. Nucleosome clusters and their structural fluctuations were visualized in confined nanochannels. Force spectroscopy revealed not only similar mechanical properties of the TAD-sized fibers as shorter fibers but also large rupture events, consistent with breaking the interactions between distant clutches of nucleosomes. Though the arrays of native human DNA, lambda-phage and ‘601’ DNA featured minor differences in reconstitution yield and nucleosome stability, the fibers’ global structural and mechanical properties were similar, including the interactions between nucleosome clutches. These single-molecule experiments quantify the mechanical forces that stabilize large TAD-sized chromatin domains consisting of disordered, dynamically interacting nucleosome clutches and their effect on the condensation of large chromatin domains.

Stochastic models of nucleosome dynamics reveal regulatory rules of stimulus-induced epigenome remodeling

The genomic positions of nucleosomes are a defining feature of the cell's epigenomic state, but signal-dependent transcription factors (SDTFs), upon activation, bind to specific genomic locations and modify nucleosome positioning. Here we leverage SDTFs as perturbation probes to learn about nucleosome dynamics in living cells. We develop Markov models of nucleosome dynamics and fit them to time course sequencing data of DNA accessibility. We find that (1) the dynamics of DNA unwrapping are significantly slower in cells than reported from cell-free experiments, (2) only models with cooperativity in wrapping and unwrapping fit the available data, (3) SDTF activity produces the highest eviction probability when its binding site is adjacent to but not on the nucleosome dyad, and (4) oscillatory SDTF activity results in high location variability. Our work uncovers the regulatory rules governing SDTF-induced nucleosome dynamics in live cells, which can predict chromatin accessibility alterations during inflammation at single-nucleosome resolution.

Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools—software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

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.

Asymmetry in histone rotation in forced unwrapping and force quench rewrapping in a nucleosome

Single molecule pulling experiments have shown that DNA in the nucleosomes unwraps in two stages from the histone protein core (HPC). The first stage, attributed to the rupture of the outer DNA turn, occurs between 3 and 5 pNs, and is reversible. The inner DNA turn ruptures irreversibly at forces between 9 and 15 pNs (or higher) in the second stage. Molecular simulations using the Self-Organized Polymer model capture the experimental findings. The unwrapping of the outer DNA turn is independent of the pulling direction. The rupture of the DNA inner turn depends on the pulling direction and involves overcoming substantial energetic (most likely electrostatic in origin) and kinetic barriers. They arise because the mechanical force has to generate sufficient torque to rotate the HPC by 180°. On the other hand, during the rewrapping process, HPC rotation is stochastic, with force playing no role. The assembly of the outer DNA wrap upon force quench nearly coincides with the unwrapping process, confirming the reversibility of the outer turn rupture. The asymmetry in HPC rotation during unwrapping and rewrapping explains the observed hysteresis in the stretch-release cycles in experiments. We propose experiments to test the prediction that HPC rotation produces kinetic barriers in the unwrapping process.

Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions

Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS.

A critical role for linker DNA in higher-order folding of chromatin fibers

Nucleosome-nucleosome interactions drive the folding of nucleosomal arrays into dense chromatin fibers. A better physical account of the folding of chromatin fibers is necessary to understand the role of chromatin in regulating DNA transactions. Here, we studied the unfolding pathway of regular chromatin fibers as a function of single base pair increments in linker length, using both rigid base-pair Monte Carlo simulations and single-molecule force spectroscopy. Both computational and experimental results reveal a periodic variation of the folding energies due to the limited flexibility of the linker DNA. We show that twist is more restrictive for nucleosome stacking than bend, and find the most stable stacking interactions for linker lengths of multiples of 10 bp. We analyzed nucleosomes stacking in both 1- and 2-start topologies and show that stacking preferences are determined by the length of the linker DNA. Moreover, we present evidence that the sequence of the linker DNA also modulates nucleosome stacking and that the effect of the deletion of the H4 tail depends on the linker length. Importantly, these results imply that nucleosome positioning in vivo not only affects the phasing of nucleosomes relative to DNA but also directs the higher-order structure of chromatin.

Twist dynamics and buckling instability of ring DNA: the effect of groove asymmetry and anisotropic bending

By combining analytical theory and Molecular Dynamics simulations we study the relaxation dynamics of DNA circular plasmids that initially undergo a local twist perturbation. In this process, the twist-bend coupling arising from the groove asymmetry in the DNA double helix clearly shows up. In the two scenarios explored, with/without this coupling, the initial perturbation relaxes diffusively. However, there are some marked differences in the value of the diffusion coefficient and the dynamics in both cases. These differences can be explained by assuming the existence of three distinctive time scales; a rapid relaxation of local bending, the slow twist spreading, and the buckling transition taking place in a much longer time scale. In particular, the separation of time scales allows deducing an effective diffusion equation in stage , with a diffusion coefficient influenced by the twist-bend coupling. We also discuss the mapping of the realistic DNA model to the simpler isotropic twistable worm-like chain using the renormalized bending and twist moduli; although useful in many cases, it fails to make a quantitative prediction on the instability mode of buckling transition.

CENP-A Nucleosome is a Sensitive Allosteric Scaffold for DNA and Chromatin Factors

Centromeric loci of chromosomes are defined by nucleosomes containing the histone H3 variant CENP-A, which bind their DNA termini more permissively than their canonical counterpart, a feature that is critical for the mitotic fidelity. A recent cryo-EM study demonstrated that the DNA termini of CENP-A nucleosomes, reconstituted with the Widom 601 DNA sequence, are asymmetrically flexible, meaning one terminus is more clearly resolved than the other. However, an earlier work claimed that both ends could be resolved in the presence of two stabilizing single chain variable fragment (scFv) antibodies per nucleosome, and thus are likely permanently bound to the histone octamer. This suggests that the binding of scFv antibodies to the histone octamer surface would be associated with CENP-A nucleosome conformational changes, including stable binding of the DNA termini. Here, we present computational evidence that allows to explain at atomistic level the structural rearrangements of CENP-A nucleosomes resulting from the antibody binding. The antibodies, while they only bind the octamer façades, are capable of altering the dynamics of the nucleosomal core, and indirectly also the surrounding DNA. This effect has more drastic implications for the structure and the dynamics of the CENP-A nucleosome in comparison to its canonical counterpart. Furthermore, we find evidence that the antibodies bind the left and the right octamer façades at different affinities, another manifestation of the DNA sequence. We speculate that the cells could use induction of similar allosteric effects to control centromere function.

Dynamic networks observed in the nucleosome core particles couple the histone globular domains with DNA

The dynamics of eukaryotic nucleosomes are essential in gene activity and well regulated by various factors. Here, we elucidated the internal dynamics at multiple timescales for the human histones hH3 and hH4 in the Widom 601 nucleosome core particles (NCP), suggesting that four dynamic networks are formed by the residues exhibiting larger-scale μs-ms motions that extend from the NCP core to the histone tails and DNA. Furthermore, despite possessing highly conserved structural features, histones in the telomeric NCP exhibit enhanced μs-ms dynamics in the globular sites residing at the identified dynamic networks and in a neighboring region. In addition, higher mobility was observed for the N-terminal tails of hH3 and hH4 in the telomeric NCP. The results demonstrate the existence of dynamic networks in nucleosomes, through which the center of the core regions could interactively communicate with histone tails and DNA to potentially propagate epigenetic changes.

Shi et al. use solid-state nuclear magnetic resonance spectroscopy to reveal the internal dynamics of human histones hH3 and hH4 in the Widom 601 and the telomeric nucleosome core particles. This work has implications for the propagation of epigenetic changes via the center of the nucleosome core communicating with histone tails and DNA.

Multiplexed Nanometric 3D Tracking of Microbeads Using an FFT-Phasor Algorithm

Many single-molecule biophysical techniques rely on nanometric tracking of microbeads to obtain quantitative information about the mechanical properties of biomolecules such as chromatin fibers. Their three-dimensional (3D) position can be resolved by holographic analysis of the diffraction pattern in wide-field imaging. Fitting this diffraction pattern to Lorenz-Mie scattering theory yields the bead's position with nanometer accuracy in three dimensions but is computationally expensive. Real-time multiplexed bead tracking therefore requires a more efficient tracking method, such as comparison with previously measured diffraction patterns, known as look-up tables. Here, we introduce an alternative 3D phasor algorithm that provides robust bead tracking with nanometric localization accuracy in a z range of over 10 μm under nonoptimal imaging conditions. The algorithm is based on a two-dimensional cross correlation using fast Fourier transforms with computer-generated reference images, yielding a processing rate of up to 10,000 regions of interest per second. We implemented the technique in magnetic tweezers and tracked the 3D position of over 100 beads in real time on a generic CPU. The accuracy of 3D phasor tracking was extensively tested and compared to a look-up table approach using Lorenz-Mie simulations, avoiding experimental uncertainties. Its easy implementation, efficiency, and robustness can improve multiplexed biophysical bead-tracking applications, especially when high throughput is required and image artifacts are difficult to avoid.

Bottom-Up Meets Top-Down: The Crossroads of Multiscale Chromatin Modeling

Chromatin can be viewed as a hierarchically structured fiber that regulates gene expression. It consists of a complex network of DNA and proteins whose characteristic dynamical modes facilitate compaction and rearrangement in the cell nucleus. These modes stem from chromatin's fundamental unit, the nucleosome, and their effects are propagated across length scales. Understanding the effects of nucleosome dynamics on the chromatin fiber, primarily through post-translational modifications that occur on the histones, is of central importance to epigenetics. Within the last decade, imaging and chromosome conformation capture techniques have revealed a number of structural and statistical features of the packaged chromatin fiber at a hitherto unavailable level of resolution. Such experiments have led to increased efforts to develop polymer models that aim to reproduce, explain, and predict the contact probability scaling and density heterogeneity. At nanometer scales, available models have focused on the role of the nucleosome and epigenetic marks on local chromatin structure. At micrometer scales, existing models have sought to explain scaling laws and density heterogeneity. Less work, however, has been done to reconcile these two approaches: bottom-up and top-down models of chromatin. In this perspective, we highlight the multiscale simulation models that are driving toward an understanding of chromatin structure and function, from the nanometer to the micron scale, and we highlight areas of opportunity and some of the prospects for new frameworks that bridge these two scales. Taken together, experimental and modeling advances over the last few years have established a robust platform for the study of chromatin fiber structure and dynamics, which will be of considerable use to the chromatin community in developing an understanding of the interplay between epigenomic regulation and molecular structure.

Sensitive Automated Measurement of Histone-DNA Affinities in Nucleosomes

The DNA of eukaryotes is wrapped around histone octamers to form nucleosomes. Although it is well established that the DNA sequence significantly influences nucleosome formation, its precise contribution has remained controversial, partially owing to the lack of quantitative affinity data. Here, we present a method to measure DNA-histone binding free energies at medium throughput and with high sensitivity. Competitive nucleosome formation is achieved through automation, and a modified epifluorescence microscope is used to rapidly and accurately measure the fractions of bound/unbound DNA based on fluorescence anisotropy. The procedure allows us to obtain full titration curves with high reproducibility. We applied this technique to measure the histone-DNA affinities for 47 DNA sequences and analyzed how the affinities correlate with relevant DNA sequence features. We found that the GC content has a significant impact on nucleosome-forming preferences, but 10 bp dinucleotide periodicities and the presence of poly(dA:dT) stretches do not.

•

Robotics permits full titration series to measure histone-DNA binding affinities

•

Fluorescence anisotropy used as a fast, sensitive readout of bound/unbound DNA

•

Free energies span three orders of magnitude, less for naturally occurring sequences

•

GC content is a major determinant of measured binding free energies

Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding

Transcription factor (TF) proteins bind to DNA to regulate gene expression. Normally, accessibility to DNA is required for their function. However, in the nucleus, the DNA is often inaccessible, wrapped around histone proteins in nucleosomes forming the chromatin. Pioneer TFs are thought to induce chromatin opening by recognizing their DNA binding sites on nucleosomes. For example, Oct4, a master regulator and inducer of stem cell pluripotency, binds to DNA in nucleosomes in a sequence-specific manner. Here, we reveal the structural dynamics of nucleosomes that mediate Oct4 binding from molecular dynamics simulations. Nucleosome flexibility and the amplitude of nucleosome motions such as breathing and twisting are enhanced in nucleosomes with multiple TF binding sites. Moreover, the regions around the binding sites display higher local structural flexibility. Probing different structures of Oct4-nucleosome complexes, we show that alternative configurations in which Oct4 recognizes partial binding sites display stable TF-DNA interactions similar to those observed in complexes with free DNA and compatible with the DNA curvature and DNA-histone interactions. Therefore, we propose a structural basis for nucleosome recognition by a pioneer TF that is essential for understanding how chromatin is unraveled during cell fate conversions.

Chromatin fibers stabilize nucleosomes under torsional stress

Torsional stress generated during DNA replication and transcription has been suggested to facilitate nucleosome unwrapping and thereby the progression of polymerases. However, the propagation of twist in condensed chromatin remains yet unresolved. Here, we measure how force and torque impact chromatin fibers with a nucleosome repeat length of 167 and 197. We find that both types of fibers fold into a left-handed superhelix that can be stabilized by positive torsion. We observe that the structural changes induced by twist were reversible, indicating that chromatin has a large degree of elasticity. Our direct measurements of torque confirmed the hypothesis of chromatin fibers as a twist buffer. Using a statistical mechanics-based torsional spring model, we extracted values of the chromatin twist modulus and the linking number per stacked nucleosome that were in good agreement with values measured here experimentally. Overall, our findings indicate that the supercoiling generated by DNA-processing enzymes, predicted by the twin-supercoiled domain model, can be largely accommodated by the higher-order structure of chromatin.

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.

DNA Sequence Is a Major Determinant of Tetrasome Dynamics

Eukaryotic genomes are hierarchically organized into protein-DNA assemblies for compaction into the nucleus. Nucleosomes, with the (H3-H4)2 tetrasome as a likely intermediate, are highly dynamic in nature by way of several different mechanisms. We have recently shown that tetrasomes spontaneously change the direction of their DNA wrapping between left- and right-handed conformations, which may prevent torque buildup in chromatin during active transcription or replication. DNA sequence has been shown to strongly affect nucleosome positioning throughout chromatin. It is not known, however, whether DNA sequence also impacts the dynamic properties of tetrasomes. To address this question, we examined tetrasomes assembled on a high-affinity DNA sequence using freely orbiting magnetic tweezers. In this context, we also studied the effects of mono- and divalent salts on the flipping dynamics. We found that neither DNA sequence nor altered buffer conditions affect overall tetrasome structure. In contrast, tetrasomes bound to high-affinity DNA sequences showed significantly altered flipping kinetics, predominantly via a reduction in the lifetime of the canonical state of left-handed wrapping. Increased mono- and divalent salt concentrations counteracted this behavior. Thus, our study indicates that high-affinity DNA sequences impact not only the positioning of the nucleosome but that they also endow the subnucleosomal tetrasome with enhanced conformational plasticity. This may provide a means to prevent histone loss upon exposure to torsional stress, thereby contributing to the integrity of chromatin at high-affinity sites.

Sequence-dependent nucleosome nanoscale structure characterized by atomic force microscopy

DNA sequence plays an important role in the assembly of nucleosomes, and DNA motifs with high specificity to nucleosomes have been identified. At the same time, important questions such as how the DNA sequence changes DNA wrapping and how the DNA sequence contributes to the interaction between the nucleosomes remain unclear. Here, we addressed these questions by comparing nanoscale properties of nucleosomes assembled on the highest nucleosome positioning sequence, the 601 motif, and essentially random DNA sequences. We used atomic force microscopy to measure the nucleosome positions and the DNA wrapping. The studies showed that nucleosomes assemble on the nonspecific sequence without any preference for position on DNA, but they wrap the same DNA length as DNA of the same length with the 601 sequence. Experiments with longer DNA containing the 601 motif along with nonspecific DNA, capable of forming dinucleosomes, revealed that dinucleosomes assembled on hybrid sequences show a preference for positioning near each other, with one always assembling on the 601 motif. These findings point to the interaction between the nucleosomes and suggest that internucleosomal interactions may play a large role in nucleosome positioning.-Stormberg, T., Stumme-Diers, M., Lyubchenko, Y. L. Sequence-dependent nucleosome nanoscale structure characterized by atomic force microscopy.

Replication-dependent histone isoforms: a new source of complexity in chromatin structure and function

Replication-dependent histones are expressed in a cell cycle regulated manner and supply the histones necessary to support DNA replication. In mammals, the replication-dependent histones are encoded by a family of genes that are located in several clusters. In humans, these include 16 genes for histone H2A, 22 genes for histone H2B, 14 genes for histone H3, 14 genes for histone H4 and 6 genes for histone H1. While the proteins encoded by these genes are highly similar, they are not identical. For many years, these genes were thought to encode functionally equivalent histone proteins. However, several lines of evidence have emerged that suggest that the replication-dependent histone genes can have specific functions and may constitute a novel layer of chromatin regulation. This Survey and Summary reviews the literature on replication-dependent histone isoforms and discusses potential mechanisms by which the small variations in primary sequence between the isoforms can alter chromatin function. In addition, we summarize the wealth of data implicating altered regulation of histone isoform expression in cancer.

Biophysics of Chromatin Dynamics

Nucleosomes and chromatin control eukaryotic genome accessibility and thereby regulate DNA processes, including transcription, replication, and repair. Conformational dynamics within the nucleosome and chromatin structure play a key role in this regulatory function. Structural fluctuations continuously expose internal DNA sequences and nucleosome surfaces, thereby providing transient access for the nuclear machinery. Progress in structural studies of nucleosomes and chromatin has provided detailed insight into local chromatin organization and has set the stage for recent in-depth investigations of the structural dynamics of nucleosomes and chromatin fibers. Here, we discuss the dynamic processes observed in chromatin over different length scales and timescales and review current knowledge about the biophysics of distinct structural transitions.

Overtwisting induces polygonal shapes in bent DNA

By combining analytical results and simulations of various coarse-grained models, we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternation of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, rather than the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants.

Structural interpretation of DNA-protein hydroxyl-radical footprinting experiments with high resolution using HYDROID

Hydroxyl-radical footprinting (HRF) is a powerful method for probing structures of nucleic acid-protein complexes with single-nucleotide resolution in solution. To tap the full quantitative potential of HRF, we describe a protocol, hydroxyl-radical footprinting interpretation for DNA (HYDROID), to quantify HRF data and integrate them with atomistic structural models. The stages of the HYDROID protocol are extraction of the lane profiles from gel images, quantification of the DNA cleavage frequency at each nucleotide and theoretical estimation of the DNA cleavage frequency from atomistic structural models, followed by comparison of experimental and theoretical results. Example scripts for each step of HRF data analysis and interpretation are provided for several nucleosome systems; they can be easily adapted to analyze user data. As input, HYDROID requires polyacrylamide gel electrophoresis (PAGE) images of HRF products and optionally can use a molecular model of the DNA-protein complex. The HYDROID protocol can be used to quantify HRF over DNA regions of up to 100 nucleotides per gel image. In addition, it can be applied to the analysis of RNA-protein complexes and free RNA or DNA molecules in solution. Compared with other methods reported to date, HYDROID is unique in its ability to simultaneously integrate HRF data with the analysis of atomistic structural models. HYDROID is freely available. The complete protocol takes ~3 h. Users should be familiar with the command-line interface, the Python scripting language and Protein Data Bank (PDB) file formats. A graphical user interface (GUI) with basic functionality (HYDROID_GUI) is also available.

Bend-Induced Twist Waves and the Structure of Nucleosomal DNA

Recent work indicates that twist-bend coupling plays an important role in DNA micromechanics. Here we investigate its effect on bent DNA. We provide an analytical solution of the minimum-energy shape of circular DNA, showing that twist-bend coupling induces sinusoidal twist waves. This solution is in excellent agreement with both coarse-grained simulations of minicircles and nucleosomal DNA data, which is bent and wrapped around histone proteins in a superhelical conformation. Our analysis shows that the observed twist oscillation in nucleosomal DNA, so far attributed to the interaction with the histone proteins, is an intrinsic feature of free bent DNA, and should be observable in other protein-DNA complexes.

Native DNA electronics: is it a matter of nanoscale assembly?

The genomic DNA is enveloped by nanotubes formed by the nuclear aggregates of polyamines (NAPs) that induce DNA conformational changes and provide protection and increased interaction abilities for the double strands. In a physiological environment, the nanotube arrangement is initiated by spontaneous interaction among the terminal amino groups of the polyamines and the phosphate ions, with the consequent formation of cyclic monomers that hook at the DNA grooves. The polymer thus formed has the morphological features of an organic semiconductor and therefore, it can be considered to be able to conduct electric charges. Phosphate ions positioned on the NAP external surface could regulate, as in a physical electric circuit, both linear and rotational (histones) protein motion, in accordance with the basilar principles of the electronics. A model of a carrier system for protein motion along the polymer wrapping the DNA strands, based on the phosphate-phosphate complexation, is proposed.

The Latest Twists in Chromatin Remodeling

In its most restrictive interpretation, the notion of chromatin remodeling refers to the action of chromatin-remodeling enzymes on nucleosomes with the aim of displacing and removing them from the chromatin fiber (the effective polymer formed by a DNA molecule and proteins). This local modification of the fiber structure can have consequences for the initiation and repression of the transcription process, and when the remodeling process spreads along the fiber, it also results in long-range effects essential for fiber condensation. There are three regulatory levels of relevance that can be distinguished for this process: the intrinsic sequence preference of the histone octamer, which rules the positioning of the nucleosome along the DNA, notably in relation to the genetic information coded in DNA; the recognition or selection of nucleosomal substrates by remodeling complexes; and, finally, the motor action on the nucleosome exerted by the chromatin remodeler. Recent work has been able to provide crucial insights at each of these three levels that add new twists to this exciting and unfinished story, which we highlight in this perspective.

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.

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.

Evolutionary analysis of nucleosome positioning sequences based on New Symmetric Relative Entropy

New Symmetric Relative Entropy (NSRE) was applied innovatively to analyze the nucleosome sequences in S. cerevisiae, S. pombe and Drosophila. NSRE distributions could well reflect the characteristic differences of nucleosome sequences among three organisms, and the differences indicate a concerted evolution in the sequence usage of nucleosome. Further analysis about the nucleosomes around TSS shows that the constitutive property of +1/-1 nucleosomes in S. cerevisiae is different from that in S. pombe and Drosophila, which indicates that S. cerevisiae has a different transcription regulation mechanism based on nucleosome. However, in either case, the nucleosome dyad region is conserved and always has a higher NSRE. Base composition analysis shows that this conservative property in nucleosome dyad region is mainly determined by base A and T, and the dependence degrees on base A and T are consistent in three organisms.

Structural transitions in torsionally constrained DNA and their dependence on solution electrostatics

Experimental studies on single molecules of DNA have reported a rich variety of structural transitions, including coexistence of three phases, in a torsionally constrained molecule. A comprehensive knowledge of these structural transitions is useful for unraveling the in vivo and in vitro behavior of DNA. Our objective is to understand the structural transitions in a torsionally constrained DNA molecule when it is pulled using optical or magnetic tweezers. We use foundational concepts from the Zimm-Bragg helix-coil transition theory and merge them with ideas from the theory of fluctuating elastic rods to model the mechanics of DNA. We also account for the electrostatic interactions between the ions and the negatively charged phosphate backbone of DNA. Using our model, we calculate the force and torque corresponding to the overstretching transition characterized by a 70% jump in the contour length of the molecule and examine the effect of salt concentration on this transition. We also deduce conditions under which the co-existence of B-, S- and P-DNA is possible. We examine how the cooperativity parameter for each transition affects the force-extension curve or torque-rotation curve. We attempt to rationalize the non-monotonic dependence of external work done on the ion concentration by connecting it to the electrostatic dependence of the interfacial energy between two phases of DNA. Our theoretical results are in agreement with multiple experiments documented in the literature and they generate falsifiable predictions that can be tested in new experiments.

STATEMENT OF SIGNIFICANCE

The overarching objective of this paper is to explore the implications of variation in ion concentration on the structural transitions driven by external forces in a torsionally constrained DNA molecule. A comprehensive understanding of the phase behavior of torsionally constrained DNA is useful because DNA in cells is tightly packaged and is acted upon by molecular machines in different ionic environments. We examine the mechanics of the overstretching transition, characterized by a 70% jump in contour length, wherein a mixture of B- and S-DNA converts into a mixture of S- and P-DNA through a triple point in the phase diagram. Our results are corroborated by experimental data at every step and we make predictions that are experimentally verifiable.

Designing nucleosomal force sensors

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Nucleosome mobility and the regulation of gene expression: Insights from single-molecule studies

Nucleosomes at the promoters of genes regulate the accessibility of the transcription machinery to DNA, and function as a basic layer in the complex regulation of gene expression. Our understanding of the role of the nucleosome's spontaneous, thermally driven position changes in modulating expression is lacking. This is the result of the paucity of experimental data on these dynamics, at high-resolution, and for DNA sequences that belong to real, transcribed genes. We have developed an assay that uses partial, reversible unzipping of nucleosomes with optical tweezers to repeatedly probe a nucleosome's position over time. Using the nucleosomes at the promoters of two model genes, Cga and Lhb, we show that the mobility of nucleosomes is modulated by the sequence of DNA and by the use of alternative histone variants, and describe how the mobility can affect transcription, at the initiation and elongation phases.

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

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.

Unique yeast histone sequences influence octamer and nucleosome stability

Yeast nucleosomes are known to be intrinsically less stable than those from higher eukaryotes. This difference presents significant challenges for the production of yeast nucleosome core particles (NCPs) and chromatin for in vitro analyses. Using recombinant yeast, human, and chimeric histone proteins, we demonstrate that three divergent amino acids in histone H3 (Q120 K121 K125 ) are responsible for the poor reconstitution of yeast histones into octamers. This QKK motif is only found in Fungi, and is located at the nucleosome dyad axis. Yeast-to-human changes at these positions render yeast histones amenable to well-established octamer reconstitution and salt dialysis methods for generating nucleosomal and longer chromatin templates. By contrast, the most divergent yeast core histones, H2A and H2B, affect the biophysical properties of NCP but not their stability. An evolutionary analysis of H3 sequences shows that a gradual divergence in H3 sequences occurred in Fungi to yield QKK in budding yeast. This likely facilitates the highly euchromatic nature of yeast genomes. Our results provide an explanation for the long recognized difference in yeast nucleosome stability and they offer a simple method to generate yeast chromatin templates for in vitro studies.