Evidence for heteromorphic chromatin fibers from analysis of nucleosome interactions

All-atom MD simulations of the HPV16 capsid

Human papillomavirus (HPV) type 16, one of over 200 known genotypes, is a major cause of anogenital and oropharyngeal cancers worldwide. The L1 protein capsid - which represents the exterior shell of the virus - plays essential roles in host cell adhesion and entry, and is the major antigen used in vaccines that protect against HPV. Here we report an all-atom molecular dynamics (MD) simulation of the intact capsid in explicit solvent, a calculation comprising 16.8 million atoms. The detailed motions and emergent biophysical properties revealed by this simulation provide critical insights into capsid function and immunogenicity, establishing a new platform for rational drug design and antibody engineering.

Compromised two-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-electron microscopy (cryo-EM) tomography, AI-assisted 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 becomes globally condensed via tight zigzag folding. This mechanism may be broadly utilized to introduce higher chromatin folding entropy for epigenomic plasticity.

Regulation of Genome Architecture in Huntington's Disease

Huntington's disease (HD) is a neurological condition caused by an excessive expansion of CAG repeats in the Huntingtin (HTT) gene. Although experiments have shown an altered epigenetic landscape and chromatin architecture upon HD development, the structural consequences on the HTT gene remain elusive. Structural data are only available for model nucleosome systems and yeast systems with human nucleosomes. Here, we use our experimentally validated nucleosome-resolution mesoscale chromatin model to investigate folding changes of the HTT gene associated with HD. We investigate how the histone fold domain of the variant macroH2A1, a biomarker of HD, affects the genome structure by modeling HD-like systems that contain (i) 100% canonical, (ii) 100% macroH2A1, (iii) 50% canonical and 50% macroH2A1, and (iv) 100% hybrid cores (one canonical H2A and one macroH2A1 per nucleosome). Then, we model the mouse HTT gene in healthy and HD conditions by incorporating the CAG expansion and macroH2A1 cores, reducing the linker histone density and tail acetylation levels, and incorporating genomic contacts. Overall, our results show that the histone fold domain of macroH2A1 affects chromatin compaction in a fiber-dependent manner (i.e., nucleosome distribution dependent) and can thus both enhance or repress HTT gene expression. Our modeling of the HTT gene shows that HTT is less compact in the diseased condition, which could accelerate the production of the mutated protein. By suggesting the structural biophysical consequences of the HTT gene under HD conditions, our findings may help in the development of diagnostic and therapeutic treatments for HD.

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.

Selective signal enhancement in Fourier space as a tool for discovering ultrastructural organization of macromolecules from in situ TEM

We present a Fourier transform (FT) based analytical method that allows to obtain of ultrastructural details from TEM images at sub-nanometer scale applying a selective filtering for singular macromolecule electron microscopy density information. It can be applied to high-pressure frozen, frozen hydrated and epoxy freeze substituted and embedded biological species. Both 2D projections and orthoslices from reconstructed tomograms can be used as a source of structural information. The key to the method is to select the macromolecule or organelle of interest with an accuracy of ≥ 7 - 3 nm (depending on pixel size of initial tilt series or singular image acquisition) and explore both the central low frequency FT intensity and diffraction regions to obtain the spatial structural organization and its dimensional characteristics, respectively. We also introduce a structure-specific selective mask FT filtering approach that can significantly improve image information even in poorly contrasted TEM of resin sections without heavy metal been used. The described method elucidates chromatin architecture without the need of averaging. A zigzag symmetry of 30 nm diameter chromatin fibers which in general is a controversial topic of research has been identified for C. elegans cells in vivo with sub-nanometer details being preserved in the images.

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.

Regulation of chromatin architecture by protein binding: insights from molecular modeling

Histone and non-histone proteins play key roles in the activation and repression of genes. In addition to experimental studies of their regulation of gene expression, molecular modeling at the nucleosome, chromatin, and chromosome levels can contribute insights into the molecular mechanisms involved. In this review, we provide an overview for protein-bound chromatin modeling, and describe how our group has integrated protein binding into genome systems across the scales, from all-atom to coarse-grained models, using explicit to implicit descriptions. We describe the associated applications to protein binding effects and biological mechanisms of genome folding and gene regulation. We end by illustrating the application of machine learning tools like AlphaFold2 to proteins relevant to chromatin systems.

Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences

Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.

Decoding silkworm spinning programmed by pH and metal ions

