Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

Thermodynamic characterization of single mismatches found in naturally occurring RNA

Many naturally occurring RNA structures contain single mismatches. However, the algorithms currently used to predict RNA structure from sequence rely on a minimal set of data for single mismatches, most of which occur rather infrequently in nature. As a result, several approximations and assumptions are used to predict the stability of RNA duplexes containing the most common single mismatches. Therefore, the relative frequency of single mismatches was determined by compiling and searching a database of 955 RNA secondary structures. Thermodynamic parameters for duplex formation, derived from optical melting experiments, are reported for 28 oligoribonucleotides containing frequently occurring single mismatches. These data were then combined with previous data to construct a dataset of 64 single mismatches, including the 30 most common in the database. Because of this increase in experimental thermodynamic parameters for single mismatches that occur frequently in nature, more accurate free energy calculations have resulted. To improve the prediction of the thermodynamic parameters for duplexes containing single mismatches that have not been experimentally measured, single mismatch-specific nearest neighbor parameters were derived. The free energy of an RNA duplex containing a single mismatch that has not been thermodynamically characterized can be calculated by: DeltaG degrees 37,single mismatch = DeltaG degrees 37,mismatch nt + DeltaG degrees 37,mismatch-NN interaction + DeltaG degrees 37,AU/GU. Here, DeltaG degrees 37,mismatch is -0.4, -2.1, and -0.3 kcal/mol for A.G, G.G, and U.U mismatches, respectively; DeltaG degrees 37,mismatch-NN interaction is 0.7, -0.5, 0.4, -0.4, and -1.0 kcal/mol for 5'YRR3'/3'RRY5', 5'RYY3'/3'YYR5', 5'YYR3'/3'RYY5', 5'YRY3'/3'RYR5', and 5'RRY3'/3'YYR5' mismatch-nearest neighbor combinations, respectively, when A and G are categorized as purines (R) and C and U are categorized as pyrimidines (Y); and DeltaG degrees 37,AU/GU is a penalty of 1.2 kcal/mol for replacing a G-C base pair with either an A-U or G-U base pair. Similar predictive models were also derived for DeltaH degrees single mismatch and DeltaS degrees single mismatch. These new predictive models, in conjunction with the reported thermodynamics for frequently occurring single mismatches, should allow for more accurate calculations of the free energy of RNA duplexes containing single mismatches and, furthermore, allow for improved prediction of secondary structure from sequence.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Base pairing geometry in GA mismatches depends entirely on the neighboring sequence.

We have synthesized nine self-complementary DNA oligomers containing different flanking sequences adjacent to a pair of contiguous GA mismatches, and have used high resolution nuclear magnetic resonance (n.m.r.) to investigate the GpA phosphodiester backbone conformation and mismatch pairing schemes in these duplexes. We found dramatic effects of the flanking base pair on the hydrogen bonding and backbone conformation, which appear to be coupled. Thus the Ganti-Aanti base pairing scheme in a NAGATN sequence switches to a more stable sheared GA base pairing scheme in a NCGAGN or NTGAAN context, while no duplex is formed (or only GA bulges occur) when NAGATN is changed to NGGACN. Furthermore, the more stable sheared GA pairing in NPyGAPuN sequences is associated with a BII rather than BI backbone conformation for the phosphodiester between the adjacent mismatched GA pairs. The overall stability of these adjacent GA mismatches as measured by imino proton n.m.r. studies is Py-GA-Pu > A-GA-T > G-GA-C.     

Base pairing involving deoxyinosine: implications for probe design.

The thermal stability of oligodeoxyribonucleotide duplexes containing deoxyinosine (I) residues matched with each of the four normal DNA bases were determined by optical melting techniques. The duplexes containing at least one I were obtained by mixing equimolar amounts of an oligonucleotide of sequence dCA3XA3G with one of sequence dCT3YT3G where X and Y were A, C, G, T, or I. Comparison of optical melting curves yielded relative stabilities for the I-containing standard base pairs in an otherwise identical base-pair sequence. I:C pairs were found to be less stable than A:T pairs in these duplexes. Large neighboring-base effects upon stability were observed. For example, when (X,Y) = (I,A), the duplex is eight-fold more stable than when (X,Y) = (A,I). Independent of sequence effects the order of stabilities is: I:C greater than I:A greater than I:T congruent to I:G. This order differs from that of deoxyguanosine which pairs less strongly with dA; otherwise each deoxyinosine base pair is less stable than its deoxyguanosine counterpart in the same sequence environment. Implications of these results for design of DNA oligonucleotide probes are discussed.     

