Partial Unwrapping and Histone Tail Dynamics in Nucleosome Revealed by Coarse-Grained Molecular Simulations

DNA spontaneously wrapping around a histone core prefers negative supercoiling: A Brownian dynamics study

In eukaryotes, DNA achieves a highly compact structure primarily due to its winding around the histone cores. The nature wrapping of DNA around histone core form a 1.7 left-handed superhelical turns, contributing to negative supercoiling in chromatin. During transcription, negative supercoils generated behind the polymerase during transcription may play a role in triggering nucleosome reassembly. To elucidate how supercoils influence the dynamics of wrapping of DNA around the histone cores, we developed a novel model to simulate the intricate interplay between DNA and histone. Our simulations reveal that both positively and negatively supercoiled DNAs are capable of wrapping around histone cores to adopt the nucleosome conformation. Notably, our findings confirm a strong preference for negative supercoiled DNA during nucleosome wrapping, and reveal that the both of the negative writhe and twist are beneficial to the formation of the DNA wrapping around histone. Additionally, the simulations of the multiple nucleosomes on the same DNA template indicate that the nucleosome tends to assemble in proximity to the original nucleosome. This advancement in understanding the spontaneous formation of nucleosomes may offer insights into the complex dynamics of chromatin assembly and the fundamental mechanisms governing the structure and function of chromatin.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, gene expression, to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse grained force field for energy estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the Self-Organized Polymer model with Side Chains for proteins and the Three Interaction Site model for DNA in a modular fashion, without re-calibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , the crystallographic B-factors for DNA-protein complexes, with varying sizes and topologies, are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts the scattering profiles that are in quantitative agreement with SAXS experiments as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which does not alter the balance of electrostatic interactions, but affects chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Role of the histone tails in histone octamer transfer

The exceptionally high positive charge of the histones, concentrated in the N- and C-terminal tails, is believed to contribute to the stability of the nucleosome by neutralizing the negative charge of the nucleosomal DNA. We find, on the contrary, that the high positive charge contributes to instability, performing an essential function in chromatin remodeling. We show that the tails are required for removal of the histone octamer by the RSC chromatin remodeling complex, and this function is not due to direct RSC-tail interaction. We also show that the tails are required for histone octamer transfer from nucleosomes to DNA, and this activity of the tails is a consequence of their positive charge. Thus, the histone tails, intrinsically disordered protein regions, perform a critical role in chromatin structure and transcription, unrelated to their well-known role in regulation through posttranscriptional modification.

Single-Molecule Human Nucleosome Spontaneously Ruptures under the Stress of Compressive Force: A New Perspective on Gene Stability and Epigenetic Pathways

Force manipulation on the biological entities from living cells to protein molecules has revealed many mechanical details of cell biology from resolving folding and unfolding pathways to finding molecular interaction forces. A nucleosome is the basic repeating unit of chromatin where the histone octamer is wrapped by DNA, important for gene stability and regulation. How the inner side of the DNA gets accessed by other DNA binding molecules has been a puzzle that has been intensively studied and debated, important to epigenetics, gene stability, and regulations. Here we report our observation of spontaneous ruptures of human nucleosomes under pico-Newton (pN) compressive force. The amplitude of the compressive force, a squeezing rather than pulling force, involved in our experiment is tens of pN, which can be thermally available by biological force fluctuation at room temperature and under physiological conditions. This kind of structural rupture can loosen up the DNA around the histone, which in turn makes the DNA accessible to transcription and epigenetic modifications.

Nucleosome breathing facilitates cooperative binding of pluripotency factors Sox2 and Oct4 to DNA

Critical lineage commitment events are staged by multiple transcription factors (TFs) binding to their cognate motifs, often positioned at nucleosome-enriched regions of chromatin. The underlying mechanism remains elusive due to difficulty in disentangling the heterogeneity in chromatin states. Using a novel coarse-grained model and molecular dynamics simulations, here we probe the association of Sox2 and Oct4 proteins that show clustered binding at the entry-exit region of a nucleosome. The model captures the conformational heterogeneity of nucleosome breathing dynamics that features repeated wrap-unwrap transitions of a DNA segment from one end of the nucleosome. During the dynamics, DNA forms bulges that diffuse stochastically and may regulate the target search dynamics of a protein by nonspecifically interacting with it. The overall search kinetics of the TF pair follows a "dissociation-compensated-association" mechanism, where Oct4 binding is facilitated by the association of Sox2. The cooperativity stems from a change in entropy caused by an alteration in the nucleosome dynamics upon TF binding. The binding pattern is consistent with a live-cell single-particle tracking experiment, suggesting the mechanism observed for clustered binding of a TF pair, which is a hallmark of cis-regulatory elements, has broader implications in understanding gene regulation in a complex chromatin environment.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

