Competition between curls and plectonemes near the buckling transition of stretched supercoiled DNA

DNA packaging by molecular motors: from bacteriophage to human chromosomes

Dense packaging of genomic DNA is crucial for organismal survival, as DNA length always far exceeds the dimensions of the cells that contain it. Organisms, therefore, use sophisticated machineries to package their genomes. These systems range across kingdoms from a single ultra-powerful rotary motor that spools the DNA into a bacteriophage head, to hundreds of thousands of relatively weak molecular motors that coordinate the compaction of mitotic chromosomes in eukaryotic cells. Recent technological advances, such as DNA proximity-based sequencing approaches, polymer modelling and in vitro reconstitution of DNA loop extrusion, have shed light on the biological mechanisms driving DNA organization in different systems. Here, we discuss DNA packaging in bacteriophage, bacteria and eukaryotic cells, which, despite their extreme variation in size, structure and genomic content, all rely on the action of molecular motors to package their genomes.

The SKMT Algorithm: A method for assessing and comparing underlying protein entanglement

We present fast and simple-to-implement measures of the entanglement of protein tertiary structures which are appropriate for highly flexible structure comparison. These are performed using the SKMT algorithm, a novel method of smoothing the Cα backbone to achieve a minimal complexity curve representation of the manner in which the protein's secondary structure elements fold to form its tertiary structure. Its subsequent complexity is characterised using measures based on the writhe and crossing number quantities heavily utilised in DNA topology studies, and which have shown promising results when applied to proteins recently. The SKMT smoothing is used to derive empirical bounds on a protein's entanglement relative to its number of secondary structure elements. We show that large scale helical geometries dominantly account for the maximum growth in entanglement of protein monomers, and further that this large scale helical geometry is present in a large array of proteins, consistent across a number of different protein structure types and sequences. We also show how these bounds can be used to constrain the search space of protein structure prediction from small angle x-ray scattering experiments, a method highly suited to determining the likely structure of proteins in solution where crystal structure or machine learning based predictions often fail to match experimental data. Finally we develop a structural comparison metric based on the SKMT smoothing which is used in one specific case to demonstrate significant structural similarity between Rossmann fold and TIM Barrel proteins, a link which is potentially significant as attempts to engineer the latter have in the past produced the former. We provide the SWRITHE interactive python notebook to calculate these metrics.

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.

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.

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.

Condensin-driven loop extrusion on supercoiled DNA

Condensin, a structural maintenance of chromosomes (SMC) complex, has been shown to be a molecular motor protein that organizes chromosomes by extruding loops of DNA. In cells, such loop extrusion is challenged by many potential conflicts, for example, the torsional stresses that are generated by other DNA-processing enzymes. It has so far remained unclear how DNA supercoiling affects loop extrusion. Here, we use time-lapse single-molecule imaging to study condensin-driven DNA loop extrusion on supercoiled DNA. We find that condensin binding and DNA looping are stimulated by positively supercoiled DNA, and condensin preferentially binds near the tips of supercoiled plectonemes. Upon loop extrusion, condensin collects nearby plectonemes into a single supercoiled loop that is highly stable. Atomic force microscopy imaging shows that condensin generates supercoils in the presence of ATP. Our findings provide insight into the topology-regulated loading and formation of supercoiled loops by SMC complexes and clarify the interplay of loop extrusion and supercoiling.

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.

DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases

Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or ‘push’ between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate by modifying the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the DNA torsional response. Using stochastic simulations, we demonstrate DNA supercoiling-mediated long-range cooperation between co-transcribing RNAPs. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern transcriptional dynamics.

Writhing and hockling instabilities in twisted elastic fibers

The buckling and twisting of slender, elastic fibers is a deep and well-studied field. A slender elastic rod that is twisted with respect to a fixed end will spontaneously form a loop, or hockle, to relieve the torsional stress that builds. Further twisting results in the formation of plectonemes-a helical excursion in the fiber that extends with additional twisting. Here we use an idealized, micron-scale experiment to investigate the energy stored, and subsequently released, by hockles and plectonemes as they are pulled apart, in analogy with force spectroscopy studies of DNA and protein folding. Hysteresis loops in the snapping and unsnapping inform the stored energy in the twisted fiber structures.

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.

