Elucidating the influence of linker histone variants on chromatosome dynamics and energetics

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.

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.

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.

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.

Post-Translation Modifications and Mutations of Human Linker Histone Subtypes: Their Manifestation in Disease

Linker histones (LH) are a critical component of chromatin in addition to the canonical histones (H2A, H2B, H3, and H4). In humans, 11 subtypes (7 somatic and 4 germinal) of linker histones have been identified, and their diverse cellular functions in chromatin structure, DNA replication, DNA repair, transcription, and apoptosis have been explored, especially for the somatic subtypes. Delineating the unique role of human linker histone (hLH) and their subtypes is highly tedious given their high homology and overlapping expression patterns. However, recent advancements in mass spectrometry combined with HPLC have helped in identifying the post-translational modifications (PTMs) found on the different LH subtypes. However, while a number of PTMs have been identified and their potential nuclear and non-nuclear functions explored in cellular processes, there are very few studies delineating the direct relevance of these PTMs in diseases. In addition, recent whole-genome sequencing of clinical samples from cancer patients and individuals afflicted with Rahman syndrome have identified high-frequency mutations and therefore broadened the perspective of the linker histone mutations in diseases. In this review, we compile the identified PTMs of hLH subtypes, current knowledge of the relevance of hLH PTMs in human diseases, and the correlation of PTMs coinciding with mutations mapped in diseases.

The Role of PARP1 and PAR in ATP-Independent Nucleosome Reorganisation during the DNA Damage Response

The functioning of the eukaryotic cell genome is mediated by sophisticated protein-nucleic-acid complexes, whose minimal structural unit is the nucleosome. After the damage to genomic DNA, repair proteins need to gain access directly to the lesion; therefore, the initiation of the DNA damage response inevitably leads to local chromatin reorganisation. This review focuses on the possible involvement of PARP1, as well as proteins acting nucleosome compaction, linker histone H1 and non-histone chromatin protein HMGB1. The polymer of ADP-ribose is considered the main regulator during the development of the DNA damage response and in the course of assembly of the correct repair complex.

Genome-wide characterization and evolutionary analysis of linker histones in castor bean (Ricinus communis)

H1s, or linker histones, are ubiquitous proteins in eukaryotic cells, consisting of a globular GH1 domain flanked by two unstructured tails. Whilst it is known that numerous non-allelic variants exist within the same species, the degree of interspecific and intraspecific variation and divergence of linker histones remain unknown. The conserved basic binding sites in GH1 and evenly distributed strong positive charges on the C-terminal domain (CTD) are key structural characters for linker histones to bind chromatin. Based on these features, we identified five linker histones from 13 GH1-containing proteins in castor bean (Ricinus communis), which were named as RcH1.1, RcH1.2a, RcH1.2b, RcH1.3, and RcH1.4 based on their phylogenetic relationships with the H1s from five other economically important Euphorbiaceae species (Hevea brasiliensis Jatropha curcas, Manihot esculenta Mercurialis annua, and Vernicia fordii) and Arabidopsis thaliana. The expression profiles of RcH1 genes in a variety of tissues and stresses were determined from RNA-seq data. We found three RcH1 genes (RcH1.1, RcH1.2a, and RcH1.3) were broadly expressed in all tissues, suggesting a conserved role in stabilizing and organizing the nuclear DNA. RcH1.2a and RcH1.4 was preferentially expressed in floral tissues, indicating potential involvement in floral development in castor bean. Lack of non-coding region and no expression detected in any tissue tested suggest that RcH1.2b is a pseudogene. RcH1.3 was salt stress inducible, but not induced by cold, heat and drought in our investigation. Structural comparison confirmed that GH1 domain was highly evolutionarily conserved and revealed that N- and C-terminal domains of linker histones are divergent between variants, but highly conserved between species for a given variant. Although the number of H1 genes varies between species, the number of H1 variants is relatively conserved in more closely related species (such as within the same family). Through comparison of nucleotide diversity of linker histone genes and oil-related genes, we found similar mutation rate of these two groups of genes. Using Tajima's D and ML-HKA tests, we found RcH1.1 and RcH1.3 may be under balancing selection.

