An NMR study of the polymorphous behavior of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) in solution. The B, Z, and hairpin forms

Effects of interrupting residues on DNA dumbbell structures formed by CCTG tetranucleotide repeats associated with myotonic dystrophy type 2

Myotonic dystrophy type 2 (DM2) is a neurogenerative disease caused by caprylic/capric triglyceride (CCTG) tetranucleotide repeat expansions in intron 1 of the cellular nucleic acid-binding protein (CNBP) gene. Non-B DNA structures formed by CCTG repeats can promote genetic instability, whereas interrupting motifs of NCTG (N = A/T/G) within CCTG repeats help to maintain genomic stability. However, whether the interrupting motifs can affect DNA structures of CCTG repeats remains unclear. Here, we report that four CCTG repeats with an interrupting 3'-A/T/G residue formed dumbbell structures, whereas a non-interrupting 3'-C residue resulted in a multi-loop structure exhibiting conformational dynamics that may contribute to a higher tendency of escaping from DNA mismatch repair and causing repeat expansions. The results provide new structural insights into the genetic instability of CCTG repeats in DM2.

High-Resolution Structures of DNA Minidumbbells Comprising Type II Tetraloops with a Purine Minor Groove Residue

Minidumbbell (MDB) is a newly discovered DNA structure formed by native sequences, which serves as a possible structural intermediate causing repeat expansion mutations in the genome and also a functional structural motif in constructing DNA-based molecular switches. Until now, all the reported MDBs containing two adjacent type II tetraloops were formed by pyrimidine-rich sequences 5'-YYYR YYYR-3' (Y and R represent pyrimidine and purine, respectively), wherein the second and sixth residues folded into the minor groove and interacted with each other. In this study, we have conducted a high-resolution nuclear magnetic resonance (NMR) spectroscopic investigation on alternative MDB-forming sequences and discovered that an MDB could also be formed stably with a purine in the minor groove, which has never been observed in any previously reported DNA type II tetraloops. Our refined NMR solution structures of the two MDBs formed by 5'-CTTG CATG-3' and 5'-CTTG CGTG-3' reveal that the sixth purine residue was driven into the minor groove via base-base stacking with the second thymine residue and adenine stacked better than guanine. The results of our present research work expand the sequence criteria for the formation of MDBs and shed light to explore the significance of MDBs.

In vitro biosynthesis and chemical identification of UDP-N-acetyl-d-quinovosamine (UDP-d-QuiNAc)

N-acetyl-d-quinovosamine (2-acetamido-2,6-dideoxy-d-glucose, QuiNAc) occurs in the polysaccharide structures of many Gram-negative bacteria. In the biosynthesis of QuiNAc-containing polysaccharides, UDP-QuiNAc is the hypothetical donor of the QuiNAc residue. Biosynthesis of UDP-QuiNAc has been proposed to occur by 4,6-dehydration of UDP-N-acetyl-d-glucosamine (UDP-GlcNAc) to UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose followed by reduction of this 4-keto intermediate to UDP-QuiNAc. Several specific dehydratases are known to catalyze the first proposed step. A specific reductase for the last step has not been demonstrated in vitro, but previous mutant analysis suggested that Rhizobium etli gene wreQ might encode this reductase. Therefore, this gene was cloned and expressed in Escherichia coli, and the resulting His6-tagged WreQ protein was purified. It was tested for 4-reductase activity by adding it and NAD(P)H to reaction mixtures in which 4,6-dehydratase WbpM had acted on the precursor substrate UDP-GlcNAc. Thin layer chromatography of the nucleotide sugars in the mixture at various stages of the reaction showed that WbpM converted UDP-GlcNAc completely to what was shown to be its 4-keto-6-deoxy derivative by NMR and that addition of WreQ and NADH led to formation of a third compound. Combined gas chromatography-mass spectrometry analysis of acid hydrolysates of the final reaction mixture showed that a quinovosamine moiety had been synthesized after WreQ addition. The two-step reaction progress also was monitored in real time by NMR. The final UDP-sugar product after WreQ addition was purified and determined to be UDP-d-QuiNAc by one-dimensional and two-dimensional NMR experiments. These results confirmed that WreQ has UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose 4-reductase activity, completing a pathway for UDP-d-QuiNAc synthesis in vitro.

