Thermodynamic stability of DNA tandem mismatches

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.     

NMR spectroscopic study of single-stranded DNA fragments of d(CGGCGAAAGCCG) and d(CGGCAAAAGCCG).

We have investigated the structures of two kinds of single-stranded DNA fragments, d(CGGCGAAAGCCG) and d(CGGCAAAAGCCG), by use of NMR spectroscopy. It was found that the former takes a hair-pin like structure stabilized by hydrogen bonding of G/C base pairs in the stem region, while the latter takes a rather extended helical structure. The stable conformation of the former DNA is considered to originate from the stability of the sequence-specific conformation of the loop region rather than the stem region.     

Conformation of DNA modified at a d(GG) or a d(AG) site by the antitumor drug cis-diamminedichloroplatinum(II)

The purpose of this work was the comparison of the conformational changes induced in the double helix by the adducts formed at d(GG) and d(AG) sites in the reaction between the antitumor drug cis-diamminedichloroplatinum(II) (cis-DDP) and DNA. Two duplexes (20-mer) containing either a single d(A*G*) or a single d(G*G) adduct were studied by means of gel electrophoresis and artificial nuclease and chemical probes. It is shown that the d(G*G*) and the d(A*G*) adducts bend DNA similarly, but at the nucleotide level they distort differently the double helix. We suggest that the weaker interactions between platinated A residues and the other nucleotides, as compared to the interactions between platinated G residues and the other nucleotides, are largely responsible for the differences in the distortions induced in DNA by the d(A*G) and d(G*G*) adducts. This suggestion is supported by the study of the distortions induced in duplexes by the d(G*G*) adducts, one of the platinated G residues being paired with a T residue.     

T-T base mismatches enhance drug binding at the branch site in a four-arm DNA junction.

Base mismatches--non Watson-Crick pairing between bases--can arise in duplex DNA as a consequence of mutational events or by recombination. In a duplex, the sequence of the two bases involved, and those flanking the site of mismatch, determines the local structure and extent of destabilization of the helix. Base mismatches can arise also in recombination of nonhomologous strands, and their occurrence in Holliday recombination intermediates can influence the outcome of general or specialized recombination events. We have previously reported that the branch site in a DNA junction can interact selectively with a variety of ligands. Here we describe the thermodynamics of junctions containing T-T mismatches flanking the branch and show that these structures bind methidium and other intercalators with higher affinity than junctions lacking mismatches.     

NMR studies of A.C mismatches in DNA dodecanucleotides at acidic pH. Wobble A(anti).C(anti) pair formation

Two-dimensional proton NMR studies were undertaken on the d(C-G-A-G-A-A-T-T-C-C-C-G) duplex (designated A.C 12-mer) where the A at the mismatch site is flanked by G residues and the d(C-G-C-G-A-A-T-T-C-A-C-G) duplex (designated C.A 12-mer) where the A at the mismatch site is flanked by C residues in an attempt to elucidate the role of flanking base pairs on the structure of the A.C mismatch. The exchangeable and nonexchangeable proton spectra of these two dodecanucleotides have been completely characterized by two-dimensional nuclear Overhauser enhancement (NOE) experiments in H2O and D2O solution at acidic pH. The NOE distance connectivities demonstrate that both A and C at the mismatch site are stacked into a right-handed helix between flanking G.C base pairs and exhibit anti-glycosidic torsion angles. The proton chemical shifts and NOE patterns are consistent with Wobble A.C pairing for the A.C 12-mer and C.A 12-mer duplexes in solution and demonstrate that the A.C mismatches introduce local conformational perturbations that do not extend to the central AATT segment. We detect that amino protons of adenosine (approximately 9.2 ppm) but not of cytidine at the A.C mismatch site in both duplexes on lowering the pH below 6.     

N.m.r. and c.d. studies of the DNA fragments d(TATATATA) and d(TATATA) in solution

