Depletion effects massively change chromatin properties and influence genome folding

Chromatin transitions triggered by LH density as epigenetic regulators of the genome

Motivated by experiments connecting linker histone (LH) deficiency to lymphoma progression and retinal disorders, we study by mesoscale chromatin modeling how LH density (ρ) induces gradual, as well sudden, changes in chromatin architecture and how the process depends on DNA linker length, LH binding dynamics and binding mode, salt concentration, tail modifications, and combinations of ρ and linker DNA length. We show that ρ tightly regulates the overall shape and compaction of the fiber, triggering a transition from an irregular disordered state to a compact and ordered structure. Such a structural transition, resembling B to A compartment transition connected with lymphoma of B cells, appears to occur around ρ = 0.5. The associated mechanism is DNA stem formation by LH binding, which is optimal when the lengths of the DNA linker and LH C-terminal domain are similar. Chromatin internal and external parameters are key regulators, promoting or impeding the transition. The LH density thus emerges as a critical tunable variable in controlling cellular functions through structural transitions of the genome.

Chromatin Compaction Multiscale Modeling: A Complex Synergy Between Theory, Simulation, and Experiment

Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.

Dynamic chromatin organization in the cell

The organization and regulation of genomic DNA as nuclear chromatin is necessary for proper DNA function inside living eukaryotic cells. While this has been extensively explored, no true consensus is currently reached regarding the exact mechanism of chromatin organization. The traditional view has assumed that the DNA is packaged into a hierarchy of structures inside the nucleus based on the regular 30-nm chromatin fiber. This is currently being challenged by the fluid-like model of the chromatin which views the chromatin as a dynamic structure based on the irregular 10-nm fiber. In this review, we focus on the recent progress in chromatin structure elucidation highlighting the paradigm shift in chromatin folding mechanism from the classical textbook perspective of the regularly folded chromatin to the more dynamic fluid-like perspective.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

Compaction of Single-Molecule Megabase-Long Chromatin under the Influence of Macromolecular Crowding

The megabase-sized length of chromatin is highly relevant to the state of chromatin in vivo, where it is subject to a highly crowded environment and is organized in topologically associating domains of similar dimension. We developed an in vitro experimental chromatin model system reconstituted from T4 DNA (approximately 166 kbp) and histone octamers and studied the monomolecular compaction of this megabase-sized chromatin fiber under the influence of macromolecular crowding. We used single-molecule fluorescence microscopy and observed compaction in aqueous solutions containing poly(ethylene glycol) in the presence of monovalent (Na+ and K+) and divalent (Mg2+) cations. Both DNA and chromatin demonstrated compaction under comparable conditions in the presence of poly(ethylene glycol) and Na+ or Mg2+ salt. However, the mechanism of the compaction changed from a first-order phase transition for DNA to a continuous folding for megabase-sized chromatin fibers. A more efficient and pronounced chromatin compaction was observed in the presence of Na+ compared to K+. A flow-stretching technique to unfold DNA and chromatin coils was used to gain further insight into the morphology of partially folded chromatin fibers. The results revealed a distribution of partially folded chromatin fibers. This variability is likely the result of the heterogeneous distribution of nucleosomes on the DNA chain. The packaging of DNA in the form of chromatin in the crowded nuclear environment appears essential to ensure gradual conformational changes of DNA.

Entropic effect of macromolecular crowding enhances binding between nucleosome clutches in heterochromatin, but not in euchromatin

Sharp increase in macromolecular crowding induces abnormal chromatin compaction in the cell nucleus, suggesting its non-negligible impact on chromatin structure and function. However, the details of the crowding-induced chromatin compaction remain poorly understood. In this work, we present a computer simulation study on the entropic effect of macromolecular crowding on the interaction between chromatin structural units called nucleosome clutches. Nucleosome clutches were modeled by a chain of nucleosomes collapsed by harmonic restraints implicitly mimicking the nucleosome association mediated by histone tails and linker histones. The nucleosome density of the clutches was set close to either that of high-density heterochromatin or that of low-density euchromatin. The effective interactions between these nucleosome clutches were calculated in various crowding conditions, and it was found that the increase in the degree of macromolecular crowding induced attractive interaction between two clutches with high nucleosome density. Interestingly, the increased degree of macromolecular crowding did not induce any attraction between two clutches with low nucleosome density. Our results suggest that the entropic effect of macromolecular crowding can enhance binding between nucleosome clutches in heterochromatin, but not in euchromatin, as a result of the difference in nucleosome packing degrees.

