Kinetic control of nucleosome displacement by ISWI/ACF chromatin remodelers

Chromatin remodelers: a concise introduction for biophysicists

Chromatin remodelers are molecular motors that act on nucleosomes: they move them along DNA or (dis-)assemble them. Despite the fact that they perform essential regulatory functions in cells-their deregulation can contribute to the development of cancers and lead to cell death-chromatin remodelers have only received meager attention in the biophysics community so far. In this short text, we attempt to present the key features of this interesting class of enzymes obtained with different experimental and theoretical methods, thereby providing a concise introduction for biophysicists to further stimulate interest in their properties.

Nucleosome sliding can influence the spreading of histone modifications

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

Nucleosome Assembly and Disassembly in vitro Are Governed by Chemical Kinetic Principles

As the elementary unit of eukaryotic chromatin, nucleosomes in vivo are highly dynamic in many biological processes, such as DNA replication, repair, recombination, or transcription, to allow the necessary factors to gain access to their substrate. The dynamic mechanism of nucleosome assembly and disassembly has not been well described thus far. We proposed a chemical kinetic model of nucleosome assembly and disassembly in vitro. In the model, the efficiency of nucleosome assembly was positively correlated with the total concentration of histone octamer, reaction rate constant and reaction time. All the corollaries of the model were well verified for the Widom 601 sequence and the six artificially synthesized DNA sequences, named CS1–CS6, by using the salt dialysis method in vitro. The reaction rate constant in the model may be used as a new parameter to evaluate the nucleosome reconstitution ability with DNAs. Nucleosome disassembly experiments for the Widom 601 sequence detected by Förster resonance energy transfer (FRET) and fluorescence thermal shift (FTS) assays demonstrated that nucleosome disassembly is the inverse process of assembly and can be described as three distinct stages: opening phase of the (H2A–H2B) dimer/(H3–H4)₂ tetramer interface, release phase of the H2A–H2B dimers from (H3–H4)₂ tetramer/DNA and removal phase of the (H3–H4)₂ tetramer from DNA. Our kinetic model of nucleosome assembly and disassembly allows to confirm that nucleosome assembly and disassembly in vitro are governed by chemical kinetic principles.

On the role of transcription in positioning nucleosomes

Nucleosome positioning is crucial for the genome’s function. Though the role of DNA sequence in positioning nucleosomes is well understood, a detailed mechanistic understanding on the impact of transcription remains lacking. Using numerical simulations, we investigated the dependence of nucleosome density profiles on transcription level across multiple species. We found that the low nucleosome affinity of yeast, but not mouse, promoters contributes to the formation of phased nucleosomes arrays for inactive genes. For the active genes, a heterogeneous distribution of +1 nucleosomes, caused by a tug-of-war between two types of remodeling enzymes, is essential for reproducing their density profiles. In particular, while positioning enzymes are known to remodel the +1 nucleosome and align it toward the transcription start site (TSS), spacer enzymes that use a pair of nucleosomes as their substrate can shift the nucleosome array away from the TSS. Competition between these enzymes results in two types of nucleosome density profiles with well- and ill-positioned +1 nucleosome. Finally, we showed that Pol II assisted histone exchange, if occurring at a fast speed, can abolish the impact of remodeling enzymes. By elucidating the role of individual factors, our study reconciles the seemingly conflicting results on the overall impact of transcription in positioning nucleosomes across species.

Nucleosome positioning plays a key role in the genome’s function by regulating the accessibility of protein binding sites as well as higher-order chromatin organization. Though significant progress has been made towards studying the role of DNA sequence in positioning the nucleosomes, our understanding on the impact of transcription lags behind. Our study uses kinetic simulations to explore the role of DNA sequence specificity, transcription factor binding, enzyme remodeling, and Pol II elongation in positioning nucleosomes. It suggests that the differences in nucleosome density profiles observed at various transcription levels in yeast and mouse embryonic stem cells can be understood from a tug-of-war between two types of remodeling enzymes.

Histone mark recognition controls nucleosome translocation via a kinetic proofreading mechanism: Confronting theory and high-throughput experiments

Chromatin remodelers are multidomain enzymatic motor complexes that displace nucleosomes along DNA and hence "remodel chromatin structure," i.e., they dynamically reorganize nucleosome positions in both gene activation and gene repression. Recently, experimental insights from structural biology methods and remodeling assays have substantially advanced the understanding of these key chromatin components. Here we confront the kinetic proofreading scenario of chromatin remodeling, which proposes a mechanical link between histone residue modifications and the ATP-dependent action of remodelers, with recent experiments. We show that recent high-throughput data on nucleosome libraries assayed with remodelers from the Imitation Switch family are in accord with our earlier predictions of the kinetic proofreading scenario. We make suggestions for experimentally verifiable predictions of the kinetic proofreading scenarios for remodelers from other families.

Theory of Active Chromatin Remodeling

