Left-handed Z form in superhelical DNA: a theoretical study

This is a comprehensive statistical mechanical treatment of the Z form formation in purinepyrimidine stretches of different length inserted into superhelical DNA. The B-Z transition for short inserts is shown to follow the "all-or-none" principle. Over some critical value of the insert length n, the B-Z transition in the insert proceeds in two stages. The flipping of m base pairs into the Z form is followed by a gradual growth of the Z-form stretch until it occupies the whole insert. By fitting the theoretical transition curves to experimental ones the fundamental thermodynamic parameters of the B-Z transition have been determined: the B-Z junction energy Fj = 4-5kcal.mol-1 and the free energy change delta FBZ = 0.5-7.0 kcal.mol-1 under standard salt conditions. Calculations show that the B-Z transition in short purinepyrimidine inserts may be seriously affected by cruciform formation in the carrier DNA.     

Preliminary spectroscopic characterization of a synthetic DNA oligomer containing a B-Z junction at high salt.

It is well known that the local conformation of a segment of DNA is dependent upon both the sequence of the segment and the conditions under which the DNA is prepared. In extreme cases, the DNA may contain regions of both right and left-handed conformations, mandating the existence of a conformational junction between the two. These B-Z junctions have been observed in plasmids but, to date, no model systems have been characterized to determine the molecular nature of these junctions. Preliminary CD, UV, and NMR studies on such a model are presented here. A 16 base pair oligonucleotide, containing a potential B-Z junction, has been synthesized and characterized by the above techniques. The results suggest that this molecule contains both right and left-handed conformations under condition of high salt, and thus a B-Z junction.     

Effects of 5 cytosine methylation on the B-Z transition in DNA restriction fragments and recombinant plasmids

Alternating (dC-dG)n regions in DNA restriction fragments and recombinant plasmids were methylated at the 5 position of the cytosine residues by the HhaI methylase. Methylation lowers the concentration of NaCl or MgCl2 necessary to cause the B-Z conformational transition in these sequences. Ionic strengths higher than physiological conditions are required to form the Z conformation when the methylated (dC-dG)n tract is contiguous with regions that do not form Z structures, in contrast to the results with the DNA polymer poly(m5dC-dG) . poly(m5dC-dG). In supercoiled plasmids containing (dC-dG)n sequences, methylation reduces the number of negative supercoils necessary to stabilize the Z conformation. Calculations of the observed free energy contributions of the B-Z junction and cytosine methylation suggest that two junctions offset the favorable effect of methylation on the Z conformation in (dC-dG)n sequences (about 29 base-pairs in length). Studies with individual methylated topoisomers demonstrate that increasing Na+ concentration up to approximately 0.2 M inhibits the formation of the Z conformation in the (m5dC-dG)n region of supercoiled plasmids. The results suggest that methylation may serve as a triggering mechanism for Z DNA formation in supercoiled DNAs.     

Left-handed Z Form in Superhelical DNA: A Theoretical Study

Abstract

This is a comprehensive statistical mechanical treatment of the Z form formation in purine- pyrimidine stretches of different length inserted into superhelical DNA. The B-Z transition for short inserts is shown to follow the “all-or-none” principle. Over some critical value of the insert length n, the B-Z transition in the insert proceeds in two stages. The flipping of m base pairs into the Z form is followed by a gradual growth of the Z-form stretch until it occupies the whole insert. By fitting the theoretical transition curves to experimental ones the fundamental thermodynamic parameters of the B-Z transition have been determined: the B-Z junction energy Fj = 4–5kcal?mol^?1 and the free energy change ΔF_B-Z = 0.5–0.7 kcal?mol^?1 under standard salt conditions. Calculations show that the B-Z transition in short purine-pyrimidine inserts may be seriously affected by cruciform formation in the carrier DNA.     

Temperature-dependent CD and NMR studies on a synthetic oligonucleotide containing a B-Z junction at high salt

It is now accepted that two or more conformations may exist within the same DNA molecule, thereby generating conformational junctions. The presence of B-Z junctions between right- and left-handed DNA conformations has been detected in plasmids by a number of techniques. Preliminary characterization of the first example of a B-Z junction is a short DNA oligonucleotide has recently been reported [Sheardy, R. D. (1988) Nucleic Acids Res. 16, 1153-1167]. We report additional CD and NMR data that support the existence of the junction in this model oligomer. These studies indicate that only three base pairs are involved in the junction and only one of these is dramatically distorted. Furthermore, the NMR saturation-transfer experiments suggest the junction's internal motion is temperature dependent.     