WASP: a software package for correctly characterizing the topological development of ribbon structures

We introduce the Writhe Application Software Package (WASP) which can be used to characterisze the topology of ribbon structures, the underlying mathematical model of DNA, Biopolymers, superfluid vorticies, elastic ropes and magnetic flux ropes. This characterization is achieved by the general twist–writhe decomposition of both open and closed ribbons, in particular through a quantity termed the polar writhe. We demonstrate how this decomposition is far more natural and straightforward than artificial closure methods commonly utilized in DNA modelling. In particular, we demonstrate how the decomposition of the polar writhe into local and non-local components distinctly characterizes the local helical structure and knotting/linking of the ribbon. This decomposition provides additional information not given by alternative approaches. As example applications, the WASP routines are used to characterise the evolving topology (writhe) of DNA minicircle and open ended plectoneme formation magnetic/optical tweezer simulations, and it is shown that the decomponsition into local and non-local components is particularly important for the detection of plectonemes. Finally it is demonstrated that a number of well known alternative writhe expressions are actually simplifications of the polar writhe measure.

Dynamics of the Buckling Transition in Double-Stranded DNA and RNA

DNA under torsional strain undergoes a buckling transition that is the fundamental step in plectoneme nucleation and supercoil dynamics, which are critical for the processing of genomic information. Despite its importance, quantitative models of the buckling transition, in particular to also explain the surprising two-orders-of-magnitude difference between the buckling times for RNA and DNA revealed by single-molecule tweezers experiments, are currently lacking. Additionally, little is known about the configurations of the DNA during the buckling transition because they are not directly observable experimentally. Here, we use a discrete worm-like chain model and Brownian dynamics to simulate the DNA/RNA buckling transition. Our simulations are in good agreement with experimentally determined parameters of the buckling transition. The simulations show that the buckling time strongly and exponentially depends on the bending stiffness, which accounts for more than half the measured difference between DNA and RNA. Analyzing the microscopic conformations of the chain revealed by our simulations, we find clear evidence for a solenoid-shaped transition state and a curl intermediate. The curl intermediate features a single loop and becomes increasingly populated at low forces. Taken together, the simulations suggest that the worm-like chain model can account semiquantitatively for the buckling dynamics of both DNA and RNA.

Topology, Geometry, and Mechanics of Strongly Stretched and Twisted Filaments: Solenoids, Plectonemes, and Artificial Muscle Fibers

Soft elastic filaments that can be stretched, bent, and twisted exhibit a range of topologically and geometrically complex morphologies. Recently, a number of experiments have shown how to use these building blocks to create filament-based artificial muscles that use the conversion of writhe to extension or contraction, exposing the connection between topology, geometry, and mechanics. Here, we combine numerical simulations of soft elastic filaments that account for geometric nonlinearities and self-contact to map out the basic structures underlying artificial muscle fibers in a phase diagram that is a function of the extension and twist density. We then use ideas from computational topology to track the interconversion of link, twist, and writhe in these geometrically complex physical structures to explain the physical principles underlying artificial muscle fibers and provide guidelines for their design.

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.

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

Kinetic Pathway of Torsional DNA Buckling

In magnetic tweezers experiments, we observe that torsional DNA buckling rates and transition state distances are insensitive to base-pairing defects. This is surprising because defects are expected to kink DNA and lower the energy of a localized loop. Nonetheless, base-pairing defects lead to pinning of buckled structures at the defects, which may be important for DNA repair in vivo. We find that the decrease in entropy from pinning roughly balances the decrease in bending energy, explaining why defects have little effect on buckling rates. Our data are generally consistent with elastic rod theory, which predicts that the transition state structure for torsional buckling is a localized wave with a specific shape ("soliton"). The transition state soliton decays to a metastable looped intermediate ("curl") that is separated from the final, fully buckled state by a second, low energy barrier. DNAs with base mismatch defects buckle at lower torque, where elastic rod theory predicts the loop structure is more stable, and manifest an intermediate buckling structure consistent with such a loop. We estimate that, under our high force, high salt experimental conditions, the soliton barrier is approximately 10 k_B T and, to reach this transition state from the unbuckled state, the system torque instantaneously decreases by approximately 1 pN·nm for DNA with or without a small defect.