NMR Structure of the (+)-CPI-indole/d(GACTAATTGAC)-d(GTCAATTAGTC) Covalent Complex

We report the NMR solution structure of (+)-CPI-indole (CPI, 1,2,8,8a-tetrahydrocyclopropa[c]pyrrolo[3,2-e]indol-4(5H)-one), an agent belonging to the CC-1065/duocarmycin family of antitumor compounds. This (+)-CPI-indole structure is covalently bound to d(G(1)ACTAATTGTC(11))-d(G(12)TCAATTAGTC(22)), a synthetic DNA duplex containing a high-affinity binding site. The three-dimensional structure has been determined by several cycles of restrained molecular dynamics calculations with a total of 563 NMR-derived constraints, both in vacuo and by using the generalized Born solvent continuum model. In-depth analysis of the structure of this ligand-DNA complex led to a detailed knowledge of the bound state conformation of the CPI-indole, the most simplified agent related to CC-1065 and duocarmycins, the parent members of a family of extremely potent antitumor compounds. Comparison of the CPI-indole bound conformation with those previously found for (+)-duocarmycin SA (DSA), its unnatural enantiomer (-)-DSA, and the demethoxylated analogue (+)-DSI in their DNA complexes provided additional evidence of the tight correlation between the catalytic effect exerted by DNA on the alkylation reaction and the extent of angular twist between the two planar heteroaromatic subunits of these agents. Additionally, comparison of the structural features of the DNA-bound state of a "naked" ligand, such as CPI-indole, with those of various other duocarmycin agents provided useful information for the interpretation of the observed effects on chemical reactivity of the different substitution patterns at the hemispheres of these types of complex.

Protonation of deoxycytidine residues in dC4 tetraloops: UV spectrophotometric study of dC10 and d(A14C4T14)

It is shown that component analysis could be applied to study the UV difference spectra of cytidine oligomers and hairpin oligonucleotides with cytidines in the loop region in order to account for the melting and titration results in terms of cytidine stacking and protonation. Upon acid titration, the dC(10) oligomer undergoes cooperative conformational transition at pH 6.3 accompanied by protonation and formation of the i-structure with half of the residues protonated. The stability of the hemiprotonated structure increases with decreasing pH, the i-structure persisting still in the region of pH<pK of cytidine. An UV difference spectrum that reflects the stacking/unstacking of hemiprotonated cytidine residues was acquired from the melting and titration experiments of the dC(10) oligomer and used to describe the behavior of the dC(4) loop of the hairpin oligonucleotide d(A(14)C(4)T(14)). It is shown that upon titration, the 50% level of protonation of the deoxycytidine tetraloop is attained at pH 5.0. Simultaneously, the stacking interactions of cytidine residues reach the maximum at this pH with two residues stacked, and thereafter decline again. Only marginal stabilization of the oligomer hairpin (DeltaT(m)=1.5 degrees C) is found to accompany the formation of this single hemiprotonated dC.dC(+) base pair. We propose that at pH 5 the cytidines of the dC(4) loop form a hemiprotonated dC.dC(+) pair stacked with the last dA.dT base pair of the hairpin stem.

Solution structures of the Huntington's disease DNA triplets, (CAG)n

Highly polymorphic DNA triplet repeats, (CAG)n, are located inside the first exon of the Huntington's disease gene. Inordinate expansion of this repeat is correlated with the onset and progression of the disease. NMR spectroscopy, gel electrophoresis, digestion by single-strand specific P1 enzyme, and in vitro replication assay have been used to investigate the structural basis of (CAG)n expansion. Nondenaturing gel electrophoresis and 1D 1H NMR studies of (CAG)5 and (CAG)6 reveal the presence of hairpins and mismatched duplexes as the major and minor populations respectively. However, at high DNA concentrations (i.e., 1.0-2.0 mM that is typically required for 2D NMR experiments) both (CAG)5 and (CAG)6 exist predominantly in mismatched duplex forms. Mismatched duplex structures of (CAG)5 and (CAG)6 are useful, because they adequately model the stem of the biologically relevant hairpins formed by (CAG)n. We, therefore, performed detailed NMR spectroscopic studies on the duplexes of (CAG)5 and (CAG)6. We also studied a model duplex, (CGCAGCG)2 that contains the underlined building block of the duplex. This duplex shows the following structural characteristics: (i) all the nucleotides are in (C2'-endo, anti) conformations, (ii) mismatched A x A base pairs are flanked by two Watson-Crick G x C base pairs and (iii) A x A base pairs are stably stacked (and intra-helical) and are formed by a single N6-H--N1 hydrogen bond. The nature of A x A pairing is confirmed by temperature-dependent HMQC and HMQC-NOESY experiments on the [(CA*G)5]2 duplex where the adenines are 15N-labeled at N6. Temperature- and pH-dependent imino proton spectra, nondenaturing electrophoresis, and P1 digestion data demonstrate that under a wide range of solution conditions longer (CAG)n repeats (n> or =10) exist exclusively in hairpin conformation with two single-stranded loops. Finally, an in vitro replication assay with (CAG)8,21 inserts in the M13 single-stranded DNA templates shows a replication bypass for the (CAG)21 insert but not for the (CAG)8 insert in the template. This demonstrates that for a sufficiently long insert (n=21 in this case), a hairpin is formed by the (CAG)n even in presence of its complementary strand. This observation implies that the formation of hairpin by the (CAG)n may cause slippage during replication and thus may explain the observed length polymorphism.