Allosteric conversion of Z DNA to an intercalated right-handed conformation by daunomycin

Absorbance and fluorescence methods were used to measure the binding of the anticancer drug daunomycin to poly (dGdC) under ionic conditions that initially favor the left-handed Z conformation of the polymer. Drug binding was cooperative under these conditions and may be fully accounted for by an allosteric model in which the drug binds preferentially (but not exclusively) to the right-handed B conformation and shifts the polymer from the Z to an intercalated right-handed conformation. Quantitative analysis of binding isotherms in terms of the allosteric model allowed for estimation of the equilibrium constants for the conversion of a base pair at a B-Z interface from the Z to the B conformation and for the formation of a base pair in the B conformation within a stretch of helix in the Z conformation. The free energy of the Z to B conversion of a base pair was calculated from this data and ranges from +0.03 to +0.3 kcal/mol over the NaCl range of 2.4-3.5 M. The free energy for the formation of a B-Z junction was nearly constant at +4.0 kcal/mol over the same range of NaCl concentrations. The salt dependence of the free energy of the Z to B transition indicates preferential Na+ binding to the Z form and that there is a net release of Na+ upon conversion of a base pair from the Z to the B conformation. The energetically unfavorable Z to B transition was found by this analysis to be driven by coupling to the energetically favorable interaction of daunomycin with B form DNA. In 3.5 M NaCl, for example, the free energy change for the overall reaction (Z DNA base pairs) + (daunomycin) in equilibrium with (right-handed complex) is -7.0 kcal/mol, nearly all of which is contributed by the binding of drug to B DNA. Analysis using the allosteric model also shows that the number of base pairs converted from the Z to the B conformation per bound drug molecule is salt dependent and provides evidence that drug molecules partition into regions of the polymer in the right-handed conformation.     

B-Z cooperativity and kinetics of poly(dG-m5dC) are controlled by an unfavorable B-Z interface energy

Thermodynamic and kinetic properties of the B-Z transition of poly(dG-m5dC) were investigated using polynucleotide samples ranging in length from 11000 to 300 base pairs. Van't Hoff enthalpy values increase with increasing polymer length for the B-Z transition in 0.35 mM MgCl2, 50 mM NaCl, 5 mM TRIS, pH 8. Rates of the B to Z transition increase with increasing polymer length for a jump of 0 to 3 mM MgCl2 in 50 mM NaCl, 5 mM TRIS, pH 8. The activation energy of the B to Z transition equals 7.9 +/- 0.3 kcal/mol and is length independent. Thermodynamic and kinetic data were fit to a model that simulates distribution of B- and Z-form tracts at the midpoint of B-Z equilibrium as a function of polymer length. A cooperative length of 1000 +/- 200 base pairs is estimated for the B-Z transition. A direct relationship between rates of the B to Z transition and the square of the van't Hoff enthalpy values of the B-Z transition reflects a dependence of kinetics and cooperativity upon the energy of the nucleation event. Faster B to Z transition rates with increasing polymer length can be explained by a mechanism rate limited by nucleation within the polymer instead of the ends.     

Applying the stochastic difference equation to DNA conformational transitions: a study of B-Z and B-A DNA transitions

Despite the existence of numerous models to account for the B-Z DNA transition, experimenters have not yet arrived at a conclusive answer to the structural and dynamical features of the B-Z transition. By applying the stochastic difference equation to simulate the B-Z DNA transition, we have shown that the stretched intermediate model of the B-Z transition is more probable than other B-Z transition models such as the Harvey model. This is accomplished by comparing potential energy profiles of various B-Z DNA transition models and calculating relative probabilities based on the stochastic difference equation with respect to length (SDEL) formalism. The results garnered in this article allow for new approaches in determining the structural transition of B-DNA to Z-DNA experimentally. We have also simulated the B-A DNA transition using the stochastic difference equation. Unlike the B-Z DNA transition, the mechanism for the B-A DNA transition is well established. The variation in the pseudorotation angle during the transition is in good agreement with experimental results. Qualitative features of the simulated B-A transition also agree well with experimental data. The SDEL approach is thus a suitable numerical technique to compute long-time molecular dynamics trajectory for DNA molecules.     