NMR studies of G:A mismatches in oligodeoxyribonucleotide duplexes modelled after ribozymes.

The structures of two oligodeoxyribonucleotide duplexes, the base sequences of which were modelled after both a hammerhead ribozyme and a small metalloribozyme, were studied by NMR. Both duplexes contain adjacent G:A mismatches; one has a PyGAPu:PyGAPu sequence and the other a PyGAPy:PuGAPu sequence. It is concluded on the basis of many characteristic NOEs that in both duplexes G:A base pairs are formed in the unique 'sheared' form, where an amino proton instead of an imino proton of G is involved in the hydrogen bonding, and G and A bases are arranged 'side by side' instead of 'head to head'. A photo-CIDNP experiment, which gives unique and independent information on the solvent accessibility of nucleotide bases, also supports G:A base pairing rather than a bulged-out structure of G and A residues. This is the first demonstration that not only the PyGAPu:PyGAPu sequence but also the PyGAPy:PuGAPu sequence can form the unique sheared G:A base pairs. Taking the previous studies on G:A mismatches into account, the idea is suggested that a PyGA:GAPu sequence is a minimum and essential element for the formation of the sheared G:A base pairs. The sheared G:A base pairs in the PyGAPu:PyGAPu sequence are suggested to be more stable than those in the PyGAPy:PuGAPu sequence. This is explained rationally by the idea proposed above.     

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.     

Base-base mismatches. Thermodynamics of double helix formation for dCA3XA3G + dCT3YT3G (X, Y = A,C,G,T).

Thermodynamic parameters for double strand formation have been measured for the sixteen double helices of the sequence dCA3XA3G.dCT3YT3G, with each of the bases A, C, G and T at the positions labelled X and Y. The results are analyzed in terms of nearest-neighbors and are compared with thermodynamic parameters for RNA secondary structure. At room temperature the sequence (Formula: see text) is more stable than (Formula: see text) and is similar in stability to (Formula: see text) and (Formula: see text) are least stable. At higher temperatures the sequences containing a G.C base pair become more stable than those containing only A.T. All molecules containing mismatches are destabilized with respect to those with only Watson-Crick pairing, but there is a wide range of destabilization. At room temperature the most stable mismatches are those containing guanine (G.T, G.G, G.A); the least stable contain cytosine (C.A, C.C). At higher temperatures pyrimidine-pyrimidine mismatches become the least stable.     

Synthesis of oligodeoxyribonucleotides containing degenerate bases and their use as primers in the polymerase chain reaction.

Heptadecaoligodeoxyribonucleotides containing one or more of the bases, 6H,8H-3,4-dihydropyrimido[4,5-c][1,2]oxazin-7-one (P), 2-amino-6-methoxyaminopurine (K), and hypoxanthine (I) and combinations of P with K and I have been synthesised on a DNA synthesiser. The stability of duplexes containing these basemodified oligomers with P/A, P/G, K/C and K/T; P/A, P/G, I/C, I/T and I/A, I/G, I/C, I/T base pairs were compared by measuring their melting transition (Tm) values. Oligomers containing both P and K and P and I were more stable than those with I alone or with mismatches. These oligomers together with one with a P base at the 3'-end were used as primers in polymerase chain reaction (PCR) experiments. They were all effective primers except one with I alone and a triple mismatch. Thus the use of the degenerate bases P and K in primer design is established.     

Solution structure of the ActD-5'-CCGTT3GTGG-3' complex: drug interaction with tandem G.T mismatches and hairpin loop backbone

