Analysis of the structure of dimeric DNA catenanes by electron microscopy

Influence of DNA sequence on the structure of minicircles under torsional stress

The sequence dependence of the conformational distribution of DNA under various levels of torsional stress is an important unsolved problem. Combining theory and coarse-grained simulations shows that the DNA sequence and a structural correlation due to topology constraints of a circle are the main factors that dictate the 3D structure of a 336 bp DNA minicircle under torsional stress. We found that DNA minicircle topoisomers can have multiple bend locations under high torsional stress and that the positions of these sharp bends are determined by the sequence, and by a positive mechanical correlation along the sequence. We showed that simulations and theory are able to provide sequence-specific information about individual DNA minicircles observed by cryo-electron tomography (cryo-ET). We provided a sequence-specific cryo-ET tomogram fitting of DNA minicircles, registering the sequence within the geometric features. Our results indicate that the conformational distribution of minicircles under torsional stress can be designed, which has important implications for using minicircle DNA for gene therapy.

Intercalation-Based Single-Molecule Fluorescence Assay To Study DNA Supercoil Dynamics

DNA supercoiling crucially affects cellular processes such as DNA replication, gene expression, and chromatin organization. However, mechanistic understanding of DNA supercoiling and the related DNA-processing enzymes has remained limited, mainly due to the lack of convenient experimental tools to probe these phenomena. Here, we report a novel high-throughput single-molecule assay for real-time visualization of supercoiled DNA molecules, named ISD (Intercalation-induced Supercoiling of DNA). We use an intercalating dye to induce supercoiling of surface-attached DNA molecules as well as to visualize coiled-loop structures (i.e., plectonemes) formed on DNA. The technique is solely based on epifluorescence microscopy and requires no mechanical manipulation of the DNA molecules. This new assay allows to track positions and sizes of individual plectonemes and characterize their position-dependent dynamics such as nucleation, termination, and diffusion. We describe the ISD technique and demonstrate its potential by establishing that plectonemes are pinned to a local 10-nucleotide long mispaired sequence along a double-stranded DNA molecule.

Structural diversity of supercoiled DNA

By regulating access to the genetic code, DNA supercoiling strongly affects DNA metabolism. Despite its importance, however, much about supercoiled DNA (positively supercoiled DNA, in particular) remains unknown. Here we use electron cryo-tomography together with biochemical analyses to investigate structures of individual purified DNA minicircle topoisomers with defined degrees of supercoiling. Our results reveal that each topoisomer, negative or positive, adopts a unique and surprisingly wide distribution of three-dimensional conformations. Moreover, we uncover striking differences in how the topoisomers handle torsional stress. As negative supercoiling increases, bases are increasingly exposed. Beyond a sharp supercoiling threshold, we also detect exposed bases in positively supercoiled DNA. Molecular dynamics simulations independently confirm the conformational heterogeneity and provide atomistic insight into the flexibility of supercoiled DNA. Our integrated approach reveals the three-dimensional structures of DNA that are essential for its function.

DNA supercoiling strongly affects its metabolism. By electron cryo-tomography, biochemical assays and molecular dynamics simulations, here the authors show that supercoiled DNA minicircles adopt unique and wide distributions of three-dimensional conformations, many with disrupted base pairs.

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.

Effects of Knots on Ring Polymers in Solvents of Varying Quality

We employ extensive computer simulations to investigate the conformations and the interactions of ring polymers under conditions of worsening solvent quality, in comparison with those for linear polymers. We determine the dependence of the Θ-temperature on knotedness by considering ring polymers of different topologies. We establish a clear decrease of the former upon changing the topology of the polymer from linear to an unknotted ring and a further decrease of the same upon introducing trefoil- or 5-fold knots but we find no difference in the Θ-point between the two knotted molecules. Our results are based on two independent methods: one considering the scaling of the gyration radius with molecular weight and one based on the dependence of the effective interaction on solvent quality. In addition, we calculate several shape-parameters of the polymers to characterize linear, unknotted, and knotted topologies in good solvents and in the proximity of the Θ-point. The shape parameters of the knotted molecules show an interesting crossover at a degree of polymerization that depends on the degree of knottedness of the molecule.

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.

