Supercoiling-dependent DNA binding: quantitative modeling and applications to bulk and single-molecule experiments

Abstract DNA stores our genetic information and is ubiquitous in applications, where it interacts with binding partners ranging from small molecules to large macromolecular complexes. Binding is modulated by mechanical strains in the molecule and can change local DNA structure. Frequently, DNA occurs in closed topological forms where topology and supercoiling add a global constraint to the interplay of binding-induced deformations and strain-modulated binding. Here, we present a quantitative model with a straight-forward numerical implementation of how the global constraints introduced by DNA topology modulate binding. We focus on fluorescent intercalators, which unwind DNA and enable direct quantification via fluorescence detection. Our model correctly describes bulk experiments using plasmids with different starting topologies, different intercalators, and over a broad range of intercalator and DNA concentrations. We demonstrate and quantitatively model supercoiling-dependent binding in a single-molecule assay, where we directly observe the different intercalator densities going from supercoiled to nicked DNA. The single-molecule assay provides direct access to binding kinetics and DNA supercoil dynamics. Our model has broad implications for the detection and quantification of DNA, including the use of psoralen for UV-induced DNA crosslinking to quantify torsional tension in vivo, and for the modulation of DNA binding in cellular contexts.


Introduction
DNA is the carrier of genetic information in all cellular life.In vivo , double-stranded DNA is often present in circular and, therefore, topological closed form.In particular, bacte-rial chromosomes and plasmids are circular DNA molecules, whose degree of supercoiling is tightly regulated in vivo (1)(2)(3) .In eukaryotes, DNA topology and supercoiling similarly play important roles in the context of a chromatinized genome, e.g. in the compaction, processing, and regulation of genetic information (4)(5)(6)(7)(8)(9)(10)(11)(12) .
Both in its biological role and in many biotechnological applications, DNA interacts with a broad range of binding partners that range from small molecules to large proteins complexes.In particular, the detection and quantification of DNA often rely on staining with fluorescent small molecules that frequently bind in an intercalative binding mode (13)(14)(15)(16)(17)(18)(19) .Ligand binding to DNA can, in general, locally alter the DNA structure and introduce strains away from the equilibrium B-form DNA conformation (20)(21)(22)(23)(24)(25) .In turn, stretching forces and torsional strains have been shown to systematically affect binding equilibria ( 20 ,26-33 ) .Having a defined DNA topology, e.g. in a plasmid or other topological domains, imposes a global constraint on the interplay between strain-dependent binding and binding-induced conformational changes, as described in seminal work by Bauer and Vinograd ( 26 ,27 ,33 ) and subsequently shown for both small molecule ( 4 , 26 , 27 , 33 , 34 ) and protein binding to DNA ( 10 , 35 , 36 ) .
Here, we present a quantitative model for DNA-ligand binding for various degrees of supercoiling under a global topological constraint and apply it to both bulk experiments and to a single-molecule assay that enables us to study DNA supercoil dynamics in real time.Our model for ligand binding to topologically closed DNA combines and integrates aspects of previous work.In particular we take into account the constraints of binding to a linear molecule ( 37 ) , both the finite DNA and intercalator concentrations ( 19 ) , the conformational changes induced by the ligand ( 13 , 14 , 19 , 27 , 31 ) , the influence of strains on binding equilibria ( 26 ,33 ) , a physical model of plasmid mechanics ( 38 ,39 ) , and the global constraint introduced by having a defined linking number due to the defined topology ( 27 , 40 , 41 ) .Importantly, the parameters of our model have a direct physical interpretation and can be determined from other assays.Therefore, our model can provide quantitative predictions without free parameters.
We first perform experiments in bulk using native topologically constrained plasmid DNA, i.e. circular DNA with both strands fully intact ( referred to as the supercoiled species, 'sc'; Figure 1 A ) and, for comparison, topologically open DNA ( nicked, open circular, 'oc', or linearized, 'lin'; Figure 1 B, C ) .We use fluorescence detection to quantify the amount of binding and find quantitative agreement with our model.We then apply our model to a single-molecule DNA assay ( 50 ,53 ) , where DNA is supercoiled in situ by intercalation.We monitor the fluorescence change upon nicking of the DNA molecule, which induces an abrupt transition from a closed topology to a topologically unconstrained, nicked form and again find excellent agreement with our model.The single-molecule assay enables us to observe the re-adjustment of the binding equilibrium upon change in topology in real time, which enables us to probe the dynamics of torque release, intercalation, and supercoil relaxation.
Our findings have direct practical applications since gelbased assays for the discrimination and detection of topoisomers are widely used to study the properties of circular DNA and of various enzymes that alter DNA topology, including topoisomerases, gyrase, reverse gyrase, and recombinases (59)(60)(61) .An unbiased quantification of the different topoisomers using fluorescence staining, which is increasingly used to replace radiolabelling due to the hazards associated with handling, storing, and disposing of radioactive materials, must take into account the observed topology dependencies.Furthermore, our model is directly applicable to single-molecule fluorescence studies of supercoiled DNA, and we demonstrate real-time monitoring of supercoiling-dependent binding dynamics.Finally, we provide practical guidelines for unbiased detection of different topoisomers and discuss consequences of topology-dependent binding more broadly.

