Nucleosomes as allosteric scaffolds for genetic regulation

Molecular dynamics simulations for the study of chromatin biology

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

Nucleosome assembly and disassembly pathways in vitro

Structural fluctuations of nucleosomes modulate the access to internal DNA in eukaryotic cells; clearly characterisation of this fundamental process is crucial to understanding gene regulation. Here we apply PhAST (Photochemical Analysis of Structural Transitions) to monitor at a base pair level, structural alterations induced all along the DNA upon histone binding or release. By offering the first reliable, detailed comparison of nucleosome assembly and disassembly in vitro, we reveal similarities and differences between the two processes. We identify multiple, sequential intermediate states characterised by specific PhAST signals whose localisation and amplitude reflect asymmetries of DNA/histone interactions with respect to the nucleosome pseudo dyad. These asymmetries involve not only the DNA extremities but also regions close to the pseudo dyad. Localisations of asymmetries develop in a consistent manner during both assembly and disassembly processes; they primarily reflect the DNA sequence effect on the efficiency of DNA-histone binding. More unexpectedly, the amplitude component of PhAST signals not only evolves as a function of intermediate states but does so differently between assembly and disassembly pathways. Our observation of differences between assembly and disassembly opens up new avenues to define the role of the DNA sequence in processes underlying the regulation of gene expression. Overall, we provide new insights into how the intrinsic properties of DNA are integrated into a holistic mechanism that controls chromatin structure.

Quantitative Analysis of Dynamic Allostery

Dynamic allostery refers to one important class of allosteric regulation that does not involve noticeable conformational changes upon effector binding. In recent years, many "quasi"-dynamic allosteric proteins have been found to only experience subtle conformational changes during allosteric regulation. However, as enthalpic and entropic contributions are coupled to each other and even tiny conformational changes could bring in noticeable free energy changes, a quantitative description is essential to understand the contribution of pure dynamic allostery. Here, by developing a unified anisotropic elastic network model (uANM) considering both side-chain information and ligand heavy atoms, we quantitatively estimated the contribution of pure dynamic allostery in a dataset of known allosteric proteins by excluding the conformational changes upon ligand binding. We found that the contribution of pure dynamic allostery is generally small (much weaker than previously expected) and robustly exhibits an allosteric activation effect, which exponentially decays with the distance between the substrate and the allosteric ligand. We further constructed toy models to study the determinant factors of dynamic allostery in monomeric and oligomeric proteins using the uANM. Analysis of the toy models revealed that a short distance, a small angle between the two ligands, strong protein-ligand interactions, and weak protein internal interactions lead to strong dynamic allostery. Our study provides a quantitative estimation of pure dynamic allostery and facilitates the understanding of dynamic-allostery-controlled biological processes and the design of allosteric drugs and proteins.

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.

Binding of regulatory proteins to nucleosomes is modulated by dynamic histone tails

Little is known about the roles of histone tails in modulating nucleosomal DNA accessibility and its recognition by other macromolecules. Here we generate extensive atomic level conformational ensembles of histone tails in the context of the full nucleosome, totaling 65 microseconds of molecular dynamics simulations. We observe rapid conformational transitions between tail bound and unbound states, and characterize kinetic and thermodynamic properties of histone tail-DNA interactions. Different histone types exhibit distinct binding modes to specific DNA regions. Using a comprehensive set of experimental nucleosome complexes, we find that the majority of them target mutually exclusive regions with histone tails on nucleosomal/linker DNA around the super-helical locations ± 1, ± 2, and ± 7, and histone tails H3 and H4 contribute most to this process. These findings are explained within competitive binding and tail displacement models. Finally, we demonstrate the crosstalk between different histone tail post-translational modifications and mutations; those which change charge, suppress tail-DNA interactions and enhance histone tail dynamics and DNA accessibility.

The intrinsic disorder of histone tails poses challenges in their characterization. Here the authors apply extensive molecular dynamics simulations of the full nucleosome to show reversible binding to DNA with specific binding modes of different types of histone tails, where charge-altering modifications suppress tail-DNA interactions and may boost interactions between nucleosomes and nucleosome-binding proteins.

Comparative distributions of RSBN1 and methylated histone H4 Lysine 20 in the mouse spermatogenesis

During spermatogenesis, nuclear architecture of male germ cells is dynamically changed and epigenetic modifications, in particular methylation of histones, highly contribute to its regulation as well as differentiation of male germ cells. Although several methyltransferases and demethylases for histone H3 are involved in the regulation of spermatogenesis, roles of either histone H4 lysine 20 (H4K20) methyltransferases or H4K20 demethylases during spermatogenesis still remain to be elucidated. Recently, RSBN1 which is a testis-specific gene expressed in round spermatids was identified as a demethylase for dimethyl H4K20. In this study, therefore, we confirm the demethylase function of RSBN1 and compare distributions between RSBN1 and methylated H4K20 in the seminiferous tubules. Unlike previous report, expression analyses for RSBN1 reveal that RSBN1 is not a testis-specific gene and is expressed not only in round spermatids but also in elongated spermatids. In addition, RSBN1 can demethylate not only dimethyl H4K20 but also trimethyl H4K20 and could convert both dimethyl H4K20 and trimethyl H4K20 into monomethyl H4K20. When distribution pattern of RSBN1 in the seminiferous tubule is compared to that of methylated H4K20, both dimethyl H4K20 and trimethyl H4K20 but not monomethyl H4K20 are disappeared from RSBN1 positive germ cells, suggesting that testis-specific distribution patterns of methylated H4K20 might be constructed by RSBN1. Thus, novel expression and function of RSBN1 could be useful to comprehend epigenetic regulation during spermatogenesis.

Interplay among transacting factors around promoter in the initial phases of transcription

The initiation signals are raised around the promoter by one of the general transcription factors, triggering a sequence of events that lead to mRNA transcript formation from target genes. Both specific noncoding DNA regions and transacting, macromolecular assemblies are intricately involved and indispensable. The transition between the subsequent transcriptional stages is accompanied by stage-specific signals and structural changes in the macromolecular assemblies and facilitated by the recruitment/removal of other chromatin and transcription-associated elements. Here, we discuss the choreography of transacting factors around promoter in the establishment and effectuation of the initial phases of transcription such as NDR formation, +1 nucleosome positioning, promoter DNA opening, and RNAPII promoter escape from a structural viewpoint.

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.

Dynamic networks observed in the nucleosome core particles couple the histone globular domains with DNA

The dynamics of eukaryotic nucleosomes are essential in gene activity and well regulated by various factors. Here, we elucidated the internal dynamics at multiple timescales for the human histones hH3 and hH4 in the Widom 601 nucleosome core particles (NCP), suggesting that four dynamic networks are formed by the residues exhibiting larger-scale μs-ms motions that extend from the NCP core to the histone tails and DNA. Furthermore, despite possessing highly conserved structural features, histones in the telomeric NCP exhibit enhanced μs-ms dynamics in the globular sites residing at the identified dynamic networks and in a neighboring region. In addition, higher mobility was observed for the N-terminal tails of hH3 and hH4 in the telomeric NCP. The results demonstrate the existence of dynamic networks in nucleosomes, through which the center of the core regions could interactively communicate with histone tails and DNA to potentially propagate epigenetic changes.

Shi et al. use solid-state nuclear magnetic resonance spectroscopy to reveal the internal dynamics of human histones hH3 and hH4 in the Widom 601 and the telomeric nucleosome core particles. This work has implications for the propagation of epigenetic changes via the center of the nucleosome core communicating with histone tails and DNA.

Human Histone Interaction Networks: An Old Concept, New Trends

To elucidate the properties of human histone interactions on the large scale, we perform a comprehensive mapping of human histone interaction networks by using data from structural, chemical cross-linking and various high-throughput studies. Histone interactomes derived from different data sources show limited overlap and complement each other. It inspires us to integrate these data into the combined histone global interaction network which includes 5308 proteins and 10,330 interactions. The analysis of topological properties of the human histone interactome reveals its scale free behavior and high modularity. Our study of histone binding interfaces uncovers a remarkably high number of residues involved in interactions between histones and non-histone proteins, 80-90% of residues in histones H3 and H4 have at least one binding partner. Two types of histone binding modes are detected: interfaces conserved in most histone variants and variant specific interfaces. Finally, different types of chromatin factors recognize histones in nucleosomes via distinct binding modes, and many of these interfaces utilize acidic patches among other sites. Interaction networks are available at https://github.com/Panchenko-Lab/Human-histone-interactome.

Nucleosome allostery in pioneer transcription factor binding

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POU_S domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.