Effect of N4-methylcytosine and 5-methylcytosine on the stability of the DNA helix

The thermodynamic parameters (delta H, delta S) of the helix-coil transition of self-complementary oligonucleotides d(CGCGCGCG), d(CG5mCGCGCG), d(CG4mCGCGCG), d(GGACCCGGGTCC), d(GGA5mCCCGGGTCC), and d(GGA4mCCCGGGTCC) were determined. The substitution of 4mC for C was found to decrease the melting temperature of the oligonucleotides. The destabilization effect of the two substitutions is equivalent to the change of A.T for G.C pair. The free energy decrease of the helix-coil transition due to the introduction of two 4mC into an octanucleotide was estimated to be 1,24 kcal/mol.     

Synthesis and physical characterization of DNA fragments containing N4-methylcytosine and 5-methylcytosine.

The synthesis of N4-methyl-2'-deoxycytidine and its fully protected mononucleotide, suitable for the oligonucleotide synthesis by phosphotriester method is described. A set of octanucleotides - d(CGCGCGCG), d(CG5mCGCGCG), d(CG4mCGCGCG) and dodecanucleotides - d(GGACCCGGGTCC), d(GGA5mCCCGGGTCC), d(GGA4mCCCGGGTCC) has been synthesized in a solution. Physical characterization of the oligonucleotide duplexes by means of UV and CD spectrometry provides the evidence that 4mC similarly to 5mC favours the B--greater than Z transition, although both of these methylated cytosines inhibit the B--greater than A conformational change. N4-Methylcytosine in contrast to 5-methylcytosine reduces the DNA double helix thermal stability.     

1H NMR and circular dichroism studies of the B and Z conformations of the self-complementary deoxyhexanucleotide d(m5C-G-C-G-m5-C-G): mechanism of the Z-B-coil transitions

The double-helical conformations of d(m5-C-G-C-G-m5-C-G) in aqueous solution were studied by circular dichroism and 1H NMR spectroscopy. In 0.1 M NaCl, only the B form is detected whereas the Z form is strongly predominant in 3 M NaCl. In the presence of 2 M NaCl, two resonance signals corresponding to the B and Z duplexes were observed for each proton below 50 degrees C, indicating a slow exchange between B and Z. However, the B-Z exchange becomes intermediate or fast in the 55-80 degrees C temperature interval. By contrast the exchange between B helix and single-stranded (or coil) forms is much faster for the same temperature conditions. The Z form is only detectable when the coil form is practically absent. With decreasing temperature the B form decreases in favor of the Z form. From proton line-width measurements under various experimental conditions, it was also shown that Z exchanges only with B, while the latter also exchanges with the single-stranded form (S): Z in equilibrium B in equilibrium S. The enthalpy value is about 8 +/- 1 kcal/mol for the B-Z transition and about 40 +/- 2 kcal/mol for the B-S dissociation (2 M NaCl solution). The activation energy is about 47 +/- 2 kcal/mol for the Z----B and 39 +/- 2 kcal/mol for the B----Z reaction. Very good agreement between the experimental results and computed data (based on the above kinetic reaction model) was found for the B, Z, and coil proportions. The B-Z transition of methylated d(C-G)n oligomers is only possible when the Watson-Crick hydrogen bonds between the CG base pairs are firmly maintained; otherwise, the transformation from B to Z would not occur, and B-S dissociation would take place instead.     

Synthesis, structure and thermodynamic properties of 8-methylguanine-containing oligonucleotides: Z-DNA under physiological salt conditions.

