Extension of torsionally stressed DNA by external force

Torsional Mechanics of Circular DNA

Circular DNA found in the cell is actively regulated to an underwound state, with their superhelical density close to σ ~ - 0 . 06 . While this underwound state is essential to life, how it impacts the torsional mechanical properties of DNA is not fully understood. In this work, we performed simulations to understand the torsional mechanics of circular DNA and validated our results with single-molecule measurements and analytical theory. We found that the torque generated at σ ~ - 0 . 06 is near but slightly below that required to melt DNA, significantly decreasing the energy barrier for proteins that interact with melted DNA. Furthermore, supercoiled circular DNA experiences force (tension) and torque that are equally distributed through the DNA contour. We have also extended a previous analytical framework to show how the plectonemic twist persistence length depends on the intrinsic bending persistence length and twist persistence length. Our work establishes a framework for understanding DNA supercoiling and torsional dynamics of circular DNA.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Chromatinization modulates topoisomerase II processivity

Type IIA topoisomerases are essential DNA processing enzymes that must robustly and reliably relax DNA torsional stress. While cellular processes constantly create varying torsional stress, how this variation impacts type IIA topoisomerase function remains obscure. Using multiple single-molecule approaches, we examined the torsional dependence of eukaryotic topoisomerase II (topo II) activity on naked DNA and chromatin. We observed that topo II is ~50-fold more processive on buckled DNA than previously estimated. We further discovered that topo II relaxes supercoiled DNA prior to plectoneme formation, but with processivity reduced by ~100-fold. This relaxation decreases with diminishing torsion, consistent with topo II capturing transient DNA loops. Topo II retains high processivity on buckled chromatin (~10,000 turns) and becomes highly processive even on chromatin under low torsional stress (~1000 turns), consistent with chromatin's predisposition to readily form DNA crossings. This work establishes that chromatin is a major stimulant of topo II function.

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.

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.

DNA fluctuations reveal the size and dynamics of topological domains

DNA supercoiling is a key regulatory mechanism that orchestrates DNA readout, recombination, and genome maintenance. DNA-binding proteins often mediate these processes by bringing two distant DNA sites together, thereby inducing (transient) topological domains. In order to understand the dynamics and molecular architecture of protein-induced topological domains in DNA, quantitative and time-resolved approaches are required. Here, we present a methodology to determine the size and dynamics of topological domains in supercoiled DNA in real time and at the single-molecule level. Our approach is based on quantifying the extension fluctuations-in addition to the mean extension-of supercoiled DNA in magnetic tweezers (MT). Using a combination of high-speed MT experiments, Monte Carlo simulations, and analytical theory, we map out the dependence of DNA extension fluctuations as a function of supercoiling density and external force. We find that in the plectonemic regime, the extension variance increases linearly with increasing supercoiling density and show how this enables us to determine the formation and size of topological domains. In addition, we demonstrate how the transient (partial) dissociation of DNA-bridging proteins results in the dynamic sampling of different topological states, which allows us to deduce the torsional stiffness of the plectonemic state and the kinetics of protein-plectoneme interactions. We expect our results to further the understanding and optimization of magnetic tweezer measurements and to enable quantification of the dynamics and reaction pathways of DNA processing enzymes in the context of physiologically relevant forces and supercoiling densities.

Statistics and topology of fluctuating ribbons

Ribbons are a class of slender structures whose length, width, and thickness are widely separated from each other. This scale separation gives a ribbon unusual mechanical properties in athermal macroscopic settings, for example, it can bend without twisting, but cannot twist without bending. Given the ubiquity of ribbon-like biopolymers in biology and chemistry, here we study the statistical mechanics of microscopic inextensible, fluctuating ribbons loaded by forces and torques. We show that these ribbons exhibit a range of topologically and geometrically complex morphologies exemplified by three phases-a twist-dominated helical phase (HT), a writhe-dominated helical phase (HW), and an entangled phase-that arise as the applied torque and force are varied. Furthermore, the transition from HW to HT phases is characterized by the spontaneous breaking of parity symmetry and the disappearance of perversions (that correspond to chirality-reversing localized defects). This leads to a universal response curve of a topological quantity, the link, as a function of the applied torque that is similar to magnetization curves in second-order phase transitions.

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.