DNA fragments d(TATATATA) and d(TATATA) were studied in low-salt aqueous solutions and found to coexist in more than one conformer. 1H-n.m.r. demonstrates that single-stranded and double-stranded states are involved in the conformational coexistence. Circular dichroism spectroscopy indicates a global B-DNA stacking of bases in the fragments. 31P-n.m.r. resonances of the TpA and ApT phosphodiester bonds are substantially separated in the spectra of both d(TATATATA) and d(TATATA) duplexes to suggest an alternating architecture of their backbones. In fact, the oligonucleotide duplexes are much more alternating than the corresponding polynucleotide under the same solution conditions. The alternating character of the d(TATATATA) double helix is further enhanced in molar caesium fluoride solutions. The oligonucleotide isomerization into X-DNA is, however, accompanied by gel formation, which makes high resolution n.m.r. measurements impossible.     

The contrasting structures of mismatched DNA sequences containing looped-out bases (bulges) and multiple mismatches (bubbles).

We have studied the structure and reactivities of two kinds of mismatched DNA sequences--unopposed bases, or bulges, and multiple mismatched pairs of bases. These were generated in a constant sequence environment, in relatively long DNA fragments, using a technique based on heteroduplex formation between sequences cloned into single-stranded M13 phage. The mismatched sequences were studied from two points of view, viz 1. The mobility of the fragments on gel electrophoresis in polyacrylamide was studied in order to examine possible bending of the DNA due to the presence of the mismatch defect. Such bending would constitute a global effect on the conformation of the molecule. 2. Sequences in and around the mismatches were studied using enzyme and chemical probes of DNA structure. This would reveal more local structural effects of the mismatched sequences. We observed that the structures of the bulges and the multiple mismatches appear to be fundamentally different. The bulged sequences exhibited a large gel retardation, consistent with a significant bending of the DNA at the bulge, and whose magnitude depends on the number of mismatched bases. The larger bulges were sensitive to cleavage by single-strand specific nucleases, and modified by diethyl pyrocarbonate (adenines) or osmium tetroxide (thymines) in a non-uniform way, suggesting that the bulges have a precise structure that leads to exposure of some, but not all, of the bases. In contrast the multiple mismatches ('bubbles') cause very much less bending of the DNA fragment in which they occur, and uniform patterns of chemical reactivity along the length of the mismatched sequences, suggesting a less well defined, and possibly flexible, structure. The precise structure of the bulges suggests that such features may be especially significant for recognition by proteins.     

Extraordinary stable structure of short single-stranded DNA fragments containing a specific base sequence: d(GCGAAAGC).

Short single-stranded DNA fragments carrying a GCGAAAGC sequence were found to move unexpectedly faster than other fragments of the same length in electrophoresis on a polyacrylamide gel containing a denaturing agent. The fragments were noted to have a stable structure even in 7M urea solution, but the stability cannot be explained simply on the basis of base pair formation alone. Physical characterization of the GCGAAAGC fragment indicated that it takes a hairpin-like structure in spite of the short chain length with only two G-C base pairs, comprised of GCG and AAAGC subsegments, each possessing a helical configuration independent of the others. Some biological implications of this unusual structure are discussed.     

Synthesis of specific diastereomers of a DNA methylphosphonate heptamer, d(CpCpApApApCpA), and stability of base pairing with the normal DNA octamer d(TPGPTPTPTPGPGPC).

DNA methylphosphonates are candidate derivatives for use in antisense DNA therapy. Their efficacy is limited by weak hybridization. One hypothesis to explain this phenomenon holds that one configuration of the chiral methylphosphonate linkage, Rp, permits stronger base pairing than the other configuration, Sp. To test this hypothesis, four specific pairs of Rp and Sp diastereomers of the DNA methylphosphonate heptamer d(CpCpApApApCpA) were prepared by block coupling of different combinations of individual diastereomers of d(CpCpApA) and d(ApCpA). Each pair of the diastereomers of the heptamer was separated into individual diastereomes using affinity chromatography on a Lichrosorb-NH2 silica column with a covalently attached complementary normal DNA octamer, d(pTpGpTpTpTpGpGpC). The stabilities of complementary complexes of phosphodiester d(TpGpTpTpTpGpGpC) with 8 individual diastereomers of methylphosphonate d(CpCpApApApCpA) were studied by measuring their melting temperatures (Tm). A direct correlation of Tm values with the number of Rp methylphosphonate centers in the heptamer was found: the more Rp centers, the higher the stability of the complex. Tm values for the diastereomers with 6 all-Rp or all-Sp methylphosphonate centers were found to be 30.5 degrees and 12.5 degrees C, respectively, in 100 mM NaCl, 10 mM Na2HPO4, 1 mM EDTA, pH 7.0 with 15 microM of each oligomer. On the average, each substitution of one Rp-center to an Sp-center in the heptamer decreased the Tm by 3 degrees C. Under the same conditions, the Tm of the normal DNA heptamer with its complement was 21 degrees C. These results are consistent with the model that all-Rp methylphosphonate DNAs hybridize much more tightly to complementary normal DNA than do racemic methylphosphonate DNAs, and may therefore exhibit greater potency as antisense inhibitors.     

