Chromatin fiber polymorphism triggered by variations of DNA linker lengths

Nanoscale analysis of human G1 and metaphase chromatin in situ

The structure of chromatin at the nucleosome level inside cells is still incompletely understood. Here we present in situ electron cryotomography analyses of chromatin in both G1 and metaphase RPE-1 cells. G1 nucleosomes are concentrated in globular chromatin domains, and metaphase nucleosomes are concentrated in the chromatids. Classification analysis reveals that canonical mononucleosomes, and in some conditions ordered stacked dinucleosomes and mononucleosomes with a disordered gyre-proximal density, are abundant in both cell-cycle states. We do not detect class averages that have more than two stacked nucleosomes or side-by-side dinucleosomes, suggesting that groups of more than two nucleosomes are heterogeneous. Large multi-megadalton structures are abundant in G1 nucleoplasm, but not found in G1 chromatin domains and metaphase chromatin. The macromolecular phenotypes studied here represent a starting point for the comparative analysis of compaction in normal vs. unhealthy human cells, in other cell-cycle states, other organisms, and in vitro chromatin assemblies.

The shifting paradigm of chromatin structure: from the 30-nm chromatin fiber to liquid-like organization

The organization and dynamics of chromatin are critical for genome functions such as transcription and DNA replication/repair. Historically, chromatin was assumed to fold into the 30-nm fiber and progressively arrange into larger helical structures, as described in the textbook model. However, over the past 15 years, extensive evidence including our studies has dramatically transformed the view of chromatin from a static, regular structure to one that is more variable and dynamic. In higher eukaryotic cells, chromatin forms condensed yet liquid-like domains, which appear to be the basic unit of chromatin structure, replacing the 30-nm fiber. These domains maintain proper accessibility, ensuring the regulation of DNA reaction processes. During mitosis, these domains assemble to form more gel-like mitotic chromosomes, which are further constrained by condensins and other factors. Based on the available evidence, I discuss the physical properties of chromatin in live cells, emphasizing its viscoelastic nature-balancing local fluidity with global stability to support genome functions.

Unveiling nucleosome dynamics: A comparative study using all-atom and coarse-grained simulations enhanced by principal component analysis

The conformational dynamics of the DNA in the nucleosome may play a role in governing gene regulation and accessibility and impact higher-order chromatin structure. This study investigates nucleosome dynamics using both all-atom and coarse-grained (CG) molecular dynamics simulations, focusing on the SIRAH force field. Simulations are performed for two nucleosomal DNA sequences-alpha satellite palindromic and Widom-601-over 6 μs at physiological salt concentrations. A comparative analysis of structural parameters, such as groove widths and base pair geometries, reveals good agreement between atomistic and CG models, although CG simulations exhibit broader conformational sampling and greater breathing motion of DNA ends. Principal component analysis is applied to DNA structural parameters, revealing multiple free energy minima, especially in CG simulations. These findings highlight the potential of the SIRAH CG force field for studying large-scale nucleosome dynamics, offering insights into DNA repositioning and sequence-dependent behavior.

Beyond the mono-nucleosome

Nucleosomes, the building block of chromatin, are responsible for regulating access to the DNA sequence. This control is critical for essential cellular processes, including transcription and DNA replication and repair. Studying chromatin can be challenging both in vitro and in vivo, leading many to use a mono-nucleosome system to answer fundamental questions relating to chromatin regulators and binding partners. However, the mono-nucleosome fails to capture essential features of the chromatin structure, such as higher-order chromatin folding, local nucleosome-nucleosome interactions, and linker DNA trajectory and flexibility. We briefly review significant discoveries enabled by the mono-nucleosome and emphasize the need to go beyond this model system in vitro. Di-, tri-, and tetra-nucleosome arrays can answer important questions about chromatin folding, function, and dynamics. These multi-nucleosome arrays have highlighted the effects of varying linker DNA lengths, binding partners, and histone post-translational modifications in a more chromatin-like environment. We identify various chromatin regulatory mechanisms yet to be explored with multi-nucleosome arrays. Combined with in-solution biophysical techniques, studies of minimal multi-nucleosome chromatin models are feasible.

Compromised 2-start zigzag chromatin folding in immature mouse retina cells driven by irregularly spaced nucleosomes with short DNA linkers

The formation of condensed heterochromatin is critical for establishing cell-specific transcriptional programs. To reveal structural transitions underlying heterochromatin formation in maturing mouse rod photoreceptors, we apply cryo-EM tomography, AI-assisted deep denoising, and molecular modeling. We find that chromatin isolated from immature retina cells contains many closely apposed nucleosomes with extremely short or absent nucleosome linkers, which are inconsistent with the typical two-start zigzag chromatin folding. In mature retina cells, the fraction of short-linker nucleosomes is much lower, supporting stronger chromatin compaction. By Cryo-EM-assisted nucleosome interaction capture we observe that chromatin in immature retina is enriched with i±1 interactions while chromatin in mature retina contains predominantly i±2 interactions typical of the two-start zigzag. By mesoscale modeling and computational simulation, we clarify that the unusually short linkers typical of immature retina are sufficient to inhibit the two-start zigzag and chromatin compaction by the interference of very short linkers with linker DNA stems. We propose that this short linker composition renders nucleosome arrays more open in immature retina and that, as the linker DNA length increases in mature retina, chromatin fibers become globally condensed via tight zigzag folding. This mechanism may be broadly utilized to introduce higher chromatin folding entropy for epigenomic plasticity.