Recent Advances in Investigating Functional Dynamics of Chromatin

Dynamics spanning the picosecond-minute time domain and the atomic-subcellular spatial window have been observed for chromatin in vitro and in vivo. The condensed organization of chromatin in eukaryotic cells prevents regulatory factors from accessing genomic DNA, which requires dynamic stabilization and destabilization of structure to initiate downstream DNA activities. Those processes are achieved through altering conformational and dynamic properties of nucleosomes and nucleosome–protein complexes, of which delineating the atomistic pictures is essential to understand the mechanisms of chromatin regulation. In this review, we summarize recent progress in determining chromatin dynamics and their modulations by a number of factors including post-translational modifications (PTMs), incorporation of histone variants, and binding of effector proteins. We focus on experimental observations obtained using high-resolution techniques, primarily including nuclear magnetic resonance (NMR) spectroscopy, Förster (or fluorescence) resonance energy transfer (FRET) microscopy, and molecular dynamics (MD) simulations, and discuss the elucidated dynamics in the context of functional response and relevance.

Can aggressive cancers be identified by the "aggressiveness" of their chromatin?

Phenotypic plasticity is a crucial feature of aggressive cancer, providing the means for cancer progression. Stochastic changes in tumor cell transcriptional programs increase the chances of survival under any condition. I hypothesize that unstable chromatin permits stochastic transitions between transcriptional programs in aggressive cancers and supports non-genetic heterogeneity of tumor cells as a basis for their adaptability. I present a mechanistic model for unstable chromatin which includes destabilized nucleosomes, mobile chromatin fibers and random enhancer-promoter contacts, resulting in stochastic transcription. I suggest potential markers for "unsettled" chromatin in tumors associated with poor prognosis. Although many of the characteristics of unstable chromatin have been described, they were mostly used to explain changes in the transcription of individual genes. I discuss approaches to evaluate the role of unstable chromatin in non-genetic tumor cell heterogeneity and suggest using the degree of chromatin instability and transcriptional noise in tumor cells to predict cancer aggressiveness.

Characteristic H3 N-tail dynamics in the nucleosome core particle, nucleosome, and chromatosome

The nucleosome core particle (NCP) comprises a histone octamer, wrapped around by ∼146-bp DNA, while the nucleosome additionally contains linker DNA. We previously showed that, in the nucleosome, H4 N-tail acetylation enhances H3 N-tail acetylation by altering their mutual dynamics. Here, we have evaluated the roles of linker DNA and/or linker histone on H3 N-tail dynamics and acetylation by using the NCP and the chromatosome (i.e., linker histone H1.4-bound nucleosome). In contrast to the nucleosome, H3 N-tail acetylation and dynamics are greatly suppressed in the NCP regardless of H4 N-tail acetylation because the H3 N-tail is strongly bound between two DNA gyres. In the chromatosome, the asymmetric H3 N-tail adopts two conformations: one contacts two DNA gyres, as in the NCP; and one contacts linker DNA, as in the nucleosome. However, the rate of H3 N-tail acetylation is similar in the chromatosome and nucleosome. Thus, linker DNA and linker histone both regulate H3-tail dynamics and acetylation.

Structure, Dynamics and Stability of the Globular Domain of Human Linker Histone H1.0 and the Role of Positive Charges

