Torsional behavior of chromatin is modulated by rotational phasing of nucleosomes

Nick Induced Dynamics in Supercoiled DNA Facilitates the Protein Target Search Process

A DNA nick, defined as a discontinuity in a double-stranded DNA molecule where the phosphodiester bond between adjacent nucleotides of one strand is absent due to enzyme action, serves as an effective mechanism to alleviate stress in supercoiled DNA. This stress release is essential for the smooth operation of transcriptional machinery. However, the underlying mechanisms and their impact on protein search dynamics, which are crucial for initiating transcription, remain unclear. Through extensive computer simulations, we unravel the molecular picture, demonstrating that intramolecular stress release due to a DNA nick is driven by a combination of writhing and twisting motions, depending on the nick's position. This stress release is quantitatively manifested as a step-like increase in the linking number. Furthermore, we elucidate that the nicked supercoiled minicircles exhibit enhanced torsional dynamics, promoting rapid conformational changes and frequent shifts in the identities of juxtaposed DNA sites on the plectoneme. The dynamics of the juxtaposition sites facilitates communication between protein and DNA, resulting in faster protein diffusion compared with native DNA with the same topology. Our findings highlight the mechanistic intricacies and underscore the importance of DNA nicks in facilitating transcription elongation by actively managing torsional stress during DNA unwinding by the RNA polymerase.

Nucleosomes play a dual role in regulating transcription dynamics

Transcription has a mechanical component, as the translocation of the transcription machinery or RNA polymerase (RNAP) on DNA or chromatin is dynamically coupled to the chromatin torsion. This posits chromatin mechanics as a possible regulator of eukaryotic transcription, however, the modes and mechanisms of this regulation are elusive. Here, we first take a statistical mechanics approach to model the torsional response of topology-constrained chromatin. Our model recapitulates the experimentally observed weaker torsional stiffness of chromatin compared to bare DNA and proposes structural transitions of nucleosomes into chirally distinct states as the driver of the contrasting torsional mechanics. Coupling chromatin mechanics with RNAP translocation in stochastic simulations, we reveal a complex interplay of DNA supercoiling and nucleosome dynamics in governing RNAP velocity. Nucleosomes play a dual role in controlling the transcription dynamics. The steric barrier aspect of nucleosomes in the gene body counteracts transcription via hindering RNAP motion, whereas the chiral transitions facilitate RNAP motion via driving a low restoring torque upon twisting the DNA. While nucleosomes with low dissociation rates are typically transcriptionally repressive, highly dynamic nucleosomes offer less of a steric barrier and enhance the transcription elongation dynamics of weakly transcribed genes via buffering DNA twist. We use the model to predict transcription-dependent levels of DNA supercoiling in segments of the budding yeast genome that are in accord with available experimental data. The model unveils a paradigm of DNA supercoiling-mediated interaction between genes and makes testable predictions that will guide experimental design.

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

Computational approaches for inferring 3D conformations of chromatin from chromosome conformation capture data

Chromosome conformation capture (3C) and its variants are powerful experimental techniques for probing intra- and inter-chromosomal interactions within cell nuclei at high resolution and in a high-throughput, quantitative manner. The contact maps derived from such experiments provide an avenue for inferring the 3D spatial organization of the genome. This review provides an overview of the various computational methods developed in the past decade for addressing the very important but challenging problem of deducing the detailed 3D structure or structure population of chromosomal domains, chromosomes, and even entire genomes from 3C contact maps.

Insight into the unwrapping of the dinucleosome

Dynamics of nucleosomes, the building blocks of chromatin, has crucial effects on the expression, replication and repair of genomes in eukaryotes. Beside the constant movements of nucleosomes by thermal fluctuations, ATP-dependent chromatin remodelling complexes cause their active displacements. Here we propose a theoretical analysis of dinucleosome wrapping and unwrapping dynamics in the presence of an external force. We explore the energy landscape and configurations of a dinucleosome in different unwrapped states. Moreover, using a dynamical Monte-Carlo simulation algorithm, we demonstrate the dynamical features of the system such as the unwrapping force for partial and full wrapping processes. Furthermore, we show that in the short length of linker DNA (∼10-90 bp), asymmetric unwrapping occurs. These findings could shed some light on chromatin dynamics and gene accessibility.

H3K4 methylation at active genes mitigates transcription-replication conflicts during replication stress