Defect-facilitated buckling in supercoiled double-helix DNA

We present a statistical-mechanical model for stretched twisted double-helix DNA, where thermal fluctuations are treated explicitly from a Hamiltonian without using any scaling hypotheses. Our model applied to defect-free supercoiled DNA describes the coexistence of multiple plectoneme domains in long DNA molecules at physiological salt concentrations (≈0.1M Na^{+}) and stretching forces (≈1pN). We find a higher (lower) number of domains at lower (higher) ionic strengths and stretching forces, in accord with experimental observations. We use our model to study the effect of an immobile point defect on the DNA contour that allows a localized kink. The degree of the kink is controlled by the defect size, such that a larger defect further reduces the bending energy of the defect-facilitated kinked end loop. We find that a defect can spatially pin a plectoneme domain via nucleation of a kinked end loop, in accord with experiments and simulations. Our model explains previously reported magnetic tweezer experiments [A. Dittmore et al., Phys. Rev. Lett. 119, 147801 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.147801] showing two buckling signatures: buckling and "rebuckling" in supercoiled DNA with a base-unpaired region. Comparing with experiments, we find that under 1 pN force, a kinked end loop nucleated at a base-mismatched site reduces the bending energy by ≈0.7 k_{B}T per unpaired base. Our model predicts the coexistence of three states at the buckling and rebuckling transitions, which warrants new experiments.

Nucleation of Multiple Buckled Structures in Intertwined DNA Double Helices

We study the statistical-mechanical properties of intertwined double-helical DNAs (DNA braids). In magnetic tweezers experiments, we find that torsionally stressed stretched braids supercoil via an abrupt buckling transition, which is associated with the nucleation of a braid end loop, and that the buckled braid is characterized by a proliferation of multiple domains. Differences between the mechanics of DNA braids and supercoiled single DNAs can be understood as an effect of the increased bulkiness in the structure of the former. The experimental results are in accord with the predictions of a statistical-mechanical model.

Activities of gyrase and topoisomerase IV on positively supercoiled DNA

Although bacterial gyrase and topoisomerase IV have critical interactions with positively supercoiled DNA, little is known about the actions of these enzymes on overwound substrates. Therefore, the abilities of Bacillus anthracis and Escherichia coli gyrase and topoisomerase IV to relax and cleave positively supercoiled DNA were analyzed. Gyrase removed positive supercoils ∼10-fold more rapidly and more processively than it introduced negative supercoils into relaxed DNA. In time-resolved single-molecule measurements, gyrase relaxed overwound DNA with burst rates of ∼100 supercoils per second (average burst size was 6.2 supercoils). Efficient positive supercoil removal required the GyrA-box, which is necessary for DNA wrapping. Topoisomerase IV also was able to distinguish DNA geometry during strand passage and relaxed positively supercoiled substrates ∼3-fold faster than negatively supercoiled molecules. Gyrase maintained lower levels of cleavage complexes with positively supercoiled (compared with negatively supercoiled) DNA, whereas topoisomerase IV generated similar levels with both substrates. Results indicate that gyrase is better suited than topoisomerase IV to safely remove positive supercoils that accumulate ahead of replication forks. They also suggest that the wrapping mechanism of gyrase may have evolved to promote rapid removal of positive supercoils, rather than induction of negative supercoils.

Supercoiling DNA Locates Mismatches