Silk is one of the toughest fibrous materials known despite spun at ambient temperature and pressure with water as a solvent. It is a great challenge to reproduce high-performance artificial fibers comparable to natural silk by bionic for the incomplete understanding of silkworm spinning in vivo. Here, we found that amphipol and digitonin stabilized the structure of natural silk fibroin (NSF) by a large-scale screening in vitro, and then studied the close-to-native ultrastructure and hierarchical assembly of NSF in the silk gland lumen. Our study showed that NSF formed reversible flexible nanofibrils mainly composed of random coils with a sedimentation coefficient of 5.8 S and a diameter of about 4 nm, rather than a micellar or rod-like structure assembled by the aggregation of globular NSF molecules. Metal ions were required for NSF nanofibril formation. The successive pH decrease from posterior silk gland (PSG) to anterior silk gland (ASG) resulted in a gradual increase in NSF hydrophobicity, thus inducing the sol-gelation transition of NSF nanofibrils. NSF nanofibrils were randomly dispersed from PSG to ASG-1, and self-assembled into anisotropic herringbone patterns at ASG-2 near the spinneret ready for silkworm spinning. Our findings reveal the controlled self-assembly mechanism of the multi-scale hierarchical architecture of NSF from nanofibrils to herringbone patterns programmed by metal ions and pH gradient, which provides novel insights into the spinning mechanism of silk-secreting animals and bioinspired design of high-performance fibers.

Prediction and Control in DNA Nanotechnology

DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.

Explicit ion modeling predicts physicochemical interactions for chromatin organization

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono- and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

The role of cryptic ancestral symmetry in histone folding mechanisms across Eukarya and Archaea

Histones compact and store DNA in both Eukarya and Archaea, forming heterodimers in Eukarya and homodimers in Archaea. Despite this, the folding mechanism of histones across species remains unclear. Our study addresses this gap by investigating 11 types of histone and histone-like proteins across humans, Drosophila, and Archaea through multiscale molecular dynamics (MD) simulations, complemented by NMR and circular dichroism experiments. We confirm and elaborate on the widely applied "folding upon binding" mechanism of histone dimeric proteins and report a new alternative conformation, namely, the inverted non-native dimer, which may be a thermodynamically metastable configuration. Protein sequence analysis indicated that the inverted conformation arises from the hidden ancestral head-tail sequence symmetry underlying all histone proteins, which is congruent with the previously proposed histone evolution hypotheses. Finally, to explore the potential formations of homodimers in Eukarya, we utilized MD-based AWSEM and AI-based AlphaFold-Multimer models to predict their structures and conducted extensive all-atom MD simulations to examine their respective structural stabilities. Our results suggest that eukaryotic histones may also form stable homodimers, whereas their disordered tails bring significant structural asymmetry and tip the balance towards the formation of commonly observed heterotypic dimers.

Explicit Ion Modeling Predicts Physicochemical Interactions for Chromatin Organization

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

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.

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.

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.

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.

Nanopore Sensing of DNA-Histone Complexes on Nucleosome Arrays

The location of nucleosomes in DNA and their structural stability are critical in regulating DNA compaction, site accessibility, and epigenetic gene regulation. Here, we combine the nanopore platform-based fast and label-free single-molecule detection technique with a voltage-dependent force rupture assay to detect distinct structures on nucleosomal arrays and then to induce breakdown of individual nucleosome complexes. Specifically, we demonstrate direct measurement of distinct nucleosome structures present on individual 12-mer arrays. A detailed event analysis showed that nucleosomes are present as a combination of complete and partial structures, during translocation through the pore. By comparing with the voltage-dependent translocation of the mononucleosomes, we find that the partial nucleosomes result from voltage-dependent structural disintegration of nucleosomes. High signal-to-noise detection of heterogeneous levels in translocation of 12-mer array molecules quantifies the heterogeneity and nucleosomal substructure sizes on the arrays. These results facilitate the understanding of electrostatic interactions responsible for the integrity of the nucleosome structure and possible mechanisms of its unraveling by chromatin remodeling enzymes. This study also has potential applications in chromatin profiling.

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

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

Molecular organization of the early stages of nucleosome phase separation visualized by cryo-electron tomography