Disordered regions tune order in chromatin organization and function

Intrinsically disordered proteins or hybrid proteins with ordered domains and disordered regions (both collectively designated as IDP(R)s) defy the well-established structure-function paradigm due to their ability to perform multiple biological functions even in the absence of a well-defined 3D structure. IDP(R)s have a unique ability to exist as a functional heterogeneous ensemble, where they adopt multiple thermodynamically stable conformations with low energy barriers between states. The resultant structural plasticity or conformational adaptability provides them with a high functional diversity and ease of regulation. Hence, IDP(R)s are highly efficient biological machinery to mediate intricate cellular functions such as signaling, gene expression, and assembly of complex structures. One such structure is the nucleoprotein complex known as Chromatin. Interestingly, the proteins involved in shaping up the structure and function of chromatin are abundant in disordered regions, which serve more than just as mere flexible linkers. The disordered regions are involved in crucial processes such as gene expression regulation, chromatin architecture maintenance, and liquid-liquid phase separation initiation. This review is an attempt to explore the advantages and the functional and regulatory roles of intrinsic disorder in several Chromatin Associated Proteins from a mechanistic standpoint.

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.

Potential of Mean Force for DNA Wrapping Around a Cationic Nanoparticle

Sharp bending and wrapping of DNA around proteins and nanoparticles (NPs) has been of extensive research interest. Here, we present the potential of mean force (PMF) for wrapping a DNA double helix around a cationic NP using coarse-grained models of a double-stranded DNA and a cationic NP. Starting from a NP wrapped around by DNA, the PMF was calculated along the distance between the center of the NP and one end of the DNA molecule. A relationship between the distance and the extent of DNA wrapping is used to calculate the PMF as a function of DNA wrapping around a NP. In particular, the PMF was compared for two DNA sequences of (AT)₂₅/(AT)₂₅ and (AC)₂₅/(GT)₂₅, for which the persistence lengths are different by ∼10 nm. The simulation results provide solid evidence of the thermodynamic preference for complex formation of a cationic NP with more flexible DNA over the less flexible DNA. Furthermore, we estimated the elastic energy of DNA bending, which was in good order-of-magnitude agreement with the theoretical prediction of elastic rods. This work suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnologies, in which the position and dynamics of NPs are regulated on large-scale DNA structures, or the structural transformation of DNA is triggered by the sequence-dependent binding of NPs.

Molecular dynamics simulations reveal how H3K56 acetylation impacts nucleosome structure to promote DNA exposure for lesion sensing

The first order of DNA packaging is the nucleosome with the DNA wrapped around the histone octamer. This leaves the nucleosomal DNA with access restrictions, which impose a significant barrier to repair of damaged DNA. The efficiency of DNA repair has been related to nucleosome structure and chromatin status, which is modulated in part by post-translational modifications (PTMs) of histones. Numerous studies have suggested a role for acetylation of lysine at position 56 of the H3 histone (H3K56ac) in various DNA transactions, including the response to DNA damage and its association with human cancer. Biophysical studies have revealed that H3K56ac increases DNA accessibility by facilitating spontaneous and transient unwrapping motions of the DNA ends. However, how this acetylation mark modulates nucleosome structure and dynamics to promote accessibility to the damaged DNA for repair factors and other proteins is still poorly understood. Here, we utilize approximately 5-6 microseconds of atomistic molecular dynamics simulations to delineate the impact of H3K56 acetylation on the nucleosome structure and dynamics, and to elucidate how these nucleosome properties are further impacted when a bulky benzo[a]pyrene-derived DNA lesion is placed near the acetylation site. Our findings reveal that H3K56ac alone induces considerable disturbance to the histone-DNA/histone-histone interactions, and amplifies the distortions imposed by the presence of the lesion. Our work highlights the important role of H3K56 acetylation in response to DNA damage and depicts how access to DNA lesions by the repair machinery can be facilitated within the nucleosome via a key acetylation event.

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.

Effect of DNA Flexibility on Complex Formation of a Cationic Nanoparticle with Double-Stranded DNA

