Partially Assembled Nucleosome Structures at Atomic Detail

Angle between DNA linker and nucleosome core particle regulates array compaction revealed by individual-particle cryo-electron tomography

The conformational dynamics of nucleosome arrays generate a diverse spectrum of microscopic states, posing challenges to their structural determination. Leveraging cryogenic electron tomography (cryo-ET), we determine the three-dimensional (3D) structures of individual mononucleosomes and arrays comprising di-, tri-, and tetranucleosomes. By slowing the rate of condensation through a reduction in ionic strength, we probe the intra-array structural transitions that precede inter-array interactions and liquid droplet formation. Under these conditions, the arrays exhibite irregular zig-zag conformations with loose packing. Increasing the ionic strength promoted intra-array compaction, yet we do not observe the previously reported regular 30-nanometer fibers. Interestingly, the presence of H1 do not induce array compaction; instead, one-third of the arrays display nucleosomes invaded by foreign DNA, suggesting an alternative role for H1 in chromatin network construction. We also find that the crucial parameter determining the structure adopted by chromatin arrays is the angle between the entry and exit of the DNA and the corresponding tangents to the nucleosomal disc. Our results provide insights into the initial stages of intra-array compaction, a critical precursor to condensation in the regulation of chromatin organization.

Sequence Dependence in Nucleosome Dynamics

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Real-Time Multistep Asymmetrical Disassembly of Nucleosomes and Chromatosomes Visualized by High-Speed Atomic Force Microscopy

During replication, expression, and repair of the eukaryotic genome, cellular machinery must access the DNA wrapped around histone proteins forming nucleosomes. These octameric protein·DNA complexes are modular, dynamic, and flexible and unwrap or disassemble either spontaneously or by the action of molecular motors. Thus, the mechanism of formation and regulation of subnucleosomal intermediates has gained attention genome-wide because it controls DNA accessibility. Here, we imaged nucleosomes and their more compacted structure with the linker histone H1 (chromatosomes) using high-speed atomic force microscopy to visualize simultaneously the changes in the DNA and the histone core during their disassembly when deposited on mica. Furthermore, we trained a neural network and developed an automatic algorithm to track molecular structural changes in real time. Our results show that nucleosome disassembly is a sequential process involving asymmetrical stepwise dimer ejection events. The presence of H1 restricts DNA unwrapping, significantly increases the nucleosomal lifetime, and affects the pathway in which heterodimer asymmetrical dissociation occurs. We observe that tetrasomes are resilient to disassembly and that the tetramer core (H3·H4)2 can diffuse along the nucleosome positioning sequence. Tetrasome mobility might be critical to the proper assembly of nucleosomes and can be relevant during nucleosomal transcription, as tetrasomes survive RNA polymerase passage. These findings are relevant to understanding nucleosome intrinsic dynamics and their modification by DNA-processing enzymes.

Spatiotemporal kinetics of CAF-1-dependent chromatin maturation ensures transcription fidelity during S-phase

Proper maintenance of epigenetic information after replication is dependent on the rapid assembly and maturation of chromatin. Chromatin Assembly Complex 1 (CAF-1) is a conserved histone chaperone that deposits (H3-H4)2 tetramers as part of the replication-dependent chromatin assembly process. Loss of CAF-1 leads to a delay in chromatin maturation, albeit with minimal impact on steady-state chromatin structure. However, the mechanisms by which CAF-1 mediates the deposition of (H3-H4)2 tetramers and the phenotypic consequences of CAF-1-associated assembly defects are not well understood. We used nascent chromatin occupancy profiling to track the spatiotemporal kinetics of chromatin maturation in both wild-type (WT) and CAF-1 mutant yeast cells. Our results show that loss of CAF-1 leads to a heterogeneous rate of nucleosome assembly, with some nucleosomes maturing at near WT kinetics and others showing significantly slower maturation kinetics. The slow-to-mature nucleosomes are enriched in intergenic and poorly transcribed regions, suggesting that transcription-dependent assembly mechanisms can reset the slow-to-mature nucleosomes following replication. Nucleosomes with slow maturation kinetics are also associated with poly(dA:dT) sequences, which implies that CAF-1 deposits histones in a manner that counteracts resistance from the inflexible DNA sequence, promoting the formation of histone octamers as well as ordered nucleosome arrays. In addition, we show that the delay in chromatin maturation is accompanied by a transient and S-phase-specific loss of gene silencing and transcriptional regulation, revealing that the DNA replication program can directly shape the chromatin landscape and modulate gene expression through the process of chromatin maturation.

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.

