Topological diversity of chromatin fibers: Interplay between nucleosome repeat length, DNA linking number and the level of transcription

Histone N-tails modulate sequence-specific positioning of nucleosomes

Spatial organization of chromatin is essential for cellular functioning. However, the precise mechanisms governing sequence-dependent positioning of nucleosomes on DNA remain unknown in detail. Existing algorithms, considering the sequence-dependent deformability of DNA and its interactions with the histone globular domains, predict rotational setting of only 65% of human nucleosomes mapped in vivo. To uncover additional factors responsible for the nucleosome positioning, we analyzed potential involvement of the histone N-tails in this process. To this aim, we reconstituted the H2A/H4 N-tailless nucleosomes on human BRCA1 DNA (∼100 kb) and compared their positions and sequences with those of the wild-type nucleosomes. We found that removal of the histone N-tails promoted displacement of the predominant positions of nucleosomes, accompanied by redistribution of the AT-rich and GC-rich motifs in nucleosome sequences. Importantly, most of these sequence changes occurred at superhelical locations (SHLs) ±4, ±1, and ± 2, where the H2A and H4 N-tails interact with the DNA minor grooves. Furthermore, a substantial number of H4-tailless nucleosomes exhibit rotational settings opposite to that of the wild-type nucleosomes, the effect known to change the topological properties of chromatin fiber. Thus, the histone N-tails are operative in the selection of nucleosome positions, which may have wide-ranging implications for epigenetic modulation of chromatin states.

Histone H3 core domain in chromatin with different DNA linker lengths studied by 1H-Detected solid-state NMR spectroscopy

Chromatin, a dynamic protein-DNA complex that regulates eukaryotic genome accessibility and essential functions, is composed of nucleosomes connected by linker DNA with each nucleosome consisting of DNA wrapped around an octamer of histones H2A, H2B, H3 and H4. Magic angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy can yield unique insights into histone structure and dynamics in condensed nucleosomes and nucleosome arrays representative of chromatin at physiological concentrations. Recently we used J-coupling-based solid-state NMR methods to investigate with residue-specific resolution the conformational dynamics of histone H3 N-terminal tails in 16-mer nucleosome arrays containing 15, 30 or 60 bp DNA linkers. Here, we probe the H3 core domain in the 16-mer arrays as a function of DNA linker length via dipolar coupling-based 1H-detected solid-state NMR techniques. Specifically, we established nearly complete assignments of backbone chemical shifts for H3 core residues in arrays with 15-60 bp DNA linkers reconstituted with 2H,13C,15N-labeled H3. Overall, these chemical shifts were similar irrespective of the DNA linker length indicating no major changes in H3 core conformation. Notably, however, multiple residues at the H3-nucleosomal DNA interface in arrays with 15 bp DNA linkers exhibited relatively pronounced differences in chemical shifts and line broadening compared to arrays with 30 and 60 bp linkers. These findings are consistent with increased heterogeneity in nucleosome packing and structural strain within arrays containing short DNA linkers that likely leads to side-chains of these interfacial residues experiencing alternate conformations or shifts in their rotamer populations relative to arrays with the longer DNA linkers.

Epigenetic liquid biopsies for minimal residual disease, what's around the corner?

Liquid biopsy assays for minimal residual disease (MRD) are used to monitor and inform oncological treatment and predict the risk of relapse in cancer patients. To-date, most MRD assay development has focused on targeting somatic mutations. However, epigenetic changes are more frequent and universal than genetic alterations in cancer and circulating tumor DNA (ctDNA) retains much of these changes. Here, we review the epigenetic signals that can be used to detect MRD, including DNA methylation alterations and fragmentation patterns that differentiate ctDNA from noncancerous circulating cell-free DNA (ccfDNA). We then summarize the current state of MRD monitoring; highlight the advantages of epigenetics over genetics-based approaches; and discuss the emerging paradigm of assaying both genetic and epigenetic targets to monitor treatment response, detect disease recurrence, and inform adjuvant therapy.

Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies

Epigenetic research holds great promise to advance our understanding of biomarkers and regulatory processes in health and disease. An increasing number of new approaches, ranging from molecular to biophysical analyses, enable identifying epigenetic changes on the level of a single gene or the whole epigenome. The aim of this review is to highlight how the field is shifting from completely molecular-biology-driven solutions to multidisciplinary strategies including more reliance on biophysical analysis tools. Biophysics not only offers technical advancements in imaging or structure analysis but also helps to explore regulatory interactions. New computational methods are also being developed to meet the demand of growing data volumes and their processing. Therefore, it is important to capture these new directions in epigenetics from a biophysical perspective and discuss current challenges as well as multiple applications of biophysical methods and tools. Specifically, we gradually introduce different biophysical research methods by first considering the DNA-level information and eventually higher-order chromatin structures. Moreover, we aim to highlight that the incorporation of bioinformatics, machine learning, and artificial intelligence into biophysical analysis allows gaining new insights into complex epigenetic processes. The gained understanding has already proven useful in translational and clinical research providing better patient stratification options or new therapeutic insights. Together, this offers a better readiness to transform bench-top experiments into industrial high-throughput applications with a possibility to employ developed methods in clinical practice and diagnostics.