We present a method of detecting sequence defects by supercoiling DNA with magnetic tweezers. The method is sensitive to a single mismatched base pair in a DNA sequence of several thousand base pairs. We systematically compare DNA molecules with 0 to 16 adjacent mismatches at 1 M monovalent salt and 3.6 pN force and show that under these conditions, a single plectoneme forms and is stably pinned at the defect. We use these measurements to estimate the energy and degree of end-loop kinking at defects. From this, we calculate the relative probability of plectoneme pinning at the mismatch under physiologically relevant conditions. Based on this estimate, we propose that DNA supercoiling could contribute to mismatch and damage sensing in vivo.

Platinum-Based Drugs and DNA Interactions Studied by Single-Molecule and Bulk Measurements

Platinum-containing molecules are widely used as anticancer drugs. These molecules exert cytotoxic effects by binding to DNA through various mechanisms. The binding between DNA and platinum-based drugs hinders the opening of DNA, and therefore, DNA duplication and transcription are severely hampered. Overall, impeding the above-mentioned important DNA mechanisms results in irreversible DNA damage and the induction of apoptosis. Several molecules, including multinuclear platinum compounds, belong to the family of platinum drugs, and there is a body of research devoted to developing more efficient and less toxic versions of these compounds. In this study, we combined different biophysical methods, including single-molecule assays (magnetic tweezers) and bulk experiments (ultraviolet absorption for thermal denaturation) to analyze the differential stability of double-stranded DNA in complex with either cisplatin or multinuclear platinum agents. Specifically, we analyzed how the binding of BBR3005 and BBR3464, two representative multinuclear platinum-based compounds, to DNA affects its stability as compared with cisplatin binding. Our results suggest that single-molecule approaches can provide insights into the drug-DNA interactions that underlie drug potency and provide information that is complementary to that generated from bulk analysis; thus, single-molecule approaches have the potential to facilitate the selection and design of optimized drug compounds. In particular, relevant differences in DNA stability at the single-molecule level are demonstrated by analyzing nanomechanically induced DNA denaturation. On the basis of the comparison between the single-molecule and bulk analyses, we suggest that transplatinated drugs are able to locally destabilize small portions of the DNA chain, whereas other regions are stabilized.

Torque and buckling in stretched intertwined double-helix DNAs

We present a statistical-mechanical model for the behavior of intertwined DNAs, with a focus on their torque and extension as a function of their catenation (linking) number and applied force, as studied in magnetic tweezers experiments. Our model produces results in good agreement with available experimental data and predicts a catenation-dependent effective twist modulus distinct from what is observed for twisted individual double-helix DNAs. We find that buckling occurs near the point where experiments have observed a kink in the extension versus linking number, and that the subsequent "supercoiled braid" state corresponds to a proliferation of multiple small plectoneme structures. We predict a discontinuity in extension at the buckling transition corresponding to nucleation of the first plectoneme domain. We also find that buckling occurs for lower linking number at lower salt; the opposite trend is observed for supercoiled single DNAs.

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.

Probing the mechanical properties, conformational changes, and interactions of nucleic acids with magnetic tweezers

Nucleic acids are central to the storage and transmission of genetic information. Mechanical properties, along with their sequence, both enable and fundamentally constrain the biological functions of DNA and RNA. For small deformations from the equilibrium conformations, nucleic acids are well described by an isotropic elastic rod model. However, external forces and torsional strains can induce conformational changes, giving rise to a complex force-torque phase diagram. This review focuses on magnetic tweezers as a powerful tool to precisely determine both the elastic parameters and conformational transitions of nucleic acids under external forces and torques at the single-molecule level. We review several variations of magnetic tweezers, in particular conventional magnetic tweezers, freely orbiting magnetic tweezers and magnetic torque tweezers, and discuss their characteristic capabilities. We then describe the elastic rod model for DNA and RNA and discuss conformational changes induced by mechanical stress. The focus lies on the responses to torque and twist, which are crucial in the mechanics and interactions of nucleic acids and can directly be measured using magnetic tweezers. We conclude by highlighting several recent studies of nucleic acid-protein and nucleic acid-small-molecule interactions as further applications of magnetic tweezers and give an outlook of some exciting developments to come.

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).

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.

