FTIR and UV spectroscopy of parallel-stranded DNAs with mixed A*T/G*C sequences and their A*T/I*C analogues

The infrared spectra of parallel-stranded (ps) hairpin duplexes with mixed A*T/G*C composition and either isolated or sequential G*C pairs were studied in comparison with antiparallel-stranded (aps) duplexes and a corresponding set of molecules with hypoxanthine as a G base analogue lacking the exocyclic amino group. The ps duplexes showed the characteristic bands for the C2=O2 and C4=O4 stretching vibrations of thymine residues in trans-Watson-Crick A*T pairing at 1683 and 1668 cm-1. The latter band was superimposed on the stretching vibration of the free C6=O6 group of guanine. Substitution of guanine by hypoxanthine inhibited the formation of ps hairpin duplexes whatever the sequence, demonstrating that in the H-bonding between G and C the 2-NH2 group is necessary for stabilizing all of the investigated ps duplexes. This result is in agreement with a model of trans-Watson-Crick G*C base pairs with two H-bonds [N2H2(G)-N3(C) and N1H(G)-O2(C)]. However, trans-Watson-Crick A*T and G*C base pairs with two H-bonds are not isomorphous, which may explain the decreased stability of the ps, but not the aps, duplexes upon increasing the number of A*T/G*C steps. Molecular modeling studies performed on two of the ps duplexes reveal the existence of propeller twist for avoiding a clash between the N2(G) and N4(C) amino groups, and favorable stacking of sequential G*C base pairs. The optimized hairpin ps duplexes invariably incorporated G*C base pairs with two H-bonds, regardless of the initial structures adopted for the force field calculations.     

Thermodynamic characterization of daunomycin-DNA interactions: comparison of complete binding profiles for a series of DNA host duplexes

Using a combination of spectroscopic and calorimetric techniques, we have determined complete thermodynamic binding profiles (delta G degree, delta H degree, and delta S degree) for the complexation of daunomycin to a series of 10 polymeric DNA duplexes. We find the resulting drug binding data to be sensitive to the base composition and sequence of the host duplex, with the binding free energies ranging from -7.5 to -10.8 kcal/mol of bound drug and the binding enthalpies ranging from +4.11 to -10.76 kcal/mol of bound drug at 25 degrees C. The smaller range in the free energy term reflects the impact of large enthalpy-entropy compensations. We observe that the three synthetic duplexes which exhibit the highest daunomycin binding affinities all contain GC (or IC) base pairs as part of alternating purine/pyrimidine sequence motifs, with these high binding affinities being strongly enthalpy driven at 25 degrees C. Specific comparisons between the binding profiles for daunomycin complexation with select pairs of host duplexes lead to the following observations: (1) The presence or absence of a major-groove methyl group does not alter daunomycin binding thermodynamics. (2) The presence or absence of a minor-groove amino group does alter daunomycin binding thermodynamics. (3) Duplexes with different base compositions but identical minor-groove functionality exhibit similar daunomycin binding thermodynamics. (4) Homopolymeric duplexes composed of either AT or AU base pairs, but not GC base pairs, exhibit large enthalpy-entropy compensations in their daunomycin binding profiles. We propose interpretations of these and other features of our thermodynamic data in terms of specific daunomycin-DNA interactions deduced from available structural data.     

G/C-modified oligodeoxynucleotides with selective complementarity: synthesis and hybridization properties

