Chromatin fiber polymorphism triggered by variations of DNA linker lengths

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.

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 ¹H-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 ²H,¹³C,¹⁵N-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.

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.

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome.

Chromatosome Structure and Dynamics from Molecular Simulations

Chromatosomes are fundamental units of chromatin structure that are formed when a linker histone protein binds to a nucleosome. The positioning of the linker histone on the nucleosome influences the packing of chromatin. Recent simulations and experiments have shown that chromatosomes adopt an ensemble of structures that differ in the geometry of the linker histone-nucleosome interaction. In this article we review the application of Brownian, Monte Carlo, and molecular dynamics simulations to predict the structure of linker histone-nucleosome complexes, to study the binding mechanisms involved, and to predict how this binding affects chromatin fiber structure. These simulations have revealed the sensitivityof the chromatosome structure to variations in DNA and linker histone sequence, as well as to posttranslational modifications, thereby explaining the structural variability observed in experiments. We propose that a concerted application of experimental and computational approaches will reveal the determinants of chromatosome structural variability and how it impacts chromatin packing.

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.

Geometrical Heterogeneity Dominates Thermal Fluctuations in Facilitating Chromatin Contacts

Within a living cell, the myriad of proteins that bind DNA introduce heterogeneously spaced kinks into an otherwise semiflexible DNA double helix. To investigate the effects of heterogeneous nucleosome binding on chromatin organization, we extend the wormlike chain model to include statistically spaced, rigid kinks. On timescales where nucleosome positions are fixed, we find that the probability of chromatin loop formation can vary by up to six orders of magnitude between two sets of nucleosome positions drawn from the same distribution. On longer timescales, we show that continuous rerandomization due to nucleosome turnover results in chromatin tracing out an effective WLC with a dramatically smaller Kuhn length than bare DNA. Together, these observations demonstrate that nucleosome spacing acts as the primary source of the structural heterogeneity that dominates local and global chromatin organization.

Irregular Chromatin: Packing Density, Fiber Width, and Occurrence of Heterogeneous Clusters

How chromatin is folded on the length scale of a gene is an open question. Recent experiments have suggested that, in vivo, chromatin is folded in an irregular manner and not as an ordered fiber with a width of 30 nm that is expected from theories of higher order packaging. Using computational methods, we examine how the interplay between DNA-bending nonhistone proteins, histone tails, intrachromatin electrostatic, and other interactions decide the nature of the packaging of chromatin. We show that although the DNA-bending nonhistone proteins make the chromatin irregular, they may not alter the packing density and size of the fiber. We find that the length of the interacting region and intrachromatin electrostatic interactions influence the packing density, clustering of nucleosomes, and the width of the chromatin fiber. Our results suggest that the heterogeneity in the interaction pattern will play an important role in deciding the nature of the packaging of chromatin.

Bridging chromatin structure and function over a range of experimental spatial and temporal scales by molecular modeling

Chromatin structure, dynamics, and function are being intensely investigated by a variety of methods, including microscopy, X-ray diffraction, nuclear magnetic resonance, biochemical crosslinking, chromosome conformation capture, and computation. A range of experimental techniques combined with modeling is clearly valuable to help interpret experimental data and, importantly, generate configurations and mechanisms related to the 3D organization and function of the genome. Contact maps, in particular, as obtained by a variety of chromosome conformation capture methods, are of increasing interest due to their implications on genome structure and regulation on many levels. In this perspective, using seven examples from our group's studies, we illustrate how molecular modeling can help interpret such experimental data. Specifically, we show how computed contact maps related to experimental systems can be used to explain structures of nucleosomes, chromatin higher-order folding, domain segregation mechanisms, gene organization, and the effect on chromatin structure of external and internal fiber parameters, such as nucleosome positioning, presence of nucleosome free regions, histone posttranslational modifications, and linker histone binding. We argue that such computations on multiple spatial and temporal scales will be increasingly important for the integration of genomic, epigenomic, and biophysical data on chromatin structure and related cellular processes.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.