We present extensive molecular dynamics simulations of a cationic nanoparticle and a double-stranded DNA molecule to discuss the effect of DNA flexibility on the complex formation of a cationic nanoparticle with double-stranded DNA. Martini coarse-grained models were employed to describe double-stranded DNA molecules with two different flexibilities and cationic nanoparticles with three different electric charges. As the electric charge of a cationic nanoparticle increases, the degree of DNA bending increases, eventually leading to the wrapping of DNA around the nanoparticle at high electric charges. However, a small increase in the persistence length of DNA by 10 nm requires a cationic nanoparticle with a markedly increased electric charge to bend and wrap DNA around. Thus, a more flexible DNA molecule bends and wraps around a cationic nanoparticle with an intermediate electric charge, whereas a less flexible DNA molecule binds to a nanoparticle with the same electric charge without notable bending. This work provides solid evidence that a small difference in DNA flexibility (as small as 10 nm in persistence length) has a substantial influence on the complex formation of DNA with proteins from a biological perspective and suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnology as a new tool to manipulate the structure of DNA molecules mediated by nanoparticle binding.

Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA

Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

Surprising Twists in Nucleosomal DNA with Implication for Higher-order Folding

While nucleosomes are dynamic entities that must undergo structural deformations to perform their functions, the general view from available high-resolution structures is a largely static one. Even though numerous examples of twist defects have been documented, the DNA wrapped around the histone core is generally thought to be overtwisted. Analysis of available high-resolution structures from the Protein Data Bank reveals a heterogeneous distribution of twist along the nucleosomal DNA, with clear patterns that are consistent with the literature, and a significant fraction of structures that are undertwisted. The subtle differences in nucleosomal DNA folding, which extend beyond twist, have implications for nucleosome disassembly and modeled higher-order structures. Simulations of oligonucleosome arrays built with undertwisted models behave very differently from those constructed from overtwisted models, in terms of compaction and inter-nucleosome contacts, introducing configurational changes equivalent to those associated with 2-3 base-pair changes in nucleosome spacing. Differences in the nucleosomal DNA pathway, which underlie the way that DNA enters and exits the nucleosome, give rise to different nucleosome-decorated minicircles and affect the topological mix of configurational states.

Histone tails cooperate to control the breathing of genomic nucleosomes

Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145–147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation.

In the cell, the DNA is packed in chromatin. Chromatin is a highly dynamic fiber structure made of arrays of nucleosomes with different degrees of compaction. Each nucleosome has 145–147 basepairs of DNA wrapped around a protein octamer made of four unique histone proteins. Each histone is present twice and has a structured part and one or two disordered terminal tails. The regulation of gene expression in the cell and during cellular transitions depends on dynamic changes in chromatin structure. Chromatin dynamics are modulated by intra and inter nucleosome motions and by posttranslational chemical modifications of the histone tails. Here we reveal how histone tails control the intra nucleosome dynamics at atomic resolution. From extensive sampling of nucleosome dynamics in atomistic molecular simulations, we show that genomic nucleosomes breath more extensively than engineered ones and we describe how two histone tails cooperate to control nucleosome breathing through interactions between clusters of positively charged residues and the DNA. Nucleosome conformations with different degrees of opening are associated with different conformations, positions, and DNA interaction patterns of the tails. With this mechanism, we contribute to the understanding of chromatin dynamics at atomic resolution.

Histone tails as signaling antennas of chromatin

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

Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface

The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome. Dynamical simulations performed with an advanced sampling technique uncovered multiple pathways that connect open chromatin configurations with the zigzag crystal structure. Intermediate states along the simulated folding pathways resemble chromatin configurations reported from in situ experiments. We further determined a six-dimensional free energy surface as a function of the inter-nucleosome distances via a deep learning approach. The zigzag structure can indeed be seen as the global minimum of the surface. However, it is not favored by a significant amount relative to the partially unfolded, in situ configurations. Chemical perturbations such as histone H4 tail acetylation and thermal fluctuations can further tilt the energetic balance to stabilize intermediate states. Our study provides insight into the connection between various reported chromatin configurations and has implications on the in situ relevance of the 30 nm fiber.

The three-dimensional organization of chromatin plays critical roles in regulating genome function. Here the authors apply a near atomistic model to study the structure and dynamics of the chromatin folding unit - the tetra-nucleosome - to provide insight into how chromatin folds.

Histone Tail Conformations: A Fuzzy Affair with DNA

The core histone tails are critical in chromatin structure and signaling. Studies over the past several decades have provided a wealth of information on the histone tails and their interaction with chromatin factors. However, the conformation of the histone tails in a chromatin relevant context has remained elusive. Only recently has enough evidence emerged to start to build a structural model of the tails in the context of nucleosomes and nucleosome arrays. Here, we review these studies and propose that the histone tails adopt a high-affinity fuzzy complex with DNA, characterized by robust but dynamic association. Furthermore, we discuss how these DNA-bound conformational ensembles promote distinct chromatin structure and signaling, and that their fuzzy nature is important in transitioning between functional states.