Conformational consequences of the incorporation of arabinofuranosylcytidine in DNA. An NMR study of the DNA fragments d(CGCTAGCG) and d(CGaCTAGCG) in solution

The self-complementary octamers d(CGCTAGCG) and d(CGaCTAGCG) (aC, arabinofuranosylcytidine) were studied by means of NMR spectroscopy. It is shown that d(CGaCTAGCG), under suitable conditions of oligonucleotide concentration, ionic strength and temperature, exclusively adopts a hairpin structure. However, under the same experimental conditions (5 mM DNA, no added salt, 295 K) d(CGCTAGCG) mainly adopts a B-DNA-type duplex. At lower temperatures (less than or equal to 290 K) the hairpin form of d(CGaCTAGCG) occurs in slow exchange with an intact B-DNA-type duplex. When the DNA concentration of d(CGCTAGCG) is dramatically reduced (less than or equal to 0.5 mM) the hairpin form becomes highly favoured at the expense of the dimer. Moreover, proton-chemical-shift considerations indicate that the structural features of the hairpin structure of d(CGCTAGCG) mimic, in part, those of the modified octamer d(CGaCTAGCG), i.e. a loop comprising only the two central residues with the thymine located into the minor groove (Pieters, J. M. L., de Vroom, E., van der Marel, G. A., van Boom, J. H., Koning, T. M. G., Kaptein, R. and Altona, C. unpublished results). Thermodynamic analysis of d(CGCTAGCG) yields an average Tmd value of 342 K (1 M DNA) and a delta Hod value of -266 kJ/mol for the dimer/coil transition and an average Tmh value of 321 K and delta Hoh - 102 kJ/mol for the hairpin/coil equilibrium. For the duplex/coil equilibrium of d(CGaCTAGCG) an average Tmd value of 336 K (1 M DNA) and delta Hod value of -253 kJ/mol are deduced. The hairpin/coil transition of d(CGaCTAGCG) is characterized by a delta Hoh value of -104 kJ/mol and an average Tmh value of 331 K. It is concluded that incorporation of an arabinofuranosylcytidine in the octamer d(CGaCTAGCG) results in stabilization of the hairpin form, whereas the dimer is destablized by two aC.dG base pairs.     

Structures of DNA polymerase beta with active-site mismatches suggest a transient abasic site intermediate during misincorporation

We report the crystallographic structures of DNA polymerase beta with dG-dAMPCPP and dC-dAMPCPP mismatches in the active site. These premutagenic structures were obtained with a nonhydrolyzable incoming nucleotide analog, dAMPCPP, and Mn(2+). Substituting Mn(2+) for Mg(2+) significantly decreases the fidelity of DNA synthesis. The structures reveal that the enzyme is in a closed conformation like that observed with a matched Watson-Crick base pair. The incorrect dAMPCPP binds in a conformation identical to that observed with the correct nucleotide. To accommodate the incorrect nucleotide and closed protein conformation, the template strand in the vicinity of the active site has shifted upstream over 3 A, removing the coding base from the active site and generating an abasic templating pocket. The primer terminus rotates as its complementary template base is repositioned. This rotation moves O3' of the primer terminus away from the alpha-phosphate of the incoming nucleotide, thereby deterring misincorporation.     

The crystal structure of an RNA oligomer incorporating tandem adenosine-inosine mismatches.