Spatiotemporal kinetics of CAF-1-dependent chromatin maturation ensures transcription fidelity during S-phase

Proper maintenance of epigenetic information after replication is dependent on the rapid assembly and maturation of chromatin. Chromatin Assembly Complex 1 (CAF-1) is a conserved histone chaperone that deposits (H3-H4)2 tetramers as part of the replication-dependent chromatin assembly process. Loss of CAF-1 leads to a delay in chromatin maturation, albeit with minimal impact on steady-state chromatin structure. However, the mechanisms by which CAF-1 mediates the deposition of (H3-H4)2 tetramers and the phenotypic consequences of CAF-1-associated assembly defects are not well understood. We used nascent chromatin occupancy profiling to track the spatiotemporal kinetics of chromatin maturation in both wild-type (WT) and CAF-1 mutant yeast cells. Our results show that loss of CAF-1 leads to a heterogeneous rate of nucleosome assembly, with some nucleosomes maturing at near WT kinetics and others exhibiting significantly slower maturation kinetics. The slow-to-mature nucleosomes are enriched in intergenic and poorly transcribed regions, suggesting that transcription-dependent assembly mechanisms can reset the slow-to-mature nucleosomes following replication. Nucleosomes with slow maturation kinetics are also associated with poly(dA:dT) sequences, which implies that CAF-1 deposits histones in a manner that counteracts resistance from the inflexible DNA sequence, promoting the formation of histone octamers as well as ordered nucleosome arrays. In addition, we demonstrate that the delay in chromatin maturation is accompanied by a transient and S-phase specific loss of gene silencing and transcriptional regulation, revealing that the DNA replication program can directly shape the chromatin landscape and modulate gene expression through the process of chromatin maturation.

Conformational Change of Nucleosome Arrays prior to Phase Separation

Chromatin phase transition serves as a regulatory mechanism for eukaryotic transcription. Understanding this process requires the characterization of the nucleosome array structure in response to external stimuli prior to phase separation. However, the intrinsic flexibility and heterogeneity hinders the arrays' structure determination. Here we exploit advances in cryogenic electron tomography (cryo-ET) to determine the three-dimensional (3D) structure of each individual particle of mono-, di-, tri-, and tetranucleosome arrays. Statistical analysis reveals the ionic strength changes the angle between the DNA linker and nucleosome core particle (NCP), which regulate the overall morphology of nucleosome arrays. The finding that one-third of the arrays in the presence of H1 contain an NCP invaded by foreign DNA suggests an alternative function of H1 in constructing nucleosomal networks. The new insights into the nucleosome conformational changes prior to the intermolecular interaction stage extends our understanding of chromatin phase separation regulation.

Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome

Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG^0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H ⁺ activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG^0.3 is between 3 and 4 k _B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG^0.3 ∼ 10k _B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG^0.3 = 2.2k _B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG^0.3 = 5.2k _B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.

Nanopore Sensing of DNA-Histone Complexes on Nucleosome Arrays

The location of nucleosomes in DNA and their structural stability are critical in regulating DNA compaction, site accessibility, and epigenetic gene regulation. Here, we combine the nanopore platform-based fast and label-free single-molecule detection technique with a voltage-dependent force rupture assay to detect distinct structures on nucleosomal arrays and then to induce breakdown of individual nucleosome complexes. Specifically, we demonstrate direct measurement of distinct nucleosome structures present on individual 12-mer arrays. A detailed event analysis showed that nucleosomes are present as a combination of complete and partial structures, during translocation through the pore. By comparing with the voltage-dependent translocation of the mononucleosomes, we find that the partial nucleosomes result from voltage-dependent structural disintegration of nucleosomes. High signal-to-noise detection of heterogeneous levels in translocation of 12-mer array molecules quantifies the heterogeneity and nucleosomal substructure sizes on the arrays. These results facilitate the understanding of electrostatic interactions responsible for the integrity of the nucleosome structure and possible mechanisms of its unraveling by chromatin remodeling enzymes. This study also has potential applications in chromatin profiling.

H2A-H2B Histone Dimer Plasticity and Its Functional Implications