Nucleosome Spacing Can Fine-Tune Higher Order Chromatin Assembly

Cellular chromatin displays heterogeneous structure and dynamics, properties that control diverse nuclear processes. Models invoke phase separation of conformational ensembles of chromatin fibers as a mechanism regulating chromatin organization in vivo. Here we combine biochemistry and molecular dynamics simulations to examine, at single base-pair resolution, how nucleosome spacing controls chromatin phase separation. We show that as DNA linkers extend from 25 bp to 30 bp, as examplars of 10N+5 and 10N (integer N) bp lengths, chromatin condensates become less thermodynamically stable and nucleosome mobility increases. Simulations reveal that this is due to trade-offs between inter- and intramolecular nucleosome stacking, favored by rigid 10N+5 and 10N bp linkers, respectively. A remodeler can induce or inhibit phase separation by moving nucleosomes, changing the balance between intra- and intermolecular stacking. The intrinsic phase separation capacity of chromatin enables fine tuning of compaction and dynamics, likely contributing to heterogeneous chromatin organization in vivo.

Nucleosome condensate and linker DNA alter chromatin folding pathways and rates

Chromatin organization is essential for DNA packaging and gene regulation in eukaryotic genomes. While significant progresses have been made, the exact atomistic arrangement of nucleosomes remains controversial. Using a well-calibrated residue-level coarse-grained model and advanced dynamics modeling techniques, particularly the non-Markovian dynamics model, we map the free energy landscape of tetra-nucleosome systems, identify both metastable conformations and intermediate states in folding pathways, and quantify the folding kinetics. Our findings show that chromatin with 10 n base pairs (bp) DNA linker lengths favor zigzag fibril structures. However, longer linker lengths destabilize this conformation. When the linker length is 10 n + 5 bp , chromatin loses unique conformations, favoring a dynamic ensemble of structures resembling folding intermediates. Embedding the tetra-nucleosome in a nucleosome condensate similarly shifts stability towards folding intermediates as a result of the competition of inter-nucleosomal contacts. These results suggest that chromatin organization observed in vivo arises from the unfolding of fibril structures due to nucleosome crowding and linker length variation. This perspective aids in unifying experimental studies to develop atomistic models for chromatin.

Significance

Atomic structures of chromatin have become increasingly accessible, largely through cryo-EM techniques. Nonetheless, these approaches often face limitations in addressing how intrinsic in vivo factors influence chromatin organization. We present a structural characterization of chromatin under the combined effects of nucleosome condensate crowding and linker DNA length variation-two critical in vivo features that have remained challenging to capture experimentally. This work leverages a novel application of non-Markovian dynamical modeling, providing accurate mapping of chromatin folding kinetics and pathways. Our findings support a hypothesis that in vivo chromatin organization arises from folding intermediates advancing toward a stable fibril configuration, potentially resolving longstanding questions surrounding chromatin atomic structure.

Unveiling Nucleosome Dynamics: A Comparative Study Using All-Atom and Coarse-Grained Simulations Enhanced by Principal Component Analysis

This study investigates nucleosome dynamics using both all-atom and coarse-grained (CG) molecular dynamics simulations, focusing on the SIRAH force field. Simulations are performed for two nucleosomal DNA sequences-ASP and Widom-601-over six microseconds at physiological salt concentrations. Comparative analysis of structural parameters, such as groove widths and base pair geometries, reveals good agreement between atomistic and CG models, though CG simulations exhibit broader conformational sampling and greater breathing motion of DNA ends. Principal component analysis (PCA) is applied to DNA structural parameters, revealing multiple free energy minima, especially in CG simulations. These findings highlight the potential of the SIRAH CG force field for studying large-scale nucleosome dynamics, offering insights into DNA repositioning and sequence-dependent behavior.

The relationship between nucleosome positioning and higher-order genome folding

Eukaryotic genomes are organized into 3D structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These structures have an important role in the regulation of transcription and other nuclear processes. Despite advances in our understanding of the properties, functions, and underlying mechanisms of genome structures, there are many open questions about the interplay between these structures across scales. In particular, it is not well understood if and how 1D features of nucleosome arrays influence large-scale 3D genome folding patterns. In this review, we discuss recent studies that address these questions and summarize our current understanding of the relationship between nucleosome positioning and higher-order genome folding.

Determining mesoscale chromatin structure parameters from spatially correlated cleavage data using a coarse-grained oligonucleosome model

