Monte Carlo simulations of locally melted supercoiled DNAs in 20 mM ionic strength

Mesoscopic models of unmelted and locally melted supercoiled DNAs in 20 mM ionic strength are simulated over a range of linking difference from deltal = 0 to -26 turns, or superhelix density from sigma = 0 to -0.062. A domain containing m = 0, 28, or 56 melted basepairs (out of 4349 total) is modeled simply by a region of suitable length with substantially reduced torsion and bending elastic constants. Average structural properties are calculated from the saved configurations, and a reversible work protocol is used to calculate the supercoiling free energy, The cross-writhe between duplex and melted regions (defined herein) is found to be negligibly small. The total writhe, radius of gyration, and ordered elements of the diagonalized inertial tensor are found to be nearly universal functions of the residual linking difference (deltal(r)) associated with the duplex region, independent of m. However, deformability of the tertiary structure, as manifested by the variance of those same properties, is not a universal function of deltal(r)), but depends upon m.delta (SC) varies with deltal(r)) more strongly than deltal(r)) (2)due to the low ionic strength. The twist energy parameter, E (T) obtained from the simulated delta G(SC), deltal(r)), and net twisting strain of the melted region T (D), is found to be independent of m, hence also of the torsion and bending elastic constants of the melted region. However, E(T) increases linearly with -deltalr), which leads to 1). a small overestimation of E (T) for any given deltal(r)) when E(T) is determined from the observed deltal and deltal (r) by the protocol of Bauer and Benham; and 2). a significant enhancement of the apparent slope, -dE(T)/d(T), obtained via the protocol of Bauer and Benham, relative to the actual slope at fixed delta l(r). After taking these two effects into account, the theoretical and experimental values E(T) and -dE(T)/d(T) values agree rather well. For the larger deltal the melted regions are found preferentially in the linker domains between interwound arms, rather than in the apical regions at the ends of interwound arms.     

Monte Carlo simulations of supercoiling free energies for unknotted and trefoil knotted DNAs.

A new Monte Carlo (MC) algorithm is proposed for simulating inextensible circular chains with finite twisting and bending rigidity. This new algorithm samples the relevant Riemann volume elements in a uniform manner, when the constraining potential vanishes. Simulations are performed for filaments comprising 170 subunits, each containing approximately 28 bp, which corresponds to a DNA length of 4770 bp. The bending rigidity is chosen to yield a persistence length, P = 500 A, and the intersubunit potential is taken to be a hard-cylinder potential with diameter d = 50 A. This value of d yields the same second virial coefficient as the electrostatic potential obtained by numerical solution of the Poisson-Boltzmann equation for 150 mM salt. Simulations are performed for unknotted circles and also for trefoil knotted circles using two different values of the torsional rigidity, C = (2.0 and 3.0) x 10(-19) dyne cm2. In the case of unknotted circles, the simulated supercoiling free energy varies practically quadratically with linking difference delta l. The simulated twist energy parameter ET, its slope dET/dT, and the mean reduced writhe <w>/delta l for C = 3 x 10(-19) dyne cm2 all agree well with recent simulations for unknotted circles using the polygon-folding algorithm with identical P, d, and C. The simulated ET vs. delta l data for C = 2.0 x 10(-19) dyne cm2 agree rather well with recent experimental data for p30 delta DNA (4752 bp), for which the torsional rigidity, C = 2.07 x 10(-19) dyne cm2, was independently measured. The experimental data for p30 delta are enormously more likely to have arisen from C = 2.0 x 10(-19) than from C = 3.0 x 10(-19) dyne cm2. Serious problems with the reported experimental assessments of ET for pBR322 and their comparison with simulated data are noted. In the case of a trefoil knotted DNA, the simulated value, (ET)tre, exceeds that of the unknotted DNA, (ET)unk, by approximately equal to 1.40-fold at magnitude of delta l = 1.0, but declines to a plateau about 1.09-fold larger than (ET)unk when magnitude of delta l > or = 15. Although the predicted ratio, (ET)tre/(ET)unk approximately equal to 1.40, agrees fairly well with recent experimental measurements on a 5600-bp DNA, the individual measured ET values, like some of those reported for pBR322, are so large that they cannot be simulated using P = 500 A, d = 50 A, and any previous experimental estimate of C.     