Plasmid DNA preparation
We used the commonly used plasmid pBR322 ( NEB, catalogue number N3033S ) as the DNA substrate for bulk measurements of intercalation.We prepared different topological states of the DNA by cutting ( i.e. introducing a doublestrand break ) or nicking ( i.e. a DNA single-strand break ) the original supercoiled DNA to obtain linear and open circular DNA, respectively.Reactions were performed in NEBuffer 3.1 ( NEB ) using the enzymes EcoRV ( NEB; incubation temperature 37   C for 20 min.The products were purified with a PCR clean-up kit ( Qiagen ) .For the topologically closed and initially negatively supercoiled condition, the plasmid was only cleaned with a PCR clean-up kit ( Qiagen ) and we assume an initial supercoiling density of σ ∼ = −5%, as previously reported (62)(63)(64)(65)(66) .Additionally, we prepared topologically closed but initially relaxed plasmids by incubation with Wheat Germ Topoisomerase I.For a total volume of 100 μl, 16.2 μl assay buffer ( 50 mM Tris HCl, 1 mM EDTA, 1 mM DTT, 20% ( v / v ) glycerol, 50 mM NaCl, pH 7.9 ) , 82.12 μl RNase-free water, 2.16 μl topoisomerase I, 0.52 μl pBR322 ( c = 1000 ng / μl ) were combined in a reaction tube and incubated for 1 h at 37 • C.
Control measurements used natively negatively supercoiled plasmids of different lengths and different sequences: pUC19 ( 2.686 kbp ) and M13mp18 ( 7.249 kbp ) , both obtained similar from NEB ( catalogue numbers N3041S and N4040S, respectively ) and purified with a PCR clean-up kit ( Qiagen ) .DNA concentrations were determined using a nanodrop UV / vis photospectrometer ( ThermoFisher Scientific ) .To every 100 μl reaction volume, we added 20 μl Gel Loading Dye Purple ( 6 ×) ( NEB, catalogue number B7024S ) prior to running the gel.For the plate reader and qPCR cycler experiments, no gel loading dye was added.

DNA dilution series
For the DNA dilution series with negatively supercoiled DNA, the different topologies were combined to an equimolar mixture and diluted with Tris-acetate-EDT A ( T AE ) buffer ( 40 mM Tris, 20 mM acetic acid and 1 mM EDTA, pH 8.6 ) and gel loading dye in a serial dilution to obtain 7 different DNA concentrations ( Table S1 ) .For the experiments using topoisomerase relaxed DNA, only the linear and the relaxed topologies were combined to an equimolar mixture and diluted with TAE buffer and gel loading dye, again in a serial dilution to obtain seven different DNA concentrations ( Table S1 ) .

AFM imaging of DNA plasmids
AFM imaging of plasmid DNA was performed as described previously (67)(68)(69) .In brief, for the AFM imaging, we deposited 20 μl of DNA at different topological states in TE buffer at a final concentration of 1 ng / μl on freshly cleaved poly-l -lysine ( Sigma Aldrich, diluted to 0.01% in milliQ wa-ter; PLL ) -coated muscovite mica.The sample was incubated 30 s before washing with 20 ml MilliQ water and drying with a gentle stream of filtered argon gas.After drying, the AFM images were recorded in tapping mode at room temperature using a Nanowizard Ultraspeed 2 ( JPK, Berlin, Germany ) AFM with silicon tips ( FASTSCAN-A, drive frequency 1400 kHz, tip radius 5 nm, Bruker, Billerica, Massachusetts, USA ) .Images were scanned over different fields of view with a scanning speed of 5 Hz.The free amplitude was set to 10 nm.The amplitude setpoint was set to 80% of the free amplitude and adjusted to maintain a good image resolution.AFM image post-processing was performed in the software SPIP ( v.6.4,Image Metrology, Hørsholm, Denmark ) to flatten and line-wise correct the images ( Supplementary Figure S2 ) .

Gel electrophoresis
For gel electrophoresis we used 1%-broad-range-agarose ( Carl Roth ) gels and T AE buffer.W e used the 1 kb gene ruler ( Thermo Scientific; 5 μl ) as a size standard.The gels were run for 120 min at 75 V at 4 • C. Subsequently, the gels were removed from the gel box and placed for 20 minutes in 100 ml of 0.5 μM ( 1:100 000 dilution of the stock; 2 independent gels ) , 5 μM ( 1:10 000 dilution of the stock; 7 independent gels ) , or 50 μM ( 1:1000 dilution of the stock; 2 independent gels ) EtBr in TAE buffer, respectively, for staining.Subsequently, the gel was de-stained in TAE buffer for 15 min.The gels were visualized using a Gel Doc XR+ system ( Biorad ) .The same procedure was employed for SYBR Gold staining, except that the gels were stained for 20 minutes in 100 ml of 3 ( 1:4000 dilution of the stock; 2 gels in total ) or 6 μM ( 1:2000 dilution of the stock; 2 gels in total ) .Since the dye concentration used in staining is reduced by the agarose gel matrix and the de-staining step, we used a staining correction factor of 0.1 as determined previously ( 19 ) , which we use to correct all gel data.

Gel electrophoresis image analysis
We saved the images from the Gel Doc system in scn format to allow for quantitative fluorescence intensity analysis.The software SPIP ( v.6.4,Image Metrology, Hørsholm, Denmark ) was used to remove spikes from the image ( without changing the intensity of the bands ) and to generate average intensity profiles along each lane of the gel.In a next step, we used Origin ( OriginLab, Northampton, MA, USA ) to flatten the background of the profiles and to convert the peaks into fractions of supercoiled, open-circular, linear, and / or relaxed DNA respectively, by calculating the area under the lane profiles ( Supplemenary Figure S3 ) .

Bulk fluorescence experiments
For the SYBR Gold bulk fluorescence measurements, we used a well plate reader ( Tecan Infinite M1000 PRO; well plates: corning black polystyrene 384 well microplate with a flat bottom, Sigma-Aldrich, catalogue number: CLS3821 ) and a qPCR cycler ( CFX96 Touch Real-Time PCR Detection System, BioRad ) .In the well-plate reader, the DNA mix including various SYBR Gold concentrations was filled in the wells and the fluorescence was read out from the bottom of the wells.The excitation and emission bandwidths were set to 5 nm, the gain to 100, the flash frequency to 400 Hz, and the integration time to 20 s.We chose the excitation and emission wave-lengths -according to the excitation and emission maxima for SYBR Gold provided by Invitrogen -to be 495 and 537 nm.
In the fluorescence bulk experiments using a qPCR cycler, the DNA was filled into low-profile PCR tubes ( 19 ) , closed with flat, optical, ultra-clear caps ( Bio Rad, product ID: TCS-0803 ) since the fluorescence was read out from the top of the tubes ( at 24 • C ) .As read-out channels, channels with absorption and emission wavelengths of 494 and 518 nm, respectively, were chosen because these were the closest match to those of SYBR Gold.

Single-molecule fluorescence experiments
Supercoiled DNA were prepared and imaged at the singlemolecule level essentially as described previously ( 50 ,53 ) .Experiments were performed in custom-made flow cells built by connecting a surface-passivated glass slide and a glass coverslip using double-sided tape ( 50 ) .The surface of the glass slides was prepared as previously described ( 70 ) with slight modifications.In brief, after extensive cleaning, the surface was silanized using APTES ( 10% v / v ) and acetic acid ( 5% v / v ) methanol solution.The surface was passivated with NHS-ester PEG ( 5000 Da ) and biotinylated NHS-ester PEG ( 5000 Da ) in relation 40:1.The biotinylated NHS-ester PEG allowed us to tether the DNA molecules to the surface via biotin-streptavidin interactions.
The DNA used for intercalation with SYT O X Orange was prepared as described in ( 71 ) with the exception of introducing multiple-biotin handles at the DNA ends to allow the torsional constrain on the DNA rotation.To introduce multiple biotins on the DNA handles we performed a PCR on pBluescript SK+ ( Stratagene ) with GoTaq G2 DNA polymerase ( Promega, M7845 ) , in the presence of biotin-16-dUTP ( Jena Bioscience, NU-803-BIO16-L ) and dTTP ( Thermo Fisher Scientific, 10520651 ) in molar ratio of 1:4, respectively.The PCR was done using primers: A CCGA GA T A GGGTTGA GTG and CA GGGTCGGAA CA GGA GA GC, resulting in a 1238 bp DNA fragment that contained multiple biotins randomly incorporated due to the presence of biotin-16-dUTP modified nucleotides.The PCR products were cleaned up using a standard purification kit ( Promega, A9282 ) and we digested both the biotin handle and large 42 kb DNA plasmid with SpeI-HF ( New England Biolabs, R3133L ) for 2 h at 37 • C. The reaction was stopped by subsequent heat-inactivation for 20 min at 80 • C.This resulted in linear 42 kb DNA and ∼600-bp biotin handles.The digested products were mixed, using a 10:1 molar excess of the biotin handle to linear plasmid.We then added T4 DNA ligase ( New England Biolabs, M0202L ) in the presence of 1 mM ATP overnight at 16 • C and subsequently heat-inactivated for 20 min at 65 • C. The resulting coilable 42 kb DNA construct was cleaned up by size exclusion chromatography on an ÄKTA pure system, with a homemade gel filtration column containing approximately 46 ml of Sephacryl S-1000 SF gel filtration media ( Cytiva ) , run with TE + 150 mM NaCl buffer.The sample was run at 0.2 ml •min −1 , and we collected 0.5 ml fractions.
For immobilization of the 42 kb biotinylated-DNA, we introduced 100 μl of ∼0.5 pM DNA molecules at a flow rate of 4.2 μl / min in imaging buffer with no oxygen scavenging system components ( 40 mM Tris-HCl, 2 mM Trolox, 2.5 mM MgCl 2 , 65 mM KCl ) .The buffer during the DNA immobilization step contained either a low SYT O X Orange concentration ( 25 nM; to allow introducing positive supercoiling in a later step, see below ) , or high concentration of SYT O X Orange ( 250 nM; to allow introducing negative supercoiling in a later step, see below ) during initial immobilization.
Immediately after the introduction of DNA molecules into the flow cell, we further flowed 100 μl of the same buffer without the DNA at the same flow rate to ensure stretching and tethering of the other end of the DNA to the surface ( which introduces the topological constraint ) as well as to remove unbound DNA molecules.To introduce plectonemic supercoils into the bound DNA molecules that are torsionally constrained, we change the SYT O X Orange concentration by introducing imaging buffer with the same salt concentration but with varying SYT O X Orange concentration and an oxygen scavenging system for imaging ( 40 mM Tris-HCl, 2 mM Trolox, 2.5 mM MgCl 2 , 65 mM KCl, 2.5 mM protocatechuic acid ( PCA ) , 50 nM protocatechuate-3,4-dioxygenase ( PCD ) ) ( 72 ) .For positive supercoiling we increased the concentration from 25 nM during immobilization to 250 nM SYT O X Orange, while for the negative supercoiling we reduced the concentration from the initial binding at 250 nM to a final concentration of 50 nM for imaging.In both cases we checked for a visible presence of supercoiling on the DNA by observing moving foci appearing along the DNA molecules ( Supplementary Movies S1 and S2 ) .Thus, supercoiling was introduced by differential SYT O X Orange concentrations between the DNA-binding step ( where torsional constraint is ensured via multiple biotin molecules ) and the imaging step where SYT O X Orange concentration is increased or reduced in order to generate positive or negative supercoils, respectively.For the direct visualization of plectonemic DNA extruded to the side ( Supplementary Movies S3 and S4 ) we used a second outlet connected to a syringe.After initial imaging of supercoiled DNA without flow, we applied a strong outflow through the side outlet to extrude the plectonemic DNA sideways.We performed the side flow experiments using the same imaging buffer as the measurements without flow, to maintain the same level of supercoiling and environmental conditions.
For fluorescence imaging, we used a home-built objective-TIRF microscope.We employed continuous excitation with a 561 nm ( 15-20 mW ) laser in Highly Inclined and Laminated Optical sheet ( HiLo ) microscopy mode, to image SYT O X Orange-stained DNA as well as to introduce DNA nicking.All images were acquired with an PrimeBSI sCMOS camera at an exposure time of 20-200 ms, depending on the experiment, with a 60 × oil immersion, 1.49 NA CFI APO TIRF ( Nikon ) .For DNA visualization, and kymograph generation, we used a custom written python software published in ( 73 ) .

Numerical implementation of the binding model under topological constraint
To model the effects of the global constraint imposed by DNA topology on binding, we developed a model that takes into account changes in linking number due to intercalation and, conversely, torque-dependent binding ( see the section 'Model for ligand binding under topological constraint' in Results ) .The model comprises coupled Equations 1 -6 that are solved iteratively using a custom routine written in Matlab ( Mathworks ) .Input parameters are the temperature T and for DNA the number of base pairs N bp , the torsional stiffness C, the initial linking difference ΔLk 0 , and the DNA concentration c DNA .N bp , ΔLk 0 , and c DNA are typically known from how the DNA was prepared and we used a value C = 100 nm, unless otherwise noted (74)(75)(76) .For the single-molecule assay, the DNA concentration is poorly defined, but low, much lower than the intercalator concentration used.In this case, we used concentrations in the range of 1-100 pM •bp, which are much lower than the intercalator concentration used, but high enough to ensure numerical stability of the calculation.We found that the calculated results are insensitive to the concentration used in this regime.
Input parameters for the intercalator are the binding site size n , the dissociation constant K d , the length increase per intercalator bound z, the change in DNA helical twist per intercalator bound Δθ , and the intercalator concentration c total .Values for n , K d , Δz , and Δθ used in this work are provided in Table 1 .The outputs of the model are the number of intercalator molecules bound per DNA molecule, the torque in the DNA, and the total fluorescence intensity expected, which contains an overall scaling factor α defined in Equation ( 7 ) .Our numerical implementation uses the numerical solutions to the McGhee-von Hippel described previously ( 19 ,77 ) and a successive over-relaxation type approach ( 78 ) to speed up convergence.The code implementing out model is available in the public repository YODA at https:// doi.org/ 10.24416/ UU01-YDNOBO .

Model for ligand binding under topological constraint
We develop a model for ligand binding to topologically closed supercoiled (e.g.plasmid) DNA, where the topology imposes a global constraint (Figure 1 D-F).Binding to plasmid DNA is different from a standard bimolecular binding equilibrium for several reasons that we take into account in our model.First, the linear structure of DNA imposes local constraints for ligand binding, if bound ligands occupy a binding size of n bases (Figure 1 D).Binding of ligands with binding site size n can be modeled using the McGhee-von Hippel model ( 31 ,37 ) in cases where the DNA concentration is much lower than the ligand concentration, such that the free and total ligand concentrations are approximately equal.For bulk measurements, however, the concentration of DNA bases can be similar to or even larger than the ligand concentration and needs to be considered ( 19 ).Therefore, we use an extension of the McGheevon Hippel model that explicitly takes into account both the ligand ( c t ot al ) and DNA ( c DNA ) concentration that was derived in ( 19 ).The fractional binding γ is given by Here, c bound is the bound ligand concentration, c f ree is the free ligand concentration, c t ot al = c f ree + c bound the total ligand concentration, K d is the dissociation constant (in M), and n is the binding site size (in base pairs).Typical values of the binding site size for intercalators are n ≈ 2, corresponding to binding every other base pair.
Intercalation affects the local geometry of the DNA helix, by locally unwinding and lengthening the helix ( 14 , 16 , 27 , 30 , 31 , 49 , 79 , 80 ) (Figure 1 E).Here we assume that each intercalation event locally lengthens the DNA by Δz and unwinds it by θ .Typical values for intercalators are in the range z ≈ 0.34 nm and θ ≈ 15 • -30 • .The fact that intercalation lengthens and unwinds the DNA helix suggests, by Le Chatelier's principle, that applying a stretching force or unwinding torque, respectively, will increase intercalative binding.Conversely, overwinding the helix will hinder intercalation.We assume an Arrhenius-like exponential dependence of the binding constant ( 30 , 45 , 49 , 81 ) on applied force F and torque Γ : Here, K d,0 is the dissociation constant for the relaxed molecule, i.e. in the absence of forces or torques, k B Boltzmann's constant and T the absolute temperature.For plasmids in free solution, the force is zero (or at least small, specifically F k B T / Δz ≈ 10 pN) and the second exponential factor in Equation ( 2) can be neglected.Values for n , K d , z and θ for selected dyes are summarized in Table 1 .
For linear or open circular DNA molecules, there is no torsional constraint, and the torque will be zero in equilibrium; consequently, binding will simply be determined by Equation ( 1 ).In contrast, for topologically closed plasmids, the topology imposes a global constraint.For a closed plasmid, the linking number Lk is a topological invariant and partitions into twist Tw and writhe Wr by White's formula (also known as Calugareanu's theorem or Calugareanu-White-Fuller theorem) ( 27 , 40 , 41 , 82 , 83 ): Tw is a measure for the local winding of the helix and directly related to the torsional strain in the molecule.Conversely, Wr corresponds to the coiling of the double helix axis in 3D space and in a plasmid is related to the number of plectonemic supercoils (and other writhed configurations, such as toroidal supercoils and curls).We express the linking number balance relative to the torsionally relaxed double-stranded DNA, for which we define Lk = Lk -Lk 0 = 0, and where Lk 0 ( = Tw 0 + Wr 0 ) is the natural linking number of torsionally relaxed DNA, Tw = Tw -Tw 0 = 0, where Tw 0 is the natural twist of DNA, equal to the number of base pairs divided by the helical turn ( ≈10.5 bp per turn for bare DNA), and Wr = Wr , i.e. the torsionally relaxed conformation has an average Wr of zero.DNA is defined to be supercoiled if Lk is different from zero.Intercalation changes the intrinsic twist of the helix and, therefore, shifts the linking number difference at which the molecule is torsionally relaxed by N bound • θ / 360 • where N bound is the number of dye molecules bound.For a given plasmid, the linking number difference relative to the torsionally relaxed state is, therefore, given by where Lk 0 is the linking number difference of the plasmid in the absence of intercalation.In general, excess linking number will partition into twist and writhe (Equation 3).For plasmids, it has been shown that the partitioning is independent of the magnitude ( 38 ) and sign ( 84 ) of the linking difference and is approximately 20% Tw and 80% Wr .We assume that this partitioning between twist and writhe also holds in the presence of intercalators and for molecules with attached ends at low forces ( 85 ,86 ), such that where Lk is given by Equation 4. In order to compute the torsional strain for a given initial linking difference Lk 0 and given number of intercalated molecules N bound , we convert the excess twist (Equation 5) to torque by taking into account the torsional stiffness of DNA: where C is the torsional stiffness of DNA in nm and L C the contour length, which in turn depends on the number of molecules bound as L C = L C,0 + N bound • z , where L C,0 is the contour length in the absence of intercalation, ≈0.34 nm per bp.The torsional stiffness of DNA has been measured using single-molecule methods ( 75 , 87 , 88 ), is independent of ionic strength ( 88 ), and reported values are in the range of C ≈ 100 nm ( 74-76 ).However, it is not well known whether or how C is altered by intercalation.Previous measurements using DNA in free solution using EtBr ( 42 , 89 , 90 ) have found lower values of the torsional stiffness in the range C ≈ 50 nm, which we take as a starting point for the EtBr data.For all other conditions we use C = 100 nm.
To determine the number of intercalated molecules per plasmid as a function of total ligand concentration c total and DNA concentration c DNA (typically expressed as the base pair concentration), we numerically solve the coupled Equations 1 -6 using an iterative approach (Materials and Methods).The main output of the model is N bound .Assuming a linear relationship between the number of intercalated molecules and the fluorescence intensity, which we have previously found to hold for a large range of dye concentrations ( 19 ), the observed fluorescence intensity I is given by where α is a proportionality factor that depends on the quantum efficiency of the dye and the details of the experiments but is constant for a given intercalator and instrumental set up.We note that the proportionality factor α cancels for our predictions of the relative binding of closed vs. open topology.

