Topoisomer gel retardation: detection of anti-Z-DNA antibodies bound to Z-DNA within supercoiled DNA minicircles.

Small DNA fragments of approximately 350 bp in length, either with or without d(CG)n tracts, are ligated into underwound DNA minicircles to generate topoisomeric rings with different topological linking numbers, Lk. These minicircles, differing by an Lk of one, can be separated by acrylamide gel electrophoresis. Furthermore, electrophoresis can be used to reveal DNA double helix conformational changes that are induced by supercoiling, such as left-handed Z-DNA. When anti-Z-DNA antibodies are added to such minicircles, their binding leads to a selective retardation of the electrophoretic migration of the Z-DNA containing circles. This effect is not seen with relaxed minicircles and those with insufficient torsional stress to induce a conformational transition. Thus the technique of 'topoisomer gel retardation' presents a very sensitive assay for the identification of proteins that selectively bind to DNA conformations stabilized by negative DNA supercoiling.     

Effect of salt-dependent stiffness on the conformation of a stressed DNA loop containing initially coplanar bends

Closed DNA loops containing one or more bent regions are important structures that occur in the regulation of gene expression. We analyze the response of structures of this type to a change in applied rotation (change in linking deficiency, delta Lk). Our results apply to a closed loop formed from an elastic rod that is intrinsically bent in Nb discrete, 20 degrees steps up to a maximum of 240 degrees, the bent regions being initially coplanar with the plane of the relaxed DNA loop. We determine the effect of changing the intrinsic elastic resistance of the DNA loop to bending and torsional deformations. This relative resistance is expressed by Poisson's ratio v, which depends upon the ratio of bending stiffness to torsional rigidity. Poisson's ratio is primarily a function of salt type and concentration. We find that the tertiary structure of DNA loops changes with delta Lk, but that the geometric response can be either of two quite different types, depending upon the precise (Nb, v) pair. For combinations of Nb and v that are above a critical curve (the Fickel curve), the response to increasing delta Lk is nonmonotonic (NMT region): the distance between the loop closure point and its diametric opposite first decreases, then increases, as delta Lk increases. For combinations of Nb and v that are below the Fickel curve (NMT region), the corresponding diameter never increases, but always decreases with increasing delta Lk. In addition to these results, we define and implement a new measure of tertiary structure in closed DNA: the absolute writhe, AWr.     

Effect of nucleosome distortion on the linking deficiency in relaxed minichromosomes

The wrapping of closed circular DNA on a protein surface, followed by relaxation with a topoisomerase and removal of proteins, produces a characteristic DNA linking deficiency, delta Lk. We show that the magnitude of delta Lk depends upon the surface shape, and we calculate changes in delta Lk caused by particular distortions of the protein wrapping surface. If the DNA remains attached to the surface during distortion, the DNA winding number, phi, is not altered. The change in delta Lk is then equal to the change in the surface linking number, SLk, which is a straightforward measure of the wrapping of the DNA around the surface. For left-handed wrapping, as in a nucleosome, SLk = -n, the number of times that the DNA axis winds around the axis of the protein complex. We calculate values of SLk for the helical wrapping of a constant length of DNA on protein surfaces having the shapes of cylinders and of ellipsoids and hyperboloids of revolution. If the equatorial radius of the protein is fixed, change in shape from a cylinder to a hyperboloid increases SLk, while the corresponding change to an ellipsoid reduces SLk. We apply the general results to the interpretation of experiments in which minichromosomes are relaxed with topoisomerase at various temperatures and delta Lk is determined. The result is that a distortion of the nucleosome core by at most 5% (the change in the radius at the axial extremity relative to the equator) is sufficient to explain the observed delta Lk changes.     

DNA Topological Context Affects Access to Eukaryotic DNA Topoisomerase I

Abstract