The three-dimensional structure of chromatin has emerged as an important feature of eukaryotic gene regulation. Recent technological advances in DNA sequencing-based assays have revealed locus- and chromatin state-specific structural patterns at the length scale of a few nucleosomes (~1 kb). However, interpreting these data sets remains challenging. Radiation-induced correlated cleavage of chromatin (RICC-seq) is one such chromatin structure assay that maps DNA-DNA-contacts at base pair resolution by sequencing single-stranded DNA fragments released from irradiated cells. Here, we develop a flexible modeling and simulation framework to enable the interpretation of RICC-seq data in terms of oligonucleosome structure ensembles. Nucleosomes are modeled as rigid bodies with excluded volume and adjustable DNA wrapping, connected by linker DNA modeled as a worm-like chain. We validate the model's parameters against cryo-electron microscopy and sedimentation data. Our results show that RICC-seq is sensitive to nucleosome spacing, nucleosomal DNA wrapping, and the strength of inter-nucleosome interactions. We show that nucleosome repeat lengths consistent with orthogonal assays can be extracted from experimental RICC-seq data using a 1D convolutional neural net trained on RICC-seq signal predicted from simulated ensembles. We thus provide a suite of analysis tools that add quantitative structural interpretability to RICC-seq experiments.

Hi-BDiSCO: folding 3D mesoscale genome structures from Hi-C data using brownian dynamics

The structure and dynamics of the eukaryotic genome are intimately linked to gene regulation and transcriptional activity. Many chromosome conformation capture experiments like Hi-C have been developed to detect genome-wide contact frequencies and quantify loop/compartment structures for different cellular contexts and time-dependent processes. However, a full understanding of these events requires explicit descriptions of representative chromatin and chromosome configurations. With the exponentially growing amount of data from Hi-C experiments, many methods for deriving 3D structures from contact frequency data have been developed. Yet, most reconstruction methods use polymer models with low resolution to predict overall genome structure. Here we present a Brownian Dynamics (BD) approach termed Hi-BDiSCO for producing 3D genome structures from Hi-C and Micro-C data using our mesoscale-resolution chromatin model based on the Discrete Surface Charge Optimization (DiSCO) model. Our approach integrates reconstruction with chromatin simulations at nucleosome resolution with appropriate biophysical parameters. Following a description of our protocol, we present applications to the NXN, HOXC, HOXA and Fbn2 mouse genes ranging in size from 50 to 100 kb. Such nucleosome-resolution genome structures pave the way for pursuing many biomedical applications related to the epigenomic regulation of chromatin and control of human disease.

Regulation of chromatin architecture by transcription factor binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome to control the transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers at nucleosome resolution to determine TF effect on chromatin architecture and shed light into its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. At the same time, the level of linker histone and tail acetylation, or different linker DNA lengths, regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Regulation of Chromatin Architecture by Transcription Factor Binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome, controlling transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers to determine TF effect on chromatin architecture and shed light on its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. On the other hand, the presence of linker histone, acetylations, or different linker DNA lengths regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Structure of native chromatin fibres revealed by Cryo-ET in situ

The structure of chromatin plays pivotal roles in regulating gene transcription, DNA replication and repair, and chromosome segregation. This structure, however, remains elusive. Here, using cryo-FIB and cryo-ET, we delineate the 3D architecture of native chromatin fibres in intact interphase human T-lymphoblasts and determine the in situ structures of nucleosomes in different conformations. These chromatin fibres are not structured as uniform 30 nm one-start or two-start filaments but are composed of relaxed, variable zigzag organizations of nucleosomes connected by straight linker DNA. Nucleosomes with little H1 and linker DNA density are distributed randomly without any spatial preference. This work will inspire future high-resolution investigations on native chromatin structures in situ at both a single-nucleosome level and a population level under many different cellular conditions in health and disease.

Cryo-ET study from in vitro to in vivo revealed a general folding mode of chromatin with two-start helical architecture

The organization and dynamics of chromatin fiber play crucial roles in regulating DNA accessibility for gene expression. Here we combine cryoelectron tomography (cryo-ET), sub-volume averaging, and 3D segmentation to visualize the in vitro and in vivo chromatin fibers folding by linker histone. We discover that an increased nucleosome repeat length and prolonged fiber length do not change the two-start helical architecture in reconstituted chromatin of homogeneous composition. Additionally, an isolated chromatin fiber with heterogeneous composition was observed, which includes short-range regions compatible with two-start helix. In vivo, sub-volume averaging reveals similar subunits of two-start helical architecture in transcriptionally inactive chromatin in frog erythrocyte nuclei. Strikingly, unambiguous DNA trajectories that displayed a zigzag pattern universally between alternate N/N+2 nucleosomes were further determined by cryo-ET with voltage phase plate. Therefore, these structural similarities suggest a general folding mode of chromatin induced by linker histone, and heterogeneous compositions mainly affect local conformation rather than changing the overall architecture.

Cryoelectron tomography reveals the multiplex anatomy of condensed native chromatin and its unfolding by histone citrullination

Nucleosome chains fold and self-associate to form higher-order structures whose internal organization is unknown. Here, cryoelectron tomography (cryo-ET) of native human chromatin reveals intrinsic folding motifs such as (1) non-uniform nucleosome stacking, (2) intermittent parallel and perpendicular orientations of adjacent nucleosome planes, and (3) a regressive nucleosome chain path, which deviates from the direct zigzag topology seen in reconstituted nucleosomal arrays. By examining the self-associated structures, we observed prominent nucleosome stacking in cis and anti-parallel nucleosome interactions, which are consistent with partial nucleosome interdigitation in trans. Histone citrullination strongly inhibits nucleosome stacking and self-association with a modest effect on chromatin folding, whereas the reconstituted arrays undergo a dramatic unfolding into open zigzag chains induced by histone citrullination. This study sheds light on the internal structure of compact chromatin nanoparticles and suggests a mechanism for how epigenetic changes in chromatin folding are retained across both open and condensed forms.