Nucleosome Core Particles Lacking H2B or H3 Tails Are Altered Structurally and Have Differential Base Excision Repair Fingerprints

A recently discovered post-translational modification of histone proteins is the irreversible proteolytic clipping of the histone N-terminal tail domains. This modification is involved in the regulation of various biological processes, including the DNA damage response. In this work, we used chemical footprinting to characterize the structural alterations to nucleosome core particles (NCPs) that result from a lack of a histone H2B or H3 tail. We also examine the influence of these histone tails on excision of the mutagenic lesion 1,N⁶-ethenoadenine (εA) by the repair enzyme alkyladenine DNA glycosylase. We found that the absence of the H2B or H3 tail results in altered DNA periodicity relative to that of native NCPs. We correlated these structural alterations to εA excision by utilizing a global analysis of 21 εA sites in NCPs and unincorporated duplex DNA. In comparison to native NCPs, there is enhanced excision of εA in tailless H2B NCPs in regions that undergo DNA unwrapping. This enhanced excision is not observed for tailless H3 NCPs; rather, excision is inhibited in more static areas of the NCP not prone to unwrapping. Our results support in vivo observations of alkylation damage profiles and the potential role of tail clipping as a mechanism for overcoming physical obstructions caused by packaging in NCPs but also reveal the potential inhibition of repair by tail clipping in some locations. Taken together, these results further our understanding of how base excision repair can be facilitated or diminished by histone tail removal and contribute to our understanding of the underlying mechanism that leads to mutational hot spots.

Linker DNA Length is a Key to Tri-nucleosome Folding

The folding of a nucleosome array has long been one of the fundamental and unsolved problems in chromatin biology. In this study, we address how nucleosome array folding depends on the length of linker DNA. We performed molecular dynamics simulations of a tri-nucleosome, a minimal model of chromatin folding, with various linker lengths (LLs) ranging from 20 to 40 base pairs (bps). We found that the tri-nucleosome folding strongly depends on LLs, and classified the structure ensemble into five classes, named from trinuc-1 to trinuc-5. As a function of LL, the different classes appear, on average, every 2 bps with a period of 10 bps, and are characterized by distinct inter-nucleosome interactions. The trinuc-1 conformation corresponds to LL ~ 10n, where n is an integer, and is stabilized by the tight packing between the first and the third nucleosomes, consistent with a zigzag fiber form. Structures of the other four classes are more diverse and distributed continuously in the space of possible configurations. Histone-DNA electrostatic interactions in the tri-nucleosome are further analyzed.

Single-molecule and in silico dissection of the interaction between Polycomb repressive complex 2 and chromatin

Polycomb repressive complex 2 (PRC2) installs and spreads repressive histone methylation marks on eukaryotic chromosomes. Because of the key roles that PRC2 plays in development and disease, how this epigenetic machinery interacts with DNA and nucleosomes is of major interest. Nonetheless, the mechanism by which PRC2 engages with native-like chromatin remains incompletely understood. In this work, we employ single-molecule force spectroscopy and molecular dynamics simulations to dissect the behavior of PRC2 on polynucleosome arrays. Our results reveal an unexpectedly diverse repertoire of PRC2 binding configurations on chromatin. Besides reproducing known binding modes in which PRC2 interacts with bare DNA, mononucleosomes, and adjacent nucleosome pairs, our data also provide direct evidence that PRC2 can bridge pairs of distal nucleosomes. In particular, the "1-3" bridging mode, in which PRC2 engages two nucleosomes separated by one spacer nucleosome, is a preferred low-energy configuration. Moreover, we show that the distribution and stability of different PRC2-chromatin interaction modes are modulated by accessory subunits, oncogenic histone mutations, and the methylation state of chromatin. Overall, these findings have implications for the mechanism by which PRC2 spreads histone modifications and compacts chromatin. The experimental and computational platforms developed here provide a framework for understanding the molecular basis of epigenetic maintenance mediated by Polycomb-group proteins.

The HMGB chromatin protein Nhp6A can bypass obstacles when traveling on DNA