DNA Methylation Landscape Reflects the Spatial Organization of Chromatin in Different Cells

The relationship between DNA methylation and chromatin structure is still largely unknown. By analyzing a large set of published sequencing data, we observed a long-range power law correlation of DNA methylation with cell class-specific scaling exponents in the range of tens of kilobases. We showed that such cell class-specific scaling exponents are caused by different patchiness of DNA methylation in different cells. By modeling the chromatin structure using high-resolution chromosome conformation capture data and mapping the methylation level onto the modeled structure, we demonstrated that the patchiness of DNA methylation is related to chromatin structure. The scaling exponents of the power law correlation are thus a display of the spatial organization of chromatin. Besides the long-range correlation, we also showed that the local correlation of DNA methylation is associated with nucleosome positioning. The local correlation of partially methylated domains is different from that of nonpartially methylated domains, suggesting that their chromatin structures differ at the scale of several hundred base pairs (covering a few nucleosomes). Our study provides a novel, to our knowledge, view of the spatial organization of chromatin structure from a perspective of DNA methylation, in which both long-range and local correlations of DNA methylation along the genome reflect the spatial organization of chromatin.

Confinement-driven organization of a histone-complexed DNA molecule in a dense array of nanoposts

The first step in the controlled storage of lengthy DNA molecules is to keep DNA molecules separated while integrated in micrometer-sized space. Herein, we present hybrid Monte Carlo simulations of a histone-complexed DNA (hcDNA) molecule confined in a dense array of nanoposts. Depending on the nanopost dimension, a single, 8.7 kilobase pair hcDNA molecule was either localized and elongated in a single inter-post space surrounded by four nanoposts or spread over several inter-post spaces through passages between two neighboring nanoposts. The conformational change of a hcDNA molecule is interpreted in terms of competitive effects of confinements in the inter-post and passage spaces. We propose that, by elaborately designing nanopost arrays, the competitive confinement effects can be adjusted such that each hcDNA molecule is localized in a single inter-post space, and thereby multiple hcDNA molecules can be physically separated from each other while stored together in the nanopost array.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

RESULTS

We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

CONCLUSION

The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Elucidating internucleosome interactions and the roles of histone tails

The nucleosome is the first level of genome organization and regulation in eukaryotes where negatively charged DNA is wrapped around largely positively charged histone proteins. Interaction between nucleosomes is dominated by electrostatics at long range and guided by specific contacts at short range, particularly involving their flexible histone tails. We have thus quantified how internucleosome interactions are modulated by salts (KCl, MgCl2) and histone tail deletions (H3, H4 N-terminal), using small-angle x-ray scattering and theoretical modeling. We found that measured effective charges at low salts are ∼1/5th of the theoretically predicted renormalized charges and that H4 tail deletion suppresses the attraction at high salts to a larger extent than H3 tail deletion.

How chromatin looping and nuclear envelope attachment affect genome organization in eukaryotic cell nuclei

To understand how interphase chromatin is organized in eukaryotic cell nuclei, it is essential to understand what kind of interactions influence the nuclear architecture and to what extent. Using a mesoscale model that incorporates chromatin-chromatin interactions as well as binding of chromatin to the nuclear envelope, we can show that chromatin loops and envelope bonds are major players in genome organization because they largely affect the entropy of the chromatin fibres. The model allows us to consistently reproduce multiple characteristic chromatin parameters in agreement with experimental data. We focus on the question of how an inversion of the nuclear architecture, in the course of which the highly active euchromatin changes its preferential position from the nuclear center to the periphery, can be achieved. We find that the transition between the common and inverted organization is driven by the strength of the envelope interaction and the nuclear chromatin density.

Chromatin organization and transcriptional regulation

Cell type specific transcriptional regulation must be adhered to in order to maintain cell identity throughout the lifetime of an organism, yet it must be flexible enough to allow for responses to endogenous and exogenous stimuli. This regulation is mediated not only by molecular factors (e.g. cell type specific transcription factors, histone and DNA modifications), but also on the level of chromatin and genome organization. In this review we focus on recent findings that have contributed to our understanding of higher order chromatin structure and genome organization within the nucleus. We highlight new findings on the dynamic positioning of genes relative to each other, as well as to their chromosome territory and the nuclear lamina, and how the position of genes correlates with their transcriptional activity.

Human mitotic chromosomes consist predominantly of irregularly folded nucleosome fibres without a 30-nm chromatin structure