DNA supercoiling modulates intercalation
Intercalation into topologically closed DNA can increase or decrease binding due to the by the global constraint, relatively to a torsionally relaxed (or linearized) plasmid ( 26 , 27 , 33 ).
Starting with a negatively supercoiled plasmid, which is the form typically found in vivo , intercalation at low ligand concentration is increased relative to the nicked DNA, due to the negative torsional strain in the molecule (Figure 1 F).As more and more molecules intercalate, the negative linking difference is compensated until the plasmid becomes torsionally relaxed, at which point binding to the topologically closed and open forms is the same.Finally, as the intercalator concentration is increased further, the closed plasmid becomes overwound, and the positive torsional strain hinders further inter-calation.Therefore, at high intercalator concentration fewer molecules bind to the closed compared to the open plasmid (Figure 1 F).
To experimentally test the predictions of the model outlined in the previous section, we used plasmid DNA in both topologically constrained, negatively supercoiled form and in open circular and linear topologies (Figure 1 A-C; Materials and Methods).We prepared mixtures with equal amounts of the three different DNA topologies to facilitate direct comparison on a gel.We then separated the mixtures on a gel (Figure 2 A), imaged the gel, and quantified the band intensities to monitor the amount of intercalation (Supplementary Figure S3).While the open circular and linear topologies exhibit similar intensities, the topologically constrained species in comparison appears less bright on the gel for high EtBr concentrations (Supplementary Figure S4A).In contrast, for the lowest EtBr concentration, the supercoiled species exhibits a higher intensity than the other two species (Supplementary Figure S4B).We use our model to quantitatively account for the topology dependent intensities (Figure 2 B, Supplementary Figure S5A).We find that for the experimental parameters used here, the number of intercalated molecules per plasmid is approximately independent of DNA concentration (Supplementary Figure S5A).Consequently, the fluorescent intensity increases with DNA concentration (Figure 2 B).At the lowest EtBr concentration (Figure 2 B, black lines and symbols), more molecules bind to the supercoiled species (Figure 2 B, lines and symbols with red highlighting) as compared to open circular and linear, while at the highest concentration (Figure 2 B, light brown lines and symbols) intercalation is reduced for supercoiled compared to the other species.At the intermediate EtBr concentration (Figure 2 B, brown lines and symbols) the topologically open and closed species bind EtBr similarly.The differences in binding between the different topologies are a consequence of the negative torsional strain at the lowest EtBr concentration and the positive strain for the highest concentration (Supplementary Figure S5B).Averaging over the different DNA concentrations, we can quantitatively compare the relative enhancement or reduction of interaction for the topologically closed species compared to the open circular and linear species (Figure 2 C) and find excellent agreement between our model and the experimental data.We note that for low EtBr concentrations, the error bars are larger and the fit of our model appears less well than at higher concentrations, which might be in part due to (i) some presence of alternative (non-B form) DNA structures, (ii) a distribution of topoisomers for negatively coiled DNA and (iii) the lower fluorescence signal at lower EtBr concentrations.Nonetheless, we achieve reasonable agreement within experimental error for all conditions tested.