Various oligonucleotides containing 8-methylguanine (m8G) have been synthesized and their structures and thermodynamic properties investigated. Introduction Of M8G into DNA sequences markedly stabilizes the Z conformation under low salt conditions. The hexamer d(CGC[M8G]CG)2 exhibits a CD spectrum characteristic of the Z conformation under physiological salt conditions. The NOE-restrained refinement unequivocally demonstrated that d(CGC[m8G]CG)2 adopts a Z structure with all guanines in the syn conformation. The refined NMR structure is very similar to the Z form crystal structure of d(CGCGCG)2, with a root mean square deviation of 0.6 between the two structures. The contribution of m8G to the stabilization of Z-DNA has been estimated from the mid-point NaCl concentrations for the B-Z transition of various m8G-containing oligomers. The presence of m8G in d(CGC[m8G]CG)2 stabilizes the Z conformation by at least deltaG = -0.8 kcal/mol relative to the unmodified hexamer. The Z conformation was further stabilized by increasing the number of m8Gs incorporated and destabilized by incorporating syn-A or syn-T, found respectively in the (A,T)-containing alternating and non-alternating pyrimidine-purine sequences. The results suggest that the chemically less reactive m8G base is a useful agent for studying molecular interactions of Z-DNA or other DNA structures that incorporate syn-G conformation.     

Competing B-Z and helix-coil conformational transitions in supercoiled plasmid DNA.

The formation of melted regions from A + T-rich sequences and left-handed Z-DNA by alternating purine-pyrimidine sequences will both be facilitated by negative supercoiling, and thus if the sequences are present within the same plasmid molecule they will compete for the free energy of supercoiling. We have studied a series of plasmids that contain either (CG)8 or (TG)12 sequences in either G + C or A + T-rich contexts, by means of two-dimensional gel electrophoresis and chemical modification. We observe both B-Z and helix-coil transitions in all plasmids at elevated temperatures and low ionic strength. The plasmids fall into a number of different classes, in terms of the conformational behavior. As the superhelix density is increased, pCG8/vec ((CG)8 in G + C-rich context) undergoes an initial B-Z transition, followed by melting transitions in sequences remote from the (CG)8 sequence. The two transitions are coupled through the topology of the molecule but are otherwise independent. When the (CG)8 sequence was placed in an A + T-rich context (pCG8/col), the helix-coil transition was perturbed by the presence of the Z-DNA segment. Replacement of the (CG)8 tracts by (TG)12 sequences resulted in a further level of interaction between the transitions. Statistical mechanical modeling of the transitions suggested that at intermediate levels of negative supercoiling the Z-DNA formed by the (TG)12 sequence has a lowered probability due to the helix-coil transition in the A + T-rich sequences. These studies illustrate the complexities of competing conformational equilibria in supercoiled DNA molecules.     

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.     

Thermodynamic behaviour of the heptadecadeoxynucleotide d(CGCGCGTTTTTCGCGCG) forming B and Z hairpins in aqueous solution.

UV and CD data of the partially self-complementary heptadecadeoxynucleotide d(CGCGCGTTTTTCGCGCG), obtained as a function of temperature, salt and strand concentration, show that: at low NaCl and strand concentration the oligomer exhibits, on increasing the temperature, a biphasic thermal profile which is indicative of two structural transitions, from dimeric duplex to hairpin and from hairpin to coil; the loop stabilizes enthalpically both B and Z hairpin structures with respect to the corresponding unconstrained hexamer d(CGCGCG) by a few Kcal/mol; the oligomer undergoes a B-Z transition which appears to be complete, at 0 degree C, when induced by NaClO4; by contrast the B-Z transition induced by NaCl does not reach completeness even at salt saturation. The independence of the denaturation temperature, at high salt conditions, on the oligomer concentration indicates that the Z structure is present also in the hairpin.     

"Tie" calorimetry: theory and experiment

The results of theoretical and experimental investigations of the interaction of ligands to DNA permit to suggest the method of "tie" calorimetry by means of which the thermodynamic parameters of the cooperative transitions of DNA such as helix-coil, B-A and B-Z transitions were obtained. The prominence of the method is that the experimental samples contain such concentration of DNA that intermolecular interaction may be neglected. Moreover the smaller is the magnitude of the measured parameter the greater is the accuracy of the experimentally obtained data.     

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.     

Thermodynamics of B-Z transition of poly[d(G-C)] in a water-ethanol solution. The "tie calorimetry"

Poly[d(G-C)] in a 55% ethanol solution undergoes a transition from the Z form to the B form when the temperature is increased from 20 degrees to 50 degrees C. The enthalpy of the transition, delta HBA = -1.4 kcal/mol, has been determined with a "tie" polyamine which stabilizes the Z conformation. This value has been shown to be practically independent of ionic strength within the range of 5 X 10(-4) M - 2 X 10(-3) M NaCl.     

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.     

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.     

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.     

