A theoretical study of formation of DNA noncanonical structures under negative superhelical stress

Torque measurements reveal sequence-specific cooperative transitions in supercoiled DNA

B-DNA becomes unstable under superhelical stress and is able to adopt a wide range of alternative conformations including strand-separated DNA and Z-DNA. Localized sequence-dependent structural transitions are important for the regulation of biological processes such as DNA replication and transcription. To directly probe the effect of sequence on structural transitions driven by torque, we have measured the torsional response of a panel of DNA sequences using single molecule assays that employ nanosphere rotational probes to achieve high torque resolution. The responses of Z-forming d(pGpC)(n) sequences match our predictions based on a theoretical treatment of cooperative transitions in helical polymers. "Bubble" templates containing 50-100 bp mismatch regions show cooperative structural transitions similar to B-DNA, although less torque is required to disrupt strand-strand interactions. Our mechanical measurements, including direct characterization of the torsional rigidity of strand-separated DNA, establish a framework for quantitative predictions of the complex torsional response of arbitrary sequences in their biological context.

Inhibitory effect of a short Z-DNA forming sequence on transcription elongation by T7 RNA polymerase

DNA sequences capable of forming unusual secondary structures can be a source of genomic instability. In some cases that instability might be affected by transcription, as recently shown for the Z-DNA forming sequence (CG)₁₄, which causes genomic instability both in mammalian cells and in bacteria, and this effect increases with its transcription. We have investigated the effect of this (CG)₁₄ sequence on transcription with T7 RNA polymerase in vitro. We detected partial transcription blockage within the sequence; the blockage increased with negative supercoiling of the template DNA. This effect was not observed in a control self-complementary sequence of identical length and base composition as the (CG)₁₄ sequence, when the purine–pyrimidine alternation required for Z-DNA formation was disrupted. These findings suggest that the inhibitory effect on T7 transcription results from Z-DNA formation in the (CG)₁₄ sequence rather than from an effect of the sequence composition or from hairpin formation in either the DNA or the RNA product.

Effect of magnesium on cruciform extrusion in supercoiled DNA

Recently, it was reported that Mg2+greatly facilitates cruciform extrusion in the short palindromes of supercoiled DNA, thereby allowing the formation of cruciform structures in vivo. Because of the potential biological importance of this phenomenon, we undertook a broader study of the effect of Mg2+on a cruciform extrusion in supercoiled DNA. The method of two-dimensional gel electrophoresis was used to detect the cruciform extrusion both in the absence and in the presence of these ions. Our results show that Mg2+shifts the cruciform extrusion in the d(CCC(AT)16GGG) palindrome to a higher, rather than to a lower level of supercoiling. In order to study possible sequence-specific properties of the short palindromes for which the unusual cruciform extrusion in the presence Mg2+was reported, we constructed a plasmid with a longer palindromic region. This region bears the same sequences in the hairpin loops and four-arm junction as the short palindrome, except that the short stems of the hairpins are extended. The extension allowed us to overcome the limitation of our experimental approach which cannot be used for very short palindromes. Our results show that Mg2+also shifts the cruciform extrusion in this palindrome to a higher level of supercoiling. These data suggest that cruciform extrusion in the short palindromes at low supercoiling is highly improbable. We performed a thermodynamic analysis of the effect of Mg2+on cruciform extrusion. The treatment accounted for the effect of Mg2+on both free energy of supercoiling and the free energy of cruciform structure per se. Our analysis showed that although the level of supercoiling required for the cruciform extrusion is not reduced by Mg2+, the ions reduce the free energy of the cruciform structure.

Denaturation of supercoiled DNA: a Monte Carlo study

A theoretical investigation of the denaturation characteristics of a supercoiled DNA has been presented employing a Metropolis Monte Carlo algorithm to examine the overall melting profiles of a supercoiled plasmid as the temperature is varied. We show that in contrast to a previously presented algorithm, this much simpler method is sufficient to explain almost all the overall denaturation characteristics and it also correctly calculates the detailed denaturation probabilities of each base pair at various degrees of supercoiling. We also present for the first time a theoretical investigation of the alkaline denaturation of a supercoiled plasmid. Although one can qualitatively reproduce the denaturation profiles using the present Monte Carlo algorithm, the agreement with experiment is not as good as in the case of thermal denaturation. The possible sources of discrepancy between theory and experiment have been discussed.

