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

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.

Combining molecular dynamics simulations and scoring method to computationally model ubiquitylated linker histones in chromatosomes

The chromatin in eukaryotic cells plays a fundamental role in all processes during a cell's life cycle. This nucleoprotein is normally tightly packed but needs to be unpacked for expression and division. The linker histones are critical for such packaging processes and while most experimental and simulation works recognize their crucial importance, the focus is nearly always set on the nucleosome as the basic chromatin building block. Linker histones can undergo several modifications, but only few studies on their ubiquitylation have been conducted. Mono-ubiquitylated linker histones (HUb), while poorly understood, are expected to influence DNA compaction. The size of ubiquitin and the globular domain of the linker histone are comparable and one would expect an increased disorder upon ubiquitylation of the linker histone. However, the formation of higher order chromatin is not hindered and ubiquitylation of the linker histone may even promote gene expression. Structural data on chromatosomes is rare and HUb has never been modeled in a chromatosome so far. Descriptions of the chromatin complex with HUb would greatly benefit from computational structural data. In this study we generate molecular dynamics simulation data for six differently linked HUb variants with the help of a sampling scheme tailored to drive the exploration of phase space. We identify conformational sub-states of the six HUb variants using the sketch-map algorithm for dimensionality reduction and iterative HDBSCAN for clustering on the excessively sampled, shallow free energy landscapes. We present a highly efficient geometric scoring method to identify sub-states of HUb that fit into the nucleosome. We predict HUb conformations inside a nucleosome using on-dyad and off-dyad chromatosome structures as reference and show that unbiased simulations of HUb produce significantly more fitting than non-fitting HUb conformations. A tetranucleosome array is used to show that ubiquitylation can even occur in chromatin without too much steric clashes.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

The Effects of Histone H2B Ubiquitylations on the Nucleosome Structure and Internucleosomal Interactions

Eukaryotic gene compaction takes place at multiple levels to package DNA to chromatin and chromosomes. Two of the most fundamental levels of DNA packaging are at the nucleosome and dinucleosome stacks. The nucleosome is the basic gene-packing unit and is composed of DNA wrapped around a histone core. Nucleosomes stack with one another for further compaction of DNA. The first stacking step leads to dinucleosome formation, which is driven by internucleosomal interactions between various parts of two nucleosomes. Histone proteins are rich targets for post-translational modifications, some of which affect the structure of the nucleosome and the interactions between nucleosomes. These effects are often implicated in the regulation of various genomic transactions. In particular, histone H2B ubiquitylation has been associated with facilitated transcription and hexasome formation. Here, we employed semi-synthetically ubiquitylated histone H2B and single-molecule FRET to investigate the effects of H2B ubiquitylations at lysine 34 (H2BK34) and lysine 120 (H2BK120) on the structure of the nucleosome and the interactions between two nucleosomes. Our results suggest that H2BK34 ubiquitylation widens the DNA gyre gap in the nucleosome and stabilizes long- and short-range internucleosomal interactions while H2BK120 ubiquitylation does not affect the nucleosome structure or internucleosomal interactions. These results suggest potential roles for H2B ubiquitylations in facilitated transcription and hexasome formation while maintaining the structural integrity of chromatin.

Chromatin as self-returning walks: From population to single cell and back

With a growing understanding of the chromatin structure, many efforts remain focused on bridging the gap between what is suggested by population-averaged data and what is visualized for single cells. A popular approach to traversing these scales is to fit a polymer model to Hi-C contact data. However, Hi-C is an average of millions to billions of cells, and each cell may not contain all population-averaged contacts. Therefore, we employ a novel approach of summing individual chromosome trajectories-determined by our Self-Returning Random Walk model-to create populations of cells. We allow single cells to consist of disparate structures and reproduce a variety of experimentally relevant contact maps. We show that the amount of shared topology between cells, and their mechanism of formation, changes the population-averaged structure. Therefore, we present a modeling technique that, with few constraints and little oversight, can be used to understand which single-cell chromatin structures underlie population-averaged behavior.

Topological Constraints with Optimal Length Promote the Formation of Chromosomal Territories at Weakened Degree of Phase Separation