How a long strand of genomic DNA is compacted into a mitotic chromosome remains one of the basic questions in biology. The nucleosome fibre, in which DNA is wrapped around core histones, has long been assumed to be folded into a 30-nm chromatin fibre and further hierarchical regular structures to form mitotic chromosomes, although the actual existence of these regular structures is controversial. Here, we show that human mitotic HeLa chromosomes are mainly composed of irregularly folded nucleosome fibres rather than 30-nm chromatin fibres. Our comprehensive and quantitative study using cryo-electron microscopy and synchrotron X-ray scattering resolved the long-standing contradictions regarding the existence of 30-nm chromatin structures and detected no regular structure >11 nm. Our finding suggests that the mitotic chromosome consists of irregularly arranged nucleosome fibres, with a fractal nature, which permits a more dynamic and flexible genome organization than would be allowed by static regular structures.

Salt-modulated structure of polyelectrolyte-macroion complex fibers

The structure and stability of strongly charged complex fibers, formed by complexation of a single long semi-flexible polyelectrolyte chain and many oppositely charged spherical macroions, are investigated numerically at the ground-state level using a chain-sphere cell model. The model takes into account chain elasticity as well as electrostatic interactions between charged spheres and chain segments. Using a numerical optimization method based on a periodically repeated unit cell, we obtain fiber configurations that minimize the total energy. The optimal fiber configurations exhibit a variety of helical structures for the arrangement of macroions including zig-zag, solenoidal and beads-on-a-string patterns. These structures result from the competition between attraction between spheres and the polyelectrolyte chain (which favors chain wrapping around the spheres), chain bending rigidity and electrostatic repulsion between chain segments (which favor unwrapping of the chain), and the interactions between neighboring sphere-chain complexes which can be attractive or repulsive depending on the system parameters such as salt concentration, macroion charge and chain length per macroion (linker size). At about physiological salt concentration, dense zig-zag patterns are found to be energetically most stable when parameters appropriate for the DNA-histone system in the chromatin fiber are adopted. In fact, the predicted fiber diameter in this regime is found to be around 30 nanometers, which roughly agrees with the thickness observed in in vitro experiments on chromatin. We also find a macroion (histone) density of 5-6 per 11nm which agrees with results from the zig-zag or cross-linker models of chromatin. Since our study deals primarily with a generic chain-sphere model, these findings suggest that structures similar to those found for chromatin should also be observable for polyelectrolyte-macroion complexes formed in solutions of DNA and synthetic nano-colloids of opposite charge. In the ensemble where the mean linear density of spheres on the chain is fixed, the present model predicts a phase separation at intermediate salt concentrations into a densely packed complex phase and a dilute phase.

Hierarchies in eukaryotic genome organization: Insights from polymer theory and simulations

Eukaryotic genomes possess an elaborate and dynamic higher-order structure within the limiting confines of the cell nucleus. Knowledge of the physical principles and the molecular machinery that govern the 3D organization of this structure and its regulation are key to understanding the relationship between genome structure and function. Elegant microscopy and chromosome conformation capture techniques supported by analysis based on polymer models are important steps in this direction. Here, we review results from these efforts and provide some additional insights that elucidate the relationship between structure and function at different hierarchical levels of genome organization.

Repulsive forces between looping chromosomes induce entropy-driven segregation

One striking feature of chromatin organization is that chromosomes are compartmentalized into distinct territories during interphase, the degree of intermingling being much smaller than expected for linear chains. A growing body of evidence indicates that the formation of loops plays a dominant role in transcriptional regulation as well as the entropic organization of interphase chromosomes. Using a recently proposed model, we quantitatively determine the entropic forces between chromosomes. This Dynamic Loop Model assumes that loops form solely on the basis of diffusional motion without invoking other long-range interactions. We find that introducing loops into the structure of chromatin results in a multi-fold higher repulsion between chromosomes compared to linear chains. Strong effects are observed for the tendency of a non-random alignment; the overlap volume between chromosomes decays fast with increasing loop number. Our results suggest that the formation of chromatin loops imposes both compartmentalization as well as order on the system without requiring additional energy-consuming processes.

Histone depletion facilitates chromatin loops on the kilobasepair scale