DNA Mechanics and Topology

We review the current understanding of the mechanics of DNA and DNA-protein complexes, from scales of base pairs up to whole chromosomes. Mechanics of the double helix as revealed by single-molecule experiments will be described, with an emphasis on the role of polymer statistical mechanics. We will then discuss how topological constraints- entanglement and supercoiling-impact physical and mechanical responses. Models for protein-DNA interactions, including effects on polymer properties of DNA of DNA-bending proteins will be described, relevant to behavior of protein-DNA complexes in vivo. We also discuss control of DNA entanglement topology by DNA-lengthwise-compaction machinery acting in concert with topoisomerases. Finally, the chapter will conclude with a discussion of relevance of several aspects of physical properties of DNA and chromatin to oncology.

Twist-bend coupling and the statistical mechanics of the twistable wormlike-chain model of DNA: Perturbation theory and beyond

The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling G between twisting and bending, known to arise from the groove asymmetry of the DNA double helix. Here the effect of G on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness C_{eff} for small twist-bend coupling. We find that the "bare" G is "screened" by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with G=0 but with long-wavelength bending and twisting rigidities that are shifted by G-dependent amounts. Using results for torsional and bending rigidities for freely fluctuating DNA, we show how our perturbative results can be extended to a nonperturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo "triad" and molecular dynamics "oxDNA" models, characterized by different degrees of coarse graining, validating the perturbative and nonperturbative analyses. While our theory is in generally good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely used coarse-grained models.

Overtwisting induces polygonal shapes in bent DNA

By combining analytical results and simulations of various coarse-grained models, we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternation of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, rather than the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants.

Simulation of DNA Supercoil Relaxation

Several recent single-molecule experiments observe the response of supercoiled DNA to nicking endonucleases and topoisomerases. Typically in these experiments, indirect measurements of supercoil relaxation are obtained by observing the motion of a large micron-sized bead. The bead, which also serves to manipulate DNA, experiences significant drag and thereby obscures supercoil dynamics. Here we employ our discrete wormlike chain model to bypass experimental limitations and simulate the dynamic response of supercoiled DNA to a single strand nick. From our simulations, we make three major observations. First, extension is a poor dynamic measure of supercoil relaxation; in fact, the linking number relaxes so fast that it cannot have much impact on extension. Second, the rate of linking number relaxation depends upon its initial partitioning into twist and writhe as determined by tension. Third, the extensional response strongly depends upon the initial position of plectonemes.

Twist-Bend Coupling and the Torsional Response of Double-Stranded DNA

Recent magnetic tweezers experiments have reported systematic deviations of the twist response of double-stranded DNA from the predictions of the twistable wormlike chain model. Here we show, by means of analytical results and computer simulations, that these discrepancies can be resolved if a coupling between twist and bend is introduced. We obtain an estimate of 40±10  nm for the twist-bend coupling constant. Our simulations are in good agreement with high-resolution, magnetic-tweezers torque data. Although the existence of twist-bend coupling was predicted long ago [J. Marko and E. Siggia, Macromolecules 27, 981 (1994)MAMOBX0024-929710.1021/ma00082a015], its effects on the mechanical properties of DNA have been so far largely unexplored. We expect that this coupling plays an important role in several aspects of DNA statics and dynamics.

Theoretical Methods for Studying DNA Structural Transitions under Applied Mechanical Constraints

Recent progress in single-molecule manipulation technologies has made it possible to exert force and torque on individual DNA biopolymers to probe their mechanical stability and interaction with various DNA-binding proteins. It was revealed in these experiments that the DNA structure and formation of nucleoprotein complexes by DNA-architectural proteins can be strongly modulated by an intricate interplay between the entropic elasticity of DNA and its global topology, which is closely related to the mechanical constraints applied to the DNA. Detailed understanding of the physical processes underlying the DNA behavior observed in single-molecule experiments requires the development of a general theoretical framework, which turned out to be a rather challenging task. Here, we review recent advances in theoretical methods that can be used to interpret single-molecule manipulation experiments on DNA.

