Search

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Purpose: To put radioprobing into context as a relatively new method of determining structural detail in deoxyribonucleic acid (DNA), and to review its use since first proposed in 1997. The key feature of the method is that, by experiment or simulation, a radionuclide such as iodine-125 ((125)I) is placed near the DNA at a known point relative to the DNA base sequence, and the number of resulting strand breaks in each nucleotide is determined. As the intensity of damage declines consistently with distance from the radionuclide, relative distances between the emitter and the nucleotides can be deduced, and hence potentially the topology or structural detail of the DNA. For simulation, appropriate software includes a Molecular Dynamics package, analysis and visualization tools, and a Monte Carlo track structure program., Conclusions: A review of published work and our own recent unpublished studies have shown that radioprobing is sufficiently sensitive and consistent to determine structural detail such as internal folding topology and flexing behavior, and can be applied to DNA or a DNA-protein complex in an approximation to its normal biological environment.

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein-DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50-100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

Purpose: Certain guanine-rich DNA sequences have the capacity to fold into four-stranded structures stabilized by the stacking of square planar arrangements of four hydrogen-bonded guanine bases. However both the overall topology of folding and the more detailed three dimensional structure of these quadruplexes is difficult to determine or predict, and they can be polymorphic, altering radically depending on environmental conditions. Radioprobing experiments, in which Auger electrons emitted during the decay of a (125)I-containing base induce strand cleavage in a distance- and structure-dependent manner, have provided possible means of determining these details. Here we have used a combination of computer simulation methods to study the information obtained by one such experiment, reported in 2004., Method: Models were constructed of three quadruplex topologies considered in the experiment, and one other topology proposed more recently. Molecular Dynamics simulations were used to equilibrate these structures and monitor how they evolved over several nanoseconds in solution. Snapshots from the trajectories were then subjected to Monte Carlo track structure prediction, from which theoretical cleavage patterns have been extracted., Results: The four topologies were found to yield quite different cleavage patterns, which allow the presence of particular conformations in an experiment to be predicted., Conclusion: Radioprobing, which is usable in biologically relevant environments, is sensitive enough to distinguish with some confidence between alternative folding topologies in a DNA structure. Monte Carlo track structure simulation can reinforce or question conclusions drawn from experiment, and Molecular Dynamics used with various restraints provides a practical means of guiding a model towards one that yields cleavage patterns closer to those found experimentally.

This paper provides a brief review of recently published work on biophysical and biological aspects of Auger processes. Three specific questions have been considered. (1) Does charge neutralisation contribute to molecular damage such as DNA strand breaks? (2) How many DNA double strand breaks are produced by a single decay of DNA bound (125)I? (3) What is the correlation between number of gammaH2AX foci and number of double strand breaks (DSB)? The paper also gives preliminary reports on two new calculations: (a) calculation of the spectrum of Auger electrons released during decay of (124)I and (b) the use of Auger electrons in the decay of (125)I as a probing agent of novel DNA structures.

Molecular modelling played a central role in the discovery of the structure of DNA by Watson and Crick. Today, such modelling is done on computers: the more powerful these computers are, the more detailed and extensive can be the study of the dynamics of such biological macromolecules. To fully harness the power of modern massively parallel computers, however, we need to develop and deploy algorithms which can exploit the structure of such hardware. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) is a scalable molecular dynamics code including long-range Coulomb interactions, which has been specifically designed to function efficiently on parallel platforms. Here we describe the implementation of the AMBER98 force field in LAMMPS and its validation for molecular dynamics investigations of DNA structure and flexibility against the benchmark of results obtained with the long-established code AMBER6 (Assisted Model Building with Energy Refinement, version 6). Extended molecular dynamics simulations on the hydrated DNA dodecamer d(CTTTTGCAAAAG)(2), which has previously been the subject of extensive dynamical analysis using AMBER6, show that it is possible to obtain excellent agreement in terms of static, dynamic and thermodynamic parameters between AMBER6 and LAMMPS. In comparison with AMBER6, LAMMPS shows greatly improved scalability in massively parallel environments, opening up the possibility of efficient simulations of order-of-magnitude larger systems and/or for order-of-magnitude greater simulation times.