Topology dependent binding depends on initial topology and intercalator affinity
Having demonstrated that our model can quantitatively account for how DNA topology alters DNA binding starting with a negatively supercoiled plasmid ( Lk 0 < 0), where intercalation is increased at low intercalator concentrations and suppressed at high concentrations compared to an open topology (Figure 2 C), we turn to initially torsionally relaxed plasmids.A clear prediction of our model is that if the DNA is initially torsionally relaxed (or even positively supercoiled), intercalation should always be reduced for the closed topology compared to an open topology.We test this prediction experimentally by again preparing and separating DNA plasmids with different topologies, but now using a sample where the plasmid has been relaxed by topoisomerase treatment (Methods), such as that Lk 0 ≈ 0 (Supplementary Figures S6 and  S7).As predicted, we find that EtBr intercalation is reduced for the closed topology (Supplementary Figure S7A-C), in excellent agreement with our model.More broadly, the relative effect of DNA topology depends both in the initial linking number Lk 0 and on the intercalator concentration (Figure 2 D).
Our model predicts that topology dependent binding is most pronounced at ligand concentration below the K d (Figure 2 D, the K d value is indicated as a vertical line).At concentrations greater than the K d , binding saturates and the modulation by the torsional strain in the topologically closed plasmid is predicted to only lead to small or negligible changes in binding compared to the torsionally relaxed forms.To test this prediction, we carried out measurements using the intercalator SYBR Gold ( 19 , 47 , 48 ), which has a much lower K d (i.e. higher affinity) compared to EtBr (Table 1 , Figure 2 E, and Supplementary Figure S8).Performing measurements with initially negatively supercoiled or relaxed plasmids at different SYBR Gold concentrations around and above its K d , we find that indeed the topologically closed and open constructs bind similar amounts of SYBR Gold, in quantitative agreement with our model (Figure 2 E).To show the broad range of applications, we used three different assays to obtain fluorescence intensity data for SYBR Gold, namely gel electrophoresis, a well plate reader, and a qPCR cycler (see Methods for details).In addition, we performed control measurement with different plasmids, notably a shorter (2.7 kb) and longer (7.2 kb) construct, compared to the pBR322 plasmid with 4361 bp.We find similar results, within experimental error, for the different plasmids (Supplementary Figure S9), suggesting that the results do not strongly depend on DNA sequence or length.We note that our model predicts that there is no dependence of the topology effects on DNA length, which is born out experimentally in the length range investigated (Supplementary Figure S9).Importantly, the almost topologyindependent binding of DNA intercalators above their K d is advantageous for assays that aim to quantitatively compare different DNA topologies, e.g. to monitor the products of integration or topoisomerization reactions ( 61 , 65 , 91 , 92 ).In particular, for SYBR Gold DNA staining is essentially unbiased by topology for dye concentrations in the range of 1-2 μM, which is the concentration range that we previously identified as optimal for achieving a linear relationship between the fluorescence signal and the amount of DNA present ( 19 ).