Modeling Bacterial DNA: Simulation of Self-Avoiding Supercoiled Worm-Like Chains Including Structural Transitions of the Helix

Under supercoiling constraints, naked DNA, such as a large part of bacterial DNA, folds into braided structures called plectonemes. The double-helix can also undergo local structural transitions, leading to the formation of denaturation bubbles and other alternative structures. Various polymer models have been developed to capture these properties, with Monte-Carlo (MC) approaches dedicated to the inference of thermodynamic properties. In this chapter, we explain how to perform such Monte-Carlo simulations, following two objectives. On one hand, we present the self-avoiding supercoiled Worm-Like Chain (ssWLC) model, which is known to capture the folding properties of supercoiled DNA, and provide a detailed explanation of a standard MC simulation method. On the other hand, we explain how to extend this ssWLC model to include structural transitions of the helix.

Effect of sequence-dependent rigidity on plectoneme localization in dsDNA

We use Monte-Carlo simulations to study the effect of variable rigidity on plectoneme formation and localization in supercoiled double-stranded DNA. We show that the presence of soft sequences increases the number of plectoneme branches and that the edges of the branches tend to be localized at these sequences. We propose an experimental approach to test our results in vitro, and discuss the possible role played by plectoneme localization in the search process of transcription factors for their targets (promoter regions) on the bacterial genome.

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.

Transfer-matrix calculations of DNA polymer micromechanics under tension and torque constraints

Recent development of single-molecule manipulation technologies has made it possible to exert constant force and torque on individual DNA biopolymers to probe their elastic characteristics and structural stability. It has been previously shown that depending on the nature of applied mechanical constraints, DNA can exist in several forms including B-, L-, and P-DNA. However, there is still a lack of understanding of how structural heterogeneity of DNA, which may naturally arise due to sequence-dependent DNA properties, protein binding, or DNA damage, influences local stability of the above DNA states. To provide a more complete and detailed description of the DNA mechanics, we developed a theoretical framework based on transfer-matrix calculations and demonstrated how it can be used to predict the DNA behavior upon application of a wide range of force and torque constraints. The resulting phase diagram shows DNA structural transitions that are in good agreement with previous experimental and theoretical studies. We further discuss how the constructed formalism can be extended to include local inhomogeneities in the DNA physical properties, thus making it possible to investigate the effect of DNA sequence as well as protein binding on DNA structural stability.

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.

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

Physical origin of DNA unzipping

In DNA transcription, the base pairs are unzipped in response to the enzymatic forces, separating apart two intertwined nucleotide strands. Consequently, the double-stranded DNA (dsDNA), in which two nucleotide strands wind about each other, transits structurally to the single-stranded DNA (ssDNA) in which two nucleotide strands are completely unwound and separated. The large interstrand separation is intimately related to the softening nucleotide strands. This conceptual framework is reinforced with the flow of the bending modulus toward zero under recursion relations derived from the momentum shell renormalization group. Interestingly, the stretch modulus remains the same under recursion relations. The renormalization of the bending modulus to zero has a profound implication that ssDNA has the shorter bending persistence length than does dsDNA in accordance with experiments.

Single-molecule portrait of DNA and RNA double helices

The composition and geometry of the genetic information carriers were described as double-stranded right helices sixty years ago. The flexibility of their sugar-phosphate backbones and the chemistry of their nucleotide subunits, which give rise to the RNA and DNA polymers, were soon reported to generate two main structural duplex states with biological relevance: the so-called A and B forms. Double-stranded (ds) RNA adopts the former whereas dsDNA is stable in the latter. The presence of flexural and torsional stresses in combination with environmental conditions in the cell or in the event of specific sequences in the genome can, however, stabilize other conformations. Single-molecule manipulation, besides affording the investigation of the elastic response of these polymers, can test the stability of their structural states and transition models. This approach is uniquely suited to understanding the basic features of protein binding molecules, the dynamics of molecular motors and to shedding more light on the biological relevance of the information blocks of life. Here, we provide a comprehensive single-molecule analysis of DNA and RNA double helices in the context of their structural polymorphism to set a rigorous interpretation of their material response both inside and outside the cell. From early knowledge of static structures to current dynamic investigations, we review their phase transitions and mechanochemical behaviour and harness this fundamental knowledge not only through biological sciences, but also for Nanotechnology and Nanomedicine.