The protein core of the nucleosome is composed of an H3-H4 histone tetramer and two H2A-H2B histone dimers. The tetramer organizes the central 60 DNA bp, while H2A-H2B dimers lock the flanking DNA segments. Being positioned at the sides of the nucleosome, H2A-H2B dimers stabilize the overall structure of the nucleosome and modulate its dynamics, such as DNA unwrapping, sliding, etc. Such modulation at the epigenetic level is achieved through post-translational modifications and the incorporation of histone variants. However, the detailed connection between the sequence of H2A-H2B histones and their structure, dynamics and implications for nucleosome functioning remains elusive. In this work, we present a detailed study of H2A-H2B dimer dynamics in the free form and in the context of nucleosomes via atomistic molecular dynamics simulations (based on X. laevis histones). We supplement simulation results by comparative analysis of information in the structural databases. Particularly, we describe a major dynamical mode corresponding to the bending movement of the longest H2A and H2B α-helices. This overall bending dynamics of the H2A-H2B dimer were found to be modulated by its interactions with DNA, H3-H4 tetramer, the presence of DNA twist-defects with nucleosomal DNA and the amino acid sequence of histones. Taken together, our results shed new light on the dynamical mechanisms of nucleosome functioning, such as nucleosome sliding, DNA-unwrapping and their epigenetic modulation.

Synthesis of Oriented Hexasomes and Asymmetric Nucleosomes Using a Template Editing Process

Nucleosomes, the structural building blocks of chromatin, possess 2-fold pseudo symmetry which can be broken through differential modification or removal of one copy of a pair of sister histones. The resultant asymmetric nucleosomes and hexasomes have been implicated in gene regulation, yet the use of these noncanonical substrates in chromatin biochemistry is limited, owing to the lack of efficient methods for their preparation. Here, we report a strategy that allows the orientation of these asymmetric species to be tightly controlled relative to the underlying DNA sequence. Our approach is based on the use of truncated DNA templates to assemble oriented hexasomes followed by DNA ligation and, in the case of asymmetric nucleosomes, addition of the missing heterotypic histones. We show that this approach is compatible with multiple nucleosome positioning sequences, allowing the generation of desymmetrized mononucleosomes and oligonucleosomes with varied DNA overhangs and heterotypic histone H2A/H2B dimer compositions. Using this technology, we examine the functional consequences of asymmetry on BRG1/BRM associated factor (BAF) complex-mediated chromatin remodeling. Our results indicate that cancer-associated histone mutations can reprogram the inherent activity of BAF chromatin remodeling to induce aberrant chromatin structure.

Histone chaperone-mediated co-expression assembly of tetrasomes and nucleosomes

The nucleosome, a basic unit of chromatin found in all eukaryotes, is thought to be assembled through the orchestrated activity of several histone chaperones and chromatin assembly factors in a stepwise manner, proceeding from tetrasome assembly, to H2A/H2B deposition, and finally to formation of the mature nucleosome. In this study, we demonstrate chaperone-mediated assembly of both tetrasomes and nucleosomes on the well-defined Widom 601 positioning sequence using a co-expression/reconstitution wheat germ cell-free system. The purified tetrasomes and nucleosomes were positioned around the center of a given sequence. The heights and diameters were measured by atomic force microscopy. Together with the reported unmodified native histones produced by the wheat germ cell-free platform, our method is expected to be useful for downstream applications in the field of chromatin research.

Personalization of medical treatments in oncology: time for rethinking the disease concept to improve individual outcomes

The agenda of pharmacology discovery in the field of personalized oncology was dictated by the search of molecular targets assumed to deterministically drive tumor development. In this perspective, genes play a fundamental "causal" role while cells simply act as causal proxies, i.e., an intermediate between the molecular input and the organismal output. However, the ceaseless genomic change occurring across time within the same primary and metastatic tumor has broken the hope of a personalized treatment based only upon genomic fingerprint. Indeed, current models are unable in capturing the unfathomable complexity behind the outbreak of a disease, as they discard the contribution of non-genetic factors, environment constraints, and the interplay among different tiers of organization. Herein, we posit that a comprehensive personalized model should view at the disease as a "historical" process, in which different spatially and timely distributed factors interact with each other across multiple levels of organization, which collectively interact with a dynamic gene-expression pattern. Given that a disease is a dynamic, non-linear process - and not a static-stable condition - treatments should be tailored according to the "timing-frame" of each condition. This approach can help in detecting those critical transitions through which the system can access different attractors leading ultimately to diverse outcomes - from a pre-disease state to an overt illness or, alternatively, to recovery. Identification of such tipping points can substantiate the predictive and the preventive ambition of the Predictive, Preventive and Personalized Medicine (PPPM/3PM). However, an unusual effort is required to conjugate multi-omics approaches, data collection, and network analysis reconstruction (eventually involving innovative Artificial Intelligent tools) to recognize the critical phases and the relevant targets, which could help in patient stratification and therapy personalization.