Formation of Z-DNA in negatively supercoiled plasmids is sensitive to small changes in salt concentration within the physiological range.

Negative supercoiling of the plasmid pBR322 with or without an insert of (dG-dC)n induces the formation of Z-DNA as measured by the binding of antibodies specific for Z-DNA. Increasing the concentration of Na+ (or K+) is shown to inhibit the B to Z-DNA conversion. This may be due to the effect of the cation on the B-Z junction. Using the data for B to Z-DNA conversion of the (dG-dC)n inserts, we have estimated the free energy change per base pair as well as the energy of the B-Z junction. In pBR322, a 14-bp segment [CACGGGTGCGCATG] is believed to form Z-DNA at bacterial negative superhelical densities under salt conditions which are similar to those found in vivo.     

Torsional rigidities of weakly strained DNAs

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

Effect of anisotropy of the bending rigidity on the supercoiling free energy of small circular DNAs

In principle, the supercoiling free energy of a small circular DNA will be enhanced by increasing the anisotropy of its bending potential at constant persistence length. The magnitude of this effect is investigated by Monte Carlo simulation using an extension of a previously proposed algorithm. The supercoiling free energy at 298 K is simulated for circular DNAs containing N = 100 bp with torsion constant α = 5.8 × 10^?12 dyne cm. persistence lengths P = 500 Å and 10,000 Å and a range of anisotropies of the bending potential from p = 1.0 to 16.0. The apparent torsion constants, reckoned from these supercoiling free energies by assuming an isotropic bending potential, are found to increase by less than 3% as the input anisotropy increases from 1.0 to 16.0. When P = 500 Å, the apparent torsion constant never rises significantly above the input value over the entire range of input anisotropies. When P = 10,000 Å, the apparent torsion constant rises only about 3% above the input value for anisotropies ρ = 8.0 and 16.0. Evidently, anisotropy of the bending potential cannot account for the fact that the torsion constants reported for small circular DNAs exceed those reported for longlinear DNAs by a factor of 1.6 or more. © 1995 John Wiley & Sons, Inc.     

A biochemical analysis of the interaction of DNA gyrase with the bacteriophage Mu, pSC101 and pBR322 strong gyrase sites: the role of DNA sequence in modulating gyrase supercoiling and biological activity

Replication of bacteriophage Mu DNA, a process requiring efficient synapsis of the prophage ends, takes place within the confines of the Escherichia coli nucleoid. Critical to ensuring rapid synapsis is the function of the SGS, a strong gyrase site, located at the centre of the Mu genome. Replacement of the SGS by the strong gyrase sites from pSC101 or pBR322 fails to support efficient prophage replication. To probe the unique SGS properties we undertook a biochemical analysis of the interaction of DNA gyrase with the Mu SGS, pSC101 and pBR322 sites. In binding and cleavage assays the order of efficacy was pSC101 > Mu SGS > pBR322. However, in supercoiling assays the Mu SGS (cloned into pUC19) exhibited a strong enhancement of gyrase-catalysed supercoiling over pUC19 alone; the pSC101 site showed none and the pBR322 site gave a moderate improvement. Most striking was the Mu SGS-dependent increase in processivity of the gyrase reaction. This highly processive supercoiling coupled with efficient binding may account for the unique biological properties of the SGS. The results emphasize the importance of the DNA substrate as an active component in modulating the gyrase supercoiling reaction, and in determining the biological roles of specialized gyrase sites.     

Effect of ethidium on the torsion constants of linear and supercoiled DNAs.