DNA binding proteins rapidly locate their specific DNA targets through a combination of 3D and 1D diffusion mechanisms, with the 1D search involving bidirectional sliding along DNA. However, even in nucleosome-free regions, chromosomes are highly decorated with associated proteins that may block sliding. Here we investigate the ability of the abundant chromatin-associated HMGB protein Nhp6A from Saccharomyces cerevisiae to travel along DNA in the presence of other architectural DNA binding proteins using single-molecule fluorescence microscopy. We observed that 1D diffusion by Nhp6A molecules is retarded by increasing densities of the bacterial proteins Fis and HU and by Nhp6A, indicating these structurally diverse proteins impede Nhp6A mobility on DNA. However, the average travel distances were larger than the average distances between neighboring proteins, implying Nhp6A is able to bypass each of these obstacles. Together with molecular dynamics simulations, our analyses suggest two binding modes: mobile molecules that can bypass barriers as they seek out DNA targets, and near stationary molecules that are associated with neighboring proteins or preferred DNA structures. The ability of mobile Nhp6A molecules to bypass different obstacles on DNA suggests they do not block 1D searches by other DNA binding proteins.

Observation of Ethanol-Induced Condensation and Decondensation Processes at a Single-DNA Molecular Level in Microfluidic Devices Equipped with a Rapid Solution Exchange System

Conformational transitions from secondary (e.g., B- to A-form DNA) to higher-order (e.g., coil to globule) transitions play important roles in genome expression and maintenance. Several single-molecule approaches using microfluidic devices have been used to determine the kinetics of DNA chromatin assembly because microfluidic devices can afford stretched DNA molecules through laminar flow and rapid solution exchange. However, some issues, particularly the uncertainty of time 0 in the solution exchange process, are encountered. In such kinetic experiments, it is critical to determine when the target solution front approaches the target DNA molecules. Therefore, a new design for a microfluidic device is developed that enables the instantaneous exchange of solutions in the observation channel, allowing accurate measurements of DNA conformational transitions; stepwise, ethanol-induced conformational transitions are revealed. Although full DNA contraction from coil to globule is observed with >50% ethanol, no outstanding change is observed at concentrations <40% in 10 min. With 50% ethanol solution, the DNA conformational transition passes through two steps: (i) fast and constant-velocity contraction and (ii) relatively slow contraction from the free end. The first process is attributed to the B to A conformational transition by gradual dehydration. The second process is due to the coil-globule transition as the free end of DNA starts the contraction. This globular structure formation counteracts the shear force from the microfluids and decelerates the contraction velocity. This real-time observation system can be applied to the kinetic analysis of DNA conformational transitions such as kinetics of chromatin assembly and gene expression.

White spot syndrome viral protein VP9 alters the cellular higher-order chromatin structure

Viral protein 9 (VP9) is a non-structural protein of white spot syndrome virus (WSSV) highly expressed during the early stage of infection. The crystal structure of VP9 suggests that the polymers of VP9 dimers resemble a DNA mimic, but its function remains elusive. In this study, we demonstrated that VP9 impedes histones binding to DNA via single-molecule manipulation. We established VP9 expression in HeLa cells due to the lack of a WSSV-susceptible cell line, and observed abundant VP9 in the nucleus, which mirrors its distribution in the hemocytes of WSSV-infected shrimp. VP9 expression increased the dynamics and rotational mobility of histones in stable H3-GFP HeLa cells as revealed by fluorescent recovery after photobleaching and fluorescence anisotropy imaging, which suggested a loosened compaction of chromatin structure. Successive salt fractionation showed that a prominent population of histones was solubilized in high salt concentrations, which implies alterations of bulk chromatin structure. Southern blotting identified that VP9 alters juxtacentromeric chromatin structures to be more accessible to micrococcal nuclease digestion. RNA microarray revealed that VP9 expression also leads to significant changes of cellular gene expression. Our findings provide evidence that VP9 alters the cellular higher-order chromatin structure, uncovering a potential strategy adopted by WSSV to facilitate its replication.

Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding

Transcription factor (TF) proteins bind to DNA to regulate gene expression. Normally, accessibility to DNA is required for their function. However, in the nucleus, the DNA is often inaccessible, wrapped around histone proteins in nucleosomes forming the chromatin. Pioneer TFs are thought to induce chromatin opening by recognizing their DNA binding sites on nucleosomes. For example, Oct4, a master regulator and inducer of stem cell pluripotency, binds to DNA in nucleosomes in a sequence-specific manner. Here, we reveal the structural dynamics of nucleosomes that mediate Oct4 binding from molecular dynamics simulations. Nucleosome flexibility and the amplitude of nucleosome motions such as breathing and twisting are enhanced in nucleosomes with multiple TF binding sites. Moreover, the regions around the binding sites display higher local structural flexibility. Probing different structures of Oct4-nucleosome complexes, we show that alternative configurations in which Oct4 recognizes partial binding sites display stable TF-DNA interactions similar to those observed in complexes with free DNA and compatible with the DNA curvature and DNA-histone interactions. Therefore, we propose a structural basis for nucleosome recognition by a pioneer TF that is essential for understanding how chromatin is unraveled during cell fate conversions.