Two superhelix density-dependent DNA transitions detected by changes in DNA adsorption/desorption behavior

The adsorption behavior of covalently closed circular plasmid DNA at the mercury/water interface was studied by means of AC impedance measurements. The dependence of the differential capacitance (C) of the electrode double layer on the potential (E) was measured in the presence of adsorbed DNA. It was found that the C-E curves of supercoiled DNA at native and highly negative superhelix densities (sigma), relaxed covalently closed circular DNA, and nicked DNA differed from each other. A detailed study of topoisomer distributions ranging from -sigma of 0 to 0.11 revealed two supercoiling-dependent transitions, at about -sigma = 0.04 (transition TI) and 0.07 (transition TII). Transition TI was detected by measuring the height of the adsorption/desorption peak 1 (at about -1.2 V against the saturated calomel electrode) and the decrease of capacitance (DeltaC) at -0.35 V. This transition may be due to a sudden change in the ability of the DNA to respond to the alternating voltage, probably caused by changes in the DNA tertiary and/or secondary structure. Transition TII was detected by measuring peak 3* (at about -1.3 V), which was absent in topoisomers with -sigma less than 0.05. This transition is due to changes in the DNA adsorption/desorption behavior related to increased accessibility of bases at elevated negative superhelix density. Opening of the duplex at highly negative superhelix density was also detected by the single-strand selective probe of DNA structure, osmium tetroxide, 2, 2'-bipyridine. Our results suggest that electrochemical techniques provide sensitive experimental analysis of changes in DNA structure.

The effect of ionic conditions on the conformations of supercoiled DNA. I. Sedimentation analysis

We studied the conformations of supercoiled DNA as a function of superhelicity and ionic conditions by determining its sedimentation coefficient both experimentally and by calculation. To cancel out unknown parameters from both calculations and experiments, we determined the ratio of the sedimentation coefficient, s, to that of open circular DNA, s(oc). Calculations of the sedimentation coefficient were based on direct solution of the Burgers-Oseen problem for an equilibrium set of DNA conformations generated for each condition by the Metropolis Monte Carlo procedure. There were no adjustable parameters in the Monte Carlo simulations because all three parameters of the DNA model used, bending and torsional elasticity of DNA and DNA effective diameter specifying electrostatic interactions, were known from independent data. The good agreement between measured and calculated values of s/s(oc) allowed us to interpret the sedimentation results in terms of DNA conformations, with particular emphasis on the marked effect of ionic conditions. As NaCl concentration decreases, s/s(oc) increases because the superhelix becomes less regular and more compact. In the presence of just 10 mM MgCl(2), supercoiled DNA adopts essentially the same set of conformations as in moderate to high concentrations of NaCl. Our simulations showed that s is a strong function of the superhelix branching frequency. At near physiological ionic conditions, there are about four branches in the 7 kb DNA molecule used in this work. We found no indication of superhelix collapse in any ionic conditions even remotely approaching physiological ones. For all ionic conditions studied, we conclude that the electrostatic interaction of DNA segments specified by the DNA effective diameter is the primary determinant of supercoiled DNA conformations.

Flow of structural information between four DNA conformational levels

Closed-circular supercoiled DNA molecules have been shown to form a cholesteric assembly within bacteria as well as in vitro under physiological DNA and salt concentrations. Circular dichroism and X-ray scattering studies indicate that the macroscopic structural properties of the chiral mesophase are directly and uniquely dictated by the supercoiling parameters of the constituent molecules. Specifically, we find that the pitch of the DNA cholesteric phase derived from supercoiled DNA is determined by the superhelical density, which, in turn, is modulated by secondary conformational changes. A direct interrelationship among four DNA structural levels, namely, DNA sequence, secondary structural transitions, the tertiary superhelical conformation, and the quaternary, supramolecular organization is accordingly pointed out. Since secondary conformational changes are both sequence and environment dependent, alterations of cellular conditions may effectively modulate the properties of the packed DNA organization, through their effects on secondary structural transitions and hence on the superhelical parameters. On the basis of these results we suggest that liquid crystallinity represents an effectively regulated packaging mode of plectonemic, nucleosome-free DNA molecules in living systems.

