Bottom-Up Meets Top-Down: The Crossroads of Multiscale Chromatin Modeling

Coarse-Grained Modeling Using Neural Networks Trained on Structural Data

We propose a method of bottom-up coarse-graining, in which interactions within a coarse-grained model are determined by an artificial neural network trained on structural data obtained from multiple atomistic simulations. The method uses ideas of the inverse Monte Carlo approach, relating changes in the neural network weights with changes in average structural properties, such as radial distribution functions. As a proof of concept, we demonstrate the method on a system interacting by a Lennard-Jones potential modeled by a simple linear network and a single-site coarse-grained model of methanol-water solutions. In the latter case, we implement a nonlinear neural network with intermediate layers trained by atomistic simulations carried out at different methanol concentrations. We show that such a network acts as a transferable potential at the coarse-grained resolution for a wide range of methanol concentrations, including those not included in the training set.

Studies of the Mechanism of Nucleosome Dynamics: A Review on Multifactorial Regulation from Computational and Experimental Cases

The nucleosome, which organizes the long coil of genomic DNA in a highly condensed, polymeric way, is thought to be the basic unit of chromosomal structure. As the most important protein–DNA complex, its structural and dynamic features have been successively revealed in recent years. However, its regulatory mechanism, which is modulated by multiple factors, still requires systemic discussion. This study summarizes the regulatory factors of the nucleosome’s dynamic features from the perspective of histone modification, DNA methylation, and the nucleosome-interacting factors (transcription factors and nucleosome-remodeling proteins and cations) and focuses on the research exploring the molecular mechanism through both computational and experimental approaches. The regulatory factors that affect the dynamic features of nucleosomes are also discussed in detail, such as unwrapping, wrapping, sliding, and stacking. Due to the complexity of the high-order topological structures of nucleosomes and the comprehensive effects of regulatory factors, the research on the functional modulation mechanism of nucleosomes has encountered great challenges. The integration of computational and experimental approaches, the construction of physical modes for nucleosomes, and the application of deep learning techniques will provide promising opportunities for further exploration.

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.

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.

Reconstruct high-resolution 3D genome structures for diverse cell-types using FLAMINGO

High-resolution reconstruction of spatial chromosome organizations from chromatin contact maps is highly demanded, but is hindered by extensive pairwise constraints, substantial missing data, and limited resolution and cell-type availabilities. Here, we present FLAMINGO, a computational method that addresses these challenges by compressing inter-dependent Hi-C interactions to delineate the underlying low-rank structures in 3D space, based on the low-rank matrix completion technique. FLAMINGO successfully generates 5 kb- and 1 kb-resolution spatial conformations for all chromosomes in the human genome across multiple cell-types, the largest resources to date. Compared to other methods using various experimental metrics, FLAMINGO consistently demonstrates superior accuracy in recapitulating observed structures with raises in scalability by orders of magnitude. The reconstructed 3D structures efficiently facilitate discoveries of higher-order multi-way interactions, imply biological interpretations of long-range QTLs, reveal geometrical properties of chromatin, and provide high-resolution references to understand structural variabilities. Importantly, FLAMINGO achieves robust predictions against high rates of missing data and significantly boosts 3D structure resolutions. Moreover, FLAMINGO shows vigorous cross cell-type structure predictions that capture cell-type specific spatial configurations via integration of 1D epigenomic signals. FLAMINGO can be widely applied to large-scale chromatin contact maps and expand high-resolution spatial genome conformations for diverse cell-types.

High-resolution reconstruction of spatial chromosome organisation is in demand. Here the authors report FLAMINGO, for reconstructing high-resolution 3D Genome Organisation from HiC data which they use to generate both 5 kb and 1 kb-resolution 3D chromosomal structures for the human genome.

On the stability and layered organization of protein-DNA condensates