Conformational Dynamics of Histone H3 Tails in Chromatin

Chromatin is a supramolecular DNA-protein complex that compacts eukaryotic genomes and regulates their accessibility and functions. Dynamically disordered histone H3 N-terminal tails are among key chromatin regulatory components. Here, we used high-resolution-magic-angle-spinning NMR measurements of backbone amide ¹⁵N spin relaxation rates to investigate, with residue-specific detail, the dynamics and interactions of H3 tails in recombinant ¹³C,¹⁵N-enriched nucleosome arrays containing 15, 30, or 60 bp linker DNA between the nucleosome repeats. These measurements were compared to analogous data available for mononucleosomes devoid of linker DNA or containing two 20 bp DNA overhangs. The H3 tail dynamics in nucleosome arrays were found to be considerably attenuated compared with nucleosomes with or without linker DNA due to transient electrostatic interactions with the linker DNA segments and the structured chromatin environment. Remarkably, however, the H3 tail dynamics were not modulated by the specific linker DNA length within the 15-60 bp range investigated here.

Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface

The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome. Dynamical simulations performed with an advanced sampling technique uncovered multiple pathways that connect open chromatin configurations with the zigzag crystal structure. Intermediate states along the simulated folding pathways resemble chromatin configurations reported from in situ experiments. We further determined a six-dimensional free energy surface as a function of the inter-nucleosome distances via a deep learning approach. The zigzag structure can indeed be seen as the global minimum of the surface. However, it is not favored by a significant amount relative to the partially unfolded, in situ configurations. Chemical perturbations such as histone H4 tail acetylation and thermal fluctuations can further tilt the energetic balance to stabilize intermediate states. Our study provides insight into the connection between various reported chromatin configurations and has implications on the in situ relevance of the 30 nm fiber.

The three-dimensional organization of chromatin plays critical roles in regulating genome function. Here the authors apply a near atomistic model to study the structure and dynamics of the chromatin folding unit - the tetra-nucleosome - to provide insight into how chromatin folds.

Topological polymorphism of nucleosome fibers and folding of chromatin

We discuss recent observations of polymorphic chromatin packaging at the oligonucleosomal level and compare them with computer simulations. Our computations reveal two topologically different families of two-start 30-nm fiber conformations distinguished by the linker length L; fibers with L ≈ 10n and L ≈ 10n+5 basepairs have DNA linking numbers per nucleosome of ΔLk ≈ -1.5 and -1.0, respectively (where n is a natural number). Although fibers with ΔLk ≈ -1.5 were observed earlier, the topoisomer with ΔLk ≈ -1.0 is novel. These predictions were confirmed experimentally for circular nucleosome arrays with precisely positioned nucleosomes. We suggest that topological polymorphism of chromatin may play a role in transcription, with the {10n+5} fibers producing transcriptionally competent chromatin structures. This hypothesis is consistent with available data for yeast and, partially, for fly. We show that both fiber topoisomers (with ΔLk ≈ -1.5 and -1.0) have to be taken into account to interpret experimental data obtained using new techniques: genome-wide Micro-C, Hi-CO, and RICC-seq, as well as self-association of nucleosome arrays in vitro. The relative stability of these topoisomers is likely to depend on epigenetic histone modifications modulating the strength of internucleosome interactions. Potentially, our findings may reflect a general tendency of functionally distinct parts of the genome to retain topologically different higher-order structures.

Single-molecule and in silico dissection of the interaction between Polycomb repressive complex 2 and chromatin