It has been proposed that the intrinsic property of nucleosome arrays to undergo liquid-liquid phase separation (LLPS) in vitro is responsible for chromatin domain organization in vivo. However, understanding nucleosomal LLPS has been hindered by the challenge to characterize the structure of the resulting heterogeneous condensates. We used cryo-electron tomography and deep-learning-based 3D reconstruction/segmentation to determine the molecular organization of condensates at various stages of LLPS. We show that nucleosomal LLPS involves a two-step process: a spinodal decomposition process yielding irregular condensates, followed by their unfavorable conversion into more compact, spherical nuclei that grow into larger spherical aggregates through accretion of spinodal materials or by fusion with other spherical condensates. Histone H1 catalyzes more than 10-fold the spinodal-to-spherical conversion. We propose that this transition involves exposure of nucleosome hydrophobic surfaces causing modified inter-nucleosome interactions. These results suggest a physical mechanism by which chromatin may transition from interphase to metaphase structures.

A simulation model of heterochromatin formation at submolecular detail

Heterochromatin is a physical state of the chromatin fiber that maintains gene repression during cell development. Although evidence exists on molecular mechanisms involved in heterochromatin formation, a detailed structural mechanism of heterochromatin formation needs a better understanding. We made use of a simple Monte Carlo simulation model with explicit representation of key molecular events to observe molecular self-organization leading to heterochromatin formation. Our simulations provide a structural interpretation of several important traits of the heterochromatinization process. In particular, this study provides a depiction of how small amounts of HP1 are able to induce a highly condensed chromatin state through HP1 dimerization and bridging of sequence-remote nucleosomes. It also elucidates structural roots of a yet poorly understood phenomenon of a nondeterministic nature of heterochromatin formation and subsequent gene repression. Experimental chromatin in vivo assay provides an unbiased estimate of time scale of repressive response to a heterochromatin-triggering event.

Histone Modifications, Internucleosome Dynamics, and DNA Stresses: How They Cooperate to “Functionalize” Nucleosomes

Tight packaging of DNA in chromatin severely constrains DNA accessibility and dynamics. In contrast, nucleosomes in active chromatin state are highly flexible, can exchange their histones, and are virtually “transparent” to RNA polymerases, which transcribe through gene bodies at rates comparable to that of naked DNA. Defining mechanisms that revert nucleosome repression, in addition to their value for basic science, is of key importance for the diagnosis and treatment of genetic diseases. Chromatin activity is largely regulated by histone posttranslational modifications, ranging from small chemical groups up to the yet understudied “bulky” ubiquitylation and sumoylation. However, it is to be revealed how histone marks are “translated” to permissive or repressive changes in nucleosomes: it is a general opinion that histone modifications act primarily as “signals” for recruiting the regulatory proteins or as a “neutralizer” of electrostatic shielding of histone tails. Here, we would like to discuss recent evidence suggesting that histone ubiquitylation, in a DNA stress–dependent manner, can directly regulate the dynamics of the nucleosome and their primary structure and can promote nucleosome decomposition to hexasome particles or additionally stabilize nucleosomes against unwrapping. In addition, nucleosome repression/ derepression studies are usually performed with single mononucleosomes as a model. We would like to review and discuss recent findings showing that internucleosomal interactions could strongly modulate the dynamics and rearrangements of nucleosomes. Our hypothesis is that bulky histone modifications, nucleosome inherent dynamics, internucleosome interactions, and DNA torsions could act in cooperation to orchestrate the formation of different dynamic states of arrayed nucleosomes and thus promote chromatin functionality and diversify epigenetic programming methods.

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.

Structural features of nucleosomes in interphase and metaphase chromosomes

Structural heterogeneity of nucleosomes in functional chromosomes is unknown. Here, we devise the template-, reference- and selection-free (TRSF) cryo-EM pipeline to simultaneously reconstruct cryo-EM structures of protein complexes from interphase or metaphase chromosomes. The reconstructed interphase and metaphase nucleosome structures are on average indistinguishable from canonical nucleosome structures, despite DNA sequence heterogeneity, cell-cycle-specific posttranslational modifications, and interacting proteins. Nucleosome structures determined by a decoy-classifying method and structure variability analyses reveal the nucleosome structural variations in linker DNA, histone tails, and nucleosome core particle configurations, suggesting that the opening of linker DNA, which is correlated with H2A C-terminal tail positioning, is suppressed in chromosomes. High-resolution (3.4-3.5 Å) nucleosome structures indicate DNA-sequence-independent stabilization of superhelical locations ±0-1 and ±3.5-4.5. The linker histone H1.8 preferentially binds to metaphase chromatin, from which chromatosome cryo-EM structures with H1.8 at the on-dyad position are reconstituted. This study presents the structural characteristics of nucleosomes in chromosomes.

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.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

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.

A critical role for linker DNA in higher-order folding of chromatin fibers