Biophysics of protein-DNA interactions and chromosome organization

The function of DNA in cells depends on its interactions with protein molecules, which recognize and act on base sequence patterns along the double helix. These notes aim to introduce basic polymer physics of DNA molecules, biophysics of protein-DNA interactions and their study in single-DNA experiments, and some aspects of large-scale chromosome structure. Mechanisms for control of chromosome topology will also be discussed.

Efficient deformation algorithm for plasmid DNA simulations

BACKGROUND

Plasmid DNA molecules are closed circular molecules that are widely used in life sciences, particularly in gene therapy research. Monte Carlo methods have been used for several years to simulate the conformational behavior of DNA molecules. In each iteration these simulation methods randomly generate a new trial conformation, which is either accepted or rejected according to a criterion based on energy calculations and stochastic rules. These simulation trials are generated using a method based on crankshaft motion that, apart from some slight improvements, has remained the same for many years.

RESULTS

In this paper, we present a new algorithm for the deformation of plasmid DNA molecules for Monte Carlo simulations. The move underlying our algorithm preserves the size and connectivity of straight-line segments of the plasmid DNA skeleton. We also present the results of three experiments comparing our deformation move with the standard and biased crankshaft moves in terms of acceptance ratio of the trials, energy and temperature evolution, and average displacement of the molecule. Our algorithm can also be used as a generic geometric algorithm for the deformation of regular polygons or polylines that preserves the connections and lengths of their segments.

CONCLUSION

Compared with both crankshaft moves, our move generates simulation trials with higher acceptance ratios and smoother deformations, making it suitable for real-time visualization of plasmid DNA coiling. For that purpose, we have adopted a DNA assembly algorithm that uses nucleotides as building blocks.

Torsional behavior of chromatin is modulated by rotational phasing of nucleosomes

Torsionally stressed DNA plays a critical role in genome organization and regulation. While the effects of torsional stresses on naked DNA have been well studied, little is known about how these stresses propagate within chromatin and affect its organization. Here we investigate the torsional behavior of nucleosome arrays by means of Brownian dynamics simulations of a coarse-grained model of chromatin. Our simulations reveal a strong dependence of the torsional response on the rotational phase angle Ψ₀ between adjacent nucleosomes. Extreme values of Ψ₀ lead to asymmetric, bell-shaped extension-rotation profiles with sharp maxima shifted toward positive or negative rotations, depending on the sign of Ψ₀, and to fast, irregular propagation of DNA twist. In contrast, moderate Ψ₀ yield more symmetric profiles with broad maxima and slow, uniform propagation of twist. The observed behavior is shown to arise from an interplay between nucleosomal transitions into states with crossed and open linker DNAs and global supercoiling of arrays into left- and right-handed coils, where Ψ₀ serves to modulate the energy landscape of nucleosomal states. Our results also explain the torsional resilience of chromatin, reconcile differences between experimentally measured extension-rotation profiles, and suggest a role of torsional stresses in regulating chromatin assembly and organization.

Probing the elasticity of DNA on short length scales by modeling supercoiling under tension

The wormlike-chain (WLC) model is widely used to describe the energetics of DNA bending. Motivated by recent experiments, alternative, so-called subelastic chain models were proposed that predict a lower elastic energy of highly bent DNA conformations. Until now, no unambiguous verification of these models has been obtained because probing the elasticity of DNA on short length scales remains challenging. Here we investigate the limits of the WLC model using coarse-grained Monte Carlo simulations to model the supercoiling of linear DNA molecules under tension. At a critical supercoiling density, the DNA extension decreases abruptly due to the sudden formation of a plectonemic structure. This buckling transition is caused by the large energy required to form the tightly bent end-loop of the plectoneme and should therefore provide a sensitive benchmark for model evaluation. Although simulations based on the WLC energetics could quantitatively reproduce the buckling measured in magnetic tweezers experiments, the buckling almost disappears for the tested linear subelastic chain model. Thus, our data support the validity of a harmonic bending potential even for small bending radii down to 3.5 nm.