High resolution solution structure of a DNA duplex alkylated by the antitumor agent duocarmycin SA

The three-dimensional solution structure of duocarmycin SA in complex with d-(G1ACTAATTGAC11).d-(G12TCATTAGTC22) has been determined by restrained molecular dynamics and relaxation matrix calculations using experimental NOE distance and torsion angle constraints derived from 1H NMR spectroscopy. The final input data consisted of a total of 858 distance and 189 dihedral angle constraints, an average of 46 constraints per residue. In the ensemble of 20 final structures, there were no distance constraint violations >0.06 A or torsion angle violations >0.8 degrees. The average pairwise root mean square deviation (RMSD) over all 20 structures for the binding site region is 0.57 A (average RMSD from the mean: 0.39 A). Although the DNA is very B-like, the sugar-phosphate backbone torsion angles beta, epsilon, and zeta are distorted from standard values in the binding site region. The structure reveals site-specific bonding of duocarmycin SA at the N3 position of adenine 19 in the AT-rich minor groove of the duplex and binding stabilization via hydrophobic interactions. Comparisons have been made to the structure of a closely related complex of duocarmycin A bound to an AT-rich DNA duplex. These results provide insights into critical aspects of the alkylation site selectivity and source of catalysis of the DNA alkylating agents, and the unusual stability of the resulting adducts.

Conformational properties of Z-forming DNA oligomers bearing terminal unpaired bases

Three sets of semi-self-complementary deoxyribonucleotide decamers with the sequence XX-(5meCG)4, (5meCG)4-XX, or Y-(5meCG)4-Y, where XX = AA, CC, GG, or TT and Y = A, C, G, or T, were synthesized along with the self-complementary octamer (5meCG)4. The 8-mer duplex readily undergoes a B-to-Z conformational conversion upon increasing the NaCl concentration with a transitional midpoint of approximately 1.1 M NaCl. The 10-mers should form 8-bp duplexes a with core sequence of [(5meCG)4]2 with 5'-XX overhangs, 3'-XX overhangs, or 5',3'-Y/Y mismatches. Circular dichroism was employed to determine the conformations of all oligomers. Salt titrations were performed to measure the effect of overhangs and terminal mismatches on the B-to-Z conversion. In general, the presence of 5'-XX overhangs results in a transition midpoint equal to or slightly higher than the control, whereas the presence of 3'-XX overhangs results in a transition midpoint slightly lower than the control. The 3'-CC and 5'-GG overhangs are exceptions, with transition midpoints much higher than the control. These oligomers apparently form duplexes with 5',3'-C/C or 5',3'-G/G mismatches abutting a [(G5meC)4]2 duplex core. The presence of terminal mismatches in the third set of oligomers results in transition midpoints higher than the control. Ultraviolet absorbance methods were used to evaluate the effect of the various stacking motifs of the 10-mers on the thermodynamics of melting relative to the 8-mer for both B and Z conformations. We found that in both the B and Z conformations, the presence of an overhang stabilizes the [(5meCG)4]2 duplex, with the 5' overhangs having a greater stabilizing effect relative to the 3' overhangs. The presence of 5',3'-Y/Y mismatches also imparts a stabilizing effect on the control 8-mer in both the B and Z conformations. These results are discussed in terms of stacking interactions of the terminal unpaired bases.