The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study

The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences.

Nucleosomes are biomolecular complexes formed by DNA wrapped around histone proteins. They represent the basic units of Eukaryotic chromosomes, compacting the genome so that it fits into the small nucleus, and regulating important biological processes such as gene expression. Nucleosomes are disassembled during disruptive events such as DNA replication, and re-assembled afterwards to preserve the correct organization of chromatin. However, the molecular details of nucleosome assembly are still not well understood. In particular, experiments found that histones and DNA may associate into a variety of non-canonical complexes, but their precise conformation and role during assembly remain unclear. In this study, we addressed these problems by performing extensive molecular dynamics simulations of nucleosomes undergoing assembly and disassembly. The simulations reveal many insights into the kinetics of assembly, the structure of non-canonical nucleosome intermediates, and the influence of salt concentration and DNA sequence on the assembly process.

Nucleosome composition regulates the histone H3 tail conformational ensemble and accessibility

Hexasomes and tetrasomes are intermediates in nucleosome assembly and disassembly. Their formation is promoted by histone chaperones, ATP-dependent remodelers, and RNA polymerase II. In addition, hexasomes are maintained in transcribed genes and could be an important regulatory factor. While nucleosome composition has been shown to affect the structure and accessibility of DNA, its influence on histone tails is largely unknown. Here, we investigate the conformational dynamics of the H3 tail in the hexasome and tetrasome. Using a combination of NMR spectroscopy, MD simulations, and trypsin proteolysis, we find that the conformational ensemble of the H3 tail is regulated by nucleosome composition. As has been found for the nucleosome, the H3 tails bind robustly to DNA within the hexasome and tetrasome, but upon loss of the H2A/H2B dimer, we determined that the adjacent H3 tail has an altered conformational ensemble, increase in dynamics, and increase in accessibility. Similar to observations of DNA dynamics, this is seen to be asymmetric in the hexasome. Our results indicate that nucleosome composition has the potential to regulate chromatin signaling and ultimately help shape the chromatin landscape.

Multiple roles of H2A.Z in regulating promoter chromatin architecture in human cells

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

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.

ATRX limits the accessibility of histone H3-occupied HSV genomes during lytic infection

Histones are rapidly loaded on the HSV genome upon entry into the nucleus of human fibroblasts, but the effects of histone loading on viral replication have not been fully defined. We showed recently that ATRX is dispensable for de novo deposition of H3 to HSV genomes after nuclear entry but restricted infection through maintenance of viral heterochromatin. To further investigate the roles that ATRX and other histone H3 chaperones play in restriction of HSV, we infected human fibroblasts that were systematically depleted of nuclear H3 chaperones. We found that the ATRX/DAXX complex is unique among nuclear H3 chaperones in its capacity to restrict ICP0-null HSV infection. Only depletion of ATRX significantly alleviated restriction of viral replication. Interestingly, no individual nuclear H3 chaperone was required for deposition of H3 onto input viral genomes, suggesting that during lytic infection, H3 deposition may occur through multiple pathways. ChIP-seq for total histone H3 in control and ATRX-KO cells infected with ICP0-null HSV showed that HSV DNA is loaded with high levels of histones across the entire viral genome. Despite high levels of H3, ATAC-seq analysis revealed that HSV DNA is highly accessible, especially in regions of high GC content, and is not organized largely into ordered nucleosomes during lytic infection. ATRX reduced accessibility of viral DNA to the activity of a TN5 transposase and enhanced accumulation of viral DNA fragment sizes associated with nucleosome-like structures. Together, these findings support a model in which ATRX restricts viral infection by altering the structure of histone H3-loaded viral chromatin that reduces viral DNA accessibility for transcription. High GC rich regions of the HSV genome, especially the S component inverted repeats of the HSV-1 genome, show increased accessibility, which may lead to increased ability to transcribe the IE genes encoded in these regions during initiation of infection.

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.

The role of histone tails in nucleosome stability: An electrostatic perspective