It is generally agreed that the nuclei of eukaryotic cells at interphase are partitioned into disjointed territories, with distinct regions occupied by certain chromosomes. However, the underlying mechanism for such territorialization is still under debate. Here we model chromosomes as coarse-grained block copolymers and to investigate the effect of loop domains (LDs) on the formation of compartments and territories based on dissipative particle dynamics. A critical length of LDs, which depends sensitively on the length of polymeric blocks, is obtained to minimize the degree of phase separation. This also applies to the two-polymer system: The critical length not only maximizes the degree of territorialization but also minimizes the degree of phase separation. Interestingly, by comparing with experimental data, we find the critical length for LDs and the corresponding length of blocks to be respectively very close to the mean length of topologically associating domains (TADs) and chromosomal segments with different densities of CpG islands for human chromosomes. The results indicate that topological constraints with optimal length can contribute to the formation of territories by weakening the degree of phase separation, which likely promotes the chromosomal flexibility in response to genetic regulations.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

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.

Is the H4 histone tail intrinsically disordered or intrinsically multifunctional?

The structural versatility of histone tails is one of the key elements in the organisation of chromatin, which allows for the compact storage of genomic information. However, this structural diversity also complicates experimental and computational studies. Here, the potential and free energy landscape for the isolated and bound H4 histone tail are explored. The landscapes exhibit a set of distinct structural ensembles separated by high energy barriers, with little difference between isolated and bound tails. This consistency is a desirable feature that facilitates the formation of transient interactions, which are required for the liquid-like chromatin organisation. The existence of multiple, distinct structures on a multifunnel energy landscape is likely to be associated with multifunctionality, i.e. a set of evolved, distinct functions. Contrasting it with previously reported results for other disordered peptides, this type of landscape may be associated with a conformational selection based binding mechanism. Given the similarity to other systems exhibiting similar multifunnel energy landscapes, the disorder in histone tails might be better described in context of multifunctionality.

Three-dimensional genome organization in epigenetic regulations: cause or consequence?

The evolution of the nucleus is an evolutionary milestone. By enabling genome compartmentalization, it contributes to the fine-tuning of genome functions. The genome is partitioned into functional domains differing in spatial positioning and topological folding at different scales. The rise of '3D Genomics' embracing experimental, theoretical, and modeling approaches allowed the proposal of a multiscale model of the eukaryotic genome, capturing its organizing principles and functionalities. In these efforts, resolving causality remains an important objective. Are positioning and folding the cause or consequence of functional states? This minireview presents emerging answers to this question, borrowing examples from recent studies of the three-dimensional genome in both plants and animals.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

3D reconstruction of genomic regions from sparse interaction data

Chromosome conformation capture (3C) technologies measure the interaction frequency between pairs of chromatin regions within the nucleus in a cell or a population of cells. Some of these 3C technologies retrieve interactions involving non-contiguous sets of loci, resulting in sparse interaction matrices. One of such 3C technologies is Promoter Capture Hi-C (pcHi-C) that is tailored to probe only interactions involving gene promoters. As such, pcHi-C provides sparse interaction matrices that are suitable to characterize short- and long-range enhancer-promoter interactions. Here, we introduce a new method to reconstruct the chromatin structural (3D) organization from sparse 3C-based datasets such as pcHi-C. Our method allows for data normalization, detection of significant interactions and reconstruction of the full 3D organization of the genomic region despite of the data sparseness. Specifically, it builds, with as low as the 2-3% of the data from the matrix, reliable 3D models of similar accuracy of those based on dense interaction matrices. Furthermore, the method is sensitive enough to detect cell-type-specific 3D organizational features such as the formation of different networks of active gene communities.

Contribution of advanced fluorescence nano microscopy towards revealing mitotic chromosome structure

The organization of chromatin into higher-order structures and its condensation process represent one of the key challenges in structural biology. This is important for elucidating several disease states. To address this long-standing problem, development of advanced imaging methods has played an essential role in providing understanding into mitotic chromosome structure and compaction. Amongst these are two fast evolving fluorescence imaging technologies, specifically fluorescence lifetime imaging (FLIM) and super-resolution microscopy (SRM). FLIM in particular has been lacking in the application of chromosome research while SRM has been successfully applied although not widely. Both these techniques are capable of providing fluorescence imaging with nanometer information. SRM or "nanoscopy" is capable of generating images of DNA with less than 50 nm resolution while FLIM when coupled with energy transfer may provide less than 20 nm information. Here, we discuss the advantages and limitations of both methods followed by their contribution to mitotic chromosome studies. Furthermore, we highlight the future prospects of how advancements in new technologies can contribute in the field of chromosome science.

