Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Atomic Description of the Reciprocal Action between Supercoils and Melting Bubbles on Linear DNA

Although the mechanical response of DNA to physiological torsion and tension is well characterized, the detailed structures are not yet known. By using molecular dynamics simulations on linear DNA with 300 base-pairs, we provide, for the first time, the conformational phase diagram at atomic resolution. Our simulations also reveal the dynamics and diffusion of supercoils. We observe a new state in negative supercoiling, where denaturation bubbles form in adenine/thymine-rich regions independently of the underlying DNA topology. We thus propose sequence-dependent bubbles could position plectonemes in longer DNA.

A high throughput single molecule platform to study DNA supercoiling effect on protein-DNA interactions

DNA supercoiling significantly influences DNA metabolic pathways. To examine its impact on DNA-protein interactions at the single-molecule level, we developed a highly efficient and reliable protocol to modify plasmid DNA at specific sites, allowing us to label plasmids with fluorophores and biotin. We then induced negative and positive supercoiling in these plasmids using gyrase and reverse gyrase, respectively. Comparing supercoiled DNA with relaxed circular DNA, we assessed the effects of supercoiling on CRISPR-Cas9 and mismatch repair protein MutS. We found that negative DNA supercoiling exacerbates off-target effects in DNA unwinding by Cas9. For MutS, we observed both negative and positive DNA supercoiling enhances the binding interaction between MutS and a mismatched base pair but does not affect the rate of ATP-induced sliding clamp formation. These findings not only underscore the versatility of our protocol but also opens new avenues for exploring the intricate dynamics of protein-DNA interactions under the influences of supercoiling.

Single-molecule visualization of twin-supercoiled domains generated during transcription

Transcription-coupled supercoiling of DNA is a key factor in chromosome compaction and the regulation of genetic processes in all domains of life. It has become common knowledge that, during transcription, the DNA-dependent RNA polymerase (RNAP) induces positive supercoiling ahead of it (downstream) and negative supercoils in its wake (upstream), as rotation of RNAP around the DNA axis upon tracking its helical groove gets constrained due to drag on its RNA transcript. Here, we experimentally validate this so-called twin-supercoiled-domain model with in vitro real-time visualization at the single-molecule scale. Upon binding to the promoter site on a supercoiled DNA molecule, RNAP merges all DNA supercoils into one large pinned plectoneme with RNAP residing at its apex. Transcription by RNAP in real time demonstrates that up- and downstream supercoils are generated simultaneously and in equal portions, in agreement with the twin-supercoiled-domain model. Experiments carried out in the presence of RNases A and H, revealed that an additional viscous drag of the RNA transcript is not necessary for the RNAP to induce supercoils. The latter results contrast the current consensus and simulations on the origin of the twin-supercoiled domains, pointing at an additional mechanistic cause underlying supercoil generation by RNAP in transcription.

Correlated Hybrid DNA Structures Explored by the oxDNA Model

Thermodynamically, perfect DNA hybridization can be formed between probes and their corresponding targets due to the favorable energy. However, this is not the case dynamically. Here, we use molecular dynamics (MD) simulations based on the oxDNA model to investigate the process of DNA microarray hybridization. In general, correlated hybrid DNA structures are formed, including one probe associated with several targets as well as one target hybrid with multiple probes leading to the target-mediated hybridization. The formation of these two types of correlated structures largely depends on the surface coverage of the DNA microarray. Moreover, DNA sequence, DNA length, and spacer length have an impact on the structural formation. Our findings shed light on the dynamics of DNA hybridization, which is important for the application of DNA microarray.

Equilibrium melting probabilities of a DNA molecule with a defect: An exact solution of the Poland-Scheraga model

In this study we derive analytically the equilibrium melting probabilities for basepairs of a DNA molecule with a defect site. We assume that the defect is characterized by a change in the Watson-Crick basepair energy of the defect basepair, and in the associated two stacking energies for the defect, as compared to the remaining parts of the DNA. The defect site could, for instance, occur due to DNA basepair mismatching, cross-linking, or by the chemical modifications when attaching fluorescent labels, such as fluorescent-quencher pairs, to DNA. Our exact solution of the Poland-Scheraga model for DNA melting provides the probability that the labeled basepair, and its neighbors, are open at different temperatures. Our work is of direct importance, for instance, for studies where fluorophore-quencher pairs are used for studying single basepair fluctuations of designed DNA molecules.

Chromatinization Modulates Topoisomerase II Processivity