The torsion elastic constants (alpha) of linear pBR322 (4363 bp) and pUC8 (2717 bp) DNAs and supercoiled pBR322 and pJMSII (4375 bp) DNAs are measured in 0.1 M NaCl as a function of added ethidium/base-pair (EB/BP) ratio by studying the fluorescence polarization anisotropy (FPA) of the intercalated ethidium. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single photon counting detection. Previously developed theory for the emission anisotropy is generalized to incorporate rotations of the transition dipole due to excitation transfer. The excitation transfers are simulated by a Monte Carlo procedure (Genest et al., Biophys. Chem. 1 (1974) 266-278) and the consequent rotations of the transition dipole are superposed on the Brownian rotations. After accounting for excitation transfer, the torsion constants of the linear DNAs are found to be essentially independent of intercalated ethidium up to a binding ratio r = 0.10 dye/bp. Dynamic light scattering measurements on linear pUC8 DNA confirm that the torsion constant is independent of binding ratio up to r = 0.20 dye/bp. If alpha d denotes the torsion constant between ethidium and a base-pair, and alpha 0 that between two base-pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65 to 1.64 with a most probable value of 1.0. The torsion constants of supercoiled DNAs decrease substantially with increasing binding ratio even after accounting for excitation transfer. At the binding ratio r* = 0.064, where the superhelix density vanishes and superhelical strain is completely relaxed, the torsion constant of the supercoiled pBR322 DNA/dye complex lies below that of the corresponding linear DNA/dye complex by about 30%. This contradicts the conventional view according to which linear, nicked circular, and supercoiled DNA/dye complexes with r = r* should coexist with the same concentration of free dye, display the same distribution of bound dye, and exhibit identical secondary structures, twisting and bending rigidities, and FPA dynamics. These and other observations imply the existence of metastable secondary structure in freshly relaxed supercoiled DNAs. A tentative explanation is presented for these and other unexpected observations on supercoiled DNAs.     

DNA topoisomerase I mutants. Increased heterogeneity in linking number and other replicon-dependent changes in DNA supercoiling

Plasmid pBR322 DNA isolated from Salmonella typhimurium supX (topoisomerase I) mutants exhibits a novel supercoiling distribution characterized by extreme heterogeneity in linking number and the presence of highly negatively supercoiled topoisomers. The most negatively supercoiled topoisomers isolated from one supX mutant have more than twice the wild-type level of supercoiling; the distribution as a whole has a median superhelix density about 1.3 times that of wild type. Surprisingly, the supercoiling distribution of plasmid pUC9 DNA isolated from supX mutants differs from that of pBR322. Escherichia coli topoisomerase I mutants have been shown to acquire compensatory mutations that reduce bacterial chromosome supercoiling to below the wild-type level even in the absence of topoisomerase I. We find that such a compensatory mutation in an E. coli topoisomerase I deletion mutant does not reduce pBR322 DNA supercoiling to a level below that of wild type. Thus, the effects of topoisomerase mutations on supercoiling depend on the replicon.     

Stress-induced cruciform formation in a cloned d(CATG)10 sequence.

The synthetic alternating purine-pyrimidine sequence, d(CATG)10.d(CATG)10, has been cloned into a 2.079-kb pBR322-derived plasmid (pLN1) and its conformation studied under torsional stress. The resultant plasmid, pLNc40, is hypersensitive to cleavage by the single strand-specific nucleases, S1 nuclease and Bal31 nuclease, and to modification by the single strand-selective reagent, osmium tetroxide. The S1-hypersensitive site of this plasmid predominates over those previously mapped in pBR322. Site-specific cleavage of pLNc40 with the resolvase T4 endonuclease VII demonstrates that this alternating purine-pyrimidine tract selectively forms a cruciform structure when stably integrated into a negatively supercoiled plasmid. Quantitative measurements of the twist change (-4.3 +/- 0.2) and free energy of formation (16.2 +/- 0.5 kcal/mol) of this cruciform have been made from two-dimensional gel electrophoresis experiments, and correspond well with the predicted values of cruciform formation for this sequence. We conclude that cruciform extrusion versus the B-Z transition is the favoured conformation of this insert under torsional stress.     