Enhanced reactivity of a B-Z junction for cleavage by the restriction enzyme MboI

We have been investigating the structure, dynamics, and ligand-binding properties of the interface that exists between a right-handed conformation and a left-handed conformation (i.e., a B-Z junction) in synthetic DNA oligomers. Since exo- and endonuclease activity is known to be sensitive to the conformation of the template DNA, we have designed and synthesized a DNA oligonucleotide of 20 base pairs (designated as BZ-III) with an MboI recognition site (GATC) at the location of a potential B-Z junction. The activity of the MboI enzyme toward this molecule and DNA oligomers that contain multiple MboI sites located at B-Z junctions was monitored in the absence and presence of the Z-conformation-inducing reagent cobalt hexaammine. In all cases, the activity of the enzyme was enhanced in the presence of cobalt hexaammine. The activity of MboI toward BZ-III, in the presence and absence of cobalt hexaammine, was also examined when the DNA oligomer is also in the presence of the DNA binding drugs actinomycin D, ametantrone, or ethidium bromide. In all cases, the activity of the enzyme was inhibited in the presence of drug. The results suggest that B-Z junctions are structurally unique and that this uniqueness may alter nuclease activity at sites in or near the junction.     

Activation of DNA cleavage by dynemicin A in a B-Z conformational junction.

We report here that the DNA strand scission by dynemicin A is not only sequence-specific but also conformation-specific. The salt-induced B----Z conformational transition dramatically enhanced the cleavage by dynemicin A in a B-Z junction region. By contrast, the bleomycin-Fe(II) complex, the elsamicin A-Fe(II) complex, and esperamicin A1 did not induce any preferential DNA cutting in such a DNA structure. The characteristic hyperreactivity of dynemicin A is observed in (dC-dG)8- and (dC-dG)12-inserted DNAs, but not in (dC-dG)5-inserted DNA. These results suggest value in the use of dynemicin A as proof of the existence of a B-Z junction in vivo and also may aid in understanding the structure of B-Z junctions.     

A model of the strain-induced B-Z transition

The B-to-Z transition in supercoiled circular DNA is modeled as a strain-induced nonlinear excitation process. Using a model, in which DNA is regarded as a chain of units with a bistable energy function along the twisting coordinate together with a harmonic inter-unit interaction, we show that a Z region and the accompanying two B-Z junctions of finite width appear naturally as a solution of nonlinear equations, when the strain exceeds a critical value. We examine the B-Z transition behaviour as a function of twist under various situations. We also analyse available experimental results on B-Z transition in supercoiled plasmid with G-C insertions by this mechanistic model in order to estimate the magnitude of model parameters. The energy barrier of the B-Z transition is estimated to be of the order of 1 kcal/mole per base pair. The analysis shows that if the length of the insertion is less than a certain value, the entire insertion converts to Z form at a transition point, but if the insertion is much longer, the B-Z transition exhibits a different behavior, in which part of the insertion flips to Z form and the Z region expands linearly upon changing linking number.     

A Model of the Strain-induced B-Z Transition

Abstract

The B-to-Z transition in supercoiled circular DNA is modeled as a strain-induced nonlinear excitation process. Using a model, in which DNA is regarded as a chain of units with a bistable energy function along the twisting coordinate together with a harmonic inter-unit interaction, we show that a Z region and the accompanying two B-Z junctions of finite width appear naturally as a solution of nonlinear equations, when the strain exceeds a critical value. We examine the B-Z transition behaviour as a function of twist under various situations. We also analyse available experimental results on B-Z transition in supercoiled plasmid with G-C insertions by this mechanistic model in order to estimate the magnitude of model parameters. The energy barrier of the B-Z transition is estimated to be of the order of 1 kcal/mole per base pair. The analysis shows that if the length of the insertion is less than a certain value, the entire insertion converts to Z form at a transition point, but if the insertion is much longer, the B-Z transition exhibits a different behavior, in which part of the insertion flips to Z form and the Z region expands linearly upon changing linking number.     

The Z-Z junction: the boundary between two out-of-phase Z-DNA regions