Collapse and coexistence for a molecular braid with an attractive interaction component subject to mechanical forces

Dual mechanical braiding experiments provide a useful tool with which to investigate the nature of interactions between rod-like molecules, for instance actin and DNA. In conditions close to molecular condensation, one would expect an appearance of a local minimum in the interaction potential between the two molecules. We investigate this situation, introducing an attractive component into the interaction potential, using a model developed for describing such experiments. We consider both attractive interactions that do not depend on molecular structure and those which depend on a DNA-like helix structure. In braiding experiments, an attractive term may lead to certain effects. A local minimum may cause molecules to collapse from a loosely braided configuration into a tight one, occurring at a critical value of the moment applied about the axis of the braid. For a fixed number of braid pitches, this may lead to coexistence between the two braiding states, tight and loose. Coexistence implies certain proportions of the braid are in each state, their relative size depending on the number of braid pitches. This manifests itself as a linear dependence in numerically calculated quantities as functions of the number of braid pitches. Also, in the collapsed state, the braid radius stays roughly constant. Furthermore, if the attractive interaction is helix dependent, the left-right handed braid symmetry is broken. For a DNA like charge distribution, using the Kornyshev-Leikin interaction model, our results suggest that significant braid collapse and coexistence only occurs for left handed braids. Regardless of the interaction model, the study highlights the possible qualitative physics of braid collapse and coexistence; and the role helix specific forces might play, if important. The model could be used to connect other microscopic theories of interaction with braiding experiments.

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.

Plectoneme tip bubbles: coupled denaturation and writhing in supercoiled DNA

We predict a novel conformational regime for DNA, where denaturation bubbles form at the tips of plectonemes, and study its properties using coarse-grained simulations. For negative supercoiling, this regime lies between bubble-dominated and plectoneme-dominated phases, and explains the broad transition between the two observed in experiment. Tip bubbles cause localisation of plectonemes within thermodynamically weaker AT-rich sequences, and can greatly suppress plectoneme diffusion by a pinning mechanism. They occur for supercoiling densities and forces that are typically encountered for DNA in vivo, and may be exploited for biological control of genomic processes.

Coexistence of twisted, plectonemic, and melted DNA in small topological domains

DNA responds to small changes in force and torque by over- or undertwisting, forming plectonemes, and/or melting bubbles. Although transitions between either twisted and plectonemic conformations or twisted and melted conformations have been described as first-order phase transitions, we report here a broadening of these transitions when the size of a topological domain spans several kilobasepairs. Magnetic tweezers measurements indicate the coexistence of three conformations at subpicoNewton force and linking number densities ∼-0.06. We present a statistical physics model for DNA domains of several kilobasepairs by calculating the full partition function that describes this three-state coexistence. Real-time analysis of short DNA tethers at constant force and torque shows discrete levels of extension, representing discontinuous changes in the size of the melting bubble, which should reflect the underlying DNA sequence. Our results provide a comprehensive picture of the structure of underwound DNA at low force and torque and could have important consequences for various biological processes, in particular those that depend on local DNA melting, such as the initiation of replication and transcription.

Self-consistent treatment of electrostatics in molecular DNA braiding through external forces

In this paper we consider a physical system in which two DNA molecules braid about each other. The distance between the two molecular ends, on either side of the braid, is held at a distance much larger than supercoiling radius of the braid. The system is subjected to an external pulling force, and a moment that induces the braiding. In a model, developed for understanding such a system, we assume that each molecule can be divided into a braided and unbraided section. We also suppose that the DNA is nicked so that there is no constraint of the individual linking numbers of the molecules. Included in the model are steric and electrostatic interactions, thermal fluctuations of the braided and unbraided sections of the molecule, as well as the constraint on the braid linking (catenation) number. We compare two approximations used in estimating the free energy of the braided section. One is where the amplitude of undulations of one molecule with respect to the other is determined only by steric interactions. The other is a self-consistent determination of the mean-squared amplitude of these undulations. In this second approximation electrostatics should play an important role in determining this quantity, as suggested by physical arguments. We see that if the electrostatic interaction is sufficiently large there are indeed notable differences between the two approximations. We go on to test the self-consistent approximation-included in the full model-against experimental data for such a system, and we find good agreement. However, there seems to be a slight left-right-handed braid asymmetry in some of the experimental results. We discuss what might be the origin of this small asymmetry.