Type IIA topoisomerases are essential DNA processing enzymes that must robustly and reliably relax DNA torsional stress in vivo. While cellular processes constantly create different degrees of torsional stress, how this stress feeds back to control type IIA topoisomerase function remains obscure. Using a suite of single-molecule approaches, we examined the torsional impact on supercoiling relaxation of both naked DNA and chromatin by eukaryotic topoisomerase II (topo II). We observed that topo II was at least ~ 50-fold more processive on plectonemic DNA than previously estimated, capable of relaxing > 6000 turns. We further discovered that topo II could relax supercoiled DNA prior to plectoneme formation, but with a ~100-fold reduction in processivity; strikingly, the relaxation rate in this regime decreased with diminishing torsion in a manner consistent with the capture of transient DNA loops by topo II. Chromatinization preserved the high processivity of the enzyme under high torsional stress. Interestingly, topo II was still highly processive (~ 1000 turns) even under low torsional stress, consistent with the predisposition of chromatin to readily form DNA crossings. This work establishes that chromatin is a major stimulant of topo II function, capable of enhancing function even under low torsional stress.

Two Different Ways of Stress Release in Supercoiled DNA Minicircles under DNA Nick

It is generally believed that DNA nick is an effective way to release stress in supercoiled DNA, resulting from the twisting motion that individual strands rotate around the axis of the DNA helix. Here, we use MD simulations based on the oxDNA model to investigate the relaxation of 336 bp supercoiled minicircular DNA under DNA nick. Our simulations show that stress release, characterized by the abrupt decrease in linking number, may be induced by two types of DNA motion depending on the nick position. Except for the twisting motion, there is a writhing motion, that is, double strands collectively rotating with one plectoneme removal, which may occur in the process of DNA relaxation with the nick position in the loop region. Moreover, the writhing motion is more likely to occur in the DNA with relatively high hardness, such as C-G pairs. Our simulation results uncover the relationship between structural transformation, stress release, and DNA motion during the dynamic process under DNA nick, indicating the influence of nick position on the relaxation of the supercoiled DNA.

Equilibrium fluctuations of DNA plectonemes

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e., over- or underwound. They are ubiquitous in cellular DNA, and their physical properties have attracted significant interest both from the experimental side and from the modeling side. In this paper, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched phase and a plectonemic phase. A variety of mechanisms are found to contribute to extension fluctuations, characterized by the variance 〈Δz^{2}〉. We find the dominant contribution to 〈Δz^{2}〉 to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the twistable wormlike chain and analyze the fluctuation of various quantities, the results of which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations.

The interplay of supercoiling and thymine dimers in DNA

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

Two-phase dynamics of DNA supercoiling based on DNA polymer physics

DNA supercoils are generated in genome regulation processes such as transcription and replication and provide mechanical feedback to such processes. Under tension, a DNA supercoil can present a coexistence state of plectonemic and stretched phases. Experiments have revealed the dynamic behaviors of plectonemes, e.g., diffusion, nucleation, and hopping. To represent these dynamics with conformational changes, we demonstrated first the fast dynamics on the DNA to reach torque equilibrium within the plectonemic and stretched phases, and then identified the two-phase boundaries as collective slow variables to describe the essential dynamics. According to the timescale separation demonstrated here, we developed a two-phase model on the dynamics of DNA supercoiling, which can capture physiologically relevant events across timescales of several orders of magnitudes. In this model, we systematically characterized the slow dynamics between the two phases and compared the numerical results with those from the DNA polymer physics-based worm-like chain model. The supercoiling dynamics, including the nucleation, diffusion, and hopping of plectonemes, have been well represented and reproduced, using the two-phase dynamic model, at trivial computational costs. Our current developments, therefore, can be implemented to explore multiscale physical mechanisms of the DNA supercoiling-dependent physiological processes.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Energetics of twisted DNA topologies

Our goal is to review the main theoretical models used to calculate free energy changes associated with common, torsion-induced conformational changes in DNA and provide the resulting equations hoping to facilitate quantitative analysis of both in vitro and in vivo studies. This review begins with a summary of work regarding the energy change of the negative supercoiling-induced B- to L-DNA transition, followed by a discussion of the energetics associated with the transition to Z-form DNA. Finally, it describes the energy changes associated with the formation of DNA curls and plectonemes, which can regulate DNA-protein interactions and promote cross talk between distant DNA elements, respectively. The salient formulas and parameters for each scenario are summarized in table format to facilitate comparison and provide a concise, user-friendly resource.