Multi-component phase separation is emerging as a key mechanism for the formation of biological condensates that play essential roles in signal sensing and transcriptional regulation. The molecular factors that dictate these condensates' stability and spatial organization are not fully understood, and it remains challenging to predict their microstructures. Using a near-atomistic, chemically accurate force field, we studied the phase behavior of chromatin regulators that are crucial for heterochromatin organization and their interactions with DNA. Our computed phase diagrams recapitulated previous experimental findings on different proteins. They revealed a strong dependence of condensate stability on the protein-DNA mixing ratio as a result of balancing protein-protein interactions and charge neutralization. Notably, a layered organization was observed in condensates formed by mixing HP1, histone H1, and DNA. This layered organization may be of biological relevance, as it enables cooperative DNA packaging between the two chromatin regulators: histone H1 softens the DNA to facilitate the compaction induced by HP1 droplets. Our study supports near-atomistic models as a valuable tool for characterizing the structure and stability of biological condensates.

First-passage statistics of colloids on fractals: Theory and experimental realization

In nature and technology, particle dynamics frequently occur in complex environments, for example in restricted geometries or crowded media. These dynamics have often been modeled invoking a fractal structure of the medium although the fractal structure was only indirectly inferred through the dynamics. Moreover, systematic studies have not yet been performed. Here, colloidal particles moving in a laser speckle pattern are used as a model system. In this case, the experimental observations can be reliably traced to the fractal structure of the underlying medium with an adjustable fractal dimension. First-passage time statistics reveal that the particles explore the speckle in a self-similar, fractal manner at least over four decades in time and on length scales up to 20 times the particle radius. The requirements for fractal diffusion to be applicable are laid out, and methods to extract the fractal dimension are established.

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.

Histone acylations and chromatin dynamics: concepts, challenges, and links to metabolism

In eukaryotic cells, DNA is tightly packed with the help of histone proteins into chromatin. Chromatin architecture can be modified by various post-translational modifications of histone proteins. For almost 60 years now, studies on histone lysine acetylation have unraveled the contribution of this acylation to an open chromatin state with increased DNA accessibility, permissive for gene expression. Additional complexity emerged from the discovery of other types of histone lysine acylations. The acyl group donors are products of cellular metabolism, and distinct histone acylations can link the metabolic state of a cell with chromatin architecture and contribute to cellular adaptation through changes in gene expression. Currently, various technical challenges limit our full understanding of the actual impact of most histone acylations on chromatin dynamics and of their biological relevance. In this review, we summarize the state of the art and provide an overview of approaches to overcome these challenges. We further discuss the concept of subnuclear metabolic niches that could regulate local CoA availability and thus couple cellular metabolisms with the epigenome.

Histone acylations and chromatin dynamics: concepts, challenges, and links to metabolism

In eukaryotic cells, DNA is tightly packed with the help of histone proteins into chromatin. Chromatin architecture can be modified by various post-translational modifications of histone proteins. For almost 60 years now, studies on histone lysine acetylation have unraveled the contribution of this acylation to an open chromatin state with increased DNA accessibility, permissive for gene expression. Additional complexity emerged from the discovery of other types of histone lysine acylations. The acyl group donors are products of cellular metabolism, and distinct histone acylations can link the metabolic state of a cell with chromatin architecture and contribute to cellular adaptation through changes in gene expression. Currently, various technical challenges limit our full understanding of the actual impact of most histone acylations on chromatin dynamics and of their biological relevance. In this review, we summarize the state of the art and provide an overview of approaches to overcome these challenges. We further discuss the concept of subnuclear metabolic niches that could regulate local CoA availability and thus couple cellular metabolisms with the epigenome.

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.

Separation and Characterization of Endogenous Nucleosomes by Native Capillary Zone Electrophoresis-Top-Down Mass Spectrometry