The boundary between two segments of Z-DNA that differ in the phase of their syn-anti alternation about the glycosidic bond is termed a Z-Z junction. Using chemical probes and two-dimensional gel electrophoresis, we examined a Z-Z junction consisting of the sequence d[(CG)8C(CG)8] inserted into a plasmid and used energy minimization techniques to devise a three-dimensional model that is consistent with the available data. We show that both alternating CG segments undergo the B-Z transition together to form a Z-Z junction. The junction is very compact, displaying a distinctive reactivity signature at the two base pairs at the junction. In particular, the 5' cytosine of the CC dinucleotide at the junction is hyperreactive toward hydroxylamine, and the two guanines of the GG dinucleotide on the complementary strand are less reactive toward diethyl pyrocarbonate than are the surrounding Z-DNA guanines. Statistical mechanical treatment of the 2-D gel data yields a delta G for forming the Z-Z junction equal to 3.5 kcal, significantly less than the cost of a B-Z junction and approximately equal to the cost of a base out of alternation (i.e., a Z-DNA pyrimidine in the syn conformation). The computer-generated model shows little distortion of the Z helix outside of the central two base pairs, and the energy of the structure and the steric accessibility of the reactive groups are consistent with the data.     

B-Z Cooperativity and Kinetics of Poly(dG-m5dC) are Controlled by an Unfavorable B-Z Interface Energy

Abstract

Thermodynamic and kinetic properties of the B-Z transition of poly(dG-m⁵dC) were investigated using polynucleotide samples ranging in length from 11000 to 300 base pairs. Van't Hoff enthalpy values increase with increasing polymer length for the B-Z transition in 0.35 mM MgCl₂, 50 mM NaCl, 5 mM TRIS, pH 8. Rates of the B to Z transition increase with increasing polymer length for a jump of 0 to 3 mM MgCl, in 50 mM Nad, 5 mM TRIS, pH 8. The activation energy of the B to Z transition equals 7.9 ± 0.3 kcal/mol and is length independent Thermodynamic and kinetic data were fit to a model that simulates distribution of B- and Z-form tracts at the midpoint of B-Z equilibrium as a function of polymer length. A cooperative length of 1000 ± 200 base pairs is estimated for the B-Z transition. A direct relationship between rates of the B to Z transition and the square of the van't Hoff enthalpy values of the B-Z transition reflects a dependence of kinetics and cooperativity upon the energy of the nucleation event Faster B to Z transition rates with increasing polymer length can be explained by a mechanism rate limited by nucleation within the polymer instead of the ends.     

Unpaired bases in the B--Z junction of the supercoiled plasmid

The restriction analysis has been used to establish that O-beta-diethylaminoethylhydroxylamine (OHA) produces modification of unpaired cytidines in the polylinker region adjacent to the Z-insert (dG-dC)10. (dG-dC)10 in the negatively supercoiled plasmid pGC20. The length of the transition region between B- and Z-portions of DNA is not less than 36 bps. The reaction of OHA with the unpaired cytidines in the B-Z junction is a fixing one and produces no secondary despiralling of the neighboring regions. The reaction with DNA proceeds much slower than the one with monomers and single-strand polynucleotides. The structural nonuniformity has been observed, which is manifested in the alternating B and "non-B" form DNA in the B-Z junction. It is suggested that these junctions may contain nucleotide sequences which are stable to violation of the B structure during the change in superhelical density of DNA.     

Base protonation facilitates B-Z interconversions of poly(dG-dC) X poly(dG-dC)

Comparative studies on the salt titration and the related kinetics for poly(dG-dC) X poly(dG-dC) in pH 7.0 and 3.8 solutions clearly suggest that base protonation facilitates the kinetics of B-Z interconversion although the midpoint for such a transition in acidic solution (2.0-2.1 M NaCl) is only slightly lower than that of neutral pH. The rates for the salt-induced B to Z and the reverse actinomycin D induced Z to B transitions in pH 3.8 solutions are at least 1 order of magnitude faster than the corresponding pH 7.0 counterparts. The lowering of the B-Z transition barrier is most likely the consequence of duplex destabilization due to protonation as indicated by a striking decrease (approximately 40 degrees C) in melting temperature upon H+ binding in low salt. The thermal denaturation curve for poly(dG-dC) X poly(dG-dC) in a pH 3.8, 2.6 M NaCl solution indicates an extremely cooperative melting at 60.5 degrees C for protonated Z DNA, which is immediately followed by aggregate formation and subsequent hydrolysis to nucleotides at higher temperatures. The corresponding protonated B-form poly(dG-dC) X poly(dG-dC) in 1 M NaCl solution exhibits a melting temperature about 15 degrees C higher, suggesting further duplex destabilization upon Z formation.     