Nucleosome positioning controls the accessible regions of chromatin and plays essential roles in DNA-templated processes. ATP driven remodeling enzymes are known to be crucial for its establishment in vivo, but their nonequilibrium nature has hindered the development of a unified theoretical framework for nucleosome positioning. Using a perturbation theory, we show that the effect of these enzymes can be well approximated by effective equilibrium models with rescaled temperatures and interactions. Numerical simulations support the accuracy of the theory in predicting both kinetic and steady-state quantities, including the effective temperature and the radial distribution function, in biologically relevant regimes. The energy landscape view emerging from our study provides an intuitive understanding for the impact of remodeling enzymes in either reinforcing or overwriting intrinsic signals for nucleosome positioning, and may help improve the accuracy of computational models for its prediction in silico.

Correlation among DNA Linker Length, Linker Histone Concentration, and Histone Tails in Chromatin

Eukaryotic cells condense their genetic material in the nucleus in the form of chromatin, a macromolecular complex made of DNA and multiple proteins. The structure of chromatin is intimately connected to the regulation of all eukaryotic organisms, from amoebas to humans, but its organization remains largely unknown. The nucleosome repeat length (NRL) and the concentration of linker histones (ρLH) are two structural parameters that vary among cell types and cell cycles; the NRL is the number of DNA basepairs wound around each nucleosome core plus the number of basepairs linking successive nucleosomes. Recent studies have found a linear empirical relationship between the variation of these two properties for different cells, but its underlying mechanism remains elusive. Here we apply our established mesoscale chromatin model to explore the mechanisms responsible for this relationship, by investigating chromatin fibers as a function of NRL and ρLH combinations. We find that a threshold of linker histone concentration triggers the compaction of chromatin into well-formed 30-nm fibers; this critical value increases linearly with NRL, except for long NRLs, where the fibers remain disorganized. Remarkably, the interaction patterns between core histone tails and chromatin elements are highly sensitive to the NRL and ρLH combination, suggesting a molecular mechanism that could have a key role in regulating the structural state of the fibers in the cell. An estimate of the minimized work and volume associated with storage of chromatin fibers in the nucleus further suggests factors that could spontaneously regulate the NRL as a function of linker histone concentration. Both the tail interaction map and DNA packing considerations support the empirical NRL/ρLH relationship and offer a framework to interpret experiments for different chromatin conditions in the cell.

Theoretical estimates of exposure timescales of protein binding sites on DNA regulated by nucleosome kinetics

It is being increasingly realized that nucleosome organization on DNA crucially regulates DNA-protein interactions and the resulting gene expression. While the spatial character of the nucleosome positioning on DNA has been experimentally and theoretically studied extensively, the temporal character is poorly understood. Accounting for ATPase activity and DNA-sequence effects on nucleosome kinetics, we develop a theoretical method to estimate the time of continuous exposure of binding sites of non-histone proteins (e.g. transcription factors and TATA binding proteins) along any genome. Applying the method to Saccharomyces cerevisiae, we show that the exposure timescales are determined by cooperative dynamics of multiple nucleosomes, and their behavior is often different from expectations based on static nucleosome occupancy. Examining exposure times in the promoters of GAL1 and PHO5, we show that our theoretical predictions are consistent with known experiments. We apply our method genome-wide and discover huge gene-to-gene variability of mean exposure times of TATA boxes and patches adjacent to TSS (+1 nucleosome region); the resulting timescale distributions have non-exponential tails.

Role of transcription factor-mediated nucleosome disassembly in PHO5 gene expression

Studying nucleosome dynamics in promoter regions is crucial for understanding gene regulation. Nucleosomes regulate gene expression by sterically occluding transcription factors (TFs) and other non–histone proteins accessing genomic DNA. How the binding competition between nucleosomes and TFs leads to transcriptionally compatible promoter states is an open question. Here, we present a computational study of the nucleosome dynamics and organization in the promoter region of PHO5 gene in Saccharomyces cerevisiae. Introducing a model for nucleosome kinetics that takes into account ATP-dependent remodeling activity, DNA sequence effects, and kinetics of TFs (Pho4p), we compute the probability of obtaining different “promoter states” having different nucleosome configurations. Comparing our results with experimental data, we argue that the presence of local remodeling activity (LRA) as opposed to basal remodeling activity (BRA) is crucial in determining transcriptionally active promoter states. By modulating the LRA and Pho4p binding rate, we obtain different mRNA distributions—Poisson, bimodal, and long-tail. Through this work we explain many features of the PHO5 promoter such as sequence-dependent TF accessibility and the role of correlated dynamics between nucleosomes and TFs in opening/coverage of the TATA box. We also obtain possible ranges for TF binding rates and the magnitude of LRA.

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.

Affinity, stoichiometry and cooperativity of heterochromatin protein 1 (HP1) binding to nucleosomal arrays