We report a novel platform [native capillary zone electrophoresis-top-down mass spectrometry (nCZE-TDMS)] for the separation and characterization of whole nucleosomes, their histone subunits, and post-translational modifications (PTMs). As the repeating unit of chromatin, mononucleosomes (Nucs) are an ∼200 kDa complex of DNA and histone proteins involved in the regulation of key cellular processes central to human health and disease. Unraveling the covalent modification landscape of histones and their defined stoichiometries within Nucs helps to explain epigenetic regulatory mechanisms. In nCZE-TDMS, online Nuc separation is followed by a three-tier tandem MS approach that measures the intact mass of Nucs, ejects and detects the constituent histones, and fragments to sequence the histone. The new platform was optimized with synthetic Nucs to significantly reduce both sample requirements and cost compared to direct infusion. Limits of detection were in the low-attomole range, with linearity of over ∼3 orders of magnitude. The nCZE-TDMS platform was applied to endogenous Nucs from two cell lines distinguished by overexpression or knockout of histone methyltransferase NSD2/MMSET, where analysis of constituent histones revealed changes in histone abundances over the course of the CZE separation. We are confident the nCZE-TDMS platform will help advance nucleosome-level research in the fields of chromatin and epigenetics.

Decoding the protein composition of whole nucleosomes with Nuc-MS

Current proteomic approaches disassemble and digest nucleosome particles, blurring readouts of the 'histone code'. To preserve nucleosome-level information, we developed Nuc-MS, which displays the landscape of histone variants and their post-translational modifications (PTMs) in a single mass spectrum. Combined with immunoprecipitation, Nuc-MS quantified nucleosome co-occupancy of histone H3.3 with variant H2A.Z (sixfold over bulk) and the co-occurrence of oncogenic H3.3K27M with euchromatic marks (for example, a >15-fold enrichment of dimethylated H3K79me2). Nuc-MS is highly concordant with chromatin immunoprecipitation-sequencing (ChIP-seq) and offers a new readout of nucleosome-level biology.

Effects of H2A.B incorporation on nucleosome structures and dynamics

The H2A.B histone variant is an epigenetic regulator involved in transcriptional upregulation, DNA synthesis, and splicing that functions by replacing the canonical H2A histone in the nucleosome core particle. Introduction of H2A.B results in less compact nucleosome states with increased DNA unwinding and accessibility at the nucleosomal entry and exit sites. Despite being well characterized experimentally, the molecular mechanisms by which H2A.B incorporation alters nucleosome stability and dynamics remain poorly understood. To study the molecular mechanisms of H2A.B, we have performed a series of conventional and enhanced sampling molecular dynamics simulation of H2A.B- and canonical H2A-containing nucleosomes. Results of conventional simulations show that H2A.B weakens protein-protein and protein-DNA interactions at specific locations throughout the nucleosome. These weakened interactions result in significantly more DNA opening from both the entry and exit sites in enhanced sampling simulations. Furthermore, free energy profiles show that H2A.B-containing nucleosomes have significantly broader free wells and that H2A.B allows for sampling of states with increased DNA breathing, which are shown to be stable on the hundreds of nanoseconds timescale with further conventional simulations. Together, our results show the molecular mechanisms by which H2A.B creates less compacted nucleosome states as a means of increasing genetic accessibility and gene transcription.

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

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

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

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.

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.

Understanding and Engineering Chromatin as a Dynamical System across Length and Timescales

Connecting the molecular structure and function of chromatin across length and timescales remains a grand challenge to understanding and engineering cellular behaviors. Across five orders of magnitude, dynamic processes constantly reshape chromatin structures, driving spaciotemporal patterns of gene expression and cell fate. Through the interplay of structure and function, the genome operates as a highly dynamic feedback control system. Recent experimental techniques have provided increasingly detailed data that revise and augment the relatively static, hierarchical view of genomic architecture with an understanding of how dynamic processes drive organization. Here, we review how novel technologies from sequencing, imaging, and synthetic biology refine our understanding of chromatin structure and function and enable chromatin engineering. Finally, we discuss opportunities to use these tools to enhance understanding of the dynamic interrelationship of chromatin structure and function.