Slow conformational exchange in DNA minihairpin loops: a conformational study of the circular dumbbell d<pCGC-TT-GCG-TT>

In recent years various examples of highly stable two-residue hairpin loops (miniloops) in DNA have been encountered. As the detailed structure and stability of miniloops appear to be determined not only by the nature and sequence of the two bases in the loop, but also by the closing base pair, it is desirable to carry out in-depth studies of especially designed small model DNA compounds. Therefore, a circular DNA dumbbell-like molecule is tailored to consist of a stem of three Watson-Crick base pairs, flanked on each side by a minihairpin loop. The resulting circular DNA decamer 5'-d<pCGC-TT-GCG-TT>-3' (I) is studied in solution by means of nmr spectroscopy. At a temperature of 269 K the molecule occurs in a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4). L2L2 contains three Watson-Crick C-G base pairs and two two-residue loops (H2-family type) in opposite parts of the molecule. On raising the temperature from 269 to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop with only one remaining solvent-accessible hydrogen bond between NH alpha of the cytosine C(1) and O6 of the guanine G(8), whereas the opposite 5'-CTTG-3' loop remains stable. The disruption of the C(1)-G(8) base pair in the L2L4 form is correlated with the presence of a syn orientation for the C(1) base at the 5'-3' loop-stem junction in the 5'-GTTC-3' loop. The two conformers, L2L2 and L2L4, occur in slow equilibrium (2-20 s-1). Moderate line broadening of specific 1H, 13C, and 31P resonances of residues C(1), G(8), T(9), and T(10) at low temperatures, due to chemical exchange between L2L2 and L2L4, show that the interconversion from an anti to syn conformer in residue C(1) has a small local effect on the structure of the dumbbell. T1 relaxation measurements, chemical-shift considerations, and complete band-shape calculations of the exchange process of the G(8) imino proton reveal a possibility for the existence of multiconformational states in the anti-syn equilibrium.

Structural characterization of d(CAACCCGTTG) and d(CAACGGGTTG) mini-hairpin loops by heteronuclear NMR: the effects of purines versus pyrimidines in DNA hairpins

The DNA decamers, d(CAACCCGTTG) and d(CAACGGGTTG) were studied in solution by proton and heteronuclear NMR. Under appropriate conditions of pH, temperature, salt concentration and DNA concentration, both decamers form hairpin conformations with similar stabilities [Avizonis and Kearns (1995) Biopolymers, 35, 187-200]. Both decamers adopt mini-hairpin loops, where the first and last four nucleotides are involved in Watson-Crick hydrogen bonding and the central two nucleotides, CC or GG respectively, form the loop. Through the use of proton-proton, proton-phosphorus and natural abundance proton-carbon NMR experiments, backbone torsion angles (beta, gamma and epsilon), sugar puckers and interproton distances were measured. The nucleotides forming the loops of these decamers were found to stack upon one another in an L1 type of loop conformation. Both show gamma tr and unusual beta torsion angles in the loop-closing nucleotide G7, as expected for mini-hairpin loop formation. Our results indicate that the beta and epsilon torsion angles of the fifth and sixth nucleotides that form the loop and the loop-closing nucleotide G7 are not in the standard trans conformation as found in B-DNA. Although the loop structures calculated from NMR-derived constraints are not well defined, the stacking of the bases in the two different hairpins is different. This difference in the base stacking of the loop may provide an explanation as to why the cytosine-containing hairpin is thermodynamically more stable than the guanine-containing hairpin.

Sampling of the conformations of the d(CGCTGCGGC) hairpin in solution by two-dimensional nuclear magnetic resonance and theoretical methods

