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

CG modeling of nucleosome arrays reveals the salt-dependent chromatin fiber conformational variability

Eukaryotic DNA is packaged in the cell nucleus into chromatin, composed of arrays of DNA-histone protein octamer complexes, the nucleosomes. Over the past decade, it has become clear that chromatin structure in vivo is not a hierarchy of well-organized folded nucleosome fibers but displays considerable conformational variability and heterogeneity. In vitro and in vivo studies, as well as computational modeling, have revealed that attractive nucleosome-nucleosome interaction with an essential role of nucleosome stacking defines chromatin compaction. The internal structure of compacted nucleosome arrays is regulated by the flexible and dynamic histone N-terminal tails. Since DNA is a highly negatively charged polyelectrolyte, electrostatic forces make a decisive contribution to chromatin formation and require the histones, particularly histone tails, to carry a significant positive charge. This also results in an essential role of mobile cations of the cytoplasm (K+, Na+, Mg2+) in regulating electrostatic interactions. Building on a previously successfully established bottom-up coarse-grained (CG) nucleosome model, we have developed a CG nucleosome array (chromatin fiber) model with the explicit presence of mobile ions and studied its conformational variability as a function of Na+ and Mg2+ ion concentration. With progressively elevated ion concentrations, we identified four main conformational states of nucleosome arrays characterized as extended, flexible, nucleosome-clutched, and globular fibers.

Nanoscale Characterization of Interaction of Nucleosomes with H1 Linker Histone

In eukaryotic nuclei, DNA is wrapped around an octamer of core histones to form nucleosomes. H1 binds to the linker DNA of nucleosome to form the chromatosome, the next structural unit of chromatin. Structural features on individual chromatosomes contribute to chromatin structure, but not fully characterized. In addition to canonical nucleosomes composed of two copies each of histones H2A, H2B, H3, and H4 (H3 nucleosomes), centromeres chromatin contain nucleosomes in which H3 is replaced with its analog CENP-A, changing structural properties of CENP-A nucleosomes. Nothing is known about the interaction of H1 with CENP-A nucleosomes. Here we filled this gap and characterized the interaction of H1 histone with both types of nucleosomes. H1 does bind both types of the nucleosomes forming more compact chromosome particles with elevated affinity to H3 nucleosomes. H1 binding significantly increases the stability of chromatosomes preventing their spontaneous dissociation. In addition to binding to the entry-exit position of the DNA arms identified earlier, H1 is capable of bridging of distant DNA segments. H1 binding leads to the assembly of mononucleosomes in aggregates, stabilized by internucleosome interactions as well as bridging of the DNA arms of chromatosomes. Contribution of these finding to the chromatin structure and functions are discussed.

Nucleosome condensate and linker DNA alter chromatin folding pathways and rates

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

Significance

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

From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization

Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.

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.

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.

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.

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.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

The effects of RNA.DNA-DNA triple helices on nucleosome structures and dynamics

Noncoding RNAs (ncRNAs) are an emerging epigenetic factor and have been recognized as playing a key role in many gene expression pathways. Structurally, binding of ncRNAs to isolated DNA is strongly dependent on sequence complementary and results in the formation of an RNA.DNA-DNA (RDD) triple helix. However, in vivo DNA is not isolated but is rather packed in chromatin fibers, the fundamental unit of which is the nucleosome. Biochemical experiments have shown that ncRNA binding to nucleosomal DNA is elevated at DNA entry and exit sites and is dependent on the presence of the H3 N-terminal tails. However, the structural and dynamical bases for these mechanisms remain unknown. Here, we have examined the mechanisms and effects of RDD formation in the context of the nucleosome using a series of all-atom molecular dynamics simulations. Results highlight the importance of DNA sequence on complex stability, elucidate the effects of the H3 tails on RDD structures, show how RDD formation impacts the structure and dynamics of the H3 tails, and show how RNA alters the local and global DNA double-helical structure. Together, our results suggest ncRNAs can modify nucleosome, and potentially higher-order chromatin, structures and dynamics as a means of exerting epigenetic control.

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.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

CENP-N promotes the compaction of centromeric chromatin

The histone variant CENP-A is the epigenetic determinant for the centromere, where it is interspersed with canonical H3 to form a specialized chromatin structure that nucleates the kinetochore. How nucleosomes at the centromere arrange into higher order structures is unknown. Here we demonstrate that the human CENP-A-interacting protein CENP-N promotes the stacking of CENP-A-containing mononucleosomes and nucleosomal arrays through a previously undefined interaction between the α6 helix of CENP-N with the DNA of a neighboring nucleosome. We describe the cryo-EM structures and biophysical characterization of such CENP-N-mediated nucleosome stacks and nucleosomal arrays and demonstrate that this interaction is responsible for the formation of densely packed chromatin at the centromere in the cell. Our results provide first evidence that CENP-A, together with CENP-N, promotes specific chromatin higher order structure at the centromere.

A glitch in the snitch: the role of linker histone H1 in shaping the epigenome in normal and diseased cells

Histone H1s or the linker histones are a family of dynamic chromatin compacting proteins that are essential for higher-order chromatin organization. These highly positively charged proteins were previously thought to function solely as repressors of transcription. However, over the last decade, there is a growing interest in understanding this multi-protein family, finding that not all variants act as repressors. Indeed, the H1 family members appear to have distinct affinities for chromatin and may potentially affect distinct functions. This would suggest a more nuanced contribution of H1 to chromatin organization. The advent of new technologies to probe H1 dynamics in vivo, combined with powerful computational biology, and in vitro imaging tools have greatly enhanced our knowledge of the mechanisms by which H1 interacts with chromatin. This family of proteins can be metaphorically compared to the Golden Snitch from the Harry Potter series, buzzing on and off several regions of the chromatin, in combat with competing transcription factors and chromatin remodellers, thereby critical to the epigenetic endgame on short and long temporal scales in the life of the nucleus. Here, we summarize recent efforts spanning structural, computational, genomic and genetic experiments which examine the linker histone as an unseen architect of chromatin fibre in normal and diseased cells and explore unanswered fundamental questions in the field.

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.

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.