NMR studies have shown that the minor groove-binding ligand Hoechst 33258 binds to the two T4/A4 tracts within the duplex d(CTTTTCGAAAAG)2 in a highly cooperative manner, such that in titration experiments no intermediate 1:1 complex can be detected. The NMR-derived structures of the free DNA and the 2:1 complex have been obtained, but can shed little light on what the origins of this cooperativity may be. Here we present the results of a series of molecular dynamics simulations on the free DNA, the 1:1 complex, and the 2:1 complex, which have been designed to enable us to calculate thermodynamic parameters associated with the molecular recognition events. The results of the molecular dynamics studies confirm that structural factors alone cannot explain the cooperativity observed, indeed when enthalpic and hydration factors are looked at in isolation, the recognition process is predicted to be slightly anticooperative. However, when changes in configurational entropy are taken into account as well, the overall free energy differences are such that the calculated cooperativity is in good agreement with that observed experimentally. The results indicate the power of molecular dynamics methods to provide reasonable explanations for phenomena that are difficult to explain on the basis of static models alone, and provide a nice example of the concept of "allostery without conformational change".

We model and analyze a short section of a DNA chain with a defect, with the aim of understanding how the frequency, amplitude, and localization of breathing events depend on the strength of the bonds between base pairs, both along the chain and between the chains. Our results show that the presence of a defect in the chain permits the existence of a localized breather mode. The models we analyze are linear and hence solvable, with solvability extending to the statistical mechanics formulation of the problem. Parameter values for the interaction energy of a base with its nearest neighbors are obtained from AMBER. The results indicate good agreement with both the amplitude and the number of base pairs affected by defect-induced breathing motion.

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)(2) and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar-phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving 'lock and key' recognition.

The dynamic behaviour of DNA is of fundamental importance to many cellular processes. One principal characteristic, central to transcription and replication, is the ability of the duplex to "melt". It has recently been shown that dynamic force spectroscopy provides information about the energetics of biomolecular dissociation. We have employed this technique to investigate the unbinding of single dodecanucleotide molecules. To separate the duplex to single-stranded DNA, forces ranging from 17 to 40 pN were required over a range of loading rates. Interpretation of the dependence of melting force on loading rate revealed that the energy barrier to rupture is between 9 and 13 kcal mol-1 in height and situated 0.58 nm from an intermediate structural state. Thermal melting studies show that, prior to dissociation, the oligonucleotide underwent a transition which required between 7 and 11 kcal mol-1 in energy. Through combined dynamic force spectroscopy and thermal melting studies we show the derivation of an energy landscape to dissociate a 12-mer duplex. Until very recently, this type of information was only accessible by computational analysis. Additionally, the force spectroscopy data allow an estimation of the kinetics of duplex formation and melting.

Purpose: To present a new method for radioprobing a DNA triple helix structure by Auger electrons emitted in the decay of 125I using theoretical/computational approaches., Materials and Methods: A Monte Carlo track structure method was used to simulate the damage to a triplex resulting from Auger electrons emitted in the decay of an incorporated 125I atom in plasmid DNA. Comparison of the theoretical frequency distributions of single-strand breaks induced on the Pu and Py strands with the experimental data and a knowledge of the distances from the strand breaks to the iodine provide information on the structures otherwise difficult to obtain with X-ray crystallography., Results: In comparing theoretical frequency distributions of single-strand breaks with the experimental data it is found that the results are very sensitive to the conformation of the triplex model used. It is found that the best fit to the experimental data results from using a hybrid triplex model, in which the base-step geometry is A-like, while the sugar puckers adopt the B-like C2'-endo conformation., Conclusions: The approach and technique presented here represent a valuable new addition to the methods available for DNA structure determination since they provide information on medium-range structure otherwize difficult to obtain in the absence of X-ray crystallography. It is concluded that currently accepted models for triplex structure are not optimal, and a modified structure is proposed that fits the radioprobing results better, while maintaining agreement with the fibre diffraction and NMR data. Although the method has proved to be very useful for scoring alternative trial solutions, further studies combining experimental data from multiple iodine positions with track structure modelling are required for directing structural optimization.

In this study we investigate the possibility of using Auger electrons as a probing agent for the study of structures of nucleic acids. To this end, we present the distribution of breaks produced in strands of a DNA duplex and a triplex-forming oligonucleotide (TFO) carrying Auger emitting radionuclide 125I. The method of calculation includes use of a molecular model of plasmid DNA duplex with bound TFO carrying a labelled 125I at position C5 of a single deoxycytosine residue, a source of Auger spectra, Monte Carlo electron track structure and the ensuing chemistry codes, to simulate the distribution of breaks produced in both strands of a plasmid DNA. Frequencies of fragment length distributions were obtained for the TFO, the purine and the pyrimidine strands. The frequency of breaks in the purine strand showed good correlation with the published experimental results, while that for the pyrimidine strand is lower by a factor of 3. It is concluded that the true structure of triplex DNA may not be purely of B-form.