Torque transfer coefficient in DNA under torsional stress

In recent years, significant progress in understanding the properties of supercoiled DNA has been obtained due to nanotechniques that made stretching and twisting of single molecules possible. Quantitative interpretation of such experiments requires accurate knowledge of torques inside manipulated DNA. This paper argues that it is not possible to transfer the entire magnitudes of external torques to the twisting stress of the double helix, and that a reducing torque transfer coefficient (TTC < 1) should always be assumed. This assertion agrees with simple physical intuition and is supported by the results of all-atom molecular dynamics (MD) simulations. According to MD, the TTCs around 0.8 are observed in nearly optimal conditions. Reaching higher values requires special efforts and it should be difficult in practice. The TTC can be partially responsible for the persistent discrepancies between the twisting rigidity of DNA measured by different methods.

On binding of DNA-bending proteins to DNA minicircles

We present a theoretical study of binding of DNA-bending proteins to circular DNA, using computer simulations of the wormlike chain model of DNA. We find that the binding affinity is affected by the bending elasticity and the conformational entropy of the polymer and that while protein adsorption is identical on open and closed long DNA molecules, there is significant enhancement of binding on DNA minicircles, compared to their linear counterparts. We also find that the ratio of the radii of gyration of open and closed chains depends on protein concentration for short DNA molecules. Experimental tests of our predictions are proposed.

Biopolymer elasticity: Mechanics and thermal fluctuations

We present an analytical study of the role of thermal fluctuations in shaping molecular elastic properties of semiflexible polymers. Our study interpolates between mechanics and statistical mechanics in a controlled way and shows how thermal fluctuations modify the elastic properties of biopolymers. We present a study of the minimum-energy configurations with explicit expressions for their energy and writhe and plots of the extension versus link for these configurations and a study of fluctuations around the local minima of energy and approximate analytical formulas for the free energy of stretched twisted polymers. The central result of our study is a closed-form expression for the leading thermal fluctuation correction to the free energy around the nonperturbative writhing family solution for the configuration of a biopolymer. From the derived formulas, the predictions of the wormlike chain model for molecular elasticity can be worked out for a comparison against numerical simulations and experiments.

DNA looping by FokI: the impact of twisting and bending rigidity on protein-induced looping dynamics

Protein-induced DNA looping is crucial for many genetic processes such as transcription, gene regulation and DNA replication. Here, we use tethered-particle motion to examine the impact of DNA bending and twisting rigidity on loop capture and release, using the restriction endonuclease FokI as a test system. To cleave DNA efficiently, FokI bridges two copies of an asymmetric sequence, invariably aligning the sites in parallel. On account of the fixed alignment, the topology of the DNA loop is set by the orientation of the sites along the DNA. We show that both the separation of the FokI sites and their orientation, altering, respectively, the twisting and the bending of the DNA needed to juxtapose the sites, have profound effects on the dynamics of the looping interaction. Surprisingly, the presence of a nick within the loop does not affect the observed rigidity of the DNA. In contrast, the introduction of a 4-nt gap fully relaxes all of the torque present in the system but does not necessarily enhance loop stability. FokI therefore employs torque to stabilise its DNA-looping interaction by acting as a ‘torsional’ catch bond.

Revisiting polymer statistical physics to account for the presence of long-range-correlated structural disorder in 2D DNA chains