Most NMR studies of DNA oligomers have focused on rigid structures that show a strong preference for one or a small set of ground-state conformations. There is an increasing interest in extending NMR methods to investigate DNA systems in which this preference does not exist. A DNA hairpin is one such system where a large number of low-energy structures coexist in solution. In this article we show how 1D/2D NMR data of the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin are used to map the conformational space of this molecule. First, we characterize the gross morphology of the hairpin by monitoring the exchangeable imino signals in the molecule. Second, we extract a set of inter-proton distances (i.e., the average values and the associated dispersions) for various pairwise interactions by performing full-matrix NOESY simulation with respect to the observed NOESY data for mixing times of 250 and 100 ms. Third, we use these distances as structural constraints to perform a 300-ps molecular dynamics simulation at 500 K. Fourth, we extract 600 snapshots (one after every 0.5 ps) from the MD trajectory and perform constrained energy minimization to map local minima on the sampled energy surface (we call this the rapid temperature quenching step). Fifth, we assign 600 structures to 14 disjoint clusters such that conformationally similar hairpins belong to the same cluster while conformationally distinct hairpins belong to different clusters. Finally, we interpret the NOESY data in terms of conformationally distinct structures by recalculating NOESY contributions taken from representative structures of different clusters. Our analyses clearly demonstrate that the NMR data correspond to an ensemble of distinct structures, i.e., a set of energetically stable but conformationally distinct structures that satisfies the constraints of loop folding in the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin. Two types of loop folding consistent with NMR data are obtained: (i) a hairpin with two G.C pairs in the stem and four residues in the loop and (ii) a hairpin with two G.C pairs and a reverse wobble G.T pair in the stem plus two residues in the loop.

An NMR study of the conformation and thermodynamics of the circular dumbbell d [formula: see text] Slow exchange between two- and four-membered hairpin loops

The circular DNA decamer 5'-d [formula: see text] 3' is studied in solution by means of NMR spectroscopy. At low temperature the molecule adopts a dumbbell structure with three Watson-Crick C-G base pairs and two two-residue loops in opposite parts of the molecule. On raising the temperature another conformer appears, in which the closing C-G base pair in the 5'-GTTC-3' loop is disrupted, whereas the opposite 5'-CTTG-3' loop remains stable. The two conformers are in slow equilibrium over a limited temperature range.

A proton NMR study of a DNA dumb-bell structure with hairpin loops of only two nucleotides: d(CACGTGTGTGCGTGCA)

One- and two-dimensional nuclear magnetic resonance (NMR) experiments were used to study the conformation of the DNA hexadecanucleotide d(CACGTGTGTGCGTGCA) in aqueous solution. NMR spectra were recorded for the compound in D2O and in H2O/D2O (90/10) over the temperature range 1 degree C-60 degrees C. Assignments of imino proton resonances and of non-exchangeable proton resonances (except for some H4', H5' and H5" resonances) are given. The 1H-NMR spectra indicate that below about 20 degrees C, the compound exists as a single monomolecular species. Between 20 degrees C and 55 degrees C the oligonucleotide occurs as a mixture of structures in fast exchange on the NMR time scale, except for the temperature region 30 degrees - 34 degrees C, where substantial line broadening indicates intermediate exchange; above 60 degrees C the single strand predominates. The imino proton spectra, chemical shift values, and scalar coupling and NOE data reveal that the monomeric form, which is exclusively present below 20 degrees C, consists of a structure with a B-DNA double helix region of six base pairs, both ends of which are closed by hairpin loops of only two nucleotides, giving the molecule a dumbbell-like structure: [sequence: see text].

NMR and molecular modeling evidence for a G.A mismatch base pair in a purine-rich DNA duplex

1H NMR experiments indicate that the oligomer 5'-d(ATGAGCGAATA) forms an unusual 10-base-pair duplex with 4 G.A base pairs (underlined) and a 3' unpaired adenosine. NMR results indicate that guanosine imino protons of the G.A mismatches are not hydrogen bonded but are stacked in the helix. A G----I substitution in either G.A base pair causes a dramatic decrease in duplex stability and indicates that hydrogen bonding of the guanosine amino group is critical. Nuclear Overhauser effect spectroscopy (NOESY) and two-dimensional correlated spectroscopy (COSY) results indicate that the overall duplex conformation is in the B-family. Cross-strand NOEs in two-dimensional NOESY spectra between a mismatched AH2 and an AH1' of the other mismatched base pair and between a mismatched GH8 and GNH1 of the other mismatch establish a purine-purine stacking pattern, adenosine over adenosine and guanosine over guanosine, which strongly stabilizes the duplex. A computer graphics molecular model of the unusual duplex was constructed with G.A base pairs containing A-NH2 to GN3 and G-NH2 to AN7 hydrogen bonds and B-form base pairs on both sides of the G.A pairs [5'-d(ATGAGC)]. The energy-minimized duplex satisfies all experimental constraints from NOESY and COSY results. A hydrogen bond from G-NH2 of the mismatch to a phosphate oxygen is predicted.