We propose a methodology for the study of protein-DNA electrostatic interactions and apply it to clarify the effect of histone tails in nucleosomes. This method can be used to correlate electrostatic interactions to structural and functional features of protein-DNA systems, and can be combined with coarse-grained representations. In particular, we focus on the electrostatic field and resulting forces acting on the DNA. We investigate the electrostatic origins of effects such as different stages in DNA unwrapping, nucleosome destabilization upon histone tail truncation, and the role of specific arginines and lysines undergoing Post-Translational Modifications. We find that the positioning of the histone tails can oppose the attractive pull of the histone core, locally deform the DNA, and tune DNA unwrapping. Small conformational variations in the often overlooked H2A C-terminal tails had significant electrostatic repercussions near the DNA entry and exit sites. The H2A N-terminal tail exerts attractive electrostatic forces towards the histone core in positions where Polymerase II halts its progress. We validate our results with comparisons to previous experimental and computational observations.

Effects of Psoralen on Histone-DNA Interactions Studied by Using Atomic Force Microscopy

The investigation of the DNA-histone interactions and factors that affect such interactions in the nucleosome is essential for understanding the role of chromatin organization in all cellular processes involved in the repair, transcription, and replication of the eukaryotic genome. As a kind of photosensitive molecule, psoralen (PSO) is used in the treatment of skin disease with ultraviolet light (PSO and ultra violet light, type A). The effect of treatment is remarkable, but the side effect is also obvious. PSO can be embedded in a 5' TA sequence in double-stranded DNA (dsDNA), and dsDNA is mainly wrapped around a histone octamer to form a nucleosome structure in human cells. Therefore, it is very necessary to explore the influence of PSO on DNA-histone interactions. To this end, the binding specificity and mode of DNA and histone in the presence or absence of PSO are investigated systematically. The results show that the presence of PSO (no matter if there is ultra violet light treatment) can increase the overall probability of histone binding to dsDNA while lowering the selectivity of histone binding to the specific DNA sequence in vitro. In addition, the increase of solution ionic strength can lower the ratio of histone binding to nonspecific DNA.

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

Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.

"Direct" and "Indirect" Effects of Histone Modifications: Modulation of Sterical Bulk as a Novel Source of Functionality

The chromatin-regulatory principles of histone post-translational modifications (PTMs) are discussed with a focus on the potential alterations in chromatin functional state due to steric and mechanical constraints imposed by bulky histone modifications such as ubiquitin and SUMO. In the classical view, PTMs operate as recruitment platforms for histone "readers," and as determinants of chromatin array compaction. Alterations of histone charges by "small" chemical modifications (e.g., acetylation, phosphorylation) could regulate nucleosome spontaneous dynamics without globally affecting nucleosome structure. These fluctuations in nucleosome wrapping can be exploited by chromatin-processing machinery. In contrast, ubiquitin and SUMO are comparable in size to histones, and it seems logical that these PTMs could conflict with canonical nucleosome organization. An experimentally testable hypothesis that by adding sterical bulk these PTMs can robustly alter nucleosome primary structure is proposed. The model presented here stresses the diversity of mechanisms by which histone PTMs regulate chromatin dynamics, primary structure and, hence, functionality.

Histone Tail Dynamics in Partially Disassembled Nucleosomes During Chromatin Remodeling

Nucleosomes are structural units of the chromosome consisting of DNA wrapped around histone proteins, and play important roles in compaction and regulation of the chromatin structure. While the structure and dynamics of canonical nucleosomes have been studied extensively, those of nucleosomes in intermediate states, that occur when their structure or positioning is modulated, have been less understood. In particular, the dynamic features of partially disassembled nucleosomes have not been discussed in previous studies. Using all-atom molecular dynamics simulations, in this study, we investigated the dynamics and stability of nucleosome structures lacking a histone-dimer. DNA in nucleosomes lacking a histone H2A/H2B dimer was drastically deformed due to loss of local interactions between DNA and histones. In contrast, conformation of DNA in nucleosomes lacking H3/H4 was similar to the canonical nucleosome, as the H2A C-terminal domain infiltrated the space originally occupied by the dissociated H3/H4 histones and restricted DNA dynamics in close proximity. Our results suggest that, besides histone chaperones, the intrinsic dynamics of nucleosomes support the exchange of H2A/H2B, which is significantly more frequent than that of H3/H4.