Unlinking of supercoiled DNA catenanes by type IIA topoisomerases

It was found recently that DNA catenanes, formed during replication of circular plasmids, become positively (+) supercoiled, and the unlinking of such catenanes by type IIA topoisomerases proceeds much more efficiently than the unlinking of negatively (-) supercoiled catenanes. In an attempt to explain this striking finding we studied, by computer simulation, conformational properties of supercoiled DNA catenanes. Although the simulation showed that conformational properties of (+) and (-) supercoiled replication catenanes are very different, these properties per se do not give any advantage to (+) supercoiled over (-) supercoiled DNA catenanes for unlinking. An advantage became evident, however, when we took into account the established features of the enzymatic reaction catalyzed by the topoisomerases. The enzymes create a sharp DNA bend in the first bound DNA segment and allow for the transport of the second segment only from inside the bend to its outside. We showed that in (-) supercoiled DNA catenanes this protein-bound bent segment becomes nearly inaccessible for segments of the other linked DNA molecule, inhibiting the unlinking.

Interplay of DNA supercoiling and catenation during the segregation of sister duplexes

The discrete regulation of supercoiling, catenation and knotting by DNA topoisomerases is well documented both in vivo and in vitro, but the interplay between them is still poorly understood. Here we studied DNA catenanes of bacterial plasmids arising as a result of DNA replication in Escherichia coli cells whose topoisomerase IV activity was inhibited. We combined high-resolution two-dimensional agarose gel electrophoresis with numerical simulations in order to better understand the relationship between the negative supercoiling of DNA generated by DNA gyrase and the DNA interlinking resulting from replication of circular DNA molecules. We showed that in those replication intermediates formed in vivo, catenation and negative supercoiling compete with each other. In interlinked molecules with high catenation numbers negative supercoiling is greatly limited. However, when interlinking decreases, as required for the segregation of newly replicated sister duplexes, their negative supercoiling increases. This observation indicates that negative supercoiling plays an active role during progressive decatenation of newly replicated DNA molecules in vivo.

Spectroscopic methods for the physical characterization and formulation of nonviral gene delivery systems

Currently, with the exception of naked DNA formulations, most pharmaceutical preparations of plasmid DNA employ some type of polycationic delivery vector such as synthetic cationic polymers and lipids to enhance delivery. A number of biophysical techniques are readily available for the structural characterization of plasmid DNA within synthetic gene delivery complexes. Here we present applications of ultraviolet (UV) absorption, circular dichroism (CD), infrared (IR), and fluorescence spectroscopies as well as dynamic light scattering to the structural analysis of the oligonucleotide component of nonviral gene delivery vectors. We also illustrate this approach for the investigation of the formulation of lipoplex and polyplex-based gene delivery systems. To summarize such data, we show how the macromolecular complexes can be represented as vectors in a highly dimensional space in which the components of the vector consist of normalized values of experimental parameters measured as a function of different solution conditions such as pH and ionic strength.

Restriction endonucleases that bridge and excise two recognition sites from DNA

Most restriction endonucleases bridge two target sites before cleaving DNA: examples include all of the translocating Type I and Type III systems, and many Type II nucleases acting at their sites. A subset of Type II enzymes, the IIB systems, recognise bipartite sequences, like Type I sites, but cut specified phosphodiester bonds near their sites, like Type IIS enzymes. However, they make two double-strand breaks, one either side of the site, to release the recognition sequence on a short DNA fragment; 34 bp long in the case of the archetype, BcgI. It has been suggested that BcgI needs to interact with two recognition sites to cleave DNA but whether this is a general requirement for Type IIB enzymes had yet to be established. Ten Type IIB nucleases were tested against DNA substrates with one or two copies of the requisite sequences. With one exception, they all bridged two sites before cutting the DNA, usually in concerted reactions at both sites. The sites were ideally positioned in cis rather than in trans and were bridged through 3-D space, like Type II enzymes, rather than along the 1-D contour of the DNA, as seen with Type I enzymes. The standard mode of action for the restriction enzymes that excise their recognition sites from DNA thus involves concurrent action at two DNA sites.