Quantitative analysis of DNA secondary structure from solvent-accessible surfaces: the B- to Z-DNA transition as a model

Solvent structure and its interactions have been suggested to play a critical role in defining the conformation of polynucleotides and other macromolecules. In this work, we attempt to quantitate solvent effects on the well-studied conformational transition between right-handed B- and left-handed Z-DNA. The solvent-accessible surfaces of the hexamer sequences d(m5CG)3, d(CG)3, d(CA)3, and d(TA)3 were calculated in their B- and Z-DNA conformations. The difference in hydration free energies between the Z and the B conformations (delta delta GH(Z-B] was determined from these surfaces to be -0.494 kcal/mol for C-5 methylated d(CG), 0.228 kcal/mol for unmethylated d(CG), 0.756 kcal/mol for d(CA)-d(TG), and 0.896 kcal/mol for d(TA) dinucleotides. These delta delta GH(Z-B) values were compared to the experimental B- to Z-DNA transition energies of -0.56 kcal/mol that we measured for C-5 methylated d(CG), 0.69-1.30 kcal/mol reported for unmethylated d(CG), 1.32-1.48 kcal/mol reported for d(CA)-d(TG), and 2.3-2.4 kcal/mol for d(TA) dinucleotides. From this comparison, we found that the calculated delta delta GH(Z-B) of these dinucleotides could account for the previous observation that the dinucleotides were ordered as d(m5CG) greater than d(CG) greater than d(CA)-d(TG) greater than d(TA) in stability as Z-DNA. Furthermore, we predicted that one of the primary reasons for the inability of d(TA) sequences to form Z-DNA results from a decrease in exposed hydrophilic surfaces of adjacent base pairs due to the C-5 methyl group of thymine; thus, d(UA) dinucleotides should be more stable as Z-DNA than the analogous d(TA) dinucleotides.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

CC/GG contacts facilitate the B to A transition of DNA in solution

Self-complementary decadeoxynucleotides, CCGATATCGG, CCAGATCTGG, CCCTGCAGGG, GGGGGCCCCC, were designed and synthesized to estimate the A-philic free energy of CC/GG contacts. First, regions of temperature-stability of the double-stranded conformation were determined for each 10-mer. Then, circular dichroism spectra were recorded for the B-family forms at different temperatures, counter-ion concentrations and trifluoroethanol contents. A cooperative change typical of the B-A transition is observed in the CD spectra at a trifluoroethanol content specific for each duplex. The positions of half-transition points were functions not only of the nucleotide sequence but of the duplex length as well: the B to A transitions were hindered in these 10-mers in comparison with a lengthy DNA. The B-phility value was estimated to be 3 kcal/mol of 10-mer. The B-A transition point was shown to drop with an increase in the number of CC/GG contacts in a duplex. The designed 10-mers made it possible to estimate quantitatively the A-phility of CC/GG contact as compared with an average DNA: (FA-FB)CC = 0.2 Kcal/mol, (FA-FB)DNA = 0.7 Kcal/mol.     

Thermodynamic parameters are sequence-dependent for the supercoil-induced B to Z transition in recombinant plasmids

The entropy and enthalpy changes which contribute to the thermodynamics of the B to Z transition were determined for three recombinant plasmids containing a (dC-dG)16 tract and for a plasmid containing a pair of (dT-dG)20 regions. For each base pair which adopts a left-handed conformation in the plasmids with (dC-dG)16 sequences, the delta HBZ and delta SBZ are -2.1 kcal/mol bp and -8.8 cal/K-mol bp, respectively. In the plasmid containing the (dT-dG)20 tracts, however, the delta HBZ and delta SBZ values are 0.58 kcal/mol bp and -0.76 cal/K-mol bp, respectively. Also, these determinations show that for each B-Z junction that forms in the plasmids containing the (dC-dG), the enthalpy and entropy changes are 24 kcal/mol junction and 65 cal/K-mol junction, whereas for the (dT-dG) plasmid, the enthalpy and entropy changes are -1.8 kcal/mol junction and -22 cal/K-mol junction, respectively. Those values for the enthalpy and entropy changes for the formation of a BZ junction in (dC-dG) and (dT-dG) plasmids suggest that the properties and possibly the structures of the junctions are different. Calculations using the enthalpy and entropy changes determined in this study reveal that the B to Z transition in plasmids containing (dC-dG) blocks are more temperature-dependent than the transitions in plasmids with (dT-dG) blocks. Surprisingly, at temperatures above 60 degrees C, calculations indicate that the B to Z transitions in (dT-dG) plasmids should be energetically favored over that transition in (dC-dG) plasmids.     

Base sequence criteria and Cartesian coordinates for stable B/Z and B/Z/B junctions in relaxed DNA.

It seems increasingly evident that if the Z form of DNA exists in the genome it must exist as short sections of alternating pyrimidine-purine sequences in the midst of very long sections of B-form DNA. We have determined the minimum length of a string of alternating CG base pairs that can go into the Z form in the middle of a long B form. Self-complimentary oligomers of the form T(M)(CG)(N)A(M) were synthesized. The conformation of the resulting duplex was determined in 6M aqueous NaCl solution by Raman scattering. We have found that 12 alternating CG base pairs is the minimum length required to form a stable Z form of DNA inside of a long B form section. Only the 4 center CG base pairs go into the Z form. These 4 CG base pairs in the Z form are flanked on each side by 4 CG base pairs in a non-Z (probably B) form as well as the ..TT.. ..AA.. sequences in the B form. We propose a model of the B/Z junction in which the double helix flips directly from the B form to the Z form so that there are no base pairs in the junction. In this model the B form is nucleated in the AT base pairs on each end and is propagated into the CG base pairs in the center. This model is supported by isotopic H/D exchange experiments that shows that the H/D exchange of the non-Z form CG base pairs is highly retarded and indicates that they remain in the B form. A Thermodynamic analysis of the concentration dependence of the melting point of the duplexes in both low and high salt, supports our model and rules out the possibility of hairpin formation. The enthalpy for the formation of a B/Z junction is determined to be about +16 kcal/junction. A comparison of these results with recent results on B/Z junctions in super-coiled DNA is given. Molecular modeling calculations permit us to obtain values for the coordinates and torsional angles of the oligomers showing both B/Z and B/Z/B junctions. The Cartesian coordinates for these oligomers as well as stereo figures of these models in color are available from the authors.     

B,Z conformations and mechanism of the Z-B-coil transitions of the self-complementary deoxy-hexanucleotide d(C-G-m5C-G-C-G) by 1H-NMR and CD spectroscopy

The helical structures of d(C-G-m5C-G-C-G) were studied in aqueous solution at various salt concentrations and temperatures by CD and 1H-NMR spectroscopy. At room temperature only the B form is observed in 0.1 M NaCl whereas the B and Z forms are simultaneously present in 1.8 M NaCl. At high salt concentration (4 M NaCl) the Z form is largely predominant (greater than 95%). The Z form proton resonances were assigned by using the polarisation transfer method (between B and Z at 1.8 M NaCl) and by proton-proton decoupling (at high salt concentration). The Z-B-Coil transitions were studied as a function of temperature with the 1.8 M NaCl solution. At high temperature (95 degrees C) only the coil form (S) is present. Below 55 degrees C the coil proportion is negligible, and the B-Z exchange is slow. The disappearance of the coil gives rise at first to the B form and on lowering the temperature the Z proportion increases to the detriment of the B form. Proton linewidth, relaxation and polarisation transfer studies confirm the conclusion in the previous report on d(m5C-G-C-G-m5C-G) (Tran-Dinh et al Biochemistry 1984 in the press) that Z exchanges only with B whereas the latter also exchanges with S,Z in equilibrium B in equilibrium S. The present data show that even at high salt concentration where only the Z form of d(C-G-m5C-G-C-G) is observed the Z-S transition also passes through the B form as an intermediate stage. The B-Z transition takes place when the Watson-Crick hydrogen bonds are firmly maintained and is greatly favoured when there are three hydrogen bonds between the base-pairs.     

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.