We have analyzed the reactivity of a 217 base pair segment of the intrinsically curved Crithidia fasciculata kinetoplast DNA towards eukaryotic DNA topoisomerase I. The substrates were open [linear fragment and nicked circle] and closed minidomains [closed relaxed circle and circles with linking differences of ?1 and ?2], We interpreted the results with the aid of a model that was used to predict the structures of the topoisomers. The modelling shows that the ΔLk(?l) form is unusually compact because of the curvature in the DNA. To determine the role of sequence-directed curvature in both the experimental and modeling studies, controls were examined in which the curved Crithidia sequence was replaced by an uncurved sequence obtained from the plasmid pBR322.

Reactivity of the Crithidia DNA [as analyzed both by the cleavage and the topoisomerization reactions] markedly varied among the DNA forms: (i) the hierarchy of overall reactivity observed is: linear fragment > nicked circular, closed circular [ΔLk(O)], interwound [ΔLk(?2)] > bent interwound [ΔLk(?l)]; (ii) the intensity of several cleavage positions differs among DNA forms.

The results show that eukaryotic DNA topoisomerase I is very sensitive to the conformation of the substrates and that its reactivity is modulated by the variation of the compactness of the DNA molecule. The C. fasciculata sequence contains a highly curved segment that determines the conformation of the closed circle in a complex way.     

The elastic rod model for DNA and its application to the tertiary structure of DNA minicircles in mononucleosomes.

Explicit solutions to the equations of equilibrium in the theory of the elastic rod model for DNA are employed to develop a procedure for finding the configuration that minimizes the elastic energy of a minicircle in a mononucleosome with specified values of the minicircle size N in base pairs, the extent w of wrapping of DNA about the histone core particle, the helical repeat h(0)b of the bound DNA, and the linking number Lk of the minicircle. The procedure permits a determination of the set Y(N, w, h(0)b) of integral values of Lk for which the minimum energy configuration does not involve self-contact, and graphs of writhe versus w are presented for such values of Lk. For the range of N of interest here, 330 < N < 370, the set Y(N, w, h(0)b) is of primary importance: when Lk is not in Y(N, w, h(0)b), the configurations compatible with Lk have elastic energies high enough to preclude the occurrence of an observable concentration of topoisomer Lk in an equilibrium distribution of topoisomers. Equilibrium distributions of Lk, calculated by setting differences in the free energy of the extranucleosomal loop equal to differences in equilibrium elastic energy, are found to be very close to Gaussian when computed under the assumption that w is fixed, but far from Gaussian when it is assumed that w fluctuates between two values. The theoretical results given suggest a method by which one may calculate DNA-histone binding energies from measured equilibrium distributions of Lk.     

Modeling protein-induced configurational changes in DNA minicircles

A method is offered for obtaining minimum energy configurations of DNA minicircles constrained by one or more DNA-binding proteins. The minicircles are modeled as elastic rods, while the presence of bound protein is implied by rigidly fixing portions of these chains. The configurations of the geometrically constrained circular rods are sampled stochastically and optimized according to a simple elastic energy model of nicked DNA. The shapes of the minimum energy structures identified after a simulated annealing process are analyzed in terms of relative protein orientation and writhing number. The procedure is applied to minicircles 500 base pairs in length, bound to two evenly spaced DNA-wrapping proteins. The presence of histone octamers is suggested by rigidly fixing the two protein-bound portions of each minicircle as small superhelices similar in dimension to nucleosomal DNA. The folded minimum energy forms of sample chains with different degrees of protein wrapping are noteworthy in themselves in that they offer a new resolution to the well-known minichromosome linking number paradox and point to future minicircle simulations of possible import. © 1997 John Wiley & Sons, Inc.     

Facile Characterization of Topology of DNA Catenanes