Polycomb repressive complex 2 (PRC2) installs and spreads repressive histone methylation marks on eukaryotic chromosomes. Because of the key roles that PRC2 plays in development and disease, how this epigenetic machinery interacts with DNA and nucleosomes is of major interest. Nonetheless, the mechanism by which PRC2 engages with native-like chromatin remains incompletely understood. In this work, we employ single-molecule force spectroscopy and molecular dynamics simulations to dissect the behavior of PRC2 on polynucleosome arrays. Our results reveal an unexpectedly diverse repertoire of PRC2 binding configurations on chromatin. Besides reproducing known binding modes in which PRC2 interacts with bare DNA, mononucleosomes, and adjacent nucleosome pairs, our data also provide direct evidence that PRC2 can bridge pairs of distal nucleosomes. In particular, the "1-3" bridging mode, in which PRC2 engages two nucleosomes separated by one spacer nucleosome, is a preferred low-energy configuration. Moreover, we show that the distribution and stability of different PRC2-chromatin interaction modes are modulated by accessory subunits, oncogenic histone mutations, and the methylation state of chromatin. Overall, these findings have implications for the mechanism by which PRC2 spreads histone modifications and compacts chromatin. The experimental and computational platforms developed here provide a framework for understanding the molecular basis of epigenetic maintenance mediated by Polycomb-group proteins.

Braiding topology and the energy landscape of chromosome organization proteins

Assemblies of structural maintenance of chromosomes (SMC) proteins and kleisin subunits are essential to chromosome organization and segregation across all kingdoms of life. While structural data exist for parts of the SMC-kleisin complexes, complete structures of the entire complexes have yet to be determined, making mechanistic studies difficult. Using an integrative approach that combines crystallographic structural information about the globular subdomains, along with coevolutionary information and an energy landscape optimized force field (AWSEM), we predict atomic-scale structures for several tripartite SMC-kleisin complexes, including prokaryotic condensin, eukaryotic cohesin, and eukaryotic condensin. The molecular dynamics simulations of the SMC-kleisin protein complexes suggest that these complexes exist as a broad conformational ensemble that is made up of different topological isomers. The simulations suggest a critical role for the SMC coiled-coil regions, where the coils intertwine with various linking numbers. The twist and writhe of these braided coils are coupled with the motion of the SMC head domains, suggesting that the complexes may function as topological motors. Opening, closing, and translation along the DNA of the SMC-kleisin protein complexes would allow these motors to couple to the topology of DNA when DNA is entwined with the braided coils.

Sensitive effect of linker histone binding mode and subtype on chromatin condensation

The complex role of linker histone (LH) on chromatin compaction regulation has been highlighted by recent discoveries of the effect of LH binding variability and isoforms on genome structure and function. Here we examine the effect of two LH variants and variable binding modes on the structure of chromatin fibers. Our mesoscale modeling considers oligonucleosomes with H1C and H1E, bound in three different on and off-dyad modes, and spanning different LH densities (0.5–1.6 per nucleosome), over a wide range of physiologically relevant nucleosome repeat lengths (NRLs). Our studies reveal an LH-variant and binding-mode dependent heterogeneous ensemble of fiber structures with variable packing ratios, sedimentation coefficients, and persistence lengths. For maximal compaction, besides dominantly interacting with parental DNA, LHs must have strong interactions with nonparental DNA and promote tail/nonparental core interactions. An off-dyad binding of H1E enables both; others compromise compaction for bendability. We also find that an increase of LH density beyond 1 is best accommodated in chromatosomes with one on-dyad and one off-dyad LH. We suggest that variable LH binding modes and concentrations are advantageous, allowing tunable levels of chromatin condensation and DNA accessibility/interactions. Thus, LHs add another level of epigenetic regulation of chromatin.

Modelling and DNA topology of compact 2-start and 1-start chromatin fibres

We have investigated the structure of the most compact 30-nm chromatin fibres by modelling those with 2-start or 1-start crossed-linker organisations. Using an iterative procedure we obtained possible structural solutions for fibres of the highest possible compaction permitted by physical constraints, including the helical repeat of linker DNA. We find that this procedure predicts a quantized nucleosome repeat length (NRL) and that only fibres with longer NRLs (≥197 bp) can more likely adopt the 1-start organisation. The transition from 2-start to 1-start fibres is consistent with reported differing binding modes of the linker histone. We also calculate that in 1-start fibres the DNA constrains more torsion (as writhe) than 2-start fibres with the same NRL and that the maximum constraint obtained is in accord with previous experimental results. We posit that the coiling of the fibre is driven by overtwisting of linker DNA which, in the most compact forms - for example, in echinoderm sperm and avian erythrocytes - could adopt a helical repeat of ∼10 bp/turn. We argue that in vivo the total twist of linker DNA could be modulated by interaction with other abundant chromatin-associated proteins and by epigenetic modifications of the C-terminal tail of linker histones.

Dynamics of Chromatin Fibers: Comparison of Monte Carlo Simulations with Force Spectroscopy