The structure of the DNA decamer duplex d(GGTAATTACC)(2) has been determined using NMR distance restraints and molecular dynamics simulations of 500 ps to 1 ns in aqueous solution at 300 K. Using both canonical A and canonical B starting structures [root-mean-square deviation (RMSD) 4.6 A; 1 A=10(-10) m], with and without experimental restraints, we show that all four simulations converge to a similar envelope of final conformations with B-like helical parameters (pairwise RMSD 1.27-2.03 A between time-averaged structures). While the two restrained simulations reach a stable trajectory after 300-400 ps, the unrestrained trajectories take longer to equilibrate. We have analysed the dynamic aspects of these structures (sugar pucker, helical twist, roll, propeller twist and groove width) and show that the minor groove width in the AATT core of the duplex fluctuates significantly, sampling both wide and narrow conformations. The structure does not have the highly pre-organized narrow minor groove generally regarded as essential for recognition and binding by small molecules, suggesting that ligand binding carries with it a significant component of 'induced-fit'. Our simulations show that there are significant differences in structure between the TpA step (where p=phosphate) and the ApA and ApT steps, where a large roll into the major groove at the TpA step appears to be an important factor in widening the minor groove at this position.

A pentacyclic acridine, 1H-2,3-dihydroindolizino[7,6,5-kl]acridinium chloride (1), related in structure to tetra- and pentacyclic marine natural products, has previously been shown to induce apoptosis in breast and non-small-cell lung tumor cell lines and shows significant differences in biological potency and antitumor profile from other intercalating agents based on the acridine framework. We report on the molecular recognition of the acridinium salt with DNA, quantified by optical spectroscopic methods, and have compared these results with the clinical agent amsacrine (m-AMSA). The results point to an intercalative association between 1 and G-C-rich sequences of DNA. We have synthesized a hexamer duplex d(ACGCGT)2, presenting two potential 5'-CpG recipient sites, and have investigated in detail by NMR and molecular modeling methods the orientational preferences of 1, particularly with regard to the pyrrolidine ring system. On the basis of the intermolecular nuclear Overhauser effect (NOE) data, four possible intercalation models were considered; no single model produced a significantly better fit than any of the others. The best fit to the experimental data was obtained by considering a dynamic equilibrium between the different intercalated orientations with the drug maximizing pi-overlap with the G-C base pairs at the intercalation site. We found little evidence for any degree of groove specificity imparted by the pyrrolidine ring. If these simulations have biological relevance they suggest that, at most, the agent induces only a transitory hot spot in the DNA which, evidently, is sufficient to be sensed by damage-recognition mechanisms of the cell.

In duplex DNA, groove width and depth are salient structural features that may influence the binding of drugs and proteins. These features are affected by movement of the bases, which for example may enforce groove compression or expansion through a rolling action of the adjacent base-pairs. Moreover, the sugar-phosphate backbone can also undergo limited movement, independently of the bases, which will affect the groove shape. We have examined how the movement of the sugar-phosphate backbone may affect the minor groove width for a fixed base geometry. In agreement with earlier studies, the sugar-phosphate backbone is found to have a certain degree of conformational flexibility in A and B-like helices, and we note a comparable freedom even in the highly curved TATA element of the TATA-binding protein/DNA complex. Phosphate mobility is highly anisotropic in all cases with favoured directions that can significantly change the groove width, independent of any changes in base geometry. We describe how the movement of the sugar-phosphate backbone may affect the accommodation of drugs and proteins in the minor groove, and we present a co-ordinate scheme which emphasises the groove adjustments associated with ligand binding. The observations have implications for the related problem of how cognate molecules are accommodated in the major groove., (Copyright 1999 Academic Press.)

Hoechst 43254 (H43254), a 2,3,4,5-tetrahydropyrimidin-1-ium analogue of the bis-benzimidazole minor groove binding agent Hoechst 33258 (H33258), has been studied by NMR and restrained molecular dynamics in its complex with d(GGTAATTACC)2. We investigate the origin of the enhanced complex stability afforded by the replacement of the N-methylpiperazine ring of H33258 with the tetrahydropyrimidinium ring of H43254, the latter presenting the opportunity for specific minor groove-directed recognition through a pyrimidinium NH. A set of 25 drug-DNA NOEs define the binding site with some precision and are used as part of the structural analysis using restrained molecular dynamics simulations considering explicit solvation and the treatment of electrostatic interactions using the particle mesh Ewald method within AMBER 4.1. Starting with three different initial structures with the drug located at different sites in the groove (pairwise RMSD 4.3-12.6 A) we arrive at three very similar structures (pairwise RMSD 0.80-1.34 A) representing one converged binding site at the centre of the AATT tract. Two of the three structures show the tetrahydropyrimidinium ring to be suitably positioned for an -NH to adenine N3 hydrogen bond suggesting that electrostatic interactions may play an important role in the enhanced affinity as well as imparting additional A-T specificity. The NMR data show that the pyrimidinium NH interaction is dynamic since signal averaging from the two sides of the ring indicate rapid rotations in the bound form.