Transcription-replication conflicts (TRCs) occur when intensive transcriptional activity compromises replication fork stability, potentially leading to gene mutations. Transcription-deposited H3K4 methylation (H3K4me) is associated with regions that are susceptible to TRCs; however, the interplay between H3K4me and TRCs is unknown. Here we show that H3K4me aggravates TRC-induced replication failure in checkpoint-defective cells, and the presence of methylated H3K4 slows down ongoing replication. Both S-phase checkpoint activity and H3K4me are crucial for faithful DNA synthesis under replication stress, especially in highly transcribed regions where the presence of H3K4me is highest and TRCs most often occur. H3K4me mitigates TRCs by decelerating ongoing replication, analogous to how speed bumps slow down cars. These findings establish the concept that H3K4me defines the transcriptional status of a genomic region and defends the genome from TRC-mediated replication stress and instability.

Transcription-replication conflicts (TRC) can contribute to genome instability. Here the authors reveal that under replication stress H3K4 methylation can play a role in TRC prevention.

Properties of gene expression and chromatin structure with mechanically regulated elongation

In recent years, physical elements of transcription have emerged as central in our understanding of gene expression. Recent work has been done introducing a simple description of the basic physical elements of transcription where RNA elongation, RNA polymerase (RNAP) rotation and DNA supercoiling are coupled (1). Here we generalize this framework to accommodate the behavior of many RNAPs operating on multiple genes on a shared piece of DNA. The resulting framework is combined with well-established stochastic processes of transcription resulting in a model which characterizes the impact of the mechanical properties of transcription on gene expression and DNA structure. Transcriptional bursting readily emerges as a common phenomenon with origins in the geometric nature of the genetic system and results in the bounding of gene expression statistics. Properties of a multiple gene system are examined with special attention paid to the role that genome composition (gene orientation, size and intergenic distance) plays in the ability of genes to transcribe. The role of transcription in shaping DNA structure is examined and the possibility of transcription driven domain formation is discussed.

Transfer-matrix calculations of the effects of tension and torque constraints on DNA-protein interactions

Organization and maintenance of the chromosomal DNA in living cells strongly depends on the DNA interactions with a plethora of DNA-binding proteins. Single-molecule studies show that formation of nucleoprotein complexes on DNA by such proteins is frequently subject to force and torque constraints applied to the DNA. Although the existing experimental techniques allow to exert these type of mechanical constraints on individual DNA biopolymers, their exact effects in regulation of DNA-protein interactions are still not completely understood due to the lack of systematic theoretical methods able to efficiently interpret complex experimental observations. To fill this gap, we have developed a general theoretical framework based on the transfer-matrix calculations that can be used to accurately describe behaviour of DNA-protein interactions under force and torque constraints. Potential applications of the constructed theoretical approach are demonstrated by predicting how these constraints affect the DNA-binding properties of different types of architectural proteins. Obtained results provide important insights into potential physiological functions of mechanical forces in the chromosomal DNA organization by architectural proteins as well as into single-DNA manipulation studies of DNA-protein interactions.

Chromatin Higher-Order Folding: A Perspective with Linker DNA Angles

The mechanism by which the "beads-on-a-string" nucleosome chain folds into various higher-order chromatin structures in eukaryotic cell nuclei is still poorly understood. The various models depicting higher-order chromatin as regular helical fibers and the very opposite "polymer melt" theory imply that interactions between nucleosome "beads" make the main contribution to the chromatin compaction. Other models in which the geometry of linker DNA "strings" entering and exiting the nucleosome define the three-dimensional structure predict that small changes in the linker DNA configuration may strongly affect nucleosome chain folding and chromatin higher-order structure. Among those studies, the cross-disciplinary approach pioneered by Jörg Langowski that combines computational modeling with biophysical and biochemical experiments was most instrumental for understanding chromatin higher-order structure in vitro. Strikingly, many recent studies, including genome-wide nucleosome interaction mapping and chromatin imaging, show an excellent agreement with the results of three-dimensional computational modeling based on the primary role of linker DNA geometry in chromatin compaction. This perspective relates nucleosome array models with experimental studies of nucleosome array folding in vitro and in situ. I argue that linker DNA configuration plays a key role in determining nucleosome chain flexibility, topology, and propensity for self-association, thus providing new implications for regulation of chromatin accessibility to DNA binding factors and RNA transcription machinery as well as long-range communications between distant genomic sites.

The dynamic interplay between DNA topoisomerases and DNA topology