Long-range communications between DNA sites by the dimeric restriction endonuclease SgrAI

The SgrAI endonuclease displays its maximal activity on DNA with two copies of its recognition sequence, cleaving both sites concertedly. While most restriction enzymes that act concurrently at two sites are tetramers, SgrAI is a dimer in solution. Its reaction at two cognate sites involves the association of two DNA-bound dimers. SgrAI can also bridge cognate and secondary sites, the latter being certain sequences that differ from the cognate by one base-pair. The mechanisms for cognate-cognate and cognate-secondary communications were examined for sites in the following topological relationships: in cis, on plasmids with two sites in a single DNA molecule; on catenanes containing two interlinked rings of DNA with one site in each ring; and in trans, on oligoduplexes carrying either a single site or the DNA termini generated by SgrAI. Both cognate-cognate and cognate-secondary interactions occur through 3-D space and not by 1-D tracking along the DNA. Both sorts of communication arise more readily when the sites are tethered to each other, either in cis on the same molecule of DNA or by the interlinking of catenane rings, than when released from the tether. However, the dimer bound to an oligoduplex carrying either a cognate or a secondary site could be activated to cleave that duplex by interacting with a second dimer bound to the recognition site, provided both duplexes are at least 30 base-pairs long: the second dimer could alternatively be bound to the two duplexes that correspond to the products of DNA cleavage by SgrAI.

Directed polymers with constrained winding angle

In this paper we study from a nonperturbative point of view the entanglement of two directed polymers subjected to repulsive interactions given by a Dirac delta-function potential. An exact formula of the so-called second moment of the winding angle is derived. This result is used to provide a thorough analysis of entanglement phenomena in the classical system of two polymers subjected to repulsive interactions and related problems. No approximation is made in treating the constraint on the winding angle and the repulsive forces. In particular, we investigate how repulsive forces influence the entanglement degree of the two-polymer system. In the limit of ideal polymers, in which the interactions are switched off, we show that our results are in agreement with those of previous works.

Effective and robust plasmid topology analysis and the subsequent characterization of the plasmid isoforms thereby observed

Within the biopharmaceutical industry, recombinant plasmid DNA is used both as a raw material (e.g. in lentiviral and AAV vector production) as well as an active ingredient (e.g. in DNA vaccines). Consequently, many analytical laboratories are routinely involved with plasmid DNA topoisoform qualitative analysis and quantification. In order to reliably determine plasmid topology, one must ensure that the methodology employed can reliably, precisely and accurately measure qualitatively and quantitatively all topological isoforms. Presented here are an anion-exchange high-performance liquid chromatography (AEC) and an agarose gel electrophoresis (AGE)-based method developed for this purpose. The strategies undertaken to overcome the respective typical problems of limited linear range of quantitation (for AGE) and isoform resolution (for AEC) are described. Also presented is a subsequent direct comparison (for assay precision/accuracy) of these two methods, as well as a package of species characterization [by chloroquine-AGE, enzymatic digestion, multi-angle laser light-scattering (MALLS) and electron microscopy] undertaken to confirm the identity of a minor supercoiled dimeric concatamer observed by both approaches.

DNA communications by Type III restriction endonucleases--confirmation of 1D translocation over 3D looping

DNA cleavage by Type III restriction enzymes is governed strictly by the relative arrangement of recognition sites on a DNA substrate--endonuclease activity is usually only triggered by sequences in head-to-head orientation. Tens to thousands of base pairs can separate these sites. Long distance communication over such distances could occur by either one-dimensional (1D) DNA translocation or 3D DNA looping. To distinguish between these alternatives, we analysed the activity of EcoPI and EcoP15I on DNA catenanes in which the recognition sites were either on the same or separate rings. While substrates with a pair of sites located on the same ring were cleaved efficiently, catenanes with sites on separate rings were not cleaved. These results exclude a simple 3D DNA-looping activity. To characterize the interactions further, EcoPI was incubated with plasmids carrying two recognition sites interspersed with two 21res sites for site-specific recombination by Tn21 resolvase; inhibition of recombination would indicate the formation of stable DNA loops. No inhibition was observed, even under conditions where EcoPI translocation could also occur.