Inhibition of DNA replication by berenil in plasmids containing Poly(dA)poly(dT) sequences

The effect of berenil on plasmid DNA replication was studied on pBR322-derived plasmids containing poly(dA)poly(dT) sequences. In comparison to the parental plasmid pBR322, plasmid pKH47 harboring 100 bp of poly(dA)poly(dT) at the PvuII site showed a decrease in plasmid yield in the presence of berenil. This effect was also observed in pVL26, a related plasmid in which the location of the poly(dA)poly(dT) region had been shifted to the EcoRV site in pBR322. [(3)H]Thymidine incorporation experiments indicated that DNA synthesis may be affected in these plasmids in the presence of the drug. Bromodeoxyuridine incorporation experiments coupled to Cs(2)SO(4) equilibrium density gradient centrifugation indicated that the lower plasmid yield was due to an inhibition of DNA replication by berenil. We have also found that berenil induces DNA degradation in plasmids containing the homopolymer. Our studies strongly suggest that the effect of berenil on plasmid replication and DNA stability results from its binding to the poly(dA)poly(dT) region present in these plasmids. Moreover, we have found a correlation between the position of the poly(dA)poly(dT) region and this inhibitory effect. Thus, plasmid pKH47, containing the poly(dA)poly(dT) region most proximal to the origin of pBR322 replication, was most severely affected.     

Plasmid DNA supercoiling and survival in long-term cultures of Escherichia coli: role of NaCl

The relationship between the survival of Escherichia coli during long-term starvation in rich medium and the supercoiling of a reporter plasmid (pBR322) has been studied. In aerated continuously shaken cultures, E. coli lost the ability to form colonies earlier in rich NaCl-free Luria-Bertani medium than in NaCl-containing medium, and the negative supercoiling of plasmid pBR322 declined more rapidly in the absence of NaCl. Addition of NaCl at the 24th hour restored both viability and negative supercoiling in proportion to the concentration of added NaCl. Addition of ofloxacin, a quinolone inhibitor of gyrase, abolished rescue by added NaCl in proportion to the ofloxacin added. This observation raises the possibility that cells had the ability to recover plasmid supercoiling even if nutrients were not available and could survive during long-term starvation in a manner linked, at least in part, to the topological state of DNA and gyrase activity.     

Molecular cloning of cDNA coding for human preprourokinase

A cDNA library was constructed in pBR322 from 18S to 20S mRNA that was extracted from human kidney cells, fractionated on oligo(dT)-cellulose column and sucrose-density gradient, and confirmed for urokinase production in Xenopus laevis oocytes. The Escherichia coli RR1 transformants were hybridized to synthetic oligonucleotide probe prepared according to the known amino acid sequence, Glu 73 to Glu 77 of human urinary urokinase chain B. The entire cloned cDNA covers a 2250-bp region, wherein the 1293-bp sequence codes for preprourokinase consisting of 431 amino acids, with the first 20 residues being a signal peptide. The 5'-untranslated region is at least 80 bp long and the 3'-untranslated region is longer than 850 bp.     

Chemical Probing of the B-Z Transition in Negatively Supercoiled DNA

Abstract