Effects of base substituents on the hydration of B- and Z-DNA: correlations to the B- to Z-DNA transition

We present a study of how substituent groups of naturally occurring and modified nucleotide bases affect the degree of hydration of right-handed B-DNA and left-handed Z-DNA. A comparison of poly(dG-dC) and poly(dG-dm5C) titrations with the lipotropic salts of the Hofmeister series infers that the methyl stabilization of cytosines as Z-DNA is primarily a hydrophobic effect. The hydration free energies of various alternating pyrimidine-purine sequences in the two DNA conformations were calculated as solvent free energies from solvent accessible surfaces. Our analysis focused on the N2 amino group of purine bases that sits in the minor groove of the double helix. Removing this amino group from guanine to form inosine (I) destabilizes Z-DNA, while adding this group to adenines to form 2-aminoadenine (A') stabilizes Z-DNA. These predictions were tested by comparing the salt concentrations required to crystallize hexanucleotide sequences that incorporate d(CG), d(CI), d(TA) and d(TA') base pairs as Z-DNA. Combining the current results with our previous analysis of major groove substituents, we derived a thermodynamic cycle that relates the systematic addition, deletion, or substitution of each base substituent to the B- to Z-DNA transition free energy.

Statistical mechanical approach for predicting the transition to non-B DNA structures in supercoiled DNA

Supercoiling causes global twist of DNA structure and the supercoiled state has wide influence on conformational transition. A statistical mechanical approach was made for prediction of the transition probability to non-B DNA structures under torsional stress. A conditional partition function was defined as the sum over all possible states of the DNA sequence with basepair 1 and basepair n being in B-form helix and a recurrence formula was developed which expressed the partition function for basepair n with those for less number of pairs. This new definition permits a quick enumeration of every configuration of secondary structures. Energetic parameters of all conformations concerned, involving B-form, interior loop, cruciform and Z-form, were included in the equation. The probability of transition to each non-B conformation could be derived from these conditional partition functions. For treatment of effects of superhelicity, supercoiling energy was considered, and a twist of each conformation was determined to minimize the supercoiling energy. As the twist itself affects the transition probability, the whole scheme of equations was solved by renormalization technique. The present method permits a simultaneous treatment of several types of conformations under a common torsional stress. A set of energetic parameters of DNA secondary structures has been chosen for calculation. Some DNA sequences were submitted to the calculation, and all the sequences that we submitted gave stable convergence. Some of them have been investigated the critical supercoil density for the transition to non-B DNA structures. Even though the reliability of the set of parameters was not enough, the prediction of secondary structure transition showed good agreement with reported observation. Hence, the present algorithm can estimate the probability of local conformational change of DNA under a given supercoil density, and also be employed to predict some specific sequences in which conformational change is sensitive to superhelicity.

Parallel-stranded DNA under topological stress: rearrangement of (dA)15.(dT)15 to a d(A.A.T)n triplex

DNA oligonucleotides with appropriate sequences can form a stable duplex in which the two strands are paired in a parallel orientation instead of as the conventional antiparallel double helix of B-DNA. In parallel-stranded DNA (ps-DNA) base pairing is noncanonical with the glycosidic bonds in a trans orientation. The two grooves are equivalent. We have synthesized DNA duplexes consisting of a central parallel-stranded (dA)15.(dT)15 tract flanked by normal antiparallel regions, and ligated them into the pUC18 plasmid. The effect of negative supercoiling on the covalently closed circular molecules was studied by two-dimensional agarose gel electrophoresis and by chemical modification with OsO4-pyridine (Os,py) and diethylpyrocarbonate (DEPC). The following results were obtained: (i) The ps insert, and by inference ps-DNA in general, adopts a right handed helical form. (ii) Upon increasing the negative superhelix density (-sigma) to greater than 0.03 the 15 bp ps insert undergoes a major transition leading to a relaxation corresponding to a reduction in twist of approximately 2.5 helical turns. The transition free surgery is approximately kcal/mol. (iii) The chemical modification pattern of the resulting structure suggests that the purine strand folds back and associates with the pyrimidine strand, forming a novel intramolecular triplex structure consisting of d(A.A.T) base triplets. A model for the triplex conformation is proposed and its thermodynamic properties are analyzed by statistical mechanics.