Heterogeneous non-canonical nucleosomes predominate in yeast cells in situ

Nuclear processes depend on the organization of chromatin, whose basic units are cylinder-shaped complexes called nucleosomes. A subset of mammalian nucleosomes in situ (inside cells) resembles the canonical structure determined in vitro 25 years ago. Nucleosome structure in situ is otherwise poorly understood. Using cryo-electron tomography (cryo-ET) and 3D classification analysis of budding yeast cells, here we find that canonical nucleosomes account for less than 10% of total nucleosomes expected in situ. In a strain in which H2A-GFP is the sole source of histone H2A, class averages that resemble canonical nucleosomes both with and without GFP densities are found ex vivo (in nuclear lysates), but not in situ. These data suggest that the budding yeast intranuclear environment favors multiple non-canonical nucleosome conformations. Using the structural observations here and the results of previous genomics and biochemical studies, we propose a model in which the average budding yeast nucleosome's DNA is partially detached in situ.

Brownian dynamics simulations of mesoscale chromatin fibers

The relationship between chromatin architecture and function defines a central problem in biology and medicine. Many computational chromatin models with atomic, coarse-grained, mesoscale, and polymer resolution have been used to shed light onto the mechanisms that dictate genome folding and regulation of gene expression. The associated simulation techniques range from Monte Carlo to molecular, Brownian, and Langevin dynamics. Here, we present an efficient Compute Unified Device Architecture (CUDA) implementation of Brownian dynamics (BD) to simulate chromatin fibers at the nucleosome resolution with our chromatin mesoscale model. With the CUDA implementation for computer architectures with graphic processing units (GPUs), we significantly accelerate compute-intensive hydrodynamic tensor calculations in the BD simulations by massive parallelization, boosting the performance a hundred-fold compared with central processing unit calculations. We validate our BD simulation approach by reproducing experimental trends on fiber diffusion and structure as a function of salt, linker histone binding, and histone-tail composition, as well as Monte Carlo equilibrium sampling results. Our approach proves to be physically accurate with performance that makes feasible the study of chromatin fibers in the range of kb or hundreds of nucleosomes (small gene). Such simulations are essential to advance the study of biological processes such as gene regulation and aberrant genome-structure related diseases.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

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.

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.

The Chromatin Regulator HMGA1a Undergoes Phase Separation in the Nucleus

The protein high mobility group A1 (HMGA1) is an important regulator of chromatin organization and function. However, the mechanisms by which it exerts its biological function are not fully understood. Here, we report that the HMGA isoform, HMGA1a, nucleates into foci that display liquid-like properties in the nucleus, and that the protein readily undergoes phase separation to form liquid condensates in vitro. By bringing together machine-leaning modelling, cellular and biophysical experiments and multiscale simulations, we demonstrate that phase separation of HMGA1a is promoted by protein-DNA interactions, and has the potential to be modulated by post-transcriptional effects such as phosphorylation. We further show that the intrinsically disordered C-terminal tail of HMGA1a significantly contributes to its phase separation through electrostatic interactions via AT hooks 2 and 3. Our work sheds light on HMGA1 phase separation as an emergent biophysical factor in regulating chromatin structure.

DNA Sequence-Dependent Properties of Nucleosome Positioning in Regions of Distinct Chromatin States in Mouse Embryonic Stem Cells

Chromatin architecture is orchestrated, and plays crucial roles during the developmental process by regulating gene expression. In embryonic stem cells (ESCs), three types of chromatin states, including active, repressive and poised states, were previously identified and characterized with specific chromatin modification marks and different transcription activity, but it is largely unknown how nucleosomes are organized in these chromatin states. In this study, by using a DNA deformation energy model, we investigated the sequence-dependent nucleosome organization within the chromatin states in mouse ESCs. The results revealed that: (1) compared with poised genes, active genes are characterized with a higher level of nucleosome occupancy around their transcription start sites (TSS) and transcription termination sites (TTS), and both types of genes do not have a nucleosome-depleted region at their TTS, contrasting with the MNase-seq based result; (2) based on our previous DNA bending energy model, we developed an improved model capable of predicting both rotational positioning and nucleosome occupancy determined by a chemical mapping approach; (3) DNA bending-energy-based analyses demonstrated that the fragile nucleosomes positioned at both gene ends could be explained largely by enhanced rotational positioning signals encoded in DNA, but nucleosome phasing around the TSS of active genes was not determined by sequence preference; (4) the nucleosome occupancy landscape around the binding sites of some developmentally important transcription factors known to bind with different chromatin contexts, was also successfully predicted; (5) the difference of nucleosome occupancy around the TSS between CpG-rich and CpG-poor promoters was partly captured by our sequence-dependent model. Taken together, by developing an improved deformation-energy-based model, we revealed some sequence-dependent properties of the nucleosome arrangements in regions of distinct chromatin states in mouse ESCs.

Single-nucleosome imaging reveals steady-state motion of interphase chromatin in living human cells