Nucleosomes as allosteric scaffolds for genetic regulation

Nucleosomes are stable yet highly dynamic complexes exhibiting diverse types of motions, such as sliding, DNA unwrapping, and disassembly, encoding a landscape with a large number of metastable states. In this review, describing recent studies on these nucleosome structure changes, we propose that the nucleosome can be viewed as an ideal allosteric scaffold: regulated by effector molecules such as transcription factors and chromatin remodelers, the nucleosome controls the downstream gene activity. Binding of transcription factors to the nucleosome can enhance DNA unwrapping or slide the DNA, altering either the binding or the unbinding of other transcription factors to nearby sites. ATP-dependent chromatin remodelers induce a series of DNA deformations, which allosterically propagate throughout the nucleosome to induce DNA sliding or histone exchange.

Variable impact of conformationally distinct DNA lesions on nucleosome structure and dynamics: Implications for nucleotide excision repair

The packaging of DNA in nucleosomes presents a barrier for biological transactions including replication, transcription and repair. However, despite years of research, how the DNA is freed from the histone proteins and thereby allows the molecular machines to access the DNA remains poorly understood. We are interested in global genomic nucleotide excision repair (GG-NER). It is established that the histones are obstacles to this process, and DNA lesions are repaired less efficiently in nucleosomes than in free DNA. In the present study, we utilized molecular dynamics simulations to elucidate the nature of the distortions and dynamics imposed in the nucleosome by a set of three structually different lesions that vary in GG-NER efficiencies in free DNA, and in nucleosomes [Shafirovich, Geacintov, et. al, 2019]. Two of these are bulky lesions derived from metabolic activation of the environmental carcinogen benzo[a]pyrene, the 10R (+)-cis-anti-B[a]P-N²-dG and the stereoisomeric 10S (+)-trans-anti-B[a]P-N²-dG, which respectively adopt base-displaced/intercalated and minor groove-aligned conformations in DNA. The third is a non-bulky lesion, the 5'R-8-cyclo-2'-deoxyguanosine cross-link, produced by reactive oxygen and nitrogen species; cyclopurine lesions are highly mutagenic. These adducts are placed near the dyad axis, and rotationally with the lesion-containing strand facing towards or away from the histones. While each lesion has distinct conformational characteristics that are retained in the nucleosome, a spectrum of structural and dynamic disturbances, from slight to substantial, are displayed that depend on the lesion's structure and position in the nucleosome. We hypothesize that these intrinsic structural and dynamic distinctions provide different signals to initiate the cascade of chromatin-opening processes, including acetylation and other post translational modifications, remodeling by ATP-dependent complexes and spontaneous unwrapping that regulate the rate of access to the lesion; this may translate ultimately into varying GG-NER efficiencies, including repair resistance when signals for access are too weak.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

Chasing Tails: Cathepsin-L Improves Structural Analysis of Histones by HX-MS

The N-terminal regions (tails) of histone proteins are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications present on these regions contribute to transcriptional regulation. Considering their biological significance, relatively few structural details have been established for histone tails, mainly because of their inherently disordered nature. Although hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed in this context, presumably because of the poor N-terminal coverage provided by pepsin. Inspired from histone-clipping events, we profiled the activity of cathepsin-L under HX-MS quench conditions and characterized its specificity employing the four core histones (H2A, H2B, H3 and H4). Cathepsin-L demonstrated cleavage patterns that were substrate- and pH-dependent. Cathepsin-L generated overlapping N-terminal peptides about 20 amino acids long for H2A, H3, and H4 proving its suitability for the analysis of histone tails dynamics. We developed a comprehensive HX-MS method in combination with pepsin and obtained full sequence coverage for all histones. We employed our method to analyze histones H3 and H4. We observe rapid deuterium exchange of the N-terminal tails and cooperative unfolding (EX1 kinetics) in the histone-fold domains of histone monomers in-solution. Overall, this novel strategy opens new avenues for investigating the dynamic properties of histones that are not apparent from the crystal structures, providing insights into the structural basis of the histone code.

Critical role of histone tail entropy in nucleosome unwinding