Linker histone H1 (H1) is an abundant chromatin‐binding protein that acts as an epigenetic regulator binding to nucleosomes and altering chromatin structures and dynamics. Nonetheless, the mechanistic details of its function remain poorly understood. Recent work suggest that the number and position of charged side chains on the globular domain (GD) of H1 influence chromatin structure and hence gene repression. Here, we solved the solution structure of the unbound GD of human H1.0, revealing that the structure is almost completely unperturbed by complex formation, except for a loop connecting two antiparallel β‐strands. We further quantified the role of the many positive charges of the GD for its structure and conformational stability through the analysis of 11 charge variants. We find that modulating the number of charges has little effect on the structure, but the stability is affected, resulting in a difference in melting temperature of 26 K between GD of net charge +5 versus +13. This result suggests that the large number of positive charges on H1‐GDs have evolved for function rather than structure and high stability. The stabilization of the GD upon binding to DNA can thus be expected to have a pronounced electrostatic component, a contribution that is amenable to modulation by posttranslational modifications, especially acetylation and phosphorylation.

PDB Code(s): 6hq1;

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.

MNase Digestion Protection Patterns of the Linker DNA in Chromatosomes

The compact nucleosomal structure limits DNA accessibility and regulates DNA-dependent cellular activities. Linker histones bind to nucleosomes and compact nucleosomal arrays into a higher-order chromatin structure. Recent developments in high throughput technologies and structural computational studies provide nucleosome positioning at a high resolution and contribute to the information of linker histone location within a chromatosome. However, the precise linker histone location within the chromatin fibre remains unclear. Using monomer extension, we mapped core particle and chromatosomal positions over a core histone-reconstituted, 1.5 kb stretch of DNA from the chicken adult β-globin gene, after titration with linker histones and linker histone globular domains. Our results show that, although linker histone globular domains and linker histones display a wide variation in their binding affinity for different positioned nucleosomes, they do not alter nucleosome positions or generate new nucleosome positions. Furthermore, the extra ~20 bp of DNA protected in a chromatosome is usually symmetrically distributed at each end of the core particle, suggesting linker histones or linker histone globular domains are located close to the nucleosomal dyad axis.

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.

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.

Unraveling linker histone interactions in nucleosomes

Considerable progress has been made recently in defining the interactions of linker histones (H1s) within nucleosomes. Major advancements include atomic resolution structures of the globular domain of full-length H1s in the context of nucleosomes containing full-length linker DNA. Although these studies have led to a detailed understanding of the interactions and dynamics of H1 globular domains in the canonical on-dyad nucleosome binding pocket, more information regarding the intrinsically disordered N-terminal and C-terminal domains is needed. In this review, we highlight studies supporting our current understanding of the structures and interactions of the N-terminal, globular, and C-terminal domains of linker histones within the nucleosome.

Role of Position K+4 in the Phosphorylation and Dephosphorylation Reaction Kinetics of the CheY Response Regulator

Two-component signaling is a primary method by which microorganisms interact with their environments. A kinase detects stimuli and modulates autophosphorylation activity. The signal propagates by phosphotransfer from the kinase to a response regulator, eliciting a response. Response regulators operate over a range of time scales, corresponding to their related biological processes. Response regulator active site chemistry is highly conserved, but certain variable residues can influence phosphorylation kinetics. An Ala-to-Pro substitution (K+4, residue 113) in the Escherichia coli response regulator CheY triggers a constitutively active phenotype; however, the A113P substitution is too far from the active site to directly affect phosphochemistry. To better understand the activating mechanism(s) of the substitution, we analyzed receiver domain sequences to characterize the evolutionary role of the K+4 position. Although most featured Pro, Leu, Ile, and Val residues, chemotaxis-related proteins exhibited atypical Ala, Gly, Asp, and Glu residues at K+4. Structural and in silico analyses revealed that CheY A113P adopted a partially active configuration. Biochemical data showed that A113P shifted CheY toward a more activated state, enhancing autophosphorylation. By characterizing CheY variants, we determined that this functionality was transmitted through a hydrophobic network bounded by the β5α5 loop and the α1 helix of CheY. This region also interacts with the phosphodonor CheA_P1, suggesting that binding generates an activating perturbation similar to the A113P substitution. Atypical residues like Ala at the K+4 position likely serve two purposes. First, restricting autophosphorylation may minimize background noise generated by intracellular phosphodonors such as acetyl phosphate. Second, optimizing interactions with upstream partners may help prime the receiver domain for phosphorylation.