Protein motion from non-specific to specific DNA by three-dimensional routes aided by supercoiling

DNA-binding proteins are generally thought to locate their target sites by first associating with the DNA at random and then translocating to the specific site by one-dimensional (1D) diffusion along the DNA. We report here that non-specific DNA conveys proteins to their target sites just as well when held near the target by catenation as when co-linear with the target. Hence, contrary to the prevalent view, proteins move from random to specific sites primarily by three-dimensional (3D) rather than 1D pathways, by multiple dissociation/re-association events within a single DNA molecule. We also uncover a role for DNA supercoiling in target-site location. Proteins find their sites more readily in supercoiled than in relaxed DNA, again indicating 3D rather than 1D routes.

Investigation of viral DNA packaging using molecular mechanics models

A simple molecular mechanics model has been used to investigate optimal spool-like packing conformations of double-stranded DNA molecules in viral capsids with icosahedral symmetry. The model represents an elastic segmented chain by using one pseudoatom for each ten basepairs (roughly one turn of the DNA double helix). Force constants for the various terms in the energy function were chosen to approximate known physical properties, and a radial restraint was used to confine the DNA into a sphere with a volume corresponding to that of a typical bacteriophage capsid. When the DNA fills 90% of the spherical volume, optimal packaging is obtained for coaxially spooled models, but this result does not hold when the void volume is larger. When only 60% of the spherical volume is filled with DNA, the lowest energy structure has two layers, with a coiled core packed at an angle to an outer coaxially spooled shell. This relieves bending strain associated with tight curvature near the poles in a model with 100% coaxial spooling. Interestingly, the supercoiling density of these models is very similar to typical values observed in plasmids in bacterial cells. Potential applications of the methodology are also discussed.

DNA cleavage reactions by type II restriction enzymes that require two copies of their recognition sites

Several type II restriction endonucleases interact with two copies of their target sequence before they cleave DNA. Three such enzymes, NgoMIV, Cfr10I and NaeI, were tested on plasmids with one or two copies of their recognition sites, and on catenanes containing two interlinked rings of DNA with one site in each ring. The enzymes showed distinct patterns of behaviour. NgoMIV and NaeI cleaved the plasmid with two sites faster than that with one site and the catenanes at an intermediate rate, while Cfr10I gave similar steady-state rates on all three substrates. Both Cfr10I and NgoMIV converted the majority of the substrates with two sites directly to the products cut at both sites, while NaeI cleaved just one site at a time. All three enzymes thus synapse two DNA sites through three-dimensional space before cleaving DNA. With Cfr10I and NgoMIV, both sites are cleaved in one turnover, in a manner consistent with their tetrameric structures, while the cleavage of a single site by NaeI indicates that the second site acts not as a substrate but as an activator, as reported previously. The complexes spanning two sites have longer lifetimes on catenanes with one site in each ring than on circular DNA with two sites, which indicates that the catenanes have more freedom for site juxtaposition than plasmids with sites in cis.

Polymer models of meiotic and mitotic chromosomes

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush-like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein-chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A + T- and G + C-rich regions revealed by chromosome banding.

The effect of ionic conditions on the conformations of supercoiled DNA. II. Equilibrium catenation

We studied the equilibrium formation of DNA catenanes to assess the conformational properties of supercoiled DNA as a function of ionic conditions and supercoiling density. Catenanes were formed by cyclizing linear DNA with long cohesive ends in the presence of supercoiled molecules. The efficiency of the catenation depends on the distance between opposing segments of DNA in the interwound superhelix. The fraction of cyclizing molecules that becomes topologically linked with the supercoiled DNA is the product of the concentration of the supercoiled DNA and a proportionality constant, B, that depends on conformations of supercoiled DNA. In parallel with these experimental studies, we calculated the values of B using Monte Carlo simulations of the equilibrium distribution of DNA conformations. There were no adjustable parameters in the calculations because all three parameters of the DNA model, bending and torsional elasticity of DNA and DNA effective diameter, specifying intersegment interactions, were known from independent studies. We found very good agreement between measured and simulated values of B for all the ionic conditions and DNA superhelix densities studied; the discrepancy was less than a factor of 2 over the 200-fold variation in B. The value of B decreases nearly exponentially with increasing superhelicity, this dependence being especially strong at low salt concentration. The dependence of B on the concentration of NaCl, MgCl(2), and spermidine can be described with good accuracy in terms of changes of the DNA effective diameter. We found no indication of superhelix collapse under any ionic conditions studied. We discuss, in light of these results, the biological importance of the effect of DNA supercoiling on the unlinking of the products of DNA replication.