Modified oligodeoxyribonucleotides (ODNs) that have unique hybridization properties were designed and synthesized for the first time. These ODNs, called selective binding complementary ODNs (SBC ODNs), are unable to form stable hybrids with each other, yet are able to form stable, sequence specific hybrids with complementary unmodified strands of nucleic acid. To make SBC ODNs, deoxyguanosine (dG) and deoxycytidine (dC) were substituted with deoxyinosine (dI) and 3-(2'-deoxy-beta-D-ribofuranosyl)pyrrolo-[2,3-d]-pyrimidine-2-(3H)-one (dP), respectively. The hybridization properties of several otherwise identical complementary ODNs containing one or both of these nucleoside analogs were studied by both UV monitored thermal denaturation and non-denaturing PAGE. The data showed that while dI and dP did form base pairs with dC and dG, respectively, dI did not form a stable base pair with dP. A self-complementary ODN uniformly substituted with dI and dP acquired single-stranded character and was able to strand invade the end of a duplex DNA better than an unsubstituted ODN. This observation implies that SBC ODNs should effectively hybridize to hairpins present in single-stranded DNA or RNA.     

Effects of the introduction of L-nucleotides into DNA. Solution structure of the heterochiral duplex d(G-C-G-(L)T-G-C-G).d(C-G-C-A-C-G-C) studied by NMR spectroscopy

The effect of the substitution of a L-nucleoside for a D-nucleoside in the duplex d(G-C-G-T-G-C-G).d(C-G-C-A-C-G-C) was studied by UV and NMR spectroscopy. These unnatural oligonucleotides have potential for antisense DNA technology [Damha, M. J., Giannaris, P. A., & Marfey, P. (1994) Biochemistry (preceding paper in this issue)]. The thermal stability of such duplexes is lower than that of the natural one and is dependent on the nucleotide type and/or sequence. Interestingly, inversion of the chirality of thymidine but not adenosine coincides with a large stabilizing enthalpy change. The structure of the heterochiral duplex d(G1-C2-G3-(L)T4-G5-C6-G7).d(C8-G9-C10-A11-C12-G13- C14), where (L)T denotes the mirror image of the natural thymidine, has been determined by NMR spectroscopy. The sugar conformation was determined using the sum of coupling constants and the distances using a model free relaxation matrix approach. The torsion angles of the backbone follow from 3JHH, 3JHP, and 4JHP coupling constants. The structure of the duplex was calculated by metric matrix distance geometry followed by simulated annealing. The structure is close to that of B-DNA. The base pair formed by (L)T and A is of the Watson-Crick type. All sugars adopt an S-type pucker. The incorporation of the L-sugar in the duplex is accomplished by changes in the backbone torsion angles around the phosphates and the glycosidic torsion angle of (L)T. The modification induces changes in the natural strand as well. The structure exhibits an unusual interaction between the aromatic rings of the (L)T4.A11 and G3.C12 base pairs, which provides a plausible explanation of the unusual thermodynamic properties of the duplex.     

Synthesis and incorporation of tricyclic-imidazo[5',4':4,5]pyrido[2,3-d]pyrimidine nucleosides into oligonucleotides and their thermal stability

We designed and synthesized four novel tricyclicnucleosides 1-4, which have imidazo[5',4':4,5]pyrido[2,3-d]pyrimidine structures, to construct theramallystable, non-natural, DNA duplexes. Both 1:2 and 3:4 possess four complementary hydrogen-bond pairs and are expected to form more stable base-pairs than the naturally occurring Watson-Crick base-pairs when incorporated into oligodeoxynucleotides (ODNs). NMR studies of these base-pairs at the nucleoside level in CDCl3 identified the expected base-pair formation with four hydrogen bonds. We also synthesized several ODNs containing these nucleosides and measured their thermal stability. We found that if several of the novel base-pairs were incorporated on complementary strands, the resulting DNA duplexes had very high thermal stability.     

Refined solution structure of 8,9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 opposite CpA in the complementary strand of an oligodeoxynucleotide duplex as determined by 1H NMR

The solution structure of d(CCATCAFBGATCC).d(GGATCAGATGG), containing the 8,9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 adduct, was refined using molecular dynamics restrained by NOE data obtained from 1H NMR. The modified guanosine was positioned opposite cytosine, while the aflatoxin moiety was positioned opposite adenosine in the complementary strand. Sequential 1H NOEs were interrupted between C5 and AFBG6, but intrastrand NOEs were traced through the aflatoxin moiety, via H6a of aflatoxin and H8 of the modified guanine. Opposite the lesion, the NOE between A16 H1' and G17 H8 was weak. A total of 43 NOEs were observed between DNA protons and aflatoxin protons. Molecular dynamics calculations restrained with 259 experimental and empirical distances, and using sp2 hybridization at AFBG6 N7, refined structures with pairwise rms differences < 0.85 A, excluding terminal base pairs. Relaxation matrix calculations yielded a sixth root rms difference between refined structures and NOE intensity data of 7.3 x 10(-2). The aflatoxin moiety intercalated on the 5'-face of the modified guanine. The extra adenine A16 was inserted between base pair AFBG6.C15 and the aflatoxin moiety. A 36 degree bending between the plane of base pair AFBG6.C15 and the plane of the aflatoxin moiety was predicted. The aflatoxin moiety stacked below the top domain of the oligodeoxynucleotide, which consisted of base pairs C1.G21, C2.G20, A3.T19, T4.A18, and C5.G17. The bottom domain consisted of base pairs AFBG6.C15, A7.T14, T8.A13, C9.G12, and C10.G11. The average winding angle between base pair C5.G17, the intercalated aflatoxin moiety, A16, and base pair AFBG6.C15 was reduced to 10 degrees. The preponderance of base pair substitutions in the aflatoxin B1 mutational spectrum, particularly G-->T transversions, suggests that the stability of this modified oligodeoxynucleotide, which models a templated +1 addition mutation, does not reliably predict the frequency of frame shifts.     

Correlation between sequence-dependent glycosylase repair and the thermal stability of oligonucleotide duplexes containing 1, N6-ethenoadenine

Previous experiments on DNA sequence context reported that base modification, replication, and repair are affected by the nature of neighbor bases. We now report that repair by mammalian alkylpurine-DNA-N-glycosylases (APNG) of 15-mer oligonucleotides with a central 1,N6-ethenoadenine (epsilonA), flanked by 5' and 3' tandem bases, is also highly sequence dependent. Oligonucleotides with the central sequences -GGepsilonAGG- or -CCepsilonACC- are repaired 3-5-fold more efficiently than those containing -AAepsilonAAA- or -TTepsilonATT- when using human or mouse APNG. Melting curves of the same duplexes showed that oligomers with G.C/C. G neighbors were less denatured than those with A.T/T.A neighbors at 37 degreesC. This sequence-dependent difference in denaturation correlates with the relative thermodynamic stability of oligomers with G.C/C.G or A.T/T.A neighbors. The dependence of repair on thermal stability was confirmed by enzyme reactions performed over 0-45 degreesC. Under these conditions, repair of epsilonA flanked by G.C/C.G was dramatically increased at 37 degreesC with continuous increase up to 45 degreesC, in contrast to that with flanking A.T/T. A pairs, which was in agreement with the degree of denaturation of these duplexes. These results indicate that the thermodynamic stability conferred by base pairs flanking epsilonA plays an essential role in maintaining the integrity of the duplex structure which is necessary for repair.     

Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs

Improved thermodynamic parameters for prediction of RNA duplex formation are derived from optical melting studies of 90 oligoribonucleotide duplexes containing only Watson-Crick base pairs. To test end or base composition effects, new sets of duplexes are included that have identical nearest neighbors, but different base compositions and therefore different ends. Duplexes with terminal GC pairs are more stable than duplexes with the same nearest neighbors but terminal AU pairs. Penalizing terminal AU base pairs by 0.45 kcal/mol relative to terminal GC base pairs significantly improves predictions of DeltaG degrees37 from a nearest-neighbor model. A physical model is suggested in which the differential treatment of AU and GC ends accounts for the dependence of the total number of Watson-Crick hydrogen bonds on the base composition of a duplex. On average, the new parameters predict DeltaG degrees37, DeltaH degrees, DeltaS degrees, and TM within 3.2%, 6.0%, 6.8%, and 1.3 degreesC, respectively. These predictions are within the limit of the model, based on experimental results for duplexes predicted to have identical thermodynamic parameters.     

Variability in the stability of DNA-peptide nucleic acid (PNA) single-base mismatched duplexes: real-time hybridization during affinity electrophoresis in PNA-containing gels

The stability of all single-base mismatched pairs between a peptide nucleic acid 11-mer and its complementary DNA has been quantified in terms of their melting temperature and compared with the limited amount of data published to date. The strength of the interaction was determined by an automated affinity-electrophoretic approach permitting the visualization, in real time, of hybridization between a physically immobilized peptide nucleic acid and a complementary DNA migrating in an electric field. The dissociation constants are in the range of 10(-7) M (for mismatches) to 10(-10) M (for fully complementary DNA), which are in excellent agreement with solution studies. These and other thermodynamic constants can be accurately, rapidly, and reproducibly measured in this system at concentrations approaching dissociation conditions by using fluorescently labeled DNA in conjunction with commercial DNA sequencers. The stability of single-base mismatched peptide nucleic acid-DNA duplexes depends both on the position as well as on the chemical nature of the mismatch. The stability is at a minimum when the mutation is positioned 4 bases from either terminus (a loss of 20 degreesC or more in the melting temperature) but regains substantial stability when the mismatch is at the center of the duplex. The most stable mismatched pairs are G:T and T:T whereas destabilization is maximal for A:A and G:G. These observations are of significance in the design of probes for detecting mutations by hybridization.     

Parallel-stranded duplexes and quartet assemblies formed by oligonucleotides containing isoguanine

Oligonucleotides with isoguanine-cytosine base pairs form duplexes with parallel chain orientation. This strong base pair dictates the strand polarity when additional dA-dT or N7Ad-dT pairs are present. The same is observed for the isocytosine-guanine pair. Quartet structures are built by the self-assembly of oligonucleotides containing short runs of isoguanine or 7-deazaisoguanine. New base pairs between isoguanine or guanine and 5-aza-7-deazaguanine are presented.     

Stability of RNA hairpins closed by wobble base pairs

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXANmAYCC, where XY is the wobble base pair, GU or UG, and the underlined loop sequences are three to eight nucleotides. A nearest-neighbor analysis indicates the free energy of loop formation is dependent upon loop size and closing base pair. Hairpin loops closed by UG base pairs are on average 1.3 kcal/mol less stable than hairpins closed by GU base pairs. The hairpin loops closed by UG have approximately the same stability as hairpin loops closed by AU/UA base pairs, while the loops closed by GU are approximately 0.7 kcal/mol more stable than hairpins loops closed by GC/CG base pairs. These results, combined with the model previously developed [Serra et al. (1997) Biochemistry 36, 4844] to predict the stability for hairpin loops closed by Watson-Crick base pairs, allow for the following model to predict the stability of hairpin loops: delta G degree 37L(n) = delta G degree 37iL(n) + delta G degree 37mm + 0.6 (if closed by AU, UA, or UB) - 0.7 (if closed by GU) - 0.7 (if first mismatch is GA or UU except for loops closed by GU). Here, delta G degree 37iL(n) is the free energy increment for initiating a loop of n nucleotides with a CG or GC pair, and delta G degree 37mm is the free energy for the interaction of the first mismatch with the closing base pair. For hairpin loops of n = 4-9, delta G037iL(n) is 4.9, 5.0, 5.0, 5.0, 4.9, and 5.5 kcal/mol, respectively. For hairpin loops of n = 3, delta G degree 37L(3) = +4.8 + 0.6 (if closed by AU, UA, or UG) kcal/mol. Thermodynamic parameters for hairpin formation in 1 M NaCl for 13 naturally occurring RNA hairpin sequences closed by wobble base pairs are reported. The model provides good agreement for both TM and delta G degree 37 for most hairpins studied. Thermodynamic values for five terminal mismatches adjacent to wobble base pairs are also reported.     

Kinetics for hybridization of peptide nucleic acids (PNA) with DNA and RNA studied with the BIAcore technique

The binding of a mixed-sequence pentadecamer PNA (peptide nucleic acid) containing all four nucleobases to the fully complementary as well as various singly mismatched RNA and DNA oligonucleotides has been systematically investigated using thermal denaturation and BIAcore surface-interaction techniques. The rate constants for association (k(a)) and dissociation (k(d)) of the duplex formation as well as the thermal stability (melting temperature, T(m)) of the duplexes have been determined. Upon binding to PNA tethered via a biotin-linker to streptavidin at the dextran/gold surface, DNA and RNA sequences containing single mismatches at various positions in the center resulted in increased dissociation and decreased association rate constants. T(m) values for PNA x RNA duplexes are on average 4 degrees C higher than for PNA x DNA duplexes and follow quantitatively the same variation with mismatches as do the PNA x DNA duplexes. Also a faster k(a) and a slower k(d) are found for PNA x RNA duplexes compared to the PNA x DNA duplexes. An overall fair correlation between T(m), k(a), and k(d) is found for a series of PNA x DNA and PNA x RNA duplexes although the determination of k(a) seemed to be prone to artifacts of the method and was not considered capable of providing absolute values representing the association rate constant in bulk solution.     

The impact of an exocyclic cytosine adduct on DNA duplex properties: significant thermodynamic consequences despite modest lesion-induced structural alterations

The exocyclic base adduct 3,N4-deoxyethenocytosine (epsilonC) is a common DNA lesion that can arise from carcinogen exposure and/or as a biproduct of cellular processes. We have examined the thermal and thermodynamic impact of this lesion on DNA duplex properties, as well as the structural alterations imparted by the lesion. For these studies, we used calorimetric and spectroscopic techniques to investigate a family of 13-mer DNA duplexes of the form (5'CGCATGNGTACGC3')x(3'GCGTACNCATGCG5'), where the central NxN base pair represents the four standard Watson-Crick base pairs (corresponding to four control duplexes), and where either one of the N bases has been replaced by epsilonC, yielding eight test duplexes. Studies on these 12 duplexes permit us to assess the impact of the epsilonC lesion as a function of sequence context. Our spectroscopic and calorimetric data allow us to reach the following conclusions: (i) The epsilonC lesion imparts a large penalty on duplex stability, with sequence context only modestly modulating the extent of this lesion-induced destabilization. This result contrasts with our recent studies of duplexes with abasic sites, where sequence context was found to be the predominant determinant of thermodynamic damage. (ii) For the epsilonC-containing duplexes, sequence context effects are most often observed in the enthalpic contribution to lesion-induced duplex destabilization. However, due to compensating entropies, the free energy changes associated with this lesion-induced duplex destablization are nearly independent of sequence context. (iii) Despite significant lesion-induced changes in duplex energetics, our spectroscopic probes detect only modest lesion-induced changes in duplex structure. In fact, the overall duplex maintains a global B-form conformation, in agreement with NMR structural data. We discuss possible interpretations of the apparent disparity between the severe thermodynamic and relatively mild structural impacts of the epsilonC lesion on duplex properties. We also note and discuss the implications of empirical correlations between biophysical and biological properties of lesion-containing duplexes.     

Watson-Crick base-pairing properties of bicyclo-DNA

A series of sequences of the DNA analog bicyclo-DNA, 6-12 nucleotides in length and containing all four natural nucleobases, were prepared and their Watson-Crick pairing properties with complementary RNA and DNA, as well as in its own series, were analyzed by UV-melting curves and CD-spectroscopy. The results can be summarized as follows: bicyclo-DNA forms stable Watson-Crick duplexes with complementary RNA and DNA, the duplexes with RNA generally being more stable than those with DNA. Pyrimidine-rich bicyclo-DNA sequences form duplexes of equal or slightly increased stability with DNA or RNA, whereas purine-rich sequences show decreased affinity to complementary DNA and RNA when compared with wild-type (DNA-DNA, DNA-RNA) duplexes. In its own system, bicyclo-DNA prefers antiparallel strand alignment and strongly discriminates for base mismatches. Duplexes are always inferior in stability compared with the natural ones. A detailed analysis of the thermodynamic properties was performed with the sequence 5'-GGATGGGAG-3'x 5'-CTCCCATCC-3' in both backbone systems. Comparison of the pairing enthalpy and entropy terms shows an enthalpic advantage for DNA association (delta deltaH = -18 kcal x (mol)-1)) and an entropic advantage for bicyclo-DNA association (delta deltaS = 49 cal x K(-1) x mol(-1), leading to a delta deltaG 25 degrees C of -3.4 kcal x mol(-1) in favor of the natural duplex. The salt dependence of Tm for this sequence is more pronounced in the case of bicyclo-DNA due to increased counter ion screening from the solvent. Furthermore bicyclo-DNA sequences are more stable towards snake venom phosphodiesterase by a factor of 10-20, and show increased stability in fetal calf serum by a factor of 8 compared with DNA.     

Crystal structure of an RNA octamer duplex r(CCCIUGGG)2 incorporating tandem I.U wobbles

The crystal structure of the RNA octamer duplex r(CCCIUGGG)2has been elucidated at 2.5 A resolution. The crystals belong to the space group P21and have unit cell constants a = 33.44 A, b = 43.41 A, c = 49.39 A and beta = 104.7 degrees with three independent duplexes (duplexes 1-3) in the asymmetric unit. The structure was solved by the molecular replacement method and refined to an Rwork/Rfree of 0.185/0.243 using 3765 reflections between 8.0 and 2.5 A. This is the first report of an RNA crystal structure incorporating I.U wobbles and three molecules in the asymmetric unit. Duplex 1 displays a kink of 24 degrees between the mismatch sites, while duplexes 2 and 3 have two kinks each of 19 degrees and 27 degrees, and 24 degrees and 29 degrees, respectively, on either side of the tandem mismatches. At the I.U/U.I mismatch steps, duplex 1 has a twist angle of 33.9 degrees, close to the average for all base pair steps, but duplexes 2 and 3 are underwound, with twist angles of 24.4 degrees and 26.5 degrees, respectively. The tandem I.U wobbles show intrastrand purine-pyrimidine stacking but exhibit interstrand purine-purine stacking with the flanking C.G pairs. The three independent duplexes are stacked non-coaxially in a head-to-tail fashion to form infinite pseudo-continuous helical columns which form intercolumn hydrogen bonding interactions through the 2'-hydroxyl groups where the minor grooves come together.     

N 2-methylguanosine is iso-energetic with guanosine in RNA duplexes and GNRA tetraloops

Modified nucleotides are resource-intensive alternatives to the four nucleotides that constitute the bulk of natural RNAs. Yet, even in cases where modifications are highly conserved, their functions are difficult to identify. One possible function might be to modulate the stability of RNA structures. To investigate this possibility for N 2-methylguanosine (m2G), which is present in a wide variety of RNAs, we have determined the thermodynamic consequences of substituting m2G for G in G-C Watson-Crick pairs and G@U wobble pairs within RNA duplexes. The m2G substitution is iso-energetic with G in all cases, except for aninternal m2G@U pair, where it has a modest (0.3 kcal/mol) stabilizing effect. We have also examined theconsequences of replacing G by m2G, and A by N 6, N 6-dimethyladenosine (m26A) in the helix 45 tetraloop of 16S rRNA, which would otherwise be a standard GNRA tetraloop. This loop is a conserved, hypermethylated region of the ribosome where methylation appears to modulate activity. m26A substitution destabilizes the tetraloop, presumably because it prevents the formation of the G@A sheared pair it would otherwise contain. m2G substitution has no effect on tetraloop stability. Together, these results suggest that m2G is equally stable as either the s-cis or s-trans rotamer. The lack of a significant effect on secondary structural stability in these systems suggests that m2G is introduced into naturally occurring RNAs for reasons other than modulation of duplex stability.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

Age dependence of total cerebral blood flow volume from childhood to adulthood

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp3 hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G-A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)...A(anti) pair of the ApG step of the d(CCAAGATTGG)2 crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.     