Circular dichroism studies of an oligodeoxyribonucleotide containing a hairpin loop made of a hexaethylene glycol chain: conformation and stability

An oligodeoxyribonucleotide, d(GCTCACAAT-X-ATTGTGAGC), where X represents a hexaethylene glycol chain, was studied using circular dichroism spectroscopy. Its conformation and conformational stability were compared to those of compounds where X was replaced by four thymines and to the duplex of same sequence without loop. The compound with the hexaethylene glycol chain can form a hairpin looped structure as well as a bulged duplex structure. In both cases the duplex region of the oligodeoxyribonucleotide exhibits the same conformation. In similar conditions the oligodeoxyribonucleotide with a four thymines loop forms exclusively a hairpin structure. Comparison between the thermodynamic parameters (delta H, delta S, delta G) associated with the formation of the structure of the three compounds are presented. In the case of the compound with the hexaethylene glycol chain it is shown that the large increase in its melting temperature (by about 35 degrees C in our experimental conditions) when compared to the non looped structure is mainly due to the fact that its melting process is intramolecular (monomolecular) whereas the other one is bimolecular.

Energetics of the hairpin to mismatched duplex transition of d(GCCGCAGC) on NaCl solution

We report Potential of Mean Force studies to describe the relative thermodynamic stabilities of d(GCCGCAGC) in a mismatched duplex and a hairpin monomer conformation in NaCl solution. The PMF calculations are combined with previous molecular mechanics and normal mode analysis in order to estimate the role of different components of the free energy in determining the relative stability of the duplex and hairpin structures. The high entropy associated with the loop region and the lack of minor groove phosphate-phosphate interactions in the hairpin compete against the gain in enthalpic contribution to the free energy due to base pairing in the mismatched duplex. The combined free energy calculations show that the hairpin is the most stable conformation at low salt and that a hairpin to duplex transition takes place at approximately 0.47 M NaCl. In addition, we studied the hairpin to partially stacked single helical conformation equilibrium at low salt. We found a small variation in transition temperature in salt concentration, delta Tm/delta log10(cs) approximately 2-3 degrees K/decade, in contrast to the duplex to hairpin or duplex to partially stacked single helix transition where the transition temperature exhibited marked dependence on salt concentration. This is in qualitative agreement with experimental data. Based on the Potential of Mean Force free energy calculation, the order of relative stability of the three-conformations studied varies with salt concentration. We observed the following orders of stability: stacked single helix greater than hairpin greater than duplex for cs less than 0.77 M NaCl; single helix greater than duplex greater than hairpin for 0.77 less than Cs less than 2.1 M; and duplex greater than hairpin greater than single strand for cs greater than 2.1 M. From the calculated PMF free energy curves in the NaCl concentration range, 0.012 less than cs less than 5.0 M, we can assign upper and lower bounds for the non-ionic differences in free energy between the duplex, hairpin, and stacked single helical states (at standard conditions: cs = 1.0 M, T = 25 degrees C, and 1 M oligomer concentration). We found that for delta G duplex single helix = G duplex - 2 x G single helix less than -7.38 Kcal/mol, the single helix is the least stable state. For the duplex-to-hairpin free energy difference in the range, -1.87 less than delta G duplex-hairpin less than 0.03 Kcal/mol, there will always be a salt-induced hairpin-to-duplex transition for 0.01 less than cs less than 1.6 M NaCl. If delta G duplex-hairpin less than -1.87, the duplex is always more stable than the hairpin; and for delta G duplex-hairpin greater than Kcal/mol, the hairpin state is always more stable than the duplex, for all salt concentrations.

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.

Hairpin formation of d(CGCG-TA-CGCG), d(CGCG-TG-CGCG) and their cytosine methylated analogs