Psoralen-based DNA crosslinking to detect DNA torsional tension
Intercalators of the psoralen family can crosslink DNA upon irradiation with UV light (Figure 3 A).They have been widely used in phototherapy ( 54 ) and to detect chromatin structure and the degree of DNA supercoiling in vivo ( 8 ,55-58 ).Often it is assumed as a simplification that the amount of DNA crosslinking varies linearly with the supercoiling density σ ( 56 ,58 ).Our model for intercalation under the global constraint induced by topology accurately captures the relative binding of the psoralen compound TMP to supercoiled DNA vs.nicked DNA, determined from a radioactivity assay using 3 H-labeled TMP ( 56 ) (Figure 3 B).Similarly, our model correctly predicts the degree of crosslinking induced by TMP for different supercoiling densities ( 58 ) (Figure 3 C and Supplementary Figure S10A).Importantly, the crosslinking conditions are chosen such that at most one crosslinking event per plasmid is induced, which means that only a small fraction of the intercalated TMP molecules reacts.For example, under the conditions of the data in Figure 3 C, there are > 10 TMP molecules bound (Supplementary Figure S10B), but < 1 on average react.However, the fact that the crosslinking signal is well approximated by our binding model using a proportionality constant analogous to Equation ( 7 ), suggests that crosslinking is directly proportional to binding.While TMP binding is at least approximately linear with supercoiling density in the range investigated in Figure 3 B and C, a linear relationship is only an approximation to the intrinsic exponential dependence on torque ( 26 , 27 , 33 ) and its validity is limited to relatively small supercoiling densities (Figure 3 D).