Global force-torque phase diagram for the DNA double helix: structural transitions, triple points, and collapsed plectonemes

We present a free energy model for structural transitions of the DNA double helix driven by tensile and torsional stress. Our model is coarse grained and is based on semiflexible polymer descriptions of B-DNA, underwound L-DNA, and highly overwound P-DNA. The statistical-mechanical model of plectonemic supercoiling previously developed for B-DNA is applied to semiflexible polymer models of P- and L-DNA to obtain a model of DNA structural transitions in quantitative accord with experiment. We identify two distinct plectonemic states, one "inflated" by electrostatic repulsion and thermal fluctuations and the other "collapsed," with the two double helices inside the supercoils driven to close contact. We find that supercoiled B and L are stable only in the inflated form, while supercoiled P is always collapsed. We also predict the behavior and experimental signatures of highly underwound "Q"-DNA, the left-handed analog of P-DNA; as for P, supercoiled Q is always collapsed. Overstretched "S"-DNA and strand-separated "stress-melted" DNA are also included in our model, allowing prediction of a global phase diagram for forces up to 1000 pN and torques between ±60 pN nm, or, in terms of linking number density, from σ=-5 to +3.

Multiplectoneme phase of double-stranded DNA under torsion

We use the wormlike chain model to study supercoiling of DNA under tension and torque. The model reproduces experimental data for a broad range of forces, salt concentrations, and contour lengths. We find a plane of first-order phase transitions ending in a smeared-out line of critical points, the multiplectoneme phase, which is characterized by a fast twist-mediated diffusion of plectonemes and a torque that rises after plectoneme formation with increasing linking number. The discovery of this phase at the same time resolves the discrepancies between existing models and experiment.

Torque spectroscopy of DNA: base-pair stability, boundary effects, backbending, and breathing dynamics

Changes in global DNA linking number can be accommodated by localized changes in helical structure. We have used single-molecule torque measurements to investigate sequence-specific strand separation and Z-DNA formation. By controlling the boundary conditions at the edges of sequences of interest, we have confirmed theoretical predictions of distinctive boundary-dependent backbending patterns in torque-twist relationships. Abrupt torque jumps are associated with the formation and collapse of DNA bubbles, permitting direct observations of DNA breathing dynamics.

The SMC1-SMC3 cohesin heterodimer structures DNA through supercoiling-dependent loop formation

Cohesin plays a critical role in sister chromatid cohesion, double-stranded DNA break repair and regulation of gene expression. However, the mechanism of how cohesin directly interacts with DNA remains unclear. We report single-molecule experiments analyzing the interaction of the budding yeast cohesin Structural Maintenance of Chromosome (SMC)1-SMC3 heterodimer with naked double-helix DNA. The cohesin heterodimer is able to compact DNA molecules against applied forces of 0.45 pN, via a series of extension steps of a well-defined size ≈130 nm. This reaction does not require ATP, but is dependent on DNA supercoiling: DNA with positive torsional stress is compacted more quickly than negatively supercoiled or nicked DNAs. Un-nicked torsionally relaxed DNA is a poor substrate for the compaction reaction. Experiments with mutant proteins indicate that the dimerization hinge region is crucial to the folding reaction. We conclude that the SMC1-SMC3 heterodimer is able to restructure the DNA double helix into a series of loops, with a preference for positive writhe.

Magnetic tweezers measurements of the nanomechanical stability of DNA against denaturation at various conditions of pH and ionic strength