Nucleosome-nucleosome interactions drive the folding of nucleosomal arrays into dense chromatin fibers. A better physical account of the folding of chromatin fibers is necessary to understand the role of chromatin in regulating DNA transactions. Here, we studied the unfolding pathway of regular chromatin fibers as a function of single base pair increments in linker length, using both rigid base-pair Monte Carlo simulations and single-molecule force spectroscopy. Both computational and experimental results reveal a periodic variation of the folding energies due to the limited flexibility of the linker DNA. We show that twist is more restrictive for nucleosome stacking than bend, and find the most stable stacking interactions for linker lengths of multiples of 10 bp. We analyzed nucleosomes stacking in both 1- and 2-start topologies and show that stacking preferences are determined by the length of the linker DNA. Moreover, we present evidence that the sequence of the linker DNA also modulates nucleosome stacking and that the effect of the deletion of the H4 tail depends on the linker length. Importantly, these results imply that nucleosome positioning in vivo not only affects the phasing of nucleosomes relative to DNA but also directs the higher-order structure of chromatin.

Chromosome inner structure investigation by electron tomography and electron diffraction in a transmission electron microscope

Our understanding of the inner structure of metaphase chromosomes remains inconclusive despite intensive studies using multiple imaging techniques. Transmission electron microscopy has been extensively used to visualize chromosome ultrastructure. This review summarizes recent results obtained using two transmission electron microscopy-based techniques: electron tomography and electron diffraction. Electron tomography allows advanced three-dimensional imaging of chromosomes, while electron diffraction detects the presence of periodic structures within chromosomes. The combination of these two techniques provides results contributing to the understanding of local structural organization of chromatin fibers within chromosomes.

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.

Nucleosome Clutches are Regulated by Chromatin Internal Parameters

Nucleosomes cluster together when chromatin folds in the cell to form heterogeneous groups termed "clutches". These structural units add another level of chromatin regulation, for example during cell differentiation. Yet, the mechanisms that regulate their size and compaction remain obscure. Here, using our chromatin mesoscale model, we dissect clutch patterns in fibers with different combinations of nucleosome positions, linker histone density, and acetylation levels to investigate their role in clutch regulation. First, we isolate the effect of each chromatin parameter by studying systems with regular nucleosome spacing; second, we design systems with naturally-occurring linker lengths that fold onto specific clutch patterns; third, we model gene-encoding fibers to understand how these combined factors contribute to gene structure. Our results show how these chromatin parameters act together to produce different-sized nucleosome clutches. The length of nucleosome free regions (NFRs) profoundly affects clutch size, while the length of linker DNA has a moderate effect. In general, higher linker histone densities produce larger clutches by a chromatin compaction mechanism, while higher acetylation levels produce smaller clutches by a chromatin unfolding mechanism. We also show that it is possible to design fibers with naturally-occurring DNA linkers and NFRs that fold onto specific clutch patterns. Finally, in gene-encoding systems, a complex combination of variables dictates a gene-specific clutch pattern. Together, these results shed light into the mechanisms that regulate nucleosome clutches and suggest a new epigenetic mechanism by which chromatin parameters regulate transcriptional activity via the three-dimensional folded state of the genome at a nucleosome level.

Computational approaches for inferring 3D conformations of chromatin from chromosome conformation capture data

Chromosome conformation capture (3C) and its variants are powerful experimental techniques for probing intra- and inter-chromosomal interactions within cell nuclei at high resolution and in a high-throughput, quantitative manner. The contact maps derived from such experiments provide an avenue for inferring the 3D spatial organization of the genome. This review provides an overview of the various computational methods developed in the past decade for addressing the very important but challenging problem of deducing the detailed 3D structure or structure population of chromosomal domains, chromosomes, and even entire genomes from 3C contact maps.

Biology and Physics of Heterochromatin-Like Domains/Complexes