To elucidate conformational dynamics of chromatin fibers, we compared available force-spectroscopy measurements with extensive Monte Carlo simulations of nucleosome arrays under external force. Our coarse-grained model of chromatin includes phenomenological energy terms for the DNA-histone adhesion and the internucleosome stacking interactions. We found that the Monte Carlo fiber ensembles simulated with increasing degrees of DNA unwrapping and the stacking energy 8 kT can account for the intricate force-extension response observed experimentally. Our analysis shows that at low external forces (F < 3.0 picoNewtons), the DNA ends in nucleosomes breathe by ∼10 bp. Importantly, under these conditions, the fiber is highly dynamic, exhibiting continuous unstacking-restacking transitions, allowing accessibility of transcription factors to DNA, while maintaining a relatively compact conformation. Of note, changing the stacking interaction by a few kT, an in silico way to mimic histone modifications, is sufficient to transform an open chromatin state into a compact fiber. The fibers are mostly two-start zigzag folds with rare occurrences of three- to five-start morphologies. The internucleosome stacking is lost during the linear response regime. At the higher forces exceeding 4 picoNewtons, the nucleosome unwrapping becomes stochastic and asymmetric, with one DNA arm opened by ∼55 bp and the other arm only by ∼10 bp. Importantly, this asymmetric unwrapping occurs for any kind of sequence, including the symmetric ones. Our analysis brings new, to our knowledge, insights in dynamics of chromatin modulated by histone epigenetic modifications and molecular motors such as RNA polymerase.

Precise genome-wide mapping of single nucleosomes and linkers in vivo

We developed a chemical cleavage method that releases single nucleosome dyad-containing fragments, allowing us to precisely map both single nucleosomes and linkers with high accuracy genome-wide in yeast. Our single nucleosome positioning data reveal that nucleosomes occupy preferred positions that differ by integral multiples of the DNA helical repeat. By comparing nucleosome dyad positioning maps to existing genomic and transcriptomic data, we evaluated the contributions of sequence, transcription, and histones H1 and H2A.Z in defining the chromatin landscape. We present a biophysical model that neglects DNA sequence and shows that steric occlusion suffices to explain the salient features of nucleosome positioning.

DNA topology in chromatin is defined by nucleosome spacing

In eukaryotic nucleosomes, DNA makes ~1.7 superhelical turns around histone octamer. However, there is a long-standing discrepancy between the nucleosome core structure determined by x-ray crystallography and measurements of DNA topology in circular minichromosomes, indicating that there is only ~1.0 superhelical turn per nucleosome. Although several theoretical assumptions were put forward to explain this paradox by conformational variability of the nucleosome linker, none was tested experimentally. We analyzed topological properties of DNA in circular nucleosome arrays with precisely positioned nucleosomes. Using topological electrophoretic assays and electron microscopy, we demonstrate that the DNA linking number per nucleosome strongly depends on the nucleosome spacing and varies from -1.4 to -0.9. For the predominant {10n + 5} class of nucleosome repeats found in native chromatin, our results are consistent with the DNA topology observed earlier. Thus, we reconcile the topological properties of nucleosome arrays with nucleosome core structure and provide a simple explanation for the DNA topology in native chromatin with variable DNA linker length. Topological polymorphism of the chromatin fibers described here may reflect a more general tendency of chromosomal domains containing active or repressed genes to acquire different nucleosome spacing to retain topologically distinct higher-order structures.

Single-molecule force spectroscopy on histone H4 tail-cross-linked chromatin reveals fiber folding

The eukaryotic genome is highly compacted into a protein-DNA complex called chromatin. The cell controls access of transcriptional regulators to chromosomal DNA via several mechanisms that act on chromatin-associated proteins and provide a rich spectrum of epigenetic regulation. Elucidating the mechanisms that fold chromatin fibers into higher-order structures is therefore key to understanding the epigenetic regulation of DNA accessibility. Here, using histone H4-V21C and histone H2A-E64C mutations, we employed single-molecule force spectroscopy to measure the unfolding of individual chromatin fibers that are reversibly cross-linked through the histone H4 tail. Fibers with covalently linked nucleosomes featured the same folding characteristics as fibers containing wild-type histones but exhibited increased stability against stretching forces. By stabilizing the secondary structure of chromatin, we confirmed a nucleosome repeat length (NRL)-dependent folding. Consistent with previous crystallographic and cryo-EM studies, the obtained force-extension curves on arrays with 167-bp NRLs best supported an underlying structure consisting of zig-zag, two-start fibers. For arrays with 197-bp NRLs, we previously inferred solenoidal folding, which was further corroborated by force-extension curves of the cross-linked fibers. The different unfolding pathways exhibited by these two types of arrays and reported here extend our understanding of chromatin structure and its potential roles in gene regulation. Importantly, these findings imply that chromatin compaction by nucleosome stacking protects nucleosomal DNA from external forces up to 4 piconewtons.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.