Single-molecule assay monitors supercoiling-dependent binding in real time
To explore consequences of supercoiling dependent binding at the single-molecule level, we investigated ligand binding to DNA under a topological constraint via single-molecule fluorescence imaging.In our assay, we attached DNA via multiple biotin-streptavidin bonds at each end to a surface (Figure 4 A), ensuring that the molecule is topologically constrained.While the surface attachment stretches the molecules, the effective stretching forces are very low and are estimated to be ≤ 0.1 pN under the conditions of our experiments (Supplementary Figure S11).
Adding the intercalator SYT O X Orange enables us both to induce supercoiling in the DNA and to visualize the molecules using fluorescence imaging ( 50 ,53 ).To systematically study the effect of topology, we performed two different types of experiments (Materials and methods).In the first case, we prepared negatively supercoiled DNA by first staining with the intercalative dye SYT O X Orange at a high concentration (250 nM), then attaching the DNA to the surface to topologically constrain it, and subsequently imaging at a lower concentration (50 nM).Following the reduction of SYT O X Orange concentration, plectonemic supercoils are clearly visible as bright spots that diffuse along the length of the DNA molecules, as has been observed previously ( 50 , 53 , 93 , 94 ) (Figure 4 B, C and Supplementary Movies S1-S4).
During continuous observation and laser exposure, DNA molecule will nick at some point, likely due to radicals generated by photochemical processes ( 95 ,96 ), which is usually an undesirable feature.However, here we use nicking upon illumination to our advantage since it enables us to observe both the supercoiled (closed) and subsequently the nicked (open) and finally torsionally relaxed state of the same DNA molecule.Upon nicking, the fluorescence intensity suddenly decreases significantly (Figure 4 D) and the bright spots indicative of plectonemic supercoils disappear (Figure 4 B, C and Supplementary Movie S1).This is in line with our previous observations: the negative supercoiling helps intercalation, consequently, once the molecule nicks, less SYT O X Orange binds and the fluorescence intensity decreases.
For the second type of experiment, we prepared positively supercoiled DNA by attaching the DNA to the surface in the presence of a low SYT O X Orange concentration (25 nM; Figure 4 E).We then increased the dye concentration (to 250 nM), but since positive supercoiling hinders intercalative binding to the DNA, the fluorescence intensity stays relatively low.Again, plectonemic supercoils appear as bright spots that diffuse along the DNA molecule (Figure 4 F,G and Supplementary Movie S2).After the positively supercoiled molecule is nicked, we observe that the fluorescence intensity increases (Figure 4 H).Importantly, our assay enables us to quantify the change in fluorescence intensity upon changes in topology by integrating the intensity over the entire molecule, before and after nicking (Figure 4 D and H).We find a decrease in fluorescence intensity upon nicking for initially negatively supercoiled DNA of 1.35 ± 0.18 (mean ± std; ratio supercoiled / nicked) and an increase in fluorescence intensity upon nicking for initially positively supercoiled DNA of 0.75 ± 0.09 (mean ± std; ratio supercoiled / nicked) (Figure 4 I).
To quantitatively model the changes in fluorescence upon nicking observed in situ , we used our model with the param- eters reported by Biebricher et al. ( 49 ) for the binding site size n , elongation per dye z , and binding constant K taken in 100 mM NaCl, which approximately corresponds to the ionic strength in our experiments (40 mM Tris-HCl, 2.5 mM MgCl 2 , 65 mM KCl).The unwinding angle per intercalation event θ is not known for SYT O X Orange.We assume Δθ = 19.1º,which is the value for SYBR Gold ( 19 ) since the dyes are relatively similar and also generally values in the range of about 20º are typical ( 19 , 31 , 97 ).Importantly, the model is applied here in two stages: We first compute the supercoiling density, relatively to relaxed, bare DNA and in the absence of intercalator, induced by attaching the DNA in the presence of 25 and 250 nM SYT O X Orange, which are σ = 2.7% and σ = 9.7%, respectively.We then compute binding to DNA and re-adjustment of the supercoiling level at the new SYT O X Orange concentrations used for imaging (250 and 50 nM, for which we find σ = +5.2% and σ = -3.9%),using the levels of supercoiling obtained in the first step as an input.For comparison, we compute binding to topologically open DNA, which enables us to calculate the changes in fluorescence upon nicking (Figure 4 J, K).Here, we observe an excellent agreement between the predictions of our computed model and the ex-

High-speed fluorescence tracking reveals binding dynamics
To quantitatively study the dynamics of intercalation into DNA under topological constraint, we performed singlemolecule fluorescence imaging at a 20 ms frame rate (Figure 5 ), 10 times faster than the data shown in Figure 4 .By fitting a simple kinetic model to the fluorescence intensity traces (Figure 5 A, B and E, F) at the transition between supercoiled and nicked DNA, we are able to determine the on-and off-rate of SYT O X Orange.
Our kinetic model for the total fluorescence intensity of the initially supercoiled molecules, reads as follows: I (t ) = I initial for t < t 0 and else where the initial intensity I initial , the final intensity I final , the rate k , and the time at which the intensity begins to change t 0 are fitting parameters.For the DNA molecules that are negatively supercoiled prior to nicking, we find a reduction in fluorescence intensity with overall rate k = (2.16 ± 0.38) s −1 (mean ± sem from 32 traces; Figure 5 D), which is close to the off-rate extrapolated to zero force reported by Biebricher et al. of (3.2 ± 0.8) s −1 using a single-molecule stretching assay ( 49 ).Conversely, starting with positively supercoiled DNA prior to nicking, we find an increase in intensity with an overall rate k = (1.92 ± 0.26) s −1 (mean ± sem from 22 traces; Figure 5 H).Assuming that this increase is due to binding with a simple bimolecular association rate, we find an on-rate of (7.68 ± 1.3) •10 6 M -1 s -1 (since the SYT O X Orange concentration is constant at 250 nM in this case), again in agreement to the on-rate reported by Biebricher et al. of (7.9 ± 2.6) •10 6 M -1 s -1 ( 49 ).The very good agreement of our fitted overall rates with previously published on-and off-rates for SYT O X Orange suggests that binding and dissociation of the dyes and not the relaxation of torsional strain is rate limiting for the observed changes in fluorescence intensity.This is consistent with estimates from simulations that suggest that relaxation of torsional strain in DNA occurs on ∼μs time scales for ∼kbp DNA segments ( 98 ).The disappearance of the fluorescent spots (Figure 5 C and G) allows us to estimate the dynamics of the writhe relaxation upon nicking, i.e. the time scale over which the plectonemic supercoils are resolved.We observe the disappearance of fluorescent spots over at least 3-4 frames (Supplementary Figure S12), suggesting that writhe relaxation occurs over ≥ 80 ms in our assay.The observed time scale for plectoneme disappearance is very similar to the lifetime of plectonemes before 'hopping' events that have previously been detected by fluorescence imaging ( 93 ), suggesting that writhe relaxation occurs via 'hopping' relaxation to a nick.Conversely, our measurements suggest that writhe relaxation does not occur predominantly by diffusion of plectonemes along the DNA to the nicking site, which would take on the order of τ ≈ L 2 / D ≥ 3 s , where L is the length over which diffusion occurs ( L ≥ 1 μm) and D the diffusion coefficient, which is ≤ 0.3 μm 2 / s ( 93 ).