Solution structure of a DNA duplex containing the exocyclic lesion 3,N4-etheno-2'-deoxycytidine opposite 2'-deoxyguanosine

Vinyl chloride reacts with cellular DNA producing 3,N4-etheno-2'-deoxycytidine (epsilonC) along with other exocyclic adducts. The solution structure of an oligodeoxynucleotide duplex containing an epsilonC.dG base pair was determined by high-resolution NMR spectroscopy and molecular dynamics simulations. NMR data indicated that the duplex adopts a right-handed helical structure having all residues in anti orientation around the glycosylic torsion angle. The epsilonC adduct has a sugar pucker in the C3'-endo/C4'-exo region while the rest of the residues are in the C2'-endo/C3'-exo range. NOE interactions established Watson-Crick alignments for canonical base pairs of the duplex. The imino proton of the lesion-containing base pair resonated as a sharp signal that was resistant to water exchange, suggesting hydrogen bonding. Restrained molecular dynamics simulations generated three-dimensional models in excellent agreement with the spectroscopic data. The refined structures are slightly bent at the lesion site without major perturbations of the sugar-phosphate backbone. The adduct is displaced and shifted toward the major groove of the helix while its partner on the complementary strand remains stacked. The epsilonC(anti).dG(anti) base pair alignment is sheared and stabilized by the formation of hydrogen bonds. The biological implications of structures of epsilonC-containing DNA duplexes are discussed.     