Binding of actinomycin D (ActD) to the seemingly single-stranded DNA (ssDNA) oligomer 5'-CCGTT3 GTGG-3' has been studied in solution using high-resolution nuclear magnetic resonance (NMR) techniques. A strong binding constant (8 x 10(6) M(-1)) and high quality NMR spectra have allowed us to determine the initial DNA structure using distance geometry as well as the final ActD-5'-CCGTT3 GTGG-3' complex structure using constrained molecular dynamics calculations. The DNA oligomer 5'-CCGTT3GTGG-3' in the complex forms a hairpin structure with tandem G.T mismatches at the stem region next to a loop of three stacked thymine bases pointing toward the major groove. Bipartite T2O-GH1 and T2O-G2NH2 hydrogen bonds were detected for the G.T mismatches that further stabilize this unusual DNA hairpin. The phenoxazone chromophore of ActD intercalates nicely between the tandem G.T mismatches in essentially one major orientation. Additional hydrophobic interactions between the ActD quinoid amino acid residues with the loop T5-T6-T7 backbone protons were also observed. The hydrophobic G-phenoxazone-G interaction in the ActD-5'-CCGTT3GTGG-3' complex is more robust than that of the classical ActD- 5'-CCGCT3GCGG-3' complex, consistent with the roughly 2-fold stronger binding of ActD to the 5'-CCGTT3GTGG-3' sequence than to its 5'-CCG CT3GCGG-3' counterpart. Stabilization by ActD of a hairpin containing non-canonical stem base pairs further strengthens the notion that ActD or other related compounds may serve as a sequence- specific ssDNA-binding agent that inhibits human immunodeficiency virus (HIV) and other retroviruses replicating through ssDNA intermediates.     

Crystal structure of a 14 bp RNA duplex with non-symmetrical tandem GxU wobble base pairs

Adjacent GxU wobble base pairs are frequently found in rRNA. Atomic structures of small RNA motifs help to provide a better understanding of the effects of various tandem mismatches on duplex structure and stability, thereby providing better rules for RNA structure prediction and validation. The crystal structure of an RNA duplex containing the sequence r(GGUAUUGC-GGUACC)2 has been solved at 2.1 A resolution using experimental phases. Novel refinement strategies were needed for building the correct solvent model. At present, this is the only short RNA duplex structure containing 5'-U-U-3'/3'-G-G-5' non-symmetric tandem GxU wobble base pairs. In the 14mer duplex, the six central base pairs are all displaced away from the helix axis, yielding significant changes in local backbone conformation, helix parameters and charge distribution that may provide specific recognition sites for biologically relevant ligand binding. The greatest deviations from A-form helix occur where the guanine of a wobble base pair stacks over a purine from the opposite strand. In this vicinity, the intra-strand phosphate distances increase significantly, and the major groove width increases up to 3 A. Structural comparisons with other short duplexes containing symmetrical tandem GxU or GxT wobble base pairs show that nearest-neighbor sequence dependencies govern helical twist and the occurrence of cross-strand purine stacks.     

Influence of pH on the conformation and stability of mismatch base-pairs in DNA

A series of self-complementary dodecanucleotide duplexes containing two symmetrically disposed mismatches have been studied by pH-dependent, ultraviolet light melting techniques. The results indicate that A.C, and C.C mismatches are strongly stabilized by protonation and that the degree of stabilization of the A.C mismatch depends greatly on the flanking bases. In one case, a duplex containing two A.C mismatches is more stable than the native sequence below pH 5.5. The G.A mismatch displays conformational flexibility, with a protonated G(syn).A(anti) base-pair occurring in certain base stacking environments but not in others. The A.A and T.C mismatches are not stabilized at low pH. These solution studies correlate well with predictions based on X-ray crystallographic data.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.     

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.     

Hydration pattern of A4T4 and T4A4 DNA: a molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.     

Synthesis of 2'-O-methyl-RNAs incorporating a 3-deazaguanine, and UV melting and computational studies on its hybridization properties