Extended and dynamic linker histone-DNA Interactions control chromatosome compaction

Chromatosomes play a fundamental role in chromatin regulation, but a detailed understanding of their structure is lacking, partially due to their complex dynamics. Using single-molecule DNA unzipping with optical tweezers, we reveal that linker histone interactions with DNA are remarkably extended, with the C-terminal domain binding both DNA linkers as far as approximately ±140 bp from the dyad. In addition to a symmetrical compaction of the nucleosome core governed by globular domain contacts at the dyad, the C-terminal domain compacts the nucleosome's entry and exit. These interactions are dynamic, exhibit rapid binding and dissociation, are sensitive to phosphorylation of a specific residue, and are crucial to determining the symmetry of the chromatosome's core. Extensive unzipping of the linker DNA, which mimics its invasion by motor proteins, shifts H1 into an asymmetric, off-dyad configuration and triggers nucleosome decompaction, highlighting the plasticity of the chromatosome structure and its potential regulatory role.

Histone tails as signaling antennas of chromatin

Histone tails, representing the N-terminal or C-terminal regions flanking the histone core, play essential roles in chromatin signaling networks. Intrinsic disorder of histone tails and their propensity for post-translational modifications allow them to serve as hubs in coordination of epigenetic processes within the nucleosomal context. Deposition of histone variants with distinct histone tail properties further enriches histone tails' repertoire in epigenetic signaling. Given the advances in experimental techniques and in silico modelling, we review the most recent data on histone tails' effects on nucleosome stability and dynamics, their function in regulating chromatin accessibility and folding. Finally, we discuss different molecular mechanisms to understand how histone tails are involved in nucleosome recognition by binding partners and formation of higher-order chromatin structures.

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.

Binding Dynamics of Disordered Linker Histone H1 with a Nucleosomal Particle

Linker histone H1 is an essential regulatory protein for many critical biological processes, such as eukaryotic chromatin packaging and gene expression. Mis-regulation of H1s is commonly observed in tumor cells, where the balance between different H1 subtypes has been shown to alter the cancer phenotype. Consisting of a rigid globular domain and two highly charged terminal domains, H1 can bind to multiple sites on a nucleosomal particle to alter chromatin hierarchical condensation levels. In particular, the disordered H1 amino- and carboxyl-terminal domains (NTD/CTD) are believed to enhance this binding affinity, but their detailed dynamics and functions remain unclear. In this work, we used a coarse-grained computational model, AWSEM-DNA, to simulate the H1.0b-nucleosome complex, namely chromatosome. Our results demonstrate that H1 disordered domains restrict the dynamics and conformation of both globular H1 and linker DNA arms, resulting in a more compact and rigid chromatosome particle. Furthermore, we identified regions of H1 disordered domains that are tightly tethered to DNA near the entry-exit site. Overall, our study elucidates at near-atomic resolution the way the disordered linker histone H1 modulates nucleosome's structural preferences and conformational 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.

Distinct Structures and Dynamics of Chromatosomes with Different Human Linker Histone Isoforms

The repeating structural unit of metazoan chromatin is the chromatosome, a nucleosome bound to a linker histone, H1. There are 11 human H1 isoforms with diverse cellular functions, but how they interact with the nucleosome remains elusive. Here, we determined the cryoelectron microscopy (cryo-EM) structures of chromatosomes containing 197 bp DNA and three different human H1 isoforms, respectively. The globular domains of all three H1 isoforms bound to the nucleosome dyad. However, the flanking/linker DNAs displayed substantial distinct dynamic conformations. Nuclear magnetic resonance (NMR) and H1 tail-swapping cryo-EM experiments revealed that the C-terminal tails of the H1 isoforms mainly controlled the flanking DNA orientations. We also observed partial ordering of the core histone H2A C-terminal and H3 N-terminal tails in the chromatosomes. Our results provide insights into the structures and dynamics of the chromatosomes and have implications for the structure and function of chromatin.

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.