Determination of salicin and related compounds in botanical dietary supplements by liquid chromatography with fluorescence detection

The centromeric dodeca-satellite of Drosophila forms unusual DNA structures in which its purine-rich strand (GTACGGGACCGA)n folds into very stable intramolecular hairpins. These intramolecular hairpins contain groups of tandem 5'-GA:GA-3' mismatches that, as judged by gel electrophoresis analysis and UV-melting studies, have a determinant contribution to their stability. Duplexes of the dodeca-satellite purine-rich strand, carrying tandem 5'-GA:GA-3' mismatches, are as stable as equivalent fully Watson-Crick duplexes containing tandem 5'-TA:TA-3' Watson-Crick pairs in place of the non-Watson-Crick G.A pairs. On the other hand, duplexes carrying any of the other three possible tandem combinations of purine.purine mismatches, including G.A pairs on the opposite orientation 5'-AG:AG-3', are very unstable. The high stability of the dodeca-satellite hairplus suggests that the tandem G.A pairs are on the sheared configuration although they are found within the less favourable 5'-G-(G-A)-C-3' sequence context. Other centromeres DNA sequences, including the AAGAG satellite of Drosophila and the mammalian CENP-B box sequence, have the potential of forming intramolecular hairpins stabilised by similar purine.purine interactions.