2'-O-methyl-RNAs incorporating 3-deazaguanine (c3G) were synthesized by use of N,N-diphenylcarbamoyl and N,N-dimethylaminomethylene as its base protecting groups to suppress sheared-type 5'-GA-3'/5'-GA-3' tandem mismatched base pairing which requires the N3 atom. These modified RNAs hybridized more weakly with the complementary and single mismatch-containing RNAs than the unmodified RNAs. The T(m) experiments were performed to clarify the effects of replacement of the fifth G with c(3)G on stabilization of 2'-O-methyl-(5'-CGGCGAGGAG-3')/5'-CUCCGAGCCG-3' and 2'-O-methyl-(5'-CGGGGACGAG-3')/5'-CUCGGACCCG-3'duplexes, which form sheared-type and face-to-face type 5'-GA-3'/5'-GA-3' tandem mismatched base pairs, respectively. Consequently, this replacement led to more pronounced destabilization of the former duplex that needs the N3 atom for the sheared-type base pair than the latter that does not need it for the face-to-face type base pair. A similar tendency was observed for 2'-O-methyl-RNA/DNA duplexes. These results suggest that the N3 atom of G plays an important role in stabilization of the canonical G/C base pair as well as the base discrimination and its loss suppressed formation of the undesired sheared-type mismatched base pair. Computational studies based on ab initio calculations suggest that the weaker hydrogen bonding ability and larger dipole moment of c3G can be the origin of the lower T(m).     

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.     

Physicochemical properties of phosphorothioate oligodeoxynucleotides.

We have recently shown that phosphorothioate (PS) oligodeoxynucleotide (ODN) analogs, unlike their normal congeners, exhibit significant anti-HIV activity (Matsukura et al., (1987) Proc. Natl. Acad. Sci. USA 84, 7706-7710). We now report the syntheses, melting temperatures (Tm), and nuclease susceptibilities of a series of phosphorothioate ODN analogs. These include all-PS duplexes, duplexes with one normal chain and the other chain either all-PS, or end-capped with several PS groups at both 3' and 5' ends. The DNase susceptibilities of the S-ODNs are much less than the normal phosphodiesters, but by contrast duplexes of poly-rA with S-dT40 are much more susceptible to RNase H digestion. The Tm's for AT base pairs of S-ODNs are significantly depressed relative to normals, while GC base pairs show much less Tm depression. The Tm's of S-dT oligomers with poly-rA are reduced relative to the duplexes with normal dA oligomers. These results have significance for the biological properties of these analogs as anti-message inhibitors of gene expression, and provide a rational basis for the S-dC/G sequences as potential effective anti-AIDS agents.     

Crystal structures of two plasmid copy control related RNA duplexes: An 18 base pair duplex at 1.20 A resolution and a 19 base pair duplex at 1.55 A resolution.

The structures of two RNA duplexes, whose sequences correspond to portions of the ColE1 plasmid copy control RNA I and RNA II, have been determined. Crystals containing the 18mers 5'-CA CCGUUGGUAGCGGUGC-3' and 5'-CACCGCUACCAACGGUGC-3' diffract to 1.20 A resolution while those containing the 19mers 5'-GCACCGUUGGUAGCGGUGC-3' and 5'-GCACCGCUACCAACGGUGC-3' diffract to 1.55 A resolution. Both duplexes are standard A form, with Watson-Crick base pairing throughout. Use of anisotropic atomic displacement factors in refinement of the 1.20 A structure dramatically improved refinement statistics, resulting in a final R(free) of 15.0% and a crystallographic R-factor of 11.6%. Perhaps surprisingly, these crystals of the 18 base pair RNA exhibit a 36-fold static disorder, resulting in a structure with a single sugar-phosphate backbone conformation and an averaged base composition at each residue. Since the sugar-phosphate backbone structure is identical in the 36 different nucleotides that are superimposed, there can be no sequence-dependent variation in the structure. The average ribose pucker amplitude is 45.8 degrees for the 18 base pair structure and 46.4 degrees for the 19 base pair structure; these values are respectively 19% and 20% larger than the average pucker amplitude reported from nucleoside crystal structures. A standard RNA water structure, based on analysis of the hydration of these crystal structures and that of the TAR RNA stem [Ippolito, J. A., and Steitz, T. A. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 9819-9824], has been derived, which has allowed us to predict water positions in lower resolution RNA crystal structures. We report a new RNA packing motif, in which three pro-S(p) phosphate oxygens interact with an ammonium ion.