A triplex-forming oligopyrimidine has been attached at its 5'-end to a photoreactive psoralen derivative and used to target a sequence which forms part of the coding region of the human aromatase gene. The 20 base pair sequence is not a perfect triplex target since it contains three pyrimidine interruptions within the purine-rich strand. Despite this, we have detected triplex-directed photoadduct formation at pH 7.0 between the psoralen-linked oligonucleotide and a 30mer duplex representing the aromatase target. Photoadduct formation was found to be sensitive to pH, temperature, cation concentration and the base composition of the third strand. By varying the base sequence of the target duplex around the psoralen intercalation site, we have characterised the site and mode of psoralen intercalation. The attached psoralen has been found to intercalate at the triplex-duplex junction with a strong preference for one orientation. We have shown that the psoralen will bind at the junction even when there is a preferred TpA step at an adjacent site. We have also compared the binding affinity and photoreactivity of oligodeoxyribonucleotides linked to two different psoralen derivatives and found differences in the rate of crosslinking and the extent of crosslink formation. Finally, we have examined oligodeoxyribonucleotides which are attached to psoralen by polymethylene linkers of different lengths.

Real-time biomolecular interaction analysis (BIA) has been applied to triplex formation between oligodeoxynucleotides. 5'-Biotinylated oligonucleotides were immobilised on the streptavidin-coated surface of a biosensor chip and subsequently hybridised to their complementary strand. Sequence-specific triplex formation was observed when a suitable third-strand oligopyrimidine was injected over the surface-bound duplex. In addition, a single-stranded oligonucleotide immobilised on the chip surface was able to capture a DNA duplex by triplex recognition. The presence of spermine increases the rate of association between the third strand and immobilised duplex, but at elevated spermine concentrations non-specific association is observed. A preliminary kinetic analysis of triplex formation at pH 5.2 by an 11mer third strand containing thymine, cytosine and uracil is reported. Values for the association and dissociation rate constants were determined to be (1.9 +/- 0.2) x 10(3) M-1 s-1 and (8.1 +/- 1.9) x 10(-5) s-1, respectively.

Molecular mechanics has been used to predict the structure of the Y+.R-.R(+)-type DNA triple helix, in which a second polypurine strand binds antiparallel to the homopurine strand of a homopurine/homopyrimidine stretch of duplex DNA. From calculations on the sequence d(C)10.d(G)10.d(G)10, two likely structures emerge. One has the glycosidic torsions of the third strand bases in the anti-conformation and Hoogsteen hydrogen-bonds to the purine strand of the duplex, the other has the third strand purines in the syn orientation and uses a reverse-Hoogsteen hydrogen-bonding pattern. Despite the large structural differences between these two types of triplex, calculations performed in vacuo with a distance-dependent dielectric constant to mimic the shielding effect of solvent show them to be energetically very similar, with the latter (syn) slightly preferred. However, if explicit solvent molecules are included in the calculation, the anti conformation is found to be much preferred. This difference in the results seems to stem from an underestimation of short-range electrostatic interactions in the in vacuo simulations. When TAA or TAT base triples are substituted for the sixth CGG triple in the sequence, it is found that, for the solvated model, the third strand base of the TAA triple prefers the syn orientation while that in the TAT triple retains a preference, though reduced, for the anti conformation.

A molecular dynamics simulation of the DNA triple helix d(TC)5.d(GA)5.d(C+T)5 is described (C+ represents a protonated cytosine residue). The simulation has been performed using the program AMBER 3.1 and includes counterions and explicit solvent under periodic boundary conditions. Both the dynamic and time-averaged behaviour of the system has been analysed. Considerable deviations from the fibre-diffraction model for DNA triple helix structure are observed, including the repuckering of the purine strand sugars that has been identified in some nuclear magnetic resonance (n.m.r.) studies. The simulation suggests that this conformational change may be driven by the possibility of improved interactions between the phosphate groups of this strand and both the solvent and counterions. Several examples of a particular conformational transition are observed, involving correlated changes in the backbone angles alpha and gamma. These transitions provide a possible explanation for some unusual n.m.r. data that have been reported. The structure of the triple helix major groove also suggests an explanation for the observed stabilization of DNA triplexes by polyvalent cations, and their ability to interact with drugs that bind in the minor groove of DNA duplexes.