An analysis of the B-to-Z transition as a function of supercoiling for a natural Z-DNA- forming sequence found in plasmid pBR322 is presented at nucleotide resolution. The analysis is based on reactivity to four chemical probes which exhibit hyperreactivity in the presence of Z-DNA: hydroxylamine, osmium tetroxide, diethyl pyrocarbonate and dimethyl sulfate. We find that the initial transition occurs largely within a 14 base pair region which is mostly alternating purines and pyrimidines. With increasing negative supercoiling, Z-DNA extends into flanking regions having less and less alternating character, first in one direction and then in the other. Evidence of B-Z junctions is seen at four sites bracketing these three adjacent regions. One of these Z-forming regions contains the non-alternating sequence CTCCT, suggesting that such sequences can form Z-DNA without great difficulty if they are adjacent to alternating sequences. A plasmid containing three copies of a 61 base pair fragment bearing the entire Z-forming region shows equal reactivity of all three copies at any given superhelical density, implying that they compete equally and independently for the torsional strain energy which promotes the B-Z transition, and are unaffected by adjacent sequences more than 20–30 base pairs away.     

Interaction of chloroquine with linear and supercoiled DNAs. Effect on the torsional dynamics, rigidity, and twist energy parameter

The magnitude and uniformity of the torsion elastic constant (alpha) of linear pBR322 DNA and supercoiled pBR322 DNAs with high-twist (sigma = -0.083) and normal-twist (sigma = -0.48) are measured in 0.1 M NaCl as a function of added chloroquine/base-pair ratio (chl/bp) by studying the fluorescence polarization anisotrophy (FPA) of intercalated ethidium dye. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single-photon counting detection. A general theory is developed for the binding of ligands that unwind superhelical DNAs, and the simultaneous binding of two different intercalators is treated in detail. The equilibrium constant (K) for binding chloroquine to linear pBR322 DNA and the number (r) of bound chloroquines per base pair are determined from the relative amplitude ratio of the slow (normally intercalated) and fast (free) components in the decay of the (probe) ethidium fluorescence intensity as a function of chl/bp. For chloroquine binding to supercoiled pBR322 DNAs, the intrinsic binding constant is assumed to be the same as for the linear DNA, but the twist energy parameter ET (N times the free energy to change the linking number from 0 to 1 in units of kBT) is regarded as adjustable. Using the best-fit ET, the binding ratios r are calculated for each chl/bp ratio. Twist energy parameters are also determined for ethidium binding to these supercoiled DNAs by competitive dialysis. For chloroquine binding, we obtain ET = 360 and 460 respectively for the normal-twist and high-twist supercoiled DNAs. For ethidium binding the corresponding values are ET = 280 +/- 70 and 347 +/- 50. Like other dye-binding values, these are substantially lower than those obtained by ligation methods. In the absence of chloroquine, the torsion constants of all three DNAs are virtually identical, alpha = (5.0 +/- 0.4) x 10(-12) dyn.cm. For linear pBR322 DNA, the magnitude and uniformity of alpha remain unaltered by intercalated chloroquine up to r = 0.19. This finding argues that the FPA is not significantly relaxed by diffusion of any kinks or solitons. If alpha d denotes the torsion constant between a dye and a base pair and alpha 0 that between two base pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65-1.64, with a most probable value of 1.0.(ABSTRACT TRUNCATED AT 400 WORDS)     

DNA gyrase can supercoil DNA circles as small as 174 base pairs.

DNA gyrase introduces negative supercoils into closed-circular DNA using the free energy of ATP hydrolysis. Consideration of steric and thermodynamic aspects of the supercoiling reaction indicates that there should be a lower limit to the size of DNA circle which can be supercoiled by gyrase. We have investigated the supercoiling reaction of circles from 116-427 base pairs (bp) in size and have determined that gyrase can supercoil certain relaxed isomers of circles as small as 174 bp, dependent on the final superhelix density of the supercoiled product. Furthermore, this limiting superhelical density (-0.11) is the same as that determined for the supercoiling of plasmid pBR322. We also find that although circles in the range 116-152 bp cannot be supercoiled, they can nevertheless be relaxed by gyrase when positively supercoiled. These data suggest that the conformational changes associated with the supercoiling reaction can be carried out by gyrase in a circle as small as 116 bp. We discuss these results with respect to the thermodynamics of DNA supercoiling and steric aspects of the gyrase mechanism.