The packing of eukaryotic DNA in the nucleus is decisive for its function; for instance, contact between remote genome sites constitutes a basic feature of gene regulation. Interactions among regulatory proteins, DNA binding, and transcription activation are facilitated by looping of the intervening chromatin. Such long-range interactions depend on the bending flexibility of chromatin, i.e., the ring-closure probability is a directly measurable indicator of polymer flexibility. The applicability of a wormlike chain model to naked DNA has been widely accepted. However, whether this model also suffices to describe the flexibility of eukaryotic interphase chromatin is still a matter of discussion. Here we compare both 5C data from a gene desert and data from fluorescence in situ hybridization with the results of a Monte Carlo simulation of chromatin fibers with and without histone depletion. We then estimate the ring-closure probabilities of simulated fibers with estimates from analytical calculations and show that the wormlike chain model grossly underestimates chromatin flexibility for sharp bends. Most importantly, we find that only fibers with random depletion of linker histones or nucleosomes can explain the probability of random chromatin contacts on small length scales that play an important role in gene regulation. It is possible that missing linker histones and nucleosomes are not just simple, unavoidable, randomly occurring defects, but instead play a regulatory role in gene expression.

Diffusion-driven looping provides a consistent framework for chromatin organization

Chromatin folding inside the interphase nucleus of eukaryotic cells is done on multiple scales of length and time. Despite recent progress in understanding the folding motifs of chromatin, the higher-order structure still remains elusive. Various experimental studies reveal a tight connection between genome folding and function. Chromosomes fold into a confined subspace of the nucleus and form distinct territories. Chromatin looping seems to play a dominant role both in transcriptional regulation as well as in chromatin organization and has been assumed to be mediated by long-range interactions in many polymer models. However, it remains a crucial question which mechanisms are necessary to make two chromatin regions become co-located, i.e. have them in spatial proximity. We demonstrate that the formation of loops can be accomplished solely on the basis of diffusional motion. The probabilistic nature of temporary contacts mimics the effects of proteins, e.g. transcription factors, in the solvent. We establish testable quantitative predictions by deriving scale-independent measures for comparison to experimental data. In this Dynamic Loop (DL) model, the co-localization probability of distant elements is strongly increased compared to linear non-looping chains. The model correctly describes folding into a confined space as well as the experimentally observed cell-to-cell variation. Most importantly, at biological densities, model chromosomes occupy distinct territories showing less inter-chromosomal contacts than linear chains. Thus, dynamic diffusion-based looping, i.e. gene co-localization, provides a consistent framework for chromatin organization in eukaryotic interphase nuclei.

Condensed DNA: condensing the concepts

DNA is stored in vivo in a highly compact, so-called condensed phase, where gene regulatory processes are governed by the intricate interplay between different states of DNA compaction. These systems often have surprising properties, which one would not predict from classical concepts of dilute solutions. The mechanistic details of DNA packing are essential for its functioning, as revealed by the recent developments coming from biochemistry, electrostatics, statistical mechanics, and molecular and cell biology. Different aspects of condensed DNA behavior are linked to each other, but the links are often hidden in the bulk of experimental and theoretical details. Here we try to condense some of these concepts and provide interconnections between the different fields. After a brief description of main experimental features of DNA condensation inside viruses, bacteria, eukaryotes and the test tube, main theoretical approaches for the description of these systems are presented. We end up with an extended discussion of the role of DNA condensation in the context of gene regulation and mention potential applications of DNA condensation in gene therapy and biotechnology.

Monte Carlo Simulations indicate that Chromati: Nanostructure is accessible by Light Microscopy

A long controversy exists about the structure of chromatin. Theoretically, this structure could be resolved by scattering experiments if one determines the scattering function - or equivalently the pair distribution function - of the nucleosomes. Unfortunately, scattering experiments with live cells are very difficult and limited to only a couple of nucleosomes.Nevertheless, new techniques like the high-resolution light microscopy supply a new approach to this problem. In this work we determine the radial pair distribution function of chromatin described by our E2A model and find that the dominant peaks which characterize the chromatin structure are very robust in several ways: They can still be identified in the case of chromatin fibers with reasonable linker histone and nucleosome defect rates as well as in the 2D case after a projection like in most high-res light microscopy experiments. This might initiate new experimental approaches like optical microscopy to finally determine the nanostructure of chromatin.Furthermore, we examine the statistics of random chromatin collisions and compare it with 5C data of a gene desert. We find that only chromatin fibers with histone depletion show a significant amount of contacts on the kbp-scale which play a important role in gene regulation. Therefore, linker histone and nucleosome depletion might not only be chromatin defects but even be necessary to facilitate transcription.PACS codes: 82.35.Pq, 87.16.A-, 87.16.af.