Topological properties of DNA influence its structure and biochemical interactions. Within the cell, DNA topology is constantly in flux. Transcription and other essential processes, including DNA replication and repair, not only alter the topology of the genome but also introduce additional complications associated with DNA knotting and catenation. These topological perturbations are counteracted by the action of topoisomerases, a specialized class of highly conserved and essential enzymes that actively regulate the topological state of the genome. This dynamic interplay among DNA topology, DNA processing enzymes, and DNA topoisomerases is a pervasive factor that influences DNA metabolism in vivo. Building on the extensive structural and biochemical characterization over the past four decades that has established the fundamental mechanistic basis of topoisomerase activity, scientists have begun to explore the unique roles played by DNA topology in modulating and influencing the activity of topoisomerases. In this review we survey established and emerging DNA topology-dependent protein-DNA interactions with a focus on in vitro measurements of the dynamic interplay between DNA topology and topoisomerase activity.

The Dynamic Interplay Between DNA Topoisomerases and DNA Topology

Topological properties of DNA influence its structure and biochemical interactions. Within the cell DNA topology is constantly in flux. Transcription and other essential processes including DNA replication and repair, alter the topology of the genome, while introducing additional complications associated with DNA knotting and catenation. These topological perturbations are counteracted by the action of topoisomerases, a specialized class of highly conserved and essential enzymes that actively regulate the topological state of the genome. This dynamic interplay among DNA topology, DNA processing enzymes, and DNA topoisomerases, is a pervasive factor that influences DNA metabolism in vivo. Building on the extensive structural and biochemical characterization over the past four decades that established the fundamental mechanistic basis of topoisomerase activity, the unique roles played by DNA topology in modulating and influencing the activity of topoisomerases have begun to be explored. In this review we survey established and emerging DNA topology dependent protein-DNA interactions with a focus on in vitro measurements of the dynamic interplay between DNA topology and topoisomerase activity.

Elastic model for dinucleosome structure and energy

The equilibrium structure of a dinucleosome is studied using an elastic model that takes into account the force and torque balance conditions. Using the proper boundary conditions, it is found that the conformational energy of the problem does not depend on the length of the linker DNA. In addition it is shown that the two histone octamers are almost perpendicular to each other, and the linker DNA in short lengths is almost straight. These findings could shed some light on the role of DNA elasticity in the chromatin structure.

Free-energy landscape of mono- and dinucleosomes: Enhanced rotational flexibility of interconnected nucleosomes

The nucleosome represents the basic unit of eukaryotic genome organization, and its conformational fluctuations play a crucial role in various cellular processes. Here we provide insights into the flipping transition of a nucleosome by computing its free-energy landscape as a function of the linking number and nucleosome orientation using the density-of-states Monte Carlo approach. To investigate how the energy landscape is affected by the presence of neighboring nucleosomes in a chromatin fiber, we also compute the free-energy landscape for a dinucleosome array. We find that the mononucleosome is bistable between conformations with negatively and positively crossed linkers while the conformation with open linkers appears as a transition state. The dinucleosome exhibits a markedly different energy landscape in which the conformation with open linkers populates not only the transition state but also the global minimum. This enhanced stability of the open state is attributed to increased rotational flexibility of nucleosomes arising from their mechanical coupling with neighboring nucleosomes. Our results provide a possible mechanism by which chromatin may enhance the accessibility of its DNA and facilitate the propagation and mitigation of DNA torsional stresses.

Topological polymorphism of the two-start chromatin fiber

Specific details concerning the spatial organization of nucleosomes in 30 nm fibers remain unknown. To investigate this, we analyzed all stereochemically possible configurations of two-start nucleosome fibers with short DNA linkers L = 13-37 bp (nucleosome repeat length (NRL) = 160-184 bp). Four superhelical parameters-inclination of nucleosomes, twist, rise, and diameter-uniquely describe a regular symmetric fiber. The energy of a fiber is defined as the sum of four terms: elastic energy of the linker DNA, steric repulsion, electrostatics, and a phenomenological (H4 tail-acidic patch) interaction between two stacked nucleosomes. By optimizing the fiber energy with respect to the superhelical parameters, we found two types of topological transition in fibers (associated with the change in inclination angle): one caused by an abrupt 360° change in the linker DNA twisting (change in the DNA linking number, ΔLk = 1), and another caused by overcrossing of the linkers (ΔLk = 2). To the best of our knowledge, this topological polymorphism of the two-start fibers was not reported in the computations published earlier. Importantly, the optimal configurations of the fibers with linkers L = 10n and 10n + 5 bp are characterized by different values of the DNA linking number-that is, they are topologically different. Our results are consistent with experimental observations, such as the inclination 60° to 70° (the angle between the nucleosomal disks and the fiber axis), helical rise, diameter, and left-handedness of the fibers. In addition, we make several testable predictions, among them different degrees of DNA supercoiling in fibers with L = 10n and 10n + 5 bp, different flexibility of the two types of fibers, and a correlation between the local NRL and the level of transcription in different parts of the yeast genome.