During the preparation of single-stranded DNA catenanes, topological isomers of different linking numbers (Lk) are intrinsically produced, and they must be separated from each other to construct sophisticated nanostructures accurately. In many previous studies, however, mixtures of these isomers were directly employed to construct nanostructures without sufficient characterization. Here, we present a method that easily and clearly characterizes the isomers by polyacrylamide gel electrophoresis. To the mixtures of topological isomers of [2]catenanes, two-strut oligonucleotides, which are complementary with a part of both rings, were added to connect the rings and fix the whole conformations of isomers. As a result, the order of migration rate was always Lk3 > Lk2 > Lk1, irrespective of gel concentration. Thus, all the topological isomers were unanimously characterized by only one polyacrylamide gel electrophoresis experiment. Well-characterized DNA catenanes are obtainable by this two-strut strategy, opening the way to more advanced nanotechnology.     

Interactions of Drosophila DNA topoisomerase II with left-handed Z-DNA in supercoiled minicircles.

The native form of Drosophila melanogaster DNA topoisomerase II was purified from Schneider's S3 tissue culture cells and studied with two supercoiled minicircle preparations, mini and mini-CG, 354 bp and 370 bp in length, respectively. Mini-CG contains a d(CG)7 insert which assumes a left-handed Z-DNA conformation in negative supercoiled topoisomers with a negative linking number difference - delta Lk greater than or equal to 2. The interactions of topoisomerase II with topoisomer families of mini and mini-CG were studied by band-shift gel electrophoresis in which the individual topoisomers and their discrete or aggregated protein complexes were resolved. A monoclonal anti-Z-DNA IgG antibody (23B6) bound and aggregated only mini-CG, thereby confirming the presence of Z-DNA. Topoisomerase II bound and relaxed mini-CG more readily than mini. In both cases, there was a preference for more highly negatively supercoiled topoisomers. The topoisomerase II inhibitor VM-26 induced the formation of stable covalent DNA-protein intermediates. In addition, the non-hydrolyzable GTP analogue GTP gamma S inhibited the binding and relaxation activities. Experiments to detect topoisomerase cleavage sites failed to elicit specific loci on either minicircle preparation. We conclude that Drosophila topoisomerase II is able to bind and process small minicircles with lengths as short as 360 bp and negative superhelix densities, - sigma, which can exceed 0.1. Furthermore, the enzyme has a preferential affinity for topoisomers containing Z-DNA segments and relaxes these molecules, presumably by cleavage external to the inserts. Thus, a potentially functional relationship between topoisomerase II, an enzyme regulating the topological state of DNA-chromatin in vivo, and left-handed Z-DNA, a conformation stabilized by negative supercoiling, has been established.     

Computer simulation of DNA supercoiling

We treat supercoiled DNA within a wormlike model with excluded volume. A modified Monte Carlo approach has been used, which allowed computer statistical-mechanical simulations of moderately and highly supercoiled DNA molecules. Even highly supercoiled molecules do not have a regular shape, though with an increase in writhing the chains look more and more like branched interwound helixes. The averaged writhing (Wr) approximately 0.7 delta Lk. The superhelical free energy F is calculated as a function of the linking number. Lk. The calculations have shown that the generally accepted quadratic dependence of F on Lk is valid for a variety of conditions, though it is by no means universal. Significant deviations from the quadratic dependence are expected at high superhelical density under ionic conditions where the effective diameter of DNA is small. The results are compared with the available experimental data.     

The influence of salt on the structure and energetics of supercoiled DNA.