The hallmarks of constitutive heterochromatin, HP1 and H3K9me2/3, assemble heterochromatin-like domains/complexes outside canonical constitutively heterochromatic territories where they regulate chromatin template-dependent processes. Domains are more than 100 kb in size; complexes less than 100 kb. They are present in the genomes of organisms ranging from fission yeast to human, with an expansion in size and number in mammals. Some of the likely functions of domains/complexes include silencing of the donor mating type region in fission yeast, preservation of DNA methylation at imprinted germline differentially methylated regions (gDMRs) and regulation of the phylotypic progression during vertebrate development. Far cis- and trans-contacts between micro-phase separated domains/complexes in mammalian nuclei contribute to the emergence of epigenetic compartmental domains (ECDs) detected in Hi-C maps. A thermodynamic description of micro-phase separation of heterochromatin-like domains/complexes may require a gestalt shift away from the monomer as the "unit of incompatibility" that determines the sign and magnitude of the Flory-Huggins parameter, χ. Instead, a more dynamic structure, the oligo-nucleosomal "clutch", consisting of between 2 and 10 nucleosomes is both the long sought-after secondary structure of chromatin and its unit of incompatibility. Based on this assumption we present a simple theoretical framework that enables an estimation of χ for domains/complexes flanked by euchromatin and thereby an indication of their tendency to phase separate. The degree of phase separation is specified by χN, where N is the number of "clutches" in a domain/complex. Our approach could provide an additional tool for understanding the biophysics of the 3D genome.

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.

Mixed secondary chromatin structure revealed by modeling radiation-induced DNA fragment length distribution

Spatial chromatin structure plays fundamental roles in many vital biological processes including DNA replication, transcription, damage and repair. However, the current understanding of the secondary structure of chromatin formed by local nucleosome-nucleosome interactions remains controversial, especially for the existence and conformation of 30 nm structure. Since chromatin structure influences the fragment length distribution (FLD) of ionizing radiation-induced DNA strand breaks, a 3D chromatin model fitting FLD patterns can help to distinguish different models of chromatin structure. Here, we developed a novel "30-C" model combining 30 nm chromatin structure models with Hi-C data, which measured the spatial contact frequency between different loci in the genome. We first reconstructed the 3D coordinates of the 25 kb bins from Hi-C heatmaps. Within the 25 kb bins, lower level chromatin structures supported by recent studies were filled. Simulated FLD patterns based on the 30-C model were compared to published FLD patterns induced by heavy ion radiation to validate the models. Importantly, the 30-C model predicted that the most probable chromatin fiber structure for human interphase fibroblasts in vivo was 45% zig-zag 30 nm fibers and 55% 10 nm fibers.

Ultrastructural Details of Mammalian Chromosome Architecture

Over the past decade, 3C-related methods have provided remarkable insights into chromosome folding in vivo. To overcome the limited resolution of prior studies, we extend a recently developed Hi-C variant, Micro-C, to map chromosome architecture at nucleosome resolution in human ESCs and fibroblasts. Micro-C robustly captures known features of chromosome folding including compartment organization, topologically associating domains, and interactions between CTCF binding sites. In addition, Micro-C provides a detailed map of nucleosome positions and localizes contact domain boundaries with nucleosomal precision. Compared to Hi-C, Micro-C exhibits an order of magnitude greater dynamic range, allowing the identification of ∼20,000 additional loops in each cell type. Many newly identified peaks are localized along extrusion stripes and form transitive grids, consistent with their anchors being pause sites impeding cohesin-dependent loop extrusion. Our analyses comprise the highest-resolution maps of chromosome folding in human cells to date, providing a valuable resource for studies of chromosome organization.

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.

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.

[Opposite Effects of Histone H1 and HMGN5 Protein on Distant Interactions in Chromatin]

Transcriptional enhancers in the cell nuclei typically interact with the target promoters in cis over long stretches of chromatin, but the mechanism of this communication remains unknown. Previously we have developed a defined in vitro system for quantitative analysis of the rate of distant enhancer-promoter communication (EPC) and have shown that the chromatin fibers maintain efficient distant EPC in cis. Here we investigate the roles of linker histone H1 and HMGN5 protein in EPC. A considerable negative effect of histone H1 on EPC depending on its C- and N-tails was shown. Protein HMGN5 that affects chromatin compaction and is associated with active chromatin counteracts EPC inhibition by H1. The data suggest that the efficiency of the interaction between the enhancer and the promoter depends on the structure and dynamics of the chromatin fiber localized between them and can be regulated by proteins associated with chromatin.

Role of nucleosome positioning in 3D chromatin organization and loop formation

We present a physics-based polymer model that can investigate 3D organization of chromatin accounting for DNA elasticity, DNA-bending due to nucleosomes, and 1D organization of nucleosomes along DNA. We find that the packing density of chromatin oscillates between densities corresponding to highly folded and extended configurations as we change the nucleosome organization (length of linker DNA). We compute the looping probability of chromatin and show that the presence of nucleosomes increases the looping probability of the chain compared to that of a bare DNA. We also show that looping probability has a large variability depending on the nature of nucleosome organization and density of linker histones.

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.