H-DNA and Z-DNA in the mouse c-Ki-ras promoter

The mouse c-Ki-ras protooncogene promoter contains a homopurine-homopyrimidine domain that exhibits S1 nuclease sensitivity in vitro. We have studied the structure of this DNA region in a supercoiled state using a number of chemical probes for non-B DNA conformations including diethyl pyrocarbonate, osmium tetroxide, chloroacetaldehyde, and dimethyl sulfate. The results demonstrate that two types of unusual DNA structures formed under different environmental conditions. A 27-bp homopurine-homopyrimidine mirror repeat adopts a triple-helical H-DNA conformation under mildly acidic conditions. This H-DNA seems to account for the S1 hypersensitivity of the promoter in vitro, since the observed pattern of S1 hypersensitivity at a single base level fits well with the H-DNA formation. Under conditions of neutral pH we have detected Z-DNA created by a (CG)5-stretch, located adjacent to the homopurine-homopyrimidine mirror repeat. The ability of the promoter DNA segment to form non-B structures has implications for models of gene regulation.

Parallel-stranded DNA under topological stress: rearrangement of (dA)15.(dT)15 to a d(A.A.T)n triplex

DNA oligonucleotides with appropriate sequences can form a stable duplex in which the two strands are paired in a parallel orientation instead of as the conventional antiparallel double helix of B-DNA. In parallel-stranded DNA (ps-DNA) base pairing is noncanonical with the glycosidic bonds in a trans orientation. The two grooves are equivalent. We have synthesized DNA duplexes consisting of a central parallel-stranded (dA)15.(dT)15 tract flanked by normal antiparallel regions, and ligated them into the pUC18 plasmid. The effect of negative supercoiling on the covalently closed circular molecules was studied by two-dimensional agarose gel electrophoresis and by chemical modification with OsO4-pyridine (Os,py) and diethylpyrocarbonate (DEPC). The following results were obtained: (i) The ps insert, and by inference ps-DNA in general, adopts a right handed helical form. (ii) Upon increasing the negative superhelix density (-sigma) to greater than 0.03 the 15 bp ps insert undergoes a major transition leading to a relaxation corresponding to a reduction in twist of approximately 2.5 helical turns. The transition free surgery is approximately kcal/mol. (iii) The chemical modification pattern of the resulting structure suggests that the purine strand folds back and associates with the pyrimidine strand, forming a novel intramolecular triplex structure consisting of d(A.A.T) base triplets. A model for the triplex conformation is proposed and its thermodynamic properties are analyzed by statistical mechanics.

Paranemic structures of DNA and their role in DNA unwinding

A DNA structure is defined as paranemic if the participating strands can be separated without mutual rotation of the opposite strands. The experimental methods employed to detect paranemic, unwound, DNA regions is described, including probing by single-strand specific nucleases (SNN), conformation-specific chemical probes, topoisomer analysis, NMR, and other physical methods. The available evidence for the following paranemic structures is surveyed: single-stranded DNA, slippage structures, cruciforms, alternating B-Z regions, triplexes (H-DNA), paranemic duplexes and RNA, protein-stabilized paranemic DNA. The problem of DNA unwinding during gene copying processes is analyzed; the possibility that extended paranemic DNA regions are transiently formed during replication, transcription, and recombination is considered, and the evidence supporting the participation of paranemic DNA forms in genes committed to or undergoing copying processes is summarized.