This study analyses the influence that the nature of the dielectric constant has on the final structures obtained from in vacuo molecular mechanics calculations on a drug-DNA complex and compares these structures with the energy minimised complex including explicit solvent molecules. Minimisations have been performed on a proflavine-decanucleotide structure, where the drug was intercalated at the d(CpG) site of the d(GpApTpApCpGpApTpApC) decamer duplex, using two expressions for the dielectric constant: a distance-independent, epsilon ij = EPS, and a distance-dependent, epsilon ij = EPS*Rij, form and for values of EPS from 1 to 8. Significantly different structures are obtained for the distance-independent and the distance-dependent expressions of the dielectric constant. The use of a distance-independent dielectric constant leads to distorted structures, which are very sensitive to slight changes in the value of EPS. The use of a distance-dependent dielectric constant leads to less distorted and more stable structures. The effects on helical parameters are analysed in detail. The structures obtained for different values of EPS (within the distance-dependent formalism) seem to converge for values of EPS equal to 4 or greater. Based on these results a distance-dependent form of the dielectric with an EPS value of 4 is recommended in order to produce reliable refined nucleic acid structures by molecular mechanics. These conclusions have been supported by molecular-mechanics minimisation of the same structure with the inclusion of explicit water molecules and counter-ions.

We have developed a technique of partially-restrained molecular mechanics enthalpy minimisation which enables the sequence-dependence of the DNA binding of a non-intercalating ligand to be studied for arbitrary sequences of considerable length (greater than = 60 base-pairs). The technique has been applied to analyse the binding of berenil to the minor groove of a 60 base-pair sequence derived from the tyrT promoter; the results are compared with those obtained by DNAse I and hydroxyl radical footprinting on the same sequence. The calculated and experimentally observed patterns of binding are in good agreement. Analysis of the modelling data highlights the importance of DNA flexibility in ligand binding. Further, the electrostatic component of the interaction tends to favour binding to AT-rich regions, whilst the van der Waals interaction energy term favours GC-rich ones. The results also suggest that an important contribution to the observed preference for binding in AT-rich regions arises from lower DNA perturbation energies and is not accompanied by reduced DNA structural perturbations in such sequences. It is therefore concluded that those modes of DNA distortion favourable to binding are probably more flexible in AT-rich regions. The structure of the modelled DNA sequence has also been analysed in terms of helical parameters. For the DNA energy-minimised in the absence of berenil, certain helical parameters show marked sequence-dependence. For example, purine-pyrimidine (R-Y) base pairs show a consistent positive buckle whereas this feature is consistently negative for Y-R pairs. Further, CG steps show lower than average values of slide while GC steps show lower than average values of rise. Similar analysis of the modelling data from the calculations including berenil highlights the importance of DNA flexibility in ligand binding. We observe that the binding of berenil induces characteristic responses in different helical parameters for the base-pairs around the binding site. For example, buckle and tilt tend to become more negative to the 5'-side of the binding site and more positive to the 3'-side, while the base steps at either side of the centre of the site show increased twist and decreased roll.

Molecular mechanics modelling was carried out on the antimicrobial aromatic diamidine, pentamidine, bound to the minor groove of several AT-rich DNA octamers. The pentamidine molecule was found to span four base pairs, with its highly charged amidinium groups forming hydrogen bonds to the O2 of thymine or the N3 of adenine, but not to the backbone phosphate groups. The mean binding energy of the pentamidine-DNA complexes was -52 kcal/mol, of the same order as that of berenil-DNA complexes. There was no significant energy difference between those models which contained GC base pairs and those which did not. The presence of bound pentamidine resulted in some distortion of helix geometry involving helix opening towards the minor groove in most cases and a decrease in the number of residues per turn.

Structural transitions in the tertiary structure of plasmid DNA have been investigated using atomic force microscopy. Changes in superhelical stress were induced by ethidium bromide intercalation, and conformational effects monitored by recording topographic images from DNA complexes of various ethidium bromide : base pair stoichiometry. Significant changes in the tertiary structure of individual DNA molecules were observed with increasing ethidium bromide concentration. The first distinct conformational transition was from a predominantly relaxed structure to one consisting solely of toroidal supercoils. A further increase in ethidium bromide concentration resulted in the formation of regions of plectonemic supercoiling. The ratio of plectonemic : toroidal supercoiling gradually increased until an extremely tightly interwound structure of solely plectonemic supercoiling was finally adopted. The toroidal form of supercoiling observed in this study is unusual as both atomic force microscopy and electron microscopy techniques have previously shown that plectonemic supercoiling is the predominant form adopted by plasmid DNA. [ABSTRACT FROM AUTHOR]