The X-ray crystallographic structure of the RNA duplex [r(CGCAIGCG)]2 has been refined to 2.5 A. It shows a symmetric internal loop of two non-Watson-Crick base pairs which form in the middle of the duplex. The tandem A-I/I-A pairs are related by a crystallographic two-fold axis. Both A(anti)-I(anti) mismatches are in a head-to-head conformation forming hydrogen bonds using the Watson-Crick positions. The octamer duplexes stack above one another in the cell forming a pseudo-infinite helix throughout the crystal. A hydrated calcium ion bridges between the 3'-terminal of one molecule and the backbone of another. The tandem A-I mismatches are incorporated with only minor distortion to the backbone. This is in contrast to the large helical perturbations often produced by sheared G-A pairs in RNA oligonucleotides.     

The interaction of nemorubicin metabolite PNU-159682 with DNA fragments d(CGTACG)2, d(CGATCG)2 and d(CGCGCG)2 shows a strong but reversible binding to G:C base pairs

The antitumor anthracycline nemorubicin is converted by human liver microsomes to a major metabolite, PNU-159682 (PNU), which was found to be much more potent than its parent drug toward cultured tumor cells and in vivo tumor models. The mechanism of action of nemorubicin appears different from other anthracyclines and until now is the object of studies. In fact PNU is deemed to play a dominant, but still unclear, role in the in vivo antitumor activity of nemorubicin. The interaction of PNU with the oligonucleotides d(CGTACG)₂, d(CGATCG)₂ and d(CGCGCG)₂ was studied with a combined use of ¹H and ³¹P NMR spectroscopy and by ESI-mass experiments. The NMR studies allowed to establish that the intercalation between the base pairs of the duplex leads to very stable complexes and at the same time to exclude the formation of covalent bonds. Melting experiments monitored by NMR, allowed to observe with high accuracy the behaviour of the imine protons with temperature, and the results showed that the re-annealing occurs after melting. The formation of reversible complexes was confirmed by HPLC–tandem mass spectra, also combined with endonuclease P1digestion. The MS/MS spectra showed the loss of neutral PNU before breaking the double helix, a behaviour typical of intercalators. After digestion with the enzyme, the spectra did not show any compound with PNU bound to the bases. The evidence of a reversible process appears from both proton and phosphorus NOESY spectra of PNU bound to d(CGTACG)₂ and to d(CGATCG)₂. The dissociation rate constants (k_off) of the slow step of the intercalation process, measured by ³¹P NMR NOE-exchange experiments, showed that the kinetics of the process is slower for PNU than for doxorubicin and nemorubicin, leading to a 10- to 20-fold increase of the residence time of PNU into the intercalation sites, with respect to doxorubicin. A relevant number of NOE interactions allowed to derive a model of the complexes in solution from restrained MD calculations. The conformation of PNU bound to the oligonucleotides was also derived from the coupling constant values.     

A formula for thermal stability (Tm) prediction of PNA/DNA duplexes.

An empirical formula for thermal stability (T m) prediction of PNA/DNA duplexes has been derived. The model is based on the T m as calculated for the corresponding DNA/DNA duplex employing a nearest neighbour approach, by including terms for the pyrimidine content and length of the PNA to take into account the increased thermostability of PNA/DNA hybrids and the asymmetry of the PNA-DNA heteroduplex. The predictive power of the T m prediction formula was challenged with an independent data set not used for model building. The T m of >90% of the sequences was predicted within 5 K; 98% of the predicted T ms differ by not more than 10 K from the experimentally determined T m.     

Electronic structure of DNA by DV-X alpha cluster calculations. I. d(GG).d(CC), d(CG)2, d(GC)2 A and B conformations

The electronic structure of d(GG).d(CC), d(CG)2, d(GC)2 which are stacked base pairs in the DNA double helix, are elucidated for both A and B conformations in detail by DV-X alpha cluster calculations. These three DNA double helix fragments are constructed from the same bases, G and C, but the electronic structure of the fragments for A and B conformations differs from each other characteristically. In particular, the electronic states of the O2 and O3 in phosphates differ drastically from each other, and might play a crucial role as recognition sites in various reaction processes concerning DNA. These differences are caused by the delicate differences in the admixture of the orbital components and the intra- and inter-bases interactions. Contour maps of the wavefunction of the HOMO and LUMO are compared among the stacking isomers.