We elaborate on a generalization of the 2D wormlike chain (WLC) model that accounts for the presence of long-range correlations (LRC) in the intrinsic curvature distribution of eukaryotic DNA. This model predicts some decrease of the DNA persistence length resulting from some large-scale intrinsic curvature induced by sequence-dependent persistent random distribution of local bending sites. When assisting exact analytical calculations by numerical DNA simulations, we show that the conjugated contributions of i) the thermal curvature fluctuations characterized by the "dynamic" persistence length ℓ(p)(d) = 2A, where A is the elastic bending modulus, and ii) the intrinsic LRC curvature disorder of amplitude σ(o) and Hurst exponent H > 1/2, characterized by a "static" persistence length ℓ(p)(H) = A(1/2H)σ(o)(-1/H) Γ(1/2H + 1), can be described by a continuum of generalized WLC (GWLC) models parametrized by the LRC exponent H. We use perturbation analysis to investigate the two limiting cases of weak static disorder (w(H) << 1 and weak dynamical fluctuations (1/w (H) << 1), where w(H) = l(p)(d)/l(p)(H) is a dimensionless parameter. From a quantitative point of view, our study demonstrates that even for a small value of the LRC (H approximately equal 0.6-0.8) static disorder amplitude σ(o) ~ 10(-2), as previously reported for genomic DNA, the decrease of the persistence length from the WLC prediction l(p)(d) can be very significant, up to twofold. The implications of these results on the first steps of compaction of DNA in eukaryotic cells are discussed.

A multiscale dynamic model of DNA supercoil relaxation by topoisomerase IB

In this study, we report what we believe to be the first multiscale simulation of the dynamic relaxation of DNA supercoils by human topoisomerase IB (topo IB). We leverage our previous molecular dynamics calculations of the free energy landscape describing the interaction between a short DNA fragment and topo IB. Herein, this landscape is used to prescribe boundary conditions for a computational, elastodynamic continuum rod model of a long length of supercoiled DNA. The rod model, which accounts for the nonlinear bending, twisting, and electrostatic interaction of the (negatively charged) DNA backbone, is extended to include the hydrodynamic drag induced by the surrounding physiological buffer. Simulations for a 200-bp-long DNA supercoil in complex with topo IB reveal a relaxation timescale of ∼0.1-1.0 μs. The relaxation follows a sequence of cascading reductions in the supercoil linking number (Lk), twist (Tw), and writhe (Wr) that follow companion cascading reductions in the supercoil elastic and electrostatic energies. The novel (to our knowledge) multiscale modeling method may enable simulations of the entire experimental setup that measures DNA supercoiling and relaxation via single molecule magnetic trapping.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Scaling of Linking and Writhing Numbers for Spherically Confined and Topologically Equilibrated Flexible Polymers

Scaling laws for Gauss linking number Ca and writhing number Wr for spherically confined flexible polymers with thermally fluctuating topology are analyzed. For ideal (phantom) polymers each of N segments of length unity confined to a spherical pore of radius R there are two scaling regimes: for sufficiently weak confinement (R ⪢ N(1/3)) each chain has |Wr| ≈ N(1/2), and each pair of chains has average |Ca| ≈ N/R(3/2); alternately for sufficiently tight confinement (N(1/3) ⪢ R), |Wr| ≈ |CA| ≈ N/R(3/2). Adding segment-segment avoidance modifies this result: for n chains with excluded volume interactions |Ca| ≈ (N/n)(1/2)f(ϕ) where f is a scaling function that depends approximately linearly on the segment concentration ϕ = nN/R(3). Scaling results for writhe are used to estimate the maximum writhe of a polymer; this is demonstrated to be realizable through a writhing instability that occurs for a polymer which is able to change knotting topology and which is subject to an applied torque. Finally, scaling results for linking are used to estimate bounds on the entanglement complexity of long chromosomal DNA molecules inside cells, and to show how "lengthwise" chromosome condensation can suppress DNA entanglement.

Nucleation at the DNA supercoiling transition

Twisting DNA under a constant applied force reveals a thermally activated transition into a state with a supercoiled structure known as a plectoneme. Using transition-state theory, we predict the rate of this plectoneme nucleation to be of order 10(4) Hz. We reconcile this with experiments that have measured hopping rates of order 10 Hz by noting that the viscous drag on the bead used to manipulate the DNA limits the measured rate. We find that the intrinsic bending caused by disorder in the base-pair sequence is important for understanding the free-energy barrier that governs the transition. Both analytic and numerical methods are used in the calculations. We provide extensive details on the numerical methods for simulating the elastic rod model with and without disorder.

Effect of spontaneous twist on DNA minicircles

Monte Carlo simulations are used to study the effect of spontaneous (intrinsic) twist on the conformation of topologically equilibrated minicircles of dsDNA. The twist, writhe, and radius of gyration distributions and their moments are calculated for different spontaneous twist angles and DNA lengths. The average writhe and twist deviate in an oscillatory fashion (with the period of the double helix) from their spontaneous values, as one spans the range between two neighboring integer values of intrinsic twist. Such deviations vanish in the limit of long DNA plasmids.

Energetics at the DNA supercoiling transition

Twisting a DNA molecule held under constant tension is accompanied by a transition from a linear to a plectonemic DNA configuration, in which part of the applied twist is absorbed in a superhelical structure. Recent experiments revealed the occurrence of an abrupt extension change at the onset of this transition. To elucidate its origin we study this abrupt DNA shortening using magnetic tweezers. We find that it strongly depends on the length of the DNA molecule and the ionic strength of the solution. This behavior can be well understood in the framework of a model in which the energy per writhe for the initial plectonemic loop is larger than for subsequent turns of the superhelix. By quantitative data analysis, relevant plectoneme energies and other parameters were extracted, providing good agreement with a simple theory. As a direct confirmation of the initial-loop model, we find that for a kinked DNA molecule the abrupt extension change occurs at significantly lower twist than the subsequent superhelix formation. This should allow pinning of the plectoneme position within supercoiled DNA if a kinked substrate is used, and enable the detection of enzymes and proteins which, themselves, bend or kink DNA.

DNA-DNA interactions in tight supercoils are described by a small effective charge density

DNA-DNA interactions are important for genome compaction and transcription regulation. In studies of such complex processes, DNA is often modeled as a homogeneously charged cylinder and its electrostatic interactions are calculated within the framework of the Poisson-Boltzmann equation. Commonly, a charge adaptation factor is used to address limitations of this theoretical approach. Despite considerable theoretical and experimental efforts, a rigorous quantitative assessment of this parameter is lacking. Here, we comprehensively characterized DNA-DNA interactions in the presence of monovalent ions by analyzing the supercoiling behavior of single DNA molecules held under constant tension. Both a theoretical model and coarse-grained simulations of this process revealed a surprisingly small effective DNA charge of 40% of the nominal charge density, which was additionally supported by all-atom molecular dynamics simulations.

Linking topology of large DNA molecules

Double helix DNA molecules, the carriers of genetic instructions in cells, are strongly affected by their topological properties. Two distinct and biologically important types of linking are associated with double helix DNAs: 'internal' linking of the two strands of individual double helices, and 'external' linking of separate double helix DNAs. Constraint of internal linking gives rise to internal torsional stress and supercoiling of circular DNAs. External linking is a likely outcome of DNA replication, and must be eliminated by the cell in order to separate duplicated DNAs. I outline some of the physics and biology connected with both internal and external linking.

Discontinuities at the DNA supercoiling transition

While slowly turning the ends of a single molecule of DNA at constant applied force, a discontinuity was recently observed at the supercoiling transition when a small plectoneme is suddenly formed. This can be understood as an abrupt transition into a state in which stretched and plectonemic DNA coexist. We argue that there should be discontinuities in both the extension and the torque at the transition and provide experimental evidence for both. To predict the sizes of these discontinuities and how they change with the overall length of DNA, we organize a phenomenological theory for the coexisting plectonemic state in terms of four parameters. We also test supercoiling theories, including our own elastic rod simulation, finding discrepancies with experiment that can be understood in terms of the four coexisting state parameters.

Elasticity and electrostatics of plectonemic DNA

We present a self-contained theory for the mechanical response of DNA in single molecule experiments. Our model is based on a one-dimensional continuum description of the DNA molecule and accounts both for its elasticity and for DNA-DNA electrostatic interactions. We consider the classical loading geometry used in experiments where one end of the molecule is attached to a substrate and the other one is pulled by a tensile force and twisted by a given number of turns. We focus on configurations relevant to the limit of a large number of turns, which are made up of two phases, one with linear DNA and the other one with superhelical DNA. The model takes into account thermal fluctuations in the linear phase and electrostatic interactions in the superhelical phase. The values of the torsional stress, of the supercoiling radius and angle, and key features of the experimental extension-rotation curves, namely the slope of the linear region and thermal buckling threshold, are predicted. They are found in good agreement with experimental data.

Measurement of the torque on a single stretched and twisted DNA using magnetic tweezers

We deduced the torque applied on a single stretched and twisted DNA by integrating the change in the molecule's extension with respect to force as it is coiled. While consistent with previous direct measurements of the torque at high forces (F>1 pN), this method, which is simple and does not require a sophisticated setup, allows for lower force estimates. We used this approach to deduce the effective torsional modulus of DNA, which decreases with force, and to estimate the buckling torque of DNA as a function of force in various salt conditions.

Topoisomerase V relaxes supercoiled DNA by a constrained swiveling mechanism

Topoisomerase V is a type I topoisomerase without structural or sequence similarities to other topoisomerases. Although it belongs to the type I subfamily of topoisomerases, it is unrelated to either type IA or IB enzymes. We used real-time single-molecule micromechanical experiments to show that topoisomerase V relaxes DNA via events that release multiple DNA turns, employing a constrained swiveling mechanism similar to that for type IB enzymes. Relaxation is powered by the torque in the supercoiled DNA and is constrained by friction between the protein and the DNA. Although all type IB enzymes share a common structure and mechanism and type IA and type II enzymes show marked structural and functional similarities, topoisomerase V represents a different type of topoisomerase that relaxes DNA in a similar overall manner as type IB molecules but by using a completely different structural and mechanistic framework.

Hin-mediated DNA knotting and recombining promote replicon dysfunction and mutation

Background

The genetic code imposes a dilemma for cells. The DNA must be long enough to encode for the complexity of an organism, yet thin and flexible enough to fit within the cell. The combination of these properties greatly favors DNA collisions, which can knot and drive recombination of the DNA. Despite the well-accepted propensity of cellular DNA to collide and react with itself, it has not been established what the physiological consequences are.

Results

Here we analyze the effects of recombined and knotted plasmids in E. coli using the Hin site-specific recombination system. We show that Hin-mediated DNA knotting and recombination (i) promote replicon loss by blocking DNA replication; (ii) block gene transcription; and (iii) cause genetic rearrangements at a rate three to four orders of magnitude higher than the rate for an unknotted, unrecombined plasmid.

Conclusion

These results show that DNA reactivity leading to recombined and knotted DNA is potentially toxic and may help drive genetic evolution.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30δ DNA

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation. The effects of 20 w/v % EG on the circular dichroism (CD) spectrum of the DNA and on the emission spectrum of the free dye were also measured. At 40 degrees C, addition of 20 w/v % EG caused a substantial (1.27- to 1.35-fold) increase in alpha, a significant change in the CD spectrum, and a very small, marginally significant increase in R(H), but little or no change in the amplitude of dye wobble in its binding site or the lifetime of intercalated dye. Together with previously reported measurements of E(T), these results imply that the bending elastic constant of DNA is significantly decreased by 20 w/v % EG at 40 degrees C. At 20 degrees C, addition of 20 w/v % EG caused a marginally significant decrease in alpha and very little change in any other measured properties. Also at 20 degrees C, addition of 30 w/v % betaine caused a marginally significant increase in alpha and significant but modest change in the CD spectrum, but very little change in any other properties.

Torsional rigidities of weakly strained DNAs

Measurements on unstrained linear and weakly strained large (> or =340 bp) circular DNAs yield torsional rigidities in the range C = 170-230 fJ fm. However, larger values, in the range C = 270-420 fJ fm, are typically obtained from measurements on sufficiently small (< or =247 bp) circular DNAs, and values in the range C = 300-450 fJ fm are obtained from experiments on linear DNAs under tension. A new method is proposed to estimate torsional rigidities of weakly supercoiled circular DNAs. Monte Carlo simulations of the supercoiling free energies of solution DNAs, and also of the structures of surface-confined supercoiled plasmids, were performed using different trial values of C. The results are compared with experimental measurements of the twist energy parameter, E(T), that governs the supercoiling free energy, and also with atomic force microscopy images of surface-confined plasmids. The results clearly demonstrate that C-values in the range 170-230 fJ fm are compatible with experimental observations, whereas values in the range C > or = 269 fJ fm, are incompatible with those same measurements. These results strongly suggest that the secondary structure of DNA is altered by either sufficient coherent bending strain or sufficient tension so as to enhance its torsional rigidity.