Anomalous gel migration of DNA oligomers containing multiple conformational junctions

We have previously shown that a short 16 base pair DNA oligomer can accommodate a B-Z conformational junction [Sheardy, R. D., & Winkle, S. A. (1989) Biochemistry 28, 720-725]. Results from 1H NMR studies indicated that only three base pairs were involved in the junction and that one of these base pairs was highly distorted. Being interested in the nature of this distortion, we constructed DNA oligomers which have the potential to contain multiple B-Z junctions for polyacrylamide electrophoretic studies. We report that the mobilities displayed by these molecules through acrylamide gels in the absence and presence of cobalt suggest that these molecules run shorter than they actually are. This anomalous migration may be due to structural/dynamic properties of the DNA helix manifested by the periodic distortions of the potential B-Z junctions.     

The energetics of the B-Z transition in DNA

The paper deals with the energetics of the transition to left-handed Z form in DNA with an arbitrary base sequence. There is a brief outline of the statistical-mechanical model of the B-Z transition allowing for three possible states of each base pair. The parameters of the model can be determined by comparing the theory with experimental data for the B-Z transition in inserts with given sequences in circular DNA. The model contains six energy parameters, most of which have been determined before. In order to find the remaining parameters of the model and test its adequacy, a number of oligonucleotide sequences were synthesized and inserted into the pUC 19 plasmid. Two-dimensional gel electrophoresis was used to determine the superhelical density at which the inserts adopt the Z form. A statistical-mechanical treatment of these data yielded a complete set of six energy parameters for the B-Z transition. The theoretical assumption that the free energy of Z-form pairs does not depend on the type of adjacent pairs proved to be in agreement with the experimental data.     

The Energetics of the B-Z Transition in DNA

Abstract

The paper deals with the energetics of the transition to left-handed Z form in DNA with an arbitrary base sequence. There is a brief outline of the statistical-mechanical model of the B-Z transition allowing for three possible states of each base pair. The parameters of the model can be determined by comparing the theory with experimental data for the B-Z transition in inserts with given sequences in circular DNA The model contains six energy parameters, most of which have been determined before. In order to find the remaining parameters of the model and test its adequacy, a number of oligonucleotide sequences were synthesized and inserted into the pUC 19 plasmid. Two-dimensional gel electrophoresis was used to determine the superhelical density at which the inserts adopt the Z form. A statistical-mechanical treatment of these data yielded a complete set of six energy parameters for the B-Z transition. The theoretical assumption that the free energy of Z-form pairs does not depend on the type of adjacent pairs proved to be in agreement with the experimental data.     

B-A and B-Z transitions in deoxyoligoduplexes containing 4- and 5- methylcytosine

Self-complementary oligodeoxynucleotides: GGACCCGGGTCC, GGA4mCCCGGGTCC, GGA5mCCCGGGTCC, CGCGCGCG, CG4mCGCGCG, CG5mCGCGCG were synthetized to study the contribution of methyl groups into the energetics of the three known cooperative transitions in DNA: helix-coil, B-A and B-Z With the use of circular dichroism and absorbtion methods the profiles of the above transitions were obtained by variation of temperature (helix-coil), trifluoroethanol fraction (B-A), NaCl and trifluorethanol contents (B-Z). On the basis of the transition widths and shifts of the transition points due to the methylations the energetics of the methyl groups was estimated. 5mC stabilizes the B form relatively the A form by 0.33 kcal/mol; while 4mC by 0.5 kcal/mol. In the B-Z transition 5 mC stabilizes the Z form by 0.28 kcal/mol relatively the B form; 4mC stabilizes also the Z form although by 0.14 kcal/mol only. Thus, these naturally occurring modifications could modulate substantially the ability of a DNA piece to shift into the A or Z form.