The intrinsic stability of H2B-ubiquitylated nucleosomes and their in vitro assembly/disassembly by histone chaperone NAP1

BACKGROUND

Apart the gene-regulatory functions as docking sites for histone 'readers', some histone modifications could directly affect nucleosome structure. The H2BK34-ubiquitylation deposited by MOF-MSL complex, increases nucleosome dynamics in vitro and promotes donation of one H2A/H2B dimer to histone acceptors.

METHODS

We evaluated temperature-depended stability of H2BK34-ubiquitylated nucleosomes under 'physiological' ionic conditions in the presence or absence of histone acceptor, and examined assembly and disassembly of ubiquitylated nucleosomes in vitro by recombinant mouse NAP1.

RESULTS

H2BK34ub modification is sufficient to promote selective eviction of only one H2A/H2B dimer independently of histone-binding agents. Despite the robust H2A/H2B dimer-displacement effect of mNAP1 with the H2BK34ub (but not unmodified) nucleosomes, NAP1 could assemble symmetrically- or asymmetrically ubiquitylated nucleosomes under 'physiological' conditions in vitro.

CONCLUSIONS AND GENERAL SIGNIFICANCE

The increased mobility of one nucleosomal H2A/H2B dimer is an intrinsic nucleosome destabilizing property of H2BK34 ubiquitylation that has the intranucleosome bases. The ability of NAP to reasonably efficiently assemble H2BK34-ubiquitylated nucleosomes supposes a potential mechanism for deposition/distribution of H2BK34ub mark in the MOF-MSL independent manner (for example, during histone dimer exchange upon transcription elongation).

DNA sequence influences hexasome orientation to regulate DNA accessibility

Nucleosomes, the fundamental organizing units of eukaryotic genomes, contain ∼146 base pairs of DNA wrapped around a histone H3–H4 tetramer and two histone H2A–H2B dimers. Converting nucleosomes into hexasomes by removal of a H2A–H2B dimer is an important regulatory event, but its regulation and functional consequences are not well-understood. To investigate the influence of hexasomes on DNA accessibility, we used the property of the Widom-601 Nucleosome Positioning Sequence (NPS) to form homogeneously oriented hexasomes in vitro. We find that DNA accessibility to transcription factors (TF) on the hexasome H2A–H2B distal side is identical to naked DNA, while the accessibility on the H2A–H2B proximal side is reduced by 2-fold, which is due to a 2-fold reduction in hexasome unwrapping probability. We then determined that a 23 bp region of the Widom-601 NPS is responsible for forming homogeneously oriented hexasomes. Analysis of published ChIP-exo data of hexasome containing genes identified two DNA sequence motifs that correlate with hexasome orientation in vivo, while ExoIII mapping studies of these sequences revealed they generate homogeneously oriented hexasomes in vitro. These results indicate that hexasome orientation, which is influenced by the underlying DNA sequence in vivo, is important for modulating DNA accessibility to regulate transcription.

Quantitative Modeling of Nucleosome Unwrapping from Both Ends

In eukaryotic cells, DNA is packaged into chromatin where nucleosomes are the basic packaging unit. Important cellular processes including gene expression, DNA replication, and DNA repair require nucleosomal DNA to be unwrapped so that functional proteins can access their target sites, which otherwise are sterically occluded. A key question in this process is what the unwrapped conformations individual nucleosomes adopt within chromatin are. Here, we develop a concurrent nucleosome unwrapping model to address this question. We hypothesize that for a given end-to-end distance of the nucleosomal DNA, the nucleosomal DNA stochastically unwraps from the histone core from both ends independently and that this combination of unwrapping from both sides results in a significant increase in the average distance between the DNA extending from both sides of the nucleosomes. We test our model on recently published experiments using a DNA origami nanocaliper that quantifies nucleosome unwrapping and achieve good agreement between experiment and model prediction. We then investigate the DNA origami caliper distribution when attached to a hexasome (a nucleosome lacking an H2A/H2B dimer). A significant shift in the caliper angle distribution caused by the asymmetric structural features of the hexasome seen experimentally is consistent with the model. Our modeling approach may be more broadly useful to the interpretation of other studies of nucleosome dynamics, chromatin dynamics, and regulatory processes involving nucleosome unwrapping, as well as more generally to optimization of future DNA origami designs to probe mechanical properties of biomolecules.

DNA accessibility of chromatosomes quantified by automated image analysis of AFM data