Heterochromatin protein 1 (HP1) participates in establishing and maintaining heterochromatin via its histone-modification-dependent chromatin interactions. In recent papers HP1 binding to nucleosomal arrays was measured in vitro and interpreted in terms of nearest-neighbour cooperative binding. This mode of chromatin interaction could lead to the spreading of HP1 along the nucleosome chain. Here, we reanalysed previous data by representing the nucleosome chain as a 1D binding lattice and showed how the experimental HP1 binding isotherms can be explained by a simpler model without cooperative interactions between neighboring HP1 dimers. Based on these calculations and spatial models of dinucleosomes and nucleosome chains, we propose that binding stoichiometry depends on the nucleosome repeat length (NRL) rather than protein interactions between HP1 dimers. According to our calculations, more open nucleosome arrays with long DNA linkers are characterized by a larger number of binding sites in comparison to chains with a short NRL. Furthermore, we demonstrate by Monte Carlo simulations that the NRL dependent folding of the nucleosome chain can induce allosteric changes of HP1 binding sites. Thus, HP1 chromatin interactions can be modulated by the change of binding stoichiometry and the type of binding to condensed (methylated) and non-condensed (unmethylated) nucleosome arrays in the absence of direct interactions between HP1 dimers.

Regulation of the nucleosome repeat length in vivo by the DNA sequence, protein concentrations and long-range interactions

The nucleosome repeat length (NRL) is an integral chromatin property important for its biological functions. Recent experiments revealed several conflicting trends of the NRL dependence on the concentrations of histones and other architectural chromatin proteins, both in vitro and in vivo, but a systematic theoretical description of NRL as a function of DNA sequence and epigenetic determinants is currently lacking. To address this problem, we have performed an integrative biophysical and bioinformatics analysis in species ranging from yeast to frog to mouse where NRL was studied as a function of various parameters. We show that in simple eukaryotes such as yeast, a lower limit for the NRL value exists, determined by internucleosome interactions and remodeler action. For higher eukaryotes, also the upper limit exists since NRL is an increasing but saturating function of the linker histone concentration. Counterintuitively, smaller H1 variants or non-histone architectural proteins can initiate larger effects on the NRL due to entropic reasons. Furthermore, we demonstrate that different regimes of the NRL dependence on histone concentrations exist depending on whether DNA sequence-specific effects dominate over boundary effects or vice versa. We consider several classes of genomic regions with apparently different regimes of the NRL variation. As one extreme, our analysis reveals that the period of oscillations of the nucleosome density around bound RNA polymerase coincides with the period of oscillations of positioning sites of the corresponding DNA sequence. At another extreme, we show that although mouse major satellite repeats intrinsically encode well-defined nucleosome preferences, they have no unique nucleosome arrangement and can undergo a switch between two distinct types of nucleosome positioning.

Nucleosome positioning and kinetics near transcription-start-site barriers are controlled by interplay between active remodeling and DNA sequence

We investigate how DNA sequence, ATP-dependent chromatin remodeling and nucleosome-depleted ‘barriers’ co-operate to determine the kinetics of nucleosome organization, in a stochastic model of nucleosome positioning and dynamics. We find that ‘statistical’ positioning of nucleosomes against ‘barriers’, hypothesized to control chromatin structure near transcription start sites, requires active remodeling and therefore cannot be described using equilibrium statistical mechanics. We show that, unlike steady-state occupancy, DNA site exposure kinetics near a barrier is dominated by DNA sequence rather than by proximity to the barrier itself. The timescale for formation of positioning patterns near barriers is proportional to the timescale for active nucleosome eviction. We also show that there are strong gene-to-gene variations in nucleosome positioning near barriers, which are eliminated by averaging over many genes. Our results suggest that measurement of nucleosome kinetics can reveal information about sequence-dependent regulation that is not apparent in steady-state nucleosome occupancy.

Meeting the methodological challenges in molecular mapping of the embryonic epigenome

The past decade of life sciences research has been driven by progress in genomics. Many voices are already proclaiming the post-genomics era, in which phenomena other than sequence polymorphism influence gene expression and also explain complex phenotypes. One of these burgeoning fields is the study of the epigenome. Although the mechanisms by which chromatin structure and reorganization as well as cytosine methylation influence gene expression are not fully understood, they are being invoked to explain the now-accepted long-term impact of the environment on gene expression, which appears to be a factor in the development of numerous diseases. Such studies are particularly relevant in early embryonic development, during which waves of epigenetic reprogramming are known to have profound impacts. Since gametes and zygotes are in the process of resetting the genome in order to create embryonic stem cells that will each differentiate to create one of many specific tissue types, this phase of life is now viewed as a window of susceptibility to epigenetic reprogramming errors. Epigenetics could explain the influence of factors such as the nutritional/metabolic status of the mother or the artificial environment of assisted reproductive technologies. However, the peculiar nature of early embryos in addition to their scarcity poses numerous technological challenges that are slowly being overcome. The principal subject of this article is to review the suitability of various current and emerging technological platforms to study oocytes and early embryonic epigenome with more emphasis on studying DNA methylation. Furthermore, the constraint of samples size, inherent to the study of preimplantation embryo development, was put in perspective with the various molecular platforms described.