Dynamic chromatin behavior plays a critical role in various genome functions. However, it remains unclear how chromatin behavior changes during interphase, where the nucleus enlarges and genomic DNA doubles. While the previously reported chromatin movements varied during interphase when measured using a minute or longer time scale, we unveil that local chromatin motion captured by single-nucleosome imaging/tracking on a second time scale remained steady throughout G₁, S, and G₂ phases in live human cells. This motion mode appeared to change beyond this time scale. A defined genomic region also behaved similarly. Combined with Brownian dynamics modeling, our results suggest that this steady-state chromatin motion was mainly driven by thermal fluctuations. Steady-state motion temporarily increased following a DNA damage response. Our findings support the viscoelastic properties of chromatin. We propose that the observed steady-state chromatin motion allows cells to conduct housekeeping functions, such as transcription and DNA replication, under similar environments during interphase.

The use of vector formalism in the analysis of hydrophobic and electric driving forces in biological assemblies

Hydrophobic forces are known to have a crucial part not only in the conformation of the three-dimensional structure of proteins, but also in the build-up of DNA–protein complexes. Electric forces also play an important role both in the tertiary as well in the quaternary structure of macromolecular associations. Sometimes both hydrophobic and electric interactions add up their strengths to accomplish these structures but in most cases they act in opposite directions. This fact, together with being overall interactions with different ranges, provides a nuanced equilibrium also modulated by the need to comply with steric hindrances and geometric frustration effects. This review focuses on the utility of using the hydrophobic and electrical dipole moment vectors to describe the interactions that give rise to the structures of biological macromolecules. Although different definitions of both electric dipole and hydrophobic moments have been described in the literature, results obtained in biological assemblies demonstrate the principle of the biological membrane model. According to this model, postulated by our group, biological macromolecules tend to associate by aligning their hydrophobic moments in a similar manner to phospholipids in a membrane. Examples of both closed and open structures are used to assess the predictability of our model. We seek agreement between our results with those described in the current literature. The review ends with possible future projections using this formalism.

Differences in nanoscale organization of regulatory active and inactive human chromatin

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.

Local chromatin fiber folding represses transcription and loop extrusion in quiescent cells

A longstanding hypothesis is that chromatin fiber folding mediated by interactions between nearby nucleosomes represses transcription. However, it has been difficult to determine the relationship between local chromatin fiber compaction and transcription in cells. Further, global changes in fiber diameters have not been observed, even between interphase and mitotic chromosomes. We show that an increase in the range of local inter-nucleosomal contacts in quiescent yeast drives the compaction of chromatin fibers genome-wide. Unlike actively dividing cells, inter-nucleosomal interactions in quiescent cells require a basic patch in the histone H4 tail. This quiescence-specific fiber folding globally represses transcription and inhibits chromatin loop extrusion by condensin. These results reveal that global changes in chromatin fiber compaction can occur during cell state transitions, and establish physiological roles for local chromatin fiber folding in regulating transcription and chromatin domain formation.

Investigating site-selection mechanisms of retroviral integration in supercoiled DNA braids

We theoretically study the integration of short viral DNA in a DNA braid made up by two entwined double-stranded DNA molecules. We show that the statistics of single integration events substantially differ in the straight and buckled, or plectonemic, phase of the braid and are more likely in the latter. We further discover that integration is most likely close to plectoneme tips, where the larger bending energy helps overcome the associated energy barrier and that successive integration events are spatio-temporally correlated, suggesting a potential mechanistic explanation of clustered integration sites in host genomes. The braid geometry we consider provides a novel experimental set-up to quantify integration in a supercoiled substrate in vitro, and to better understand the role of double-stranded DNA topology during this process.

Protein intrinsic disorder on a dynamic nucleosomal landscape

The complex nucleoprotein landscape of the eukaryotic cell nucleus is rich in dynamic proteins that lack a stable three-dimensional structure. Many of these intrinsically disordered proteins operate directly on the first fundamental level of genome compaction: the nucleosome. Here we give an overview of how disordered interactions with and within nucleosomes shape the dynamics, architecture, and epigenetic regulation of the genetic material, controlling cellular transcription patterns. We highlight experimental and computational challenges in the study of protein disorder and illustrate how integrative approaches are increasingly unveiling the fine details of nuclear interaction networks. We finally dissect sequence properties encoded in disordered regions and assess common features of disordered nucleosome-binding proteins. As drivers of many critical biological processes, disordered proteins are integral to a comprehensive molecular view of the dynamic nuclear milieu.

Liquid-like chromatin in the cell: What can we learn from imaging and computational modeling?

Chromatin in eukaryotic cells is a negatively charged long polymer consisting of DNA, histones, and various associated proteins. With its highly charged and heterogeneous nature, chromatin structure varies greatly depending on various factors (e.g. chemical modifications and protein enrichment) and the surrounding environment (e.g. cations): from a 10-nm fiber, a folded 30-nm fiber, to chromatin condensates/droplets. Recent advanced imaging has observed that chromatin exhibits a dynamic liquid-like behavior and undergoes structural variations within the cell. Current computational modeling has made it possible to reconstruct the liquid-like chromatin in the cell by dealing with a number of nucleosomes on multiscale levels and has become a powerful technique to inspect the molecular mechanisms giving rise to the observed behavior, which imaging methods cannot do on their own. Based on new findings from both imaging and modeling studies, we discuss the dynamic aspect of chromatin in living cells and its functional relevance.