The opening of DNA double strands is extremely relevant to several biological functions, such as replication and transcription or binding of specific proteins. Such opening phenomenon is particularly sensitive to the aqueous solvent conditions in which the DNA molecule is dispersed, as it can be observed by considering the classical dependence of DNA melting temperature on pH and salt concentration. In the present work, we report a single-molecule study of the stability of DNA against denaturation when subjected to changes in solvent. We investigated the appearance of DNA instability under specific external applied force and imposed twist values, which was revealed by an increase in the temporal fluctuations in the DNA extension. These fluctuations occur in the presence of a continuous interval of equilibrium states, ranging from a plectonemic state to a state characterized by denaturation bubbles. In particular, we observe the fluctuations only around a characteristic force value. Moreover, this characteristic force is demonstrated to be notably sensitive to variations in the pH and ionic strength. Finally, an extension of a theoretical model of plectoneme formation is used to estimate the average denaturation energy, which is found to be linearly correlated to the melting temperature of the DNA double strands.

Denaturation of circular DNA: supercoils and overtwist

The denaturation transition of circular DNA is studied within a Poland-Scheraga-type approach, generalized to account for the fact that the total linking number (LK), which measures the number of windings of one strand around the other, is conserved. In the model the LK conservation is maintained by invoking both overtwisting and writhing (supercoiling) mechanisms. This generalizes previous studies, which considered each mechanism separately. The phase diagram of the model is analyzed as a function of the temperature and the elastic constant κ associated with the overtwisting energy for any given loop entropy exponent c. As in the case where the two mechanisms apply separately, the model exhibits no denaturation transition for c ≤ 2. For c > 2 and κ = 0 we find that the model exhibits a first-order transition. The transition becomes of higher order for any κ > 0. We also calculate the contribution of the two mechanisms separately in maintaining the conservation of the linking number and find that it is weakly dependent on the loop exponent c.

Dynamics of DNA supercoils

DNA in cells exhibits a supercoiled state in which the double helix is additionally twisted to form extended intertwined loops called plectonemes. Although supercoiling is vital to many cellular processes, its dynamics remain elusive. In this work, we directly visualize the dynamics of individual plectonemes. We observe that multiple plectonemes can be present and that their number depends on applied stretching force and ionic strength. Plectonemes moved along DNA by diffusion or, unexpectedly, by a fast hopping process that facilitated very rapid (<20 milliseconds) long-range plectoneme displacement by nucleating a new plectoneme at a distant position. These observations directly reveal the dynamics of plectonemes and identify a mode of movement that allows long-distance reorganization of the conformation of the genome on a millisecond time scale.

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.

Electromagnetic torque tweezers: a versatile approach for measurement of single-molecule twist and torque

The well-established single-molecule force-spectroscopy techniques have recently been complemented by methods that can measure torque and twist directly, notably magnetic torque tweezers and the optical torque wrench. A limitation of the current torque measurement schemes is the intrinsic coupling between the force and torque degrees of freedom. Here we present electromagnetic torque tweezers (eMTT) that combine permanent and electromagnets to enable independent control of the force and torsional trap stiffness for sensitive measurements of single molecule torque and twist. Using the eMTT, we demonstrate sensitive torque measurements on tethered DNA molecules from simple tracking of the beads' (x,y)-position, obviating the need for any angular tracking algorithms or markers. Employing the eMTT for high-resolution torque measurements, we experimentally confirm the theoretically predicted torque overshoot at the DNA buckling transition in high salt conditions. We envision that the flexibility and control afforded by the eMTT will enable a range of new torque and twist measurement schemes from single-molecules to living cells.

Recent developments in single-molecule DNA mechanics

Over the past two decades, measurements on individual stretched and twisted DNA molecules have helped define the basic elastic properties of the double helix and enabled real-time functional assays of DNA-associated molecular machines. Recently, new magnetic tweezers approaches for simultaneously measuring freely fluctuating twist and extension have begun to shed light on the structural dynamics of large nucleoprotein complexes. Related technical advances have facilitated direct measurements of DNA torque, contributing to a better understanding of abrupt structural transitions in mechanically stressed DNA. The new measurements have also been exploited in studies that hint at a developing synergistic relationship between single-molecule manipulation and structural DNA nanotechnology.