We present a detailed computational study of the influence of salt on the configurations, energies, and dynamics of supercoiled DNA. A potential function that includes both elastic and electrostatic energy components is employed. Specifically, the electrostatic term, with salt-dependent coefficients, is modeled after Stigter's pioneering work on the effective diameter of DNA as a function of salt concentration. Because an effective charge per unit length is used, the electrostatic formulation does not require explicit modeling of phosphates and can be used to study long DNAs at any desired resolution of charge. With explicit consideration of the electrostatic energy, an elastic bending constant corresponding to the nonelectrostatic part of the bending contribution to the persistence length is used. We show, for a series of salt concentrations ranging from 0.005 to 1.0 M sodium, how configurations and energies of supercoiled DNA (1000 and 3000 base pairs) change dramatically with the simulated salt environment. At high salt, the DNA adopts highly compact and bent interwound states, with the bending energy dominating over the other components, and the electrostatic energy playing a minor role in comparison to the bending and twisting terms. At low salt, the DNA supercoils are much more open and loosely interwound, and the electrostatic components are dominant. Over the range of three decades of salt examined, the electrostatic energy changes by a factor of 10. The buckling transition between the circle and figure-8 is highly sensitive to salt concentration: this transition is delayed as salt concentration decreases, with a particularly sharp increase below 0.1 M. For example, for a bending-to-twisting force constant ratio of A/C = 1.5, the linking number difference (delta LK) corresponding to equal energies for the circle and figure-8 increases from 2.1 to 3.25 as salt decreases from 1.0 to 0.005 M. We also present in detail a family of three-lobed supercoiled DNA configurations that are predicted by elasticity theory to be stable at low delta Lk. To our knowledge, such three-dimensional structures have not been previously presented in connection with DNA supercoiling. These branched forms have a higher bending energy than the corresponding interwound configurations at the same delta Lk but, especially at low salt, this bending energy difference is relatively small in comparison with the total energy, which is dominated by the electrostatic contributions. Significantly, the electrostatic energies of the three-lobed and (straight) interwound forms are comparable at each salt environment.(ABSTRACT TRUNCATED AT 400 WORDS)     

Temperature dependence of the gel electrophoretic mobility of superhelical DNA.

We have determined the gel electrophoretic behavior of closed circular plasmid pSM1 DNA (5420 bp) as a function of both temperature and of linking number (Lk). At temperatures below 37 degrees, the electrophoretic mobility first increases, then becomes constant as Lk is decreased below that of the relaxed closed DNA. As the temperature is increased above 37 degrees the electrophoretic mobility first increases as Lk decreases and then varies in a cyclic manner with further decreases in Lk. As the temperature is increased over the range 37 degrees - 65 degrees the cyclic behavior is manifested at progressively smaller decreases in Lk and the amplitude of the cycles increases. We interpret the results in terms of the early melting of superhelical DNA, in which the free energy associated with superhelix formation is progressively transferred to local denaturation. Using a two state approximation, we estimate the free energy change in the first cyclic transition to be 35 Kcal/mole DNA at 37 degrees and to decrease linearly with temperature. The free energy becomes equal to zero at a temperature of 71.6 degrees, which lies within 3 degrees of the melting temperature for the corresponding nicked circular DNA. From the slope of this relationship we estimate the apparent entropy and enthalpy of the first mobility transition to be 6.0 Kcal/mole base pair and 17.3 cal/mole base pair/degree, values consistent with duplex melting.     

DNA topology in hyperthermophilic archaea: reference states and their variation with growth phase, growth temperature, and temperature stresses

In order to address the dynamics of DNA topology in hyperthermophilic archaea, we analysed the topological state of several plasmids recently discovered in Thermococcales and Sulfolobales. All of these plasmids were from relaxed to highly positively super-coiled in vitro, i.e. they exhibited a significant linking excess compared to the negatively supercoiled plasmids from mesophilic organisms (both Archaea and Bacteria). In the two archaeai orders, plasmid linking number (Lk) decreased as growth temperature was lowered from its optimal value, i.e. positively super-coiled plasmids were relaxed whereas relaxed plasmids became negatively supercoiled. Growth temperatures above the optimum correlated with higher positive supercoiling in Sulfolobales (Lk increase) but with relaxation of positive supercoils in Thermococcus sp. GE31. The topological variation of plasmid DNA isolated from cells at different growth phases were found to be species specific in both archaeai orders. In contrast, the direction of topological variation under temperature stress was the same, i.e. a heat shock correlated with an increase in plasmid positive supercoiling, whilst a cold shock induced negative supercoiling. The kinetics of these effects were analysed in Sulfolobales. In both temperature upshift (from 80 to 85C) and downshift (from 80 to 65C), a transient sharp variation of Lk occurred first, and then DNA supercoiling progressively reached levels typical of steady-state growth at the final temperature. These results indicate that DNA topology can change with physiological states and environmental modifications in hyperthermophilic archaea.     