Extended and dynamic linker histone-DNA Interactions control chromatosome compaction

Chromatosomes play a fundamental role in chromatin regulation, but a detailed understanding of their structure is lacking, partially due to their complex dynamics. Using single-molecule DNA unzipping with optical tweezers, we reveal that linker histone interactions with DNA are remarkably extended, with the C-terminal domain binding both DNA linkers as far as approximately ±140 bp from the dyad. In addition to a symmetrical compaction of the nucleosome core governed by globular domain contacts at the dyad, the C-terminal domain compacts the nucleosome's entry and exit. These interactions are dynamic, exhibit rapid binding and dissociation, are sensitive to phosphorylation of a specific residue, and are crucial to determining the symmetry of the chromatosome's core. Extensive unzipping of the linker DNA, which mimics its invasion by motor proteins, shifts H1 into an asymmetric, off-dyad configuration and triggers nucleosome decompaction, highlighting the plasticity of the chromatosome structure and its potential regulatory role.

Structures and Functions of Chromatin Fibers

In eukaryotes, genomic DNA is packaged into chromatin in the nucleus. The accessibility of DNA is dependent on the chromatin structure and dynamics, which essentially control DNA-related processes, including transcription, DNA replication, and repair. All of the factors that affect the structure and dynamics of nucleosomes, the nucleosome-nucleosome interaction interfaces, and the binding of linker histones or other chromatin-binding proteins need to be considered to understand the organization and function of chromatin fibers. In this review, we provide a summary of recent progress on the structure of chromatin fibers in vitro and in the nucleus, highlight studies on the dynamic regulation of chromatin fibers, and discuss their related biological functions and abnormal organization in diseases.

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.

Physical Nature of Chromatin in the Nucleus

Genomic information is encoded on long strands of DNA, which are folded into chromatin and stored in a tiny nucleus. Nuclear chromatin is a negatively charged polymer composed of DNA, histones, and various nonhistone proteins. Because of its highly charged nature, chromatin structure varies greatly depending on the surrounding environment (e.g., cations, molecular crowding, etc.). New technologies to capture chromatin in living cells have been developed over the past 10 years. Our view on chromatin organization has drastically shifted from a regular and static one to a more variable and dynamic one. Chromatin forms numerous compact dynamic domains that act as functional units of the genome in higher eukaryotic cells and locally appear liquid-like. By changing DNA accessibility, these domains can govern various functions. Based on new evidences from versatile genomics and advanced imaging studies, we discuss the physical nature of chromatin in the crowded nuclear environment and how it is regulated.

The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome Array

Linker histones bind to nucleosomes and modify chromatin structure and dynamics as a means of epigenetic regulation. Biophysical studies have shown that chromatin fibers can adopt a plethora of conformations with varying levels of compaction. Linker histone condensation, and its specific binding disposition, has been associated with directly tuning this ensemble of states. However, the atomistic dynamics and quantification of this mechanism remains poorly understood. Here, we present molecular dynamics simulations of octa-nucleosome arrays, based on a cryo-EM structure of the 30-nm chromatin fiber, with and without the globular domains of the H1 linker histone to determine how they influence fiber structures and dynamics. Results show that when bound, linker histones inhibit DNA flexibility and stabilize repeating tetra-nucleosomal units, giving rise to increased chromatin compaction. Furthermore, upon the removal of H1, there is a significant destabilization of this compact structure as the fiber adopts less strained and untwisted states. Interestingly, linker DNA sampling in the octa-nucleosome is exaggerated compared to its mono-nucleosome counterparts, suggesting that chromatin architecture plays a significant role in DNA strain even in the absence of linker histones. Moreover, H1-bound states are shown to have increased stiffness within tetra-nucleosomes, but not between them. This increased stiffness leads to stronger long-range correlations within the fiber, which may result in the propagation of epigenetic signals over longer spatial ranges. These simulations highlight the effects of linker histone binding on the internal dynamics and global structure of poly-nucleosome arrays, while providing physical insight into a mechanism of chromatin compaction.

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.

Mesoscale Modeling and Single-Nucleosome Tracking Reveal Remodeling of Clutch Folding and Dynamics in Stem Cell Differentiation

Nucleosomes form heterogeneous groups in vivo, named clutches. Clutches are smaller and less dense in mouse embryonic stem cells (ESCs) compared to neural progenitor cells (NPCs). Using coarse-grained modeling of the pluripotency Pou5f1 gene, we show that the genome-wide clutch differences between ESCs and NPCs can be reproduced at a single gene locus. Larger clutch formation in NPCs is associated with changes in the compaction and internucleosome contact probability of the Pou5f1 fiber. Using single-molecule tracking (SMT), we further show that the core histone protein H2B is dynamic, and its local mobility relates to the structural features of the chromatin fiber. H2B is less stable and explores larger areas in ESCs compared to NPCs. The amount of linker histone H1 critically affects local H2B dynamics. Our results have important implications for how nucleosome organization and H2B dynamics contribute to regulate gene activity and cell identity.

