Chromatosome Structure and Dynamics from Molecular Simulations

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.

PARP3 Affects Nucleosome Compaction Regulation

Genome compaction is one of the important subject areas for understanding the mechanisms regulating genes' expression and DNA replication and repair. The basic unit of DNA compaction in the eukaryotic cell is the nucleosome. The main chromatin proteins responsible for DNA compaction have already been identified, but the regulation of chromatin architecture is still extensively studied. Several authors have shown an interaction of ARTD proteins with nucleosomes and proposed that there are changes in the nucleosomes' structure as a result. In the ARTD family, only PARP1, PARP2, and PARP3 participate in the DNA damage response. Damaged DNA stimulates activation of these PARPs, which use NAD⁺ as a substrate. DNA repair and chromatin compaction need precise regulation with close coordination between them. In this work, we studied the interactions of these three PARPs with nucleosomes by atomic force microscopy, which is a powerful method allowing for direct measurements of geometric characteristics of single molecules. Using this method, we evaluated perturbations in the structure of single nucleosomes after the binding of a PARP. We demonstrated here that PARP3 significantly alters the geometry of nucleosomes, possibly indicating a new function of PARP3 in chromatin compaction regulation.

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.

DNA Sequence-Dependent Binding of Linker Histone gH1 Regulates Nucleosome Conformations

Sequence-dependent binding between DNA and proteins in chromatin is an essential part of gene expression. Linker histone H1 is an important protein in the regulation of chromatin compartmentalization and compaction, and its binding with the nucleosome is sensitive to the DNA sequence. Although the interactions of H1 and DNA have been widely investigated, the mechanism of nucleosome conformation changes induced by the DNA-sequence-dependent binding with gH1 (globular H1.0) remains largely unclear at the atomic level. In the present molecular dynamics simulations, both linker and dyad DNAs were mutated to investigate the conformational changes of the nucleosome induced by the sequence-dependent binding of gH1 based on the on-dyad binding mode. Our results indicate that gH1 is insensitive to the DNA sequence of the dyad DNA but presents an apparent preference to linker DNA with an AT-rich sequence. Moreover, this specific binding induces the entry/exit region of a nucleosome to a tight conformation and regulates the accessibility of core histones. Considering that the entry/exit region of the nucleosome is a crucial binding site for many functional proteins related to gene expression, the conformational change at this region could represent an important gene regulation signal.

Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations

Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Release of linker histone from the nucleosome driven by polyelectrolyte competition with a disordered protein

Highly charged intrinsically disordered proteins are essential regulators of chromatin structure and transcriptional activity. Here we identify a surprising mechanism of molecular competition that relies on the pronounced dynamical disorder present in these polyelectrolytes and their complexes. The highly positively charged human linker histone H1.0 (H1) binds to nucleosomes with ultrahigh affinity, implying residence times incompatible with efficient biological regulation. However, we show that the disordered regions of H1 retain their large-amplitude dynamics when bound to the nucleosome, which enables the highly negatively charged and disordered histone chaperone prothymosin α to efficiently invade the H1-nucleosome complex and displace H1 via a competitive substitution mechanism, vastly accelerating H1 dissociation. By integrating experiments and simulations, we establish a molecular model that rationalizes the remarkable kinetics of this process structurally and dynamically. Given the abundance of polyelectrolyte sequences in the nuclear proteome, this mechanism is likely to be widespread in cellular regulation.

Pharmaceutical Applications of Quantum Dots

Nanotechnology has been utilized in developing novel drug formulations with minimal adverse effects. Nanoparticles in a lower size range with great surface area, increased potency, and easy permeability could be an approach for the treatment of cancer and other diseases. Unlike other nanoparticles, quantum dots have specific functional groups, have charges over their surface, and are extremely small in size (2-10nm), which makes them more permeable through tight junctions. Quantum dots are interesting materials that offer diagnosis and treatment concurrently. Quantum dots are reported to have several applications in pharmaceuticals as well as drug delivery, diagnosis, immunolabeling, and cell labeling tools. However, the existence of heavy metals in quantum dots such as cadmium poses a potential challenge for future medical applications, where quantum dots may be deliberately injected into the body. In this review, we are focusing on various pharmaceutical applications of quantum dots. Graphical Abstract.