DNA compaction and accessibility in eukaryotes are governed by nucleosomes and orchestrated through interactions between DNA and DNA-binding proteins. Using QuantAFM, a method for automated image analysis of atomic force microscopy (AFM) data, we performed a detailed statistical analysis of structural properties of mono-nucleosomes. QuantAFM allows fast analysis of AFM images, including image preprocessing, object segmentation, and quantification of different structural parameters to assess DNA accessibility of nucleosomes. A comparison of nucleosomes reconstituted with and without linker histone H1 quantified H1's already described ability of compacting the nucleosome. We further employed nucleosomes bearing two charge-modifying mutations at position R81 and R88 in histone H2A (H2A R81E/R88E) to characterize DNA accessibility under destabilizing conditions. Upon H2A mutation, even in presence of H1, the DNA opening angle at the entry/exit site was increased and the DNA wrapping length around the histone core was reduced. Interestingly, a distinct opening of the less bendable DNA side was observed upon H2A mutation, indicating an enhancement of the intrinsic asymmetry of the Widom-601 nucleosomes. This study validates AFM as a technique to investigate structural parameters of nucleosomes and highlights how the DNA sequence, together with nucleosome modifications, can influence the DNA accessibility.

Effects of histone H2B ubiquitylation on the nucleosome structure and dynamics

DNA in nucleosomes has restricted nucleosome dynamics and is refractory to DNA-templated processes. Histone post-translational modifications play important roles in regulating DNA accessibility in nucleosomes. Whereas most histone modifications function either by mitigating the electrostatic shielding of histone tails or by recruiting 'reader' proteins, we show that ubiquitylation of H2B K34, which is located in a tight space protected by two coils of DNA superhelix, is able to directly influence the canonical nucleosome conformation via steric hindrances by ubiquitin groups. H2B K34 ubiquitylation significantly enhances nucleosome dynamics and promotes generation of hexasomes both with symmetrically or asymmetrically modified nucleosomes. Our results indicate a direct mechanism by which a histone modification regulates the chromatin structural states.

Generalized Born Implicit Solvent Models for Biomolecules

It would often be useful in computer simulations to use an implicit description of solvation effects, instead of explicitly representing the individual solvent molecules. Continuum dielectric models often work well in describing the thermodynamic aspects of aqueous solvation and can be very efficient compared to the explicit treatment of the solvent. Here, we review a particular class of so-called fast implicit solvent models, generalized Born (GB) models, which are widely used for molecular dynamics (MD) simulations of proteins and nucleic acids. These approaches model hydration effects and provide solvent-dependent forces with efficiencies comparable to molecular-mechanics calculations on the solute alone; as such, they can be incorporated into MD or other conformational searching strategies in a straightforward manner. The foundations of the GB model are reviewed, followed by examples of newer, emerging models and examples of important applications. We discuss their strengths and weaknesses, both for fidelity to the underlying continuum model and for the ability to replace explicit consideration of solvent molecules in macromolecular simulations.

Molecular recognition of nucleosomes by binding partners

Nucleosomes represent the elementary units of chromatin packing and hubs in epigenetic signaling pathways. Across the chromatin and over the lifetime of the eukaryotic cell, nucleosomes experience a broad repertoire of alterations that affect their structure and binding with various chromatin factors. Dynamics of the histone core, nucleosomal and linker DNA, and intrinsic disorder of histone tails add further complexity to the nucleosome interaction landscape. In light of our understanding through the growing number of experimental and computational studies, we review the emerging patterns of molecular recognition of nucleosomes by their binding partners and assess the basic mechanisms of its regulation.

Analysis of histone ubiquitylation by MSL1/MSL2 proteins in vitro

Histone posttranslational modifications (PTM) control gene activity by targeting chromatin-regulatory proteins. By altering histone charges PTMs could also modulate inter- and intra-nucleosomal interactions, and thus affect chromatin high-order compaction and nucleosome stochastic folding, respectively. However, recently it has been shown that histone H2BK34- ubiquitylation (which is deposited in vivo by MOF-MSL) can destabilize one of the nucleosomal H2A-H2B dimers in symmetrically and (albeit to a lesser extend) asymmetrically modified nucleosomes, and thus promote formation of a hexasome particle. Here we have studied ubiquitylation patterns by purified MSL1/MSL2 using nucleosomes and different histone substrates. We have shown that H2B-ubiquitylation by MSL1/2 depends on substrate configuration. In addition, MSL1/2 efficiently ubiquitylate histone substrates but very poorly modify nucleosomes, which implies a requirement for nucleosome structural alteration for efficient ubiquitylation of H2BK34. Nucleosome modification by MSL1/MSL2 in vitro was analyzed directly using nucleosome gel-mobility shift assay, which suggested that MSL1/2 can deposit two ubiquitin moieties in one nucleosome.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.