Linker DNA Length is a Key to Tri-nucleosome Folding

The folding of a nucleosome array has long been one of the fundamental and unsolved problems in chromatin biology. In this study, we address how nucleosome array folding depends on the length of linker DNA. We performed molecular dynamics simulations of a tri-nucleosome, a minimal model of chromatin folding, with various linker lengths (LLs) ranging from 20 to 40 base pairs (bps). We found that the tri-nucleosome folding strongly depends on LLs, and classified the structure ensemble into five classes, named from trinuc-1 to trinuc-5. As a function of LL, the different classes appear, on average, every 2 bps with a period of 10 bps, and are characterized by distinct inter-nucleosome interactions. The trinuc-1 conformation corresponds to LL ~ 10n, where n is an integer, and is stabilized by the tight packing between the first and the third nucleosomes, consistent with a zigzag fiber form. Structures of the other four classes are more diverse and distributed continuously in the space of possible configurations. Histone-DNA electrostatic interactions in the tri-nucleosome are further analyzed.

Breaths, Twists, and Turns of Atomistic Nucleosomes

Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145-147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.

Significant compaction of H4 histone tail upon charge neutralization by acetylation and its mimics, possible effects on chromatin structure

The intrinsically disordered, positively charged H4 histone tail is important for chromatin structure and function. We have explored conformational ensembles of human H4 tail in solution, with varying levels of charge neutralization via acetylation or amino-acid substitutions such as K→Q. We have employed an explicit water model shown recently to be well suited for simulations of intrinsically disordered proteins. Upon progressive neutralization of the H4, its radius of gyration decreases linearly with the tail charge q, the trend is explained using a simple polymer model. While the wild type state (q=+8) is essentially a random coil, hyper-acetylated H4 (q=+3) is virtually as compact and stable as a globular protein of the same number of amino-acids. Conformational ensembles of acetylated H4 match the corresponding K→X substitutions only approximately: based on the ensemble similarity, we propose K→M as a possible alternative to the commonly used K→Q. Possible effects of the H4 tail compaction on chromatin structure are discussed within a qualitative model in which the chromatin is highly heterogeneous, easily inter-converting between various structural forms. We predict that upon progressive charge neutralization of the H4 tail, the least compact sub-states of chromatin de-condense first, followed by de-condensation of more compact structures, e.g. those that harbor a high fraction of stacked di-nucleosomes. The predicted hierarchy of DNA accessibility increase upon progressive acetylation of H4 might be utilized by the cell for selective DNA accessibility control.

The role of histone tails in nucleosome stability: An electrostatic perspective

We propose a methodology for the study of protein-DNA electrostatic interactions and apply it to clarify the effect of histone tails in nucleosomes. This method can be used to correlate electrostatic interactions to structural and functional features of protein-DNA systems, and can be combined with coarse-grained representations. In particular, we focus on the electrostatic field and resulting forces acting on the DNA. We investigate the electrostatic origins of effects such as different stages in DNA unwrapping, nucleosome destabilization upon histone tail truncation, and the role of specific arginines and lysines undergoing Post-Translational Modifications. We find that the positioning of the histone tails can oppose the attractive pull of the histone core, locally deform the DNA, and tune DNA unwrapping. Small conformational variations in the often overlooked H2A C-terminal tails had significant electrostatic repercussions near the DNA entry and exit sites. The H2A N-terminal tail exerts attractive electrostatic forces towards the histone core in positions where Polymerase II halts its progress. We validate our results with comparisons to previous experimental and computational observations.

The asynchronous establishment of chromatin 3D architecture between in vitro fertilized and uniparental preimplantation pig embryos

BACKGROUND

Pigs are important animals for agricultural and biomedical research, and improvement is needed for use of the assisted reproductive technologies. Determining underlying mechanisms of epigenetic reprogramming in the early stage of preimplantation embryos derived from in vitro fertilization (IVF), parthenogenesis, and androgenesis will not only contribute to assisted reproductive technologies of pigs but also will shed light into early human development. However, the reprogramming of three-dimensional architecture of chromatin in this process in pigs is poorly understood.

RESULTS

We generate three-dimensional chromatin profiles for pig somatic cells, IVF, parthenogenesis, and androgenesis preimplantation embryos. We find that the chromosomes in the pig preimplantation embryos are enriched for superdomains, which are more rare in mice. However, p(s) curves, compartments, and topologically associated domains (TADs) are largely conserved in somatic cells and are gradually established during preimplantation embryogenesis in both mammals. In the uniparental pig embryos, the establishment of chromatin architecture is highly asynchronized at all levels from IVF embryos, and a remarkably strong decompartmentalization is observed during zygotic genome activation (ZGA). Finally, chromosomes originating from oocytes always establish TADs faster than chromosomes originating from sperm, both before and during ZGA.

CONCLUSIONS

Our data highlight a potential unique 3D chromatin pattern of enriched superdomains in pig preimplantation embryos, an unusual decompartmentalization process during ZGA in the uniparental embryos, and an asynchronized TAD reprogramming between maternal and paternal genomes, implying a severe dysregulation of ZGA in the uniparental embryos in pigs.

Nucleosome positioning and chromatin organization