A Paradigm Revolution or Just Better Resolution-Will Newly Emerging Superresolution Techniques Identify Chromatin Architecture as a Key Factor in Radiation-Induced DNA Damage and Repair Regulation?

DNA double-strand breaks (DSBs) have been recognized as the most serious lesions in irradiated cells. While several biochemical pathways capable of repairing these lesions have been identified, the mechanisms by which cells select a specific pathway for activation at a given DSB site remain poorly understood. Our knowledge of DSB induction and repair has increased dramatically since the discovery of ionizing radiation-induced foci (IRIFs), initiating the possibility of spatiotemporally monitoring the assembly and disassembly of repair complexes in single cells. IRIF exploration revealed that all post-irradiation processes-DSB formation, repair and misrepair-are strongly dependent on the characteristics of DSB damage and the microarchitecture of the whole affected chromatin domain in addition to the cell status. The microscale features of IRIFs, such as their morphology, mobility, spatiotemporal distribution, and persistence kinetics, have been linked to repair mechanisms. However, the influence of various biochemical and structural factors and their specific combinations on IRIF architecture remains unknown, as does the hierarchy of these factors in the decision-making process for a particular repair mechanism at each individual DSB site. New insights into the relationship between the physical properties of the incident radiation, chromatin architecture, IRIF architecture, and DSB repair mechanisms and repair efficiency are expected from recent developments in optical superresolution microscopy (nanoscopy) techniques that have shifted our ability to analyze chromatin and IRIF architectures towards the nanoscale. In the present review, we discuss this relationship, attempt to correlate still rather isolated nanoscale studies with already better-understood aspects of DSB repair at the microscale, and consider whether newly emerging "correlated multiscale structuromics" can revolutionarily enhance our knowledge in this field.

Nano-Surveillance: Tracking Individual Molecules in a Sea of Chromatin

Chromatin is the epigenomic platform for diverse nuclear processes such as DNA repair, replication, transcription, telomere, and centromere function. In cancer cells, mutations in key processes result in DNA amplification, chromosome translocations, and chromothripsis, severely distorting the natural chromatin state. In normal and diseased states, dozens of chromatin effectors alter the physical integrity and dynamics of chromatin at the level of both single nucleosomes and arrays of nucleosomes folded into 3-dimensional shapes. Integrating these length scales, from the 10 nm sized nucleosome to mitotic chromosomes, whilst jostling within the crowded environment of the cell, cannot yet be achieved by a single technology. In this review, we discuss tools that have proven powerful in the investigation of nucleosome and chromatin fiber dynamics. We also provide a deeper focus into atomic force microscopy (AFM) applications that can bridge diverse length and time scales. Using time course AFM, we observe that chromatin condensation by H1.5 is dynamic, whereas using nano-indentation force spectroscopy we observe that both histone variants and nucleosome binding partners alter material properties of individual nucleosomes. Finally, we demonstrate how high-speed AFM can visualize plasmid DNA dynamics, intermittent nucleosome-nucleosome contacts, and changes in nucleosome phasing along a contiguous chromatin fiber. Altogether, the development of innovative technologies holds the promise of revealing the secret lives of nucleosomes, potentially bridging the gaps in our understanding of how chromatin works within living cells and tissues.

Cryo-nanoscale chromosome imaging-future prospects

The high-order structure of mitotic chromosomes remains to be fully elucidated. How nucleosomes compact at various structural levels into a condensed mitotic chromosome is unclear. Cryogenic preservation and imaging have been applied for over three decades, keeping biological structures close to the native in vivo state. Despite being extensively utilized, this field is still wide open for mitotic chromosome research. In this review, we focus specifically on cryogenic efforts for determining the mitotic nanoscale chromatin structures. We describe vitrification methods, current status, and applications of advanced cryo-microscopy including future tools required for resolving the native architecture of these fascinating structures that hold the instructions to life.