The nucleosome is the fundamental packaging unit for the genome. It must remain tightly wound to ensure genome stability while simultaneously being flexible enough to keep the DNA molecule accessible for genome function. The set of physicochemical interactions responsible for the delicate balance between these naturally opposed processes have not been determined due to challenges in resolving partially unwound nucleosome configurations at atomic resolution. Using a near atomistic protein-DNA model and advanced sampling techniques, we calculate the free energy cost of nucleosome DNA unwinding. Our simulations identify a large energetic barrier that decouples the outer and the inner DNA unwinding into two separate processes, occurring on different time scales. This dynamical decoupling allows the exposure of outer DNA at a modest cost to ensure accessibility while keeping the inner DNA and the histone core intact to maintain stability. We also reveal that this energetic barrier arises from a delayed loss of contacts between disordered histone tails and the DNA and is, surprisingly, largely offset by an entropic contribution from these tails. Implications of this enthalpy entropy compensation for the regulation of nucleosome stability and genome function are discussed.

Nucleosome Crowding in Chromatin Slows the Diffusion but Can Promote Target Search of Proteins

Dynamics of nuclear proteins in crowded chromatin has only been poorly understood. Here, we address the diffusion, target search, and structural dynamics of three proteins in a model chromatin using coarse-grained molecular simulations run on the K computer. We prepared two structures of chromatin made of 20 nucleosomes with different nucleosome densities and investigated dynamics of two transcription factors, HMGB1 and p53, and one signaling protein, ERK, embedded in the chromatin. We found fast and normal diffusion of the nuclear proteins in the low-density chromatins and slow and subdiffusional movements in the high-density chromatin. The diffusion of the largest transcription factor, p53, is slowed by high-density chromatin most markedly. The on rates and off rates for DNA binding are increased and decreased, respectively, in the high-density chromatin. To our surprise, the DNA sequence search was faster in chromatin with high nucleosome density, though the diffusion is slower. We also found that the three nuclear proteins preferred to bind on the linker DNA and the entry and exit regions of nucleosomal DNA. In addition to these regions, HMGB1 and p53 also bound to the dyad.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

AWSEM-IDP: A Coarse-Grained Force Field for Intrinsically Disordered Proteins

The associative memory, water-mediated, structure and energy model (AWSEM) has been successfully used to study protein folding, binding, and aggregation problems. In this work, we introduce AWSEM-IDP, a new AWSEM branch for simulating intrinsically disordered proteins (IDPs), where the weights of the potentials determining secondary structure formation have been finely tuned, and a novel potential is introduced that helps to precisely control both the average extent of protein chain collapse and the chain's fluctuations in size. AWSEM-IDP can efficiently sample large conformational spaces, while retaining sufficient molecular accuracy to realistically model proteins. We applied this new model to two IDPs, demonstrating that AWSEM-IDP can reasonably well reproduce higher-resolution reference data, thus providing the foundation for a transferable IDP force field. Finally, we used thermodynamic perturbation theory to show that, in general, the conformational ensembles of IDPs are highly sensitive to fine-tuning of force field parameters.

Chromatin remodelers couple inchworm motion with twist-defect formation to slide nucleosomal DNA

ATP-dependent chromatin remodelers are molecular machines that control genome organization by repositioning, ejecting, or editing nucleosomes, activities that confer them essential regulatory roles on gene expression and DNA replication. Here, we investigate the molecular mechanism of active nucleosome sliding by means of molecular dynamics simulations of the Snf2 remodeler translocase in complex with a nucleosome. During its inchworm motion driven by ATP consumption, the translocase overwrites the original nucleosome energy landscape via steric and electrostatic interactions to induce sliding of nucleosomal DNA unidirectionally. The sliding is initiated at the remodeler binding location via the generation of a pair of twist defects, which then spontaneously propagate to complete sliding throughout the entire nucleosome. We also reveal how remodeler mutations and DNA sequence control active nucleosome repositioning, explaining several past experimental observations. These results offer a detailed mechanistic picture of remodeling important for the complete understanding of these key biological processes.

Nucleosomes are the protein-DNA complexes underlying Eukaryotic genome organization, and serve as regulators of gene expression by occluding DNA to other proteins. This regulation requires the precise positioning of nucleosomes along DNA. Chromatin remodelers are the molecular machines that consume ATP to slide nucleosome at their correct locations, but the mechanisms of remodeling are still unclear. Based on the static structural information of a remodeler bound on nucleosome, we performed molecular dynamics computer simulations revealing the details of how remodelers slide nucleosomal DNA: the inchworm-like motion of remodelers create small DNA deformations called twist defects, which then spontaneously propagate throughout the nucleosome to induce sliding. These simulations explain several past experimental findings and are important for our understanding of genome organization.