Discussion
We have developed a quantitative model to describe DNA binding under supercoiling constraints that integrates aspects of previous models ( 26 , 27 , 33 ).Importantly, all parameters in our model have clear physical interpretations and can be determined independently (see e.g.Table 1 ), enabling quantitative predictions.Intriguingly, supercoiling DNA with a global topological constraint can increase or decrease binding, depending on the concentration and binding regimes ( 27 ) (Figures 1 F and 2 C).This is an important difference to applied stretching forces, which similarly can modulate binding in an Arrhenius-like exponential dependence ( 30 ,49 ), but will bias binding in only one direction.
Using the well-characterized and widely used intercalators SYBR Gold and ethidium bromide, we could show that topology-dependent binding depends on initial topology and intercalator affinity.In particular, we demonstrate that binding of DNA intercalators above their K d is almost supercoiling independent.With this, we can provide recommendations for optimal use of intercalative dyes to visualize DNA under a topological constraint: To avoid biases due to torque-dependent binding, a dye with a low K d should be used, ideally at concentrations well above the K d .Specifically, SYBR Gold at a concentration of 1-2 μM satisfies this criterion.Importantly, SYBR Gold at 1-2 μM concentration is also the optimal concentration range to obtain a high signal as well as a linear relation between DNA amount and fluorescence intensity ( 19 ).In general, the choice of intercalator type and concentration range is crucial for reliably topology-unbiased DNA staining and quantification and to obtain consistent results with other approaches of DNA quantification.
Moreover, our model can quantitatively account for observations made with the widely-used intercalator psoralen for UV -induced DNA crosslinking.Importantly , in applications where intercalation is used to detect supercoiling, a dye with a high K d (i.e.low affinity) and large unwinding angle is desirable.
We present a single-molecule assay to observe supercoiling dependent binding of SYT O X Orange to DNA in real time.We find good agreement between the theoretical predictions of our model and the experimentally observed changes in fluorescence intensity upon torsional relaxation.Highspeed tracking allows us to infer information about the binding dynamics of SYT O X Orange to DNA.We find overall on-and off-rates in good agreement with previously published rates for SYT O X Orange, suggesting that binding and dissociation of the dyes, rather than torsional relaxation, is rate limiting for the observed changes in fluorescence intensity.
Taken together, our work shows how combining theoretical modeling and multiple complementary experimental techniques can provide a highly quantitative and comprehensive view of DNA-ligand interactions under a global topological constraint.As future directions, our modeling approach could be extended to take into account topoisomer distributions, non-canonical DNA structures (e.g.flipped bases or melting bubbles), or sequence dependencies ( 53 ,99-101 ).We anticipate our approach to be broadly applicable to other DNA binding agents and allow for reliable and unbiased detection and quantification of different topological states of DNA.In addition, the interplay of DNA binding and topology has important implications for the processing and regulation of genetic information in vivo .In the cell, DNA forms topological domains and e.g.advancing polymerases introduce torsional strains of different handedness ( 1 , 2 , 102-105 ).Together, these factors are expected to modulate binding to DNA, e.g. by nucleosomes and transcription factors, also in a cellular context.Our results provide a baseline to quantitatively investigate these complex processes in the future.

Figure 1 .
Figure 1.Ov ervie w of different topological conf ormations of plasmid DNA and outline of the binding model under global constraint.( A ) Schematic and an AFM height image of linear DNA.( B ) Schematic and an AFM height image of open circular DNA, i.e. of plasmid DNA that is nicked at a single site.( C ) Schematic and an AFM height image of a negatively supercoiled DNA plasmid.AFM images in panels A-C are of pBR322 plasmid DNA ( 4361 bp; see Materials and Methods ) .( D ) Schematic of the McGhee-von Hippel binding model for ligand binding to DNA, whereby the binding site size modulates the binding equilibrium.( E ) Illustration of how intercalation into DNA lengthens and underwinds the B-form helix.Left: B-form DNA, rendered from PDB entry 4C64 ( 108 ) .Right: DNA in the presence of an intercalator, here daunom y cin rendered from PDB entry 1D11 ( 110 ) .( F ) Schematic of how the global constraint due to topology leads to increased binding as long as the topologically closed DNA is negatively supercoiled ( left ) , but will decrease binding when the DNA is positively supercoiled ( right ) , compared to the torsionally relaxed form ( center ) .Intercalation, in turn, locally underwinds the DNA and, therefore, leads to an increase in Wr with increasing intercalation in topologically closed DNA.For comparison, a topologically open DNA is shown below.