Structural interpretation of DNA-protein hydroxyl-radical footprinting experiments with high resolution using HYDROID

Hydroxyl-radical footprinting (HRF) is a powerful method for probing structures of nucleic acid-protein complexes with single-nucleotide resolution in solution. To tap the full quantitative potential of HRF, we describe a protocol, hydroxyl-radical footprinting interpretation for DNA (HYDROID), to quantify HRF data and integrate them with atomistic structural models. The stages of the HYDROID protocol are extraction of the lane profiles from gel images, quantification of the DNA cleavage frequency at each nucleotide and theoretical estimation of the DNA cleavage frequency from atomistic structural models, followed by comparison of experimental and theoretical results. Example scripts for each step of HRF data analysis and interpretation are provided for several nucleosome systems; they can be easily adapted to analyze user data. As input, HYDROID requires polyacrylamide gel electrophoresis (PAGE) images of HRF products and optionally can use a molecular model of the DNA-protein complex. The HYDROID protocol can be used to quantify HRF over DNA regions of up to 100 nucleotides per gel image. In addition, it can be applied to the analysis of RNA-protein complexes and free RNA or DNA molecules in solution. Compared with other methods reported to date, HYDROID is unique in its ability to simultaneously integrate HRF data with the analysis of atomistic structural models. HYDROID is freely available. The complete protocol takes ~3 h. Users should be familiar with the command-line interface, the Python scripting language and Protein Data Bank (PDB) file formats. A graphical user interface (GUI) with basic functionality (HYDROID_GUI) is also available.

Effects of DNA Superhelical Stress on the Stability of H2B-Ubiquitylated Nucleosomes

On the nucleosome level, histone posttranslational modifications function mainly as the regulatory signals; in addition, some posttranslational modifications can enhance nucleosome stochastic folding, which is restricted in "canonic" nucleosomes. Recently, it has been shown in vitro that symmetric or asymmetric nucleosome ubiquitylation at H2BK34 (and H2BK120, to a lesser extent) can destabilize one of the nucleosomal H2A-H2B dimers and promote nucleosome conversion to a hexasome particle [Krajewski et al. (2018). Nucleic Acids Res., 46, 7631-7642]. Such lability of H2Bub nucleosomes raises a question of whether they could accommodate transient changes in DNA torsional tensions, which are generated by virtually any process that manipulates DNA strands. Using positively or negatively supercoiled DNA minicircles and homogeneously-modified H2Bub histones, we have found that DNA topology could strongly and selectively affect nucleosome stability depending on its ubiquitylation state (here the term "nucleosome stability" means the nucleosome property to maintain its structural integrity and dynamics characteristic to "canonic" nucleosomes). The results point to a role for H2B ubiquitylation in amplifying or mitigating the effects of a DNA torque on the nucleosome stability and dynamics.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core

BACKGROUND

Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects.

RESULTS

A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects.

CONCLUSION

Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

Hydroxyl-radical footprinting combined with molecular modeling identifies unique features of DNA conformation and nucleosome positioning

Nucleosomes are the most abundant protein–DNA complexes in eukaryotes that provide compaction of genomic DNA and are implicated in regulation of transcription, DNA replication and repair. The details of DNA positioning on the nucleosome and the DNA conformation can provide key regulatory signals. Hydroxyl-radical footprinting (HRF) of protein–DNA complexes is a chemical technique that probes nucleosome organization in solution with a high precision unattainable by other methods. In this work we propose an integrative modeling method for constructing high-resolution atomistic models of nucleosomes based on HRF experiments. Our method precisely identifies DNA positioning on nucleosome by combining HRF data for both DNA strands with the pseudo-symmetry constraints. We performed high-resolution HRF for Saccharomyces cerevisiae centromeric nucleosome of unknown structure and characterized it using our integrative modeling approach. Our model provides the basis for further understanding the cooperative engagement and interplay between Cse4p protein and the A-tracts important for centromere function.