In our cells, DNA is folded and packed with the help of many proteins into chromatin whose basic unit is a nucleosome-DNA wrapped around octamer of histone proteins. The chain of nucleosomes is further folded and arranged into many layers and has a dynamic organization. How does the complex chromatin organization emerge from interactions among DNA, histones, and non-histone proteins have been a question of great interest. Here we review recent literature that investigated how nucleosome positioning and nucleosome-mediated interactions drive chromatin organization. Unlike our earlier understanding, chromatin is organized into 3D domains of various sizes having irregularly organized nucleosomes. These domains emerge due to heterogeneous nucleosome positioning and diverse inter-nucleosome interactions that vary in space and time.

Cryo-electron microscopy of the chromatin fiber

The three-dimensional (3D) organization of chromatin plays a crucial role in the regulation of gene expression. Chromatin conformation is strongly affected by the composition, structural features and dynamic properties of the nucleosome, which in turn determine the nature and geometry of interactions that can occur between neighboring nucleosomes. Understanding how chromatin is spatially organized above the nucleosome level is thus essential for understanding how gene regulation is achieved. Towards this end, great effort has been made to understand how an array of nucleosomes folds into a regular chromatin fiber. This review summarizes new insights into the 3D structure of the chromatin fiber that were made possible by recent advances in cryo-electron microscopy.

Mesoscale Modeling of Nucleosome-Binding Antibody PL2-6: Mono- versus Bivalent Chromatin Complexes

Interactions of chromatin with bivalent immunoglobin nucleosome-binding antibodies and their monovalent (papain-derived) antigen-binding fragment analogs are useful probes for examining chromatin conformational states. To help interpret antibody-chromatin interactions and explore how antibodies might compete for interactions with chromatin components, we incorporate coarse-grained PL2-6 antibody modeling into our mesoscale chromatin model. We analyze interactions and fiber structures for the antibody-chromatin complexes in open and condensed chromatin, with and without H1 linker histone (LH). Despite minimal and transient interactions at physiological salt, we capture significant differences in antibody-chromatin complex configurations in open fibers, with more intense interactions between the bivalent antibody and chromatin compared to monovalent antigen-binding fragments. For these open chromatin fiber morphologies, antibody binding to histone tails is increased and compaction is greater for bivalent compared to monovalent and antibody-free systems. Differences between monovalent and bivalent binding result from antibody competition with internal chromatin fiber components (nucleosome core and linker DNA) for histone tail (H3, H4, H2A, H2B) interactions. This antibody competition for tail contacts reduces tail-core and tail-linker interactions and increases tail-antibody interactions. Such internal structural changes in open fibers resemble mechanisms of LH condensation, driven by charge screening and entropy changes. For condensed fibers at physiological salt, the three systems are much more similar overall, but some subtle tail interaction differences can be noted. Adding LH results in less-dramatic changes for all systems, except that the bivalent complex at physiological salt shows cooperative effects between LH and the antibodies in condensing chromatin fibers. Such dynamic interactions that depend on the internal structure and complex-stabilizing interactions within the chromatin fiber have implications for gene regulation and other chromatin complexes such as with LH, remodeling proteins, and small molecular chaperones that bind and modulate chromatin structure.

Fluid-like chromatin: Toward understanding the real chromatin organization present in the cell

Eukaryotic chromatin is a negatively charged polymer consisting of genomic DNA, histones, and various nonhistone proteins. Because of its highly charged character, the structure of chromatin varies greatly depending on the surrounding environment (i.e. cations etc.): from an extended 10-nm fiber, to a folded 30-nm fiber, to chromatin condensates/liquid-droplets. Over the last ten years, newly developed technologies have drastically shifted our view on chromatin from a static regular structure to a more irregular and dynamic one, locally like a fluid. Since no single imaging (or genomics) method can tell us everything and beautiful images (or models) can fool our minds, comprehensive analyses based on many technical approaches are important to capture actual chromatin organization inside the cell. Here we critically discuss our current view on chromatin and methodology used to support the view.

Emergence of chromatin hierarchical loops from protein disorder and nucleosome asymmetry

Protein flexibility and disorder is emerging as a crucial modulator of chromatin structure. Histone tail disorder enables transient binding of different molecules to the nucleosomes, thereby promoting heterogeneous and dynamic internucleosome interactions and making possible recruitment of a wide-range of regulatory and remodeling proteins. On the basis of extensive multiscale modeling we reveal the importance of linker histone H1 protein disorder for chromatin hierarchical looping. Our multiscale approach bridges microsecond-long bias-exchange metadynamics molecular dynamics simulations of atomistic 211-bp nucleosomes with coarse-grained Monte Carlo simulations of 100-nucleosome systems. We show that the long C-terminal domain (CTD) of H1-a ubiquitous nucleosome-binding protein-remains disordered when bound to the nucleosome. Notably, such CTD disorder leads to an asymmetric and dynamical nucleosome conformation that promotes chromatin structural flexibility and establishes long-range hierarchical loops. Furthermore, the degree of condensation and flexibility of H1 can be fine-tuned, explaining chromosomal differences of interphase versus metaphase states that correspond to partial and hyperphosphorylated H1, respectively. This important role of H1 protein disorder in large-scale chromatin organization has a wide range of biological implications.