The early melting of closed duplex DNA: analysis by banding in buoyant neutral rubidium trichloroacetate.

Aqueous RbTCA permits the buoyant banding of both native and denatured DNA at room temperature and neutral pH. A unique property of this solvent is the bouyant resolution of closed circular, underwound DNA (I) from the corresponding nicked (II) species. Conditions are reported here in which PM-2 DNA I is physically resolved from native PM-2 DNA II, the buoyant separation being 1.27 mq/ml in 3.3 M RbTCA at 25 degrees C. The separation between nicked and closed DNAs increases with temperature up to 35.5 degrees C, at which PM-2 DNA II cooperatively melts and subsequently pellets. The isothermal buoyant density of a cloed DNA increases linearly as the linking number (Lk) of the closed DNA decreases. The early melting of closed DNA may be monitored with high precision by buoyant banding in RbTCA, it being possible to detect the disruption of as few as 40 base pairs in PM-2 DNA (10,000 base pairs). The constraint that the linking number be conserved in closed DNA requires that a change in duplex winding be accompanied by a compensating change in supercoiling. We estimate the linking number deficiency of PM-2 DNA I to be 0.094 turns per decibase pair. This result permits the estimation of the EtdBr unwinding angle, phi, by comparison with alternative determinations of the linking number deficiency which depend upom the value of phi. The result obtained here is that phi = 27.7 degrees +/- 0.5 degrees and is approximately independent of temperature over the range 15 degrees-35 degrees.     

The Dynamics of Growing Islets and Transmission of Schistosomiasis Japonica in the Yangtze River

We formulate and analyze a system of ordinary differential equations for the transmission of schistosomiasis japonica on the islets in the Yangtze River, China. The impact of growing islets on the spread of schistosomiasis is investigated by the bifurcation analysis. Using the projection technique developed by Hassard, Kazarinoff and Wan, the normal form of the cusp bifurcation of codimension 2 is derived to overcome the technical difficulties in studying the existence, stability, and bifurcation of the multiple endemic equilibria in high-dimensional phase space. We show that the model can also undergo transcritical bifurcations, saddle-node bifurcations, a pitchfork bifurcation, and Hopf bifurcations. The bifurcation diagrams and epidemiological interpretations are given. We conclude that when the islet reaches a critical size, the transmission cycle of the schistosomiasis japonica between wild rats Rattus norvegicus and snails Oncomelania hupensis could be established, which serves as a possible source of schistosomiasis transmission along the Yangtze River.     

The absence of supercoiling in kinetoplast DNA minicircles.

Crithidia fasciculata kinetoplast DNA is a mitochondrial DNA composed of 5000 minicircles and approximately 25 maxicircles, all catenated into a giant network. By comparing the linking number of minicircles released from the network by limited sonication with that of control minicircles, we demonstrate that not only does the elaborate catenation of the network not cause supercoiling, but that there is no minicircle supercoiling at all. The absence of catenation-induced supercoiling is explained by our finding [using electron microscopy (EM) and gel electrophoresis] that network minicircles are joined by only one interlock; single interlocking can be accommodated without helix distortion. EM revealed that propidium diiodide supertwists all the network minicircles and thereby condenses the network into a much smaller size while maintaining its planarity. At high dye concentration the network is condensed to a size comparable to that found in vivo. Nevertheless, network minicircles bind less propidium than free minicircles, indicating that catenation into a network restricts the supercoiling of individual rings. These studies show that the mitochondrion of trypanosomatids may be a unique niche in nature where a covalently-closed circular DNA is not supercoiled. This absence of supercoiling may be a major factor in promoting the formation of the network.