Contrasting enzymatic activities of topoisomerase IV and DNA gyrase from Escherichia coli

DNA gyrase and topoisomerase IV (Topo IV) have distinct roles as unlinking enzymes during DNA replication despite 40% sequence identity between them. DNA gyrase unlinks replicating DNA by introducing negative supercoils while Topo IV decatenates the two daughter molecules. For this study, we measured the rates of unlinking of various topoisomers of DNA by DNA gyrase and Topo IV. Each enzyme has marked preferences for certain strand-passage reactions. DNA gyrase is a relatively poor decatenase, catalyzing strand-passage events that result in supercoiling at rates several orders of magnitude faster than those causing decatenation. Topo IV, in contrast, decatenates linked circles 10-40 times more quickly than it removes the intramolecular crossings from supercoiled DNA. Supercoiled catenanes are unlinked at an even more increased rate by Topo IV. Thus, the supercoils augment decatenation rather than compete with catenane crossings for their removal. Knot crossings and the crossings of multiply interlinked catenanes are also preferentially removed by Topo IV. This ability of Topo IV to selectively unlink catenated molecules mirrors its key role in decatenation of replicating chromosomes in vivo.

Repercussions of DNA tracking by the type IC restriction endonuclease EcoR124I on linear, circular and catenated substrates

Type I restriction endonucleases such as EcoR124I cleave DNA at undefined loci, distant from their recognition sequences, by a mechanism that involves the enzyme tracking along the DNA between recognition and cleavage sites. This mechanism was examined on plasmids that carried recognition sites for EcoR124I and recombination sites for resolvase, the latter to create DNA catenanes. Supercoiled substrates with either one or two restriction sites were linearized by EcoR124I at similar rates, although the two-site molecule underwent further cleavage more readily than the one-site DNA. The catenane from the plasmid with one EcoR124I site, carrying the site on the smaller of the two rings, was cleaved by EcoR124I exclusively in the small ring, and this underwent multiple cleavage akin to the two-site plasmid. Linear substrates derived from the plasmids were cleaved by EcoR124I at very slow rates. The communication between recognition and cleavage sites therefore cannot stem from random looping. Instead, it must follow the DNA contour between the sites. On a circular DNA, the translocation of non-specific DNA past the specifically bound protein should increase negative supercoiling in one domain and decrease it in the other. The ensuing topological barrier may be the trigger for DNA cleavage.

Recombination by resolvase to analyse DNA communications by the SfiI restriction endonuclease

The SfiI endonuclease differs from other type II restriction enzymes by cleaving DNA concertedly at two copies of its recognition site, its optimal activity being with two sites on the same DNA molecule. The nature of this communication event between distant DNA sites was analysed on plasmids with recognition sites for SfiI interspersed with recombination sites for resolvase. These were converted by resolvase to catenanes carrying one SfiI site on each ring. The catenanes were cleaved by SfiI almost as readily as a single ring with two sites, in contrast to the slow reactions on DNA rings with one SfiI site. Interactions between SfiI sites on the same DNA therefore cannot follow the DNA contour and, instead, must stem from their physical proximity. In buffer lacking Mg2+, where SfiI is inactive while resolvase is active, the addition of SfiI to a plasmid with target sites for both proteins blocked recombination by resolvase, due to the restriction enzyme bridging its sites and thus isolating the sites for resolvase into separate loops. The extent of DNA looping by SfiI matched its extent of DNA cleavage in the presence of Mg2+.