The missing linker: emerging trends for H1 variant-specific functions

In this review, Prendergast and Reinberg discuss the likelihood that the family of histone H1 variants may be key to understanding several fundamental processes in chromatin biology and underscore their particular contributions to distinctly significant chromatin-related processes.

Major advances in the chromatin and epigenetics fields have uncovered the importance of core histones, histone variants and their post-translational modifications (PTMs) in modulating chromatin structure. However, an acutely understudied related feature of chromatin structure is the role of linker histone H1. Previous assumptions of the functional redundancy of the 11 nonallelic H1 variants are contrasted by their strong evolutionary conservation, variability in their potential PTMs, and increased reports of their disparate functions, sub-nuclear localizations and unique expression patterns in different cell types. The commonly accepted notion that histone H1 functions solely in chromatin compaction and transcription repression is now being challenged by work from multiple groups. These studies highlight histone H1 variants as underappreciated facets of chromatin dynamics that function independently in various chromatin-based processes. In this review, we present notable findings involving the individual somatic H1 variants of which there are seven, underscoring their particular contributions to distinctly significant chromatin-related processes.

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.

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 dynamics mediate DNA unwrapping and sliding in nucleosomes

Nucleosomes are elementary building blocks of chromatin in eukaryotes. They tightly wrap ∼147 DNA base pairs around an octamer of histone proteins. How nucleosome structural dynamics affect genome functioning is not completely clear. Here we report all-atom molecular dynamics simulations of nucleosome core particles at a timescale of 15 microseconds. At this timescale, functional modes of nucleosome dynamics such as spontaneous nucleosomal DNA breathing, unwrapping, twisting, and sliding were observed. We identified atomistic mechanisms of these processes by analyzing the accompanying structural rearrangements of the histone octamer and histone-DNA contacts. Octamer dynamics and plasticity were found to enable DNA unwrapping and sliding. Through multi-scale modeling, we showed that nucleosomal DNA dynamics contribute to significant conformational variability of the chromatin fiber at the supranucleosomal level. Our study further supports mechanistic coupling between fine details of histone dynamics and chromatin functioning, provides a framework for understanding the effects of various chromatin modifications.

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.

CENP-A nucleosome-a chromatin-embedded pedestal for the centromere: lessons learned from structural biology

The centromere is a chromosome locus that directs equal segregation of chromosomes during cell division. A nucleosome containing the histone H3 variant CENP-A epigenetically defines the centromere. Here, we summarize findings from recent structural biology studies, including several CryoEM structures, that contributed to elucidate specific features of the CENP-A nucleosome and molecular determinants of its interactions with CENP-C and CENP-N, the only two centromere proteins that directly bind to it. Based on those findings, we propose a role of the CENP-A nucleosome in the organization of centromeric chromatin beyond binding centromeric proteins.

Linker histone defines structure and self-association behaviour of the 177 bp human chromatosome

Linker histones play essential roles in the regulation and maintenance of the dynamic chromatin structure of higher eukaryotes. The influence of human histone H1.0 on the nucleosome structure and biophysical properties of the resulting chromatosome were investigated and compared with the 177-bp nucleosome using Cryo-EM and SAXS. The 4.5 Å Cryo-EM chromatosome structure showed that the linker histone binds at the nucleosome dyad interacting with both linker DNA arms but in a tilted manner leaning towards one of the linker sides. The chromatosome is laterally compacted and rigid in the dyad and linker DNA area, in comparison with the nucleosome where linker DNA region is more flexible and displays structural variability. In solution, the chromatosomes appear slightly larger than the nucleosomes, with the volume increase compared to the bound linker histone, according to solution SAXS measurements. SAXS X-ray diffraction characterisation of Mg-precipitated samples showed that the different shapes of the 177 chromatosome enabled the formation of a highly ordered lamello-columnar phase when precipitated by Mg²⁺, indicating the influence of linker histone on the nucleosome stacking. The biological significance of linker histone, therefore, may be affected by the change in the polyelectrolyte and DNA conformation properties of the chromatosomes, in comparison to nucleosomes.