Figure 2 .
Figure 2. DNA topology dependent intercalation.( A ) Agarose gel stained with EtBr at a final concentration of 0.5 μM.Different DNA topologies are separated on the gel.L: DNA size ladders (1 kb gene ruler, Thermo Scientific, 5 μl).Lanes 2-4 are the stock solutions of negatively supercoiled, linear, and open circular DNA, respectively.Lanes 5-11 are equimolar mixtures of the three topologies, at different total DNA concentrations.( B ) Experimentally determined fluorescence intensity for supercoiled DNA (circles with red highlight) and linear DNA (squares) as a function EtBr and DNA concentration for pBR322 DNA (4361 bp).Different colors correspond to different EtBr concentrations: from dark to light 0.05, 0.5 and 5 μM.Thin lines are for topologically open DNA (linear and open circular); thick lines with red highlights are for topologically closed DNA (supercoiled, here with supercoiling density σ≈ −5%, corresponding to Lk 0 ≈ −20 t urns).S ymbols are the mean and std from at least two gels.Lines are predictions of our binding model (same color code as in panel B), with the scale factor α as the only fitting parameter (Equation 7 ).( C ) R elativ e fluorescence intensity of a topologically closed DNA relative to the topologically open constructs.Data points are obtained by averaging the different DNA concentration at the same EtBr condition.Triangles are experimental data from at least two independent gels.Magenta symbols are the prediction of or model.( D ) R elativ e binding to topologically closed DNA with different le v els of initial supercoiling compared to topologically open constructs as a function of EtBr concentration.Colored lines are the predictions of our model.The v ertical gre y line indicates the K d of EtBr.Symbols are the data from Figure 2 E. Further analyses are shown in Supplementary Figure S5.( E ) Same as in panel D for SYBR Gold.The data shown are from three different experimental modalities: fluorescent readout using a 96-well platereader, a qPCR cycler, and gel electrophoresis.Due to the much lo w er K d f or SYBR Gold compared to EtBr, the binding is almost independent of topology for the SYBR Gold conditions investigated.

Figure 3 .
Figure 3. Topology-dependent DNA crosslinking by trimethylpsoralen.( A ) 4,5 ,8-trimethylpsoralen (TMP) intercalates into DNA and causes DNA crosslinking upon UV irradiation.DNA supercoiling-dependent crosslinking is widely used to probe DNA supercoiling and chromatin conformations in vivo .( B ) Binding of TMP to supercoiled DNA plasmids with different initial supercoiling densities relative to open circular DNA.Experimental binding data are from Ref. ( 56 ) and were determined using the radioactivity of 3 H-labeled TMP. ( C ) Binding of TMP to supercoiled DNA plasmids.Experimental data are from ( 58 ) and were obtained by quantifying the amount of DNA crosslinking after irradiation.Experiment al dat a are normalized to the dat a point at zero supercoiling density.The model in panels B and C uses the parameters in Table 1 and quantified binding to supercoiled relative to topologically open DNA.( D ) Dependence of the torque-dependent association constant (the in v erse of the dissociation constant) on supercoiling density using Equations ( 2 )-( 6 ) (red solid line).The black dashed line shows the linearization of Equation ( 2 ), i.e. the approximation exp( −x ) ≈ 1 -x , with x = Γ • θ/ k B T.

Figure 4 .
Figure 4. Single-molecule fluorescence assay to quantify topology dependent binding.A)-D): Negatively supercoiled DNA.( A ) Schematic representation of the experimental set up of the single-molecule fluorescence experiments.The SYTOX Orange-stained DNA is tethered at both its ends via multible biotin-streptavidin bonds to the surface.Top: negatively supercoiled DNA; bottom: nicked DNA.( B ) Fluorescence image snapshots at 40.6 s when the DNA is still negatively supercoiled and at 76 s after the DNA was nicked.( C ) Kymograph of SYTOX-Orange stained-DNA.When the DNA is nicked, the intensity decreases abruptly.( D ) Integrated fluorescence intensity of the kymograph shown in panel C.The dotted orange line is a fit of a t wo-st ate hidden Mark o v model to the data ( 109 ).( E-H ) Positiv ely supercoiled DNA.(E) Schematic representation of the e xperimental lineup as in panel A only now the DNA is positively supercoiled before the nicking.( F ) Fluorescence image snapshots at 19 s when the DNA is still positively supercoiled and at 82 s after the DNA was nicked.( G ) Kymograph of SYTOX-Orange stained-DNA.When the DNA is nicked, the intensity increases abruptly.(H) Fluorescence intensity of the kymograph shown in panel G.The dotted orange line is a fit (same as panel D) to the data.( I ) Intensity ratios before and af ter nic king f or originally negativ ely supercoiled DNA (blue) and f or originally positiv ely supercoiled DNA (turquoise).Av eraging giv es 1.35 ± 0.18 ( N = 51, mean ± std; ratio sc / nicked) for originally negatively supercoiled DNA and 0.75 ± 0.09 ( N = 51, mean ± std; ratio supercoiled / nicked) for originally positively supercoiled DNA.( J ) Comparison of the experimental intensity ratio to the value from theoretical modeling for originally negatively supercoiled DNA from N = 51 molecules.( K ) Comparison of the experimental intensity ratio to the value from theoretical modelling for originally positively supercoiled DNA from N = 51 molecules.

Figure 5 .
Figure 5. Probing the interplay of dye and supercoil dynamics.(A-D): Negatively supercoiled DNA.( A ) Kymograph of SYTOX-Orange stained-DNA.When the DNA is nicked, the intensity decreases abruptly.( B ) Integrated fluorescence intensity of the kymograph shown in panel A. The dotted orange line is a fit of the model shown in Equation 8 to the data.( C ) Fluorescence image snapshots at 7.2 s when the DNA is still negatively supercoiled, at 8.2 s when the nick occurs, and at 9.2 s after the DNA was nicked.Time points are indicated by matching numbers in panel B. ( D ) Experimentally determined rates k for fluorescence decrease from the fits of Equation 8 to time traces from N = 32 independent measurements.The mean ± sem are (2.16 ± 0.38) s −1 .(E-H): Positively supercoiled DNA .( E ) Kymograph of SYTOX-Orange stained-DNA.When the DNA is nicked, the intensity increases abruptly.( F ) Integrated fluorescence intensity of the kymograph shown in panel E. The dotted orange line is a fit of Equation 8 to the data.( G ) Fluorescence image snapshots at 6.7 s when the DNA is still positively supercoiled, at 7.7 s when the nick occurs, and at 8.7 s after the DNA w as nick ed.Time points are indicated by matching numbers in panel F. ( H ) Experimentally determined rates k for the increase in fluorescence after nicking from N = 22 independent measurements.The mean ± sem are (1.92 ± 0.26) s −1 .
• C ) to prepare linear DNA and Nt.BspQI ( NEB; incubation temperature 50 • C ) to create open circular DNA.The reactions were stopped after one hour by heat inactivation at

Table 1 .
Parameters of selected intercalators used in this work We used the values from ( 49 ) in 100 mM NaCl, close to the ionic strength used in this work.The value for Δθ is an estimate based on data for SYBR Gold.b Ref. ( 106 ) only determined the dissociation constant approximately.The binding site size and elongation per dye for TMP were assumed to assume the average values for intercalators indicated in the table. a