Hairpin formations of decamers d(CGCG-TA-CGCG), d(CGCG-TG-CGCG), and their m5dC analogs are evidenced by the existence of biphasic absorbance melting profiles in which the lower transition temperature increases with increasing oligomer concentration, whereas the higher melting temperature is concentration independent. The corresponding temperature dependent CD intensity at 285 nm exhibits a maximum around 55 degrees C. These observations are consistent with the interpretation that the lower temperature transition corresponds to the duplex to hairpin transformation while the melting of hairpins into single strands constitutes the higher temperature transition. The CD spectrum of the hairpin conformation appears to be characterized by a couplet with nearly equal positive and negative intensities at 285 and 255 nm, respectively, while a significantly smaller intensity at 285 nm is apparent for the duplex form. The hairpin conformation is suspected to contain a two-nucleotide loop. Titrations with NaCl further suggest that, in contrast to the TA sequence, the TG sequence with wobble base pairing favors Z formation under high salt conditions.

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.

1H NMR study of the solution structure of the self-complementary dodecanucleotide d(TGCA)3

The deoxyoligonucleotide d(TGCA)3 is a candidate for exhibiting unusual conformations. Its 1H NMR spectrum under low salt conditions has been obtained at 400 MHz and assigned using two-dimensional NMR techniques. The sugar puckers and glycosidic torsions have been determined by inspecting the relative intensities of the intranucleotide NOEs and COSY crosspeaks. At low electrolyte concentration (100 mM NaCl) the molecule exists as a right-handed duplex with twelve Watson-Crick base-pairs and deoxyribose moieties assuming the O1'-endo to C1'-exo pucker.

Effects of A:T base pairs on the B-Z conformational transitions of DNA

Effects of A:T base pairs on the propensity of B to Z conformational transitions have been investigated by the CD salt titrations on d(CG)5' d(GC)5' terminal or central A:T replaced decamers, and terminal A:T appended dodecamers. The presence of A:T at the center greatly inhibits the B to Z transition of both G:C decamers. Moderate Z inhibitions are shown by terminal A:T replacements and additions to d(CG)5' with the former exhibiting a stronger effect. In contrast, the addition and replacement with A:T at the terminals of d(GC)5 facilitate the B to Z conversion, with the replacement exhibiting a somewhat more pronounced effect. These results may be rationalized in terms of the number of contigous CG sequences present in an oligomer and the relative inhibitory effects of other dinucleotide sequences. Our results also suggest that some short oligomers with purine at the 5'-end, such as d[A(CG)nT] with n greater than or equal to 2, may likely crystallize as Z conformations.

An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

By means of one- and two-dimensional NMR spectroscopy the solution structures of the partly self-complementary octamer d(m5C-G-m5C-G-A-G-m5C-G) were investigated. It is shown that this DNA fragment, under conditions of high DNA concentration (8 mM DNA) and/or high ionic strength prefers to adopt a duplex structure. At low DNA concentration (0.4 mM DNA), the duplex exists in a 1:1 slow equilibrium with a monomeric hairpin form. Addition of salt destabilizes the hairpin structure in favour of the dimer. At high temperatures the hairpin form, as well as the dimer structure, exist in a fast equilibrium with the random-coil form. For the hairpin/random-coil equilibrium a Tm of 329 K and a delta H degree of -121 kJ.mol-1 were deduced. These thermodynamic parameters are independent of the DNA concentration, as is expected for a monomeric structure. For the dimer to coil transition a Tm of 359 K (1 M DNA) and a delta H degree of -285 kJ.mol duplex-1 were derived. The thermodynamic data of the hairpin-coil transition mutually agree with those recently reported for the hairpin to random coil equilibrium of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) [Orbons, L. P. M., van der Marel, G. A., van Boom, J. H. & Altona, C. (1987) J. Biomol. Struct. Dyns. 4, 939-963]. It is demonstrated that the dimer structure exhibits B-DNA characteristics, as is witnessed by the NOESY experiments and the analysis of the proton-proton coupling data. It is shown that the base-pair formation of the G x A mismatches is anti-anti. A comparison of 1H and 31P chemical-shift data of the title compound with those of a well-characterized B-DNA structure reveals large differences in the dm5C(3)-dG(4)-dA(5) part of the mismatched dimer structure. These differences apparently indicate some major local structural changes due to the incorporation of the G x A mismatches. Under the most extreme conditions used (i.e. up to 3 M NaCl or 75% CH3OH in the presence of 10 mM MgCl2) no Z-DNA structure was observed. It is shown that the structural features of the hairpin form of the title compound mimic those of the hairpin structure of d(m5C-G-m5C-G-T-G-m5C-G). An energy-minimized model of the hairpin form is given.

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.