Dynamics of Chromatin Fibers: Comparison of Monte Carlo Simulations with Force Spectroscopy

To elucidate conformational dynamics of chromatin fibers, we compared available force-spectroscopy measurements with extensive Monte Carlo simulations of nucleosome arrays under external force. Our coarse-grained model of chromatin includes phenomenological energy terms for the DNA-histone adhesion and the internucleosome stacking interactions. We found that the Monte Carlo fiber ensembles simulated with increasing degrees of DNA unwrapping and the stacking energy 8 kT can account for the intricate force-extension response observed experimentally. Our analysis shows that at low external forces (F < 3.0 picoNewtons), the DNA ends in nucleosomes breathe by ∼10 bp. Importantly, under these conditions, the fiber is highly dynamic, exhibiting continuous unstacking-restacking transitions, allowing accessibility of transcription factors to DNA, while maintaining a relatively compact conformation. Of note, changing the stacking interaction by a few kT, an in silico way to mimic histone modifications, is sufficient to transform an open chromatin state into a compact fiber. The fibers are mostly two-start zigzag folds with rare occurrences of three- to five-start morphologies. The internucleosome stacking is lost during the linear response regime. At the higher forces exceeding 4 picoNewtons, the nucleosome unwrapping becomes stochastic and asymmetric, with one DNA arm opened by ∼55 bp and the other arm only by ∼10 bp. Importantly, this asymmetric unwrapping occurs for any kind of sequence, including the symmetric ones. Our analysis brings new, to our knowledge, insights in dynamics of chromatin modulated by histone epigenetic modifications and molecular motors such as RNA polymerase.

Native DNA electronics: is it a matter of nanoscale assembly?

The genomic DNA is enveloped by nanotubes formed by the nuclear aggregates of polyamines (NAPs) that induce DNA conformational changes and provide protection and increased interaction abilities for the double strands. In a physiological environment, the nanotube arrangement is initiated by spontaneous interaction among the terminal amino groups of the polyamines and the phosphate ions, with the consequent formation of cyclic monomers that hook at the DNA grooves. The polymer thus formed has the morphological features of an organic semiconductor and therefore, it can be considered to be able to conduct electric charges. Phosphate ions positioned on the NAP external surface could regulate, as in a physical electric circuit, both linear and rotational (histones) protein motion, in accordance with the basilar principles of the electronics. A model of a carrier system for protein motion along the polymer wrapping the DNA strands, based on the phosphate-phosphate complexation, is proposed.

Mechanical evolution of DNA double-strand breaks in the nucleosome

Double strand breaks (DSB) in the DNA backbone are the most lethal type of defect induced in the cell nucleus by chemical and radiation treatments of cancer. However, little is known about the outcomes of damage in nucleosomal DNA, and on its effects on damage repair. We performed microsecond-long molecular dynamics computer simulations of nucleosomes including a DSB at various sites, to characterize the early stages of the evolution of this DNA lesion. The damaged structures are studied by the essential dynamics of DNA and histones, and compared to the intact nucleosome, thus exposing key features of the interactions. All DSB configurations tend to remain compact, with only the terminal bases interacting with histone proteins. Umbrella sampling calculations show that broken DNA ends at the DSB must overcome a free-energy barrier to detach from the nucleosome core. Finally, by calculating the covariant mechanical stress, we demonstrate that the coupled bending and torsional stress can force the DSB free ends to open up straight, thus making it accessible to damage signalling proteins.

Interactions of HP1 Bound to H3K9me3 Dinucleosome by Molecular Simulations and Biochemical Assays

Heterochromatin protein 1 (HP1), associated with heterochromatin formation, recognizes an epigenetically repressive marker, trimethylated lysine 9 in histone H3 (H3K9me3), and generally contributes to long-term silencing. How HP1 induces heterochromatin is not fully understood. Recent experiments suggested that not one, but two nucleosomes provide a platform for this recognition. Integrating previous and new biochemical assays with computational modeling, we provide near-atomic structural models for HP1 binding to the dinucleosomes. We found that the dimeric HP1α tends to bind two H3K9me3s that are in adjacent nucleosomes, thus bridging two nucleosomes. We identified, to our knowledge, a novel DNA binding motif in the hinge region that is specific to HP1α and is essential for recognizing the H3K9me3 sites of two nucleosomes. An HP1 isoform, HP1γ, does not easily bridge two nucleosomes in extended conformations because of the absence of the above binding motif and its shorter hinge region. We propose a molecular mechanism for chromatin structural changes caused by HP1.