Thermodynamic stability and solution conformation of tandem G.A mismatches in RNA and RNA.DNA hybrid duplexes

Multiple approaches for the investigation of bacterial small regulatory RNAs self-assembly

RNAs are flexible molecules involved in a multitude of roles in the cell. Specifically, noncoding RNAs (i.e., RNAs that do not encode a protein) have important functions in the regulation of biological processes such as RNA decay, translation, or protein translocation. In bacteria, most of those noncoding RNAs have been shown to be critical for posttranscriptional control through their binding to the untranslated regions of target mRNAs. Recent evidence shows that some of these noncoding RNAs have the propensity to self-assemble in prokaryotes. Although the function of this self-assembly is not known and may vary from one RNA to another, it offers new insights into riboregulation pathways. We present here the various approaches that can be used for the detection and analysis of bacterial small noncoding RNA self-assemblies.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

Understanding the role of base stacking in nucleic acids. MD and QM analysis of tandem GA base pairs in RNA duplexes

Preceding NMR experiments show that the conformation of tandem GA base pairs, an important recurrent non-canonical building block in RNA duplexes, is context dependent. The GA base pairs adopt "sheared" N3(G)-N6(A), N2(G)-N7(A) geometry in the r(CGAG)(2) and r(iGGAiC)(2) contexts while switching to "imino" N1(G)-N1(A), O6(G)-N6(A) geometry in the r(GGAC)(2) and r(iCGAiG)(2) contexts (iC and iG stand for isocytosine and isoguanine, respectively). As base stacking is likely to be one of the key sources of the context dependence of the conformation of GA base pairs, we calculated base stacking energies in duplexes containing such base pairs, to see if this dependence can be predicted by stacking energy calculations. When investigating the context dependence of the GA geometry two different conformations of the same duplex were compared (imino vs. sheared). The geometries were generated via explicit solvent MD simulations of the respective RNA duplexes, while the subsequent QM energy calculations focused on base stacking interactions of the four internal base pairs. Geometrical relaxation of nucleobase atoms prior to the stacking energy computations has a non-negligible effect on the results. The stacking energies were derived at the DFT-D/6-311++G(3df,3pd) level. We show a rather good correspondence between the intrinsic gas-phase stacking energies and the NMR-determined GA geometries. The conformation with more favorable gas-phase stacking is in most cases the one observed in experiments. This correlation is not improved when including solvent effects via the COSMO method. On the other side, the stacking calculations do not predict the relative thermodynamic stability of duplex formation for different sequences.

Probing hybridization parameters from microarray experiments: nearest-neighbor model and beyond

In this article, it is shown how optimized and dedicated microarray experiments can be used to study the thermodynamics of DNA hybridization for a large number of different conformations in a highly parallel fashion. In particular, free energy penalties for mismatches are obtained in two independent ways and are shown to be correlated with values from melting experiments in solution reported in the literature. The additivity principle, which is at the basis of the nearest-neighbor model, and according to which the penalty for two isolated mismatches is equal to the sum of the independent penalties, is thoroughly tested. Additivity is shown to break down for a mismatch distance below 5 nt. The behavior of mismatches in the vicinity of the helix edges, and the behavior of tandem mismatches are also investigated. Finally, some thermodynamic outlying sequences are observed and highlighted. These sequences contain combinations of GA mismatches. The analysis of the microarray data reported in this article provides new insights on the DNA hybridization parameters and can help to increase the accuracy of hybridization-based technologies.

Secondary structure and stability of the selenocysteine insertion sequences (SECIS) for human thioredoxin reductase and glutathione peroxidase

We have used high resolution NMR and thermodynamics to characterize the secondary structure and stability of the selenocysteine insertion sequences (SECIS) of human glutathione peroxidase (58 nt) and thioredoxin reductase (51 nt). These sequences are members of the two classes of SECIS recently identified with two distinct structures capable of directing selenocysteine incorporation into proteins in eukaryotes. UV melting experiments showed a single cooperative and reversible transition for each RNA, which indicates the presence of stable secondary structures. Despite their large size, the RNAs gave well resolved NMR spectra for the exchangeable protons. Using NOESY, the imino protons as well as the cytosine amino protons of all of the Watson-Crick base pairs were assigned. In addition, a number of non-canonical base pairs including the wobble G.U pairs were identified. The interbase-pair NOEs allowed definition of the hydrogen-bonded structure of the oligonucleotides, providing an experimental model of the secondary structure of these elements. The derived secondary structures are consistent with several features of the predicted models, but with some important differences, especially regarding the conserved sequence motifs.

The solution structure of [d(CGC)r(aaa)d(TTTGCG)](2): hybrid junctions flanked by DNA duplexes

The solution structure and hydration of the chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2), in which the central hybrid segment is flanked by DNA duplexes at both ends, was determined using two-dimensional NMR, simulated annealing and restrained molecular dynamics. The solution structure of this chimeric duplex differs from the previously determined X-ray structure of the analogous B-DNA duplex [d(CGCAAATTTGCG)](2)as well as NMR structure of the analogous A-RNA duplex [r(cgcaaauuugcg)](2). Long-lived water molecules with correlation time tau(c)longer than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not with the other two thymines (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation, while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations. The exchange rates for RNA C2'-OH were found to be approximately 5-20 s(-1). This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2)is wider than its B-DNA analog but narrower than that of the A-RNA analog. It was further confirmed by its titration with the minor groove binding drug distamycin. A possible 2:1 binding mode was found by the titration experiments, suggesting that this chimeric duplex contains a wider minor groove than its B-DNA analog but still narrow enough to hold two distamycin molecules. These distinct structural features and hydration patterns of this chimeric duplex provide a molecular basis for further understanding the structure and recognition of DNA. RNA hybrid and chimeric duplexes.

The 8-nucleotide-long RNA:DNA hybrid is a primary stability determinant of the RNA polymerase II elongation complex

The sliding clamp model of transcription processivity, based on extensive studies of Escherichia coli RNA polymerase, suggests that formation of a stable elongation complex requires two distinct nucleic acid components: an 8-9-nt transcript-template hybrid, and a DNA duplex immediately downstream from the hybrid. Here, we address the minimal composition of the processive elongation complex in the eukaryotes by developing a method for promoter-independent assembly of functional elongation complex of S. cerevisiae RNA polymerase II from synthetic DNA and RNA oligonucleotides. We show that only one of the nucleic acid components, the 8-nt RNA:DNA hybrid, is necessary for the formation of a stable elongation complex with RNA polymerase II. The double-strand DNA upstream and downstream of the hybrid does not affect stability of the elongation complex. This finding reveals a significant difference in processivity determinants of RNA polymerase II and E. coli RNA polymerase. In addition, using the imperfect RNA:DNA hybrid disturbed by the mismatches in the RNA, we show that nontemplate DNA strand may reduce the elongation complex stability via the reduction of the RNA:DNA hybrid length. The structure of a "minimal stable" elongation complex suggests a key role of the RNA:DNA hybrid in RNA polymerase II processivity.

Structure and recognition of sheared tandem G x A base pairs associated with human centromere DNA sequence at atomic resolution

G x A mismatched base pairs are frequently found in nucleic acids. Human centromere DNA sequences contain unusual repeating motifs, e.g. , (GAATG)n x (CATTC)n found in the human chromosome. The purine-rich strand of this repeating pentamer sequence forms duplex and hairpin structures with unusual stability. The high stability of these structures is contributed by the "sheared" G x A base pairs which present a novel recognition surface for ligands and proteins. We have solved the crystal structure, by the multiple-wavelength anomalous diffraction (MAD) method of d(CCGAATGAGG) in which the centromere core sequence motif GAATG is embedded. Three crystal forms were refined to near-atomic resolution. The structures reveal the detailed conformation of tandem G x A base pairs whose unique hydrogen-bonding surface has interesting interactions with bases, hydrated magnesium ions, cobalt(III)hexaammine, spermine, and water molecules. The results are relevant in understanding the structure associated with human centromere sequence in particular and G x A base pairs in nucleic acids (including RNA, like ribozyme) in general.

Secondary structure dimorphism and interconversion between hairpin and duplex form of oligoribonucleotides

RNA hairpins can alternatively form a dimer with a bulged loop flanked by regularly base paired regions. [1H]NMR spectroscopy and native gel electrophoresis were used to study how the sequence of nucleotides in the loop of the hairpin affect the hairpin-duplex interconversion. As a model system, a hairpin containing 7 nucleotides in the loop and 5 base pairs in the stem was used. The loop size was gradually reduced from 7 to 4 nucleotides, yielding finally the stable UNCG tetraloop. Single nucleotide mutations were performed to investigate the influence of the self-complementarity of the loop sequence on the dimerization. The results demonstrate that (1) the initial fraction of hairpin is determined by concentration of the oligonucleotide, the annealing procedure, and the relative stability of the loop, (2) the degree of self-complementarity of the loop sequence of the hairpin governs the dimerization kinetics, and (3) oligonucleotides complementary to the loop sequence decrease the dimerization rate. We propose a secondary structure-based model for the dimerization reaction of RNA hairpins in which the formation of intermolecular base pairs between self-complementary nucleotides of the loops represents the nucleation step.

The detailed structure of tandem G.A mismatched base-pair motifs in RNA duplexes is context dependent

The solution structure of the RNA duplex (rGGGCUGAAGCCCU), containing tandem G.A mismatches has been determined by NMR spectroscopy and restrained molecular dynamics. A homonuclear 3D TOCSY-NOESY was used to derive 18 to 30 distance restraints per nucleotide, as well as all gamma torsion angles and sugar puckers for the central UGAA part of the molecule. Using these constraints, together with cross-strand distances, involving exchangeable imino protons, and essentially all other torsion angles that can accurately be determined (i.e. beta, epsilon) otherwise, the structure of the UGAA domain could be determined with high precision (r.m.s.d. 0.62 A), without the aid of isotopically enriched RNA. The G.A base-pairs are of the sheared pairing type, with both nucleotides in the anti conformation, and hydrogen bonds between the guanine 2-amino and the adenine N7 and between the guanine N3 and the adenine 6-amino. Surprisingly the sugar of the guanosine of the G.A. mismatch adopts a 2'-endo sugar pucker conformation. Comparison with other RNA structures, in which two such G.A base-pairs are formed reveals that this detailed structure depends on the identity of the base 5' to the guanosine in the tandem G.A base-pairs. A geometrical model for the incorporation of sheared tandem G.A base-pairs in A-form helices is formulated, which explains the distinct different stacking properties and helical parameters in sequences containing tandem, sheared G.A base-pairs.

Conformational properties and thermodynamics of the RNA duplex r(CGCAAAUUUGCG)2: comparison with the DNA analogue d(CGCAAATTTGCG)2

The thermodynamic stability of nine dodecamers (four DNA and five RNA) of the same base composition has been compared by UV-melting. TheDeltaG of stabilisation were in the order: r(GACUGAUCAGUC)2>r(CGCAAATTTGCG)2 approximately r(CGCAUAUAUGCG)2>d(CGCAAATTTGCG)2 approximately r(CGCAAAUUUGCG)2>d(CGCATATATGCG)2 approximately d(GACTGATCAGTC)2>r(CGCUUUAAAGCG)2 approximately d(CGCTTTAAAGCG)2. Compared with the mixed sequences, both r(AAAUUU) and r(UUUAAA) are greatly destablising in RNA, whereas in DNA, d(TTTAAA) is destabilising but d(AAATTT) is stabilising, which has been attributed to the formation of a special B'structure involving large propeller twists of the A-T base pairs. The solution structure of the RNA dodecamer r(CGCAAAUUUGCG)2has been determined using NMR and restrained molecular dynamics calculations to assess the conformational reasons for its stability in comparison with d(CGCAAATTTGCG)2. The structures refined to a mean pairwise r.m.s.d. of 0.89+/-0.29 A. The nucleotide conformations are typical of the A family of structures. However, although the helix axis displacement is approximately 4.6 A into the major groove, the rise (3.0 A) and base inclination ( approximately 6 degrees ) are different from standard A form RNA. The extensive base-stacking found in the AAATTT tract of the DNA homologue that is largely responsible for the higher thermodynamic stability of the DNA duplex is reduced in the RNA structure, which may account for its low relative stability.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

Hydration of the RNA duplex r(CGCAAAUUUGCG)2 determined by NMR

The so-called spine of hydration in the minor groove of AnTn tracts in DNA is thought to stabilise the structure, and kinetically bound water detected in the minor groove of such DNA species by NMR has been attributed to a narrow minor groove [Liepinsh, E., Leupin, W. and Otting, G. (1994) Nucleic Acids Res. 22, 2249-2254]. We report here an NMR study of hydration of an RNA dodecamer which has a wide, shallow minor groove. Complete assignments of exchangeable protons, and a large number of non-exchangeable protons in r(CGCAAAUUUGCG)2 have been obtained. In addition, ribose C2'-OH resonances have been detected, which are probably involved in hydrogen bonds. Hydration at different sites in the dodecamer has been measured using ROESY and NOESY experiments at 11.75 and 14.1 T. Base protons in both the major and minor grooves are in contact with water, with effective correlation times for the interaction of approximately 0.5 ns, indicating weak hydration, in contrast to the hydration of adenine C2H in the homologous DNA sequence. NOEs to H1' in the minor groove are consistent with hydration water present that is not observed in the analogous DNA sequence. Hydration kinetics in nucleic acids may be determined by chemical factors such as hydrogen-bonding more than by simple conformational factors such as groove width.

Comparison of the thermodynamic stabilities and solution conformations of DNA.RNA hybrids containing purine-rich and pyrimidine-rich strands with DNA and RNA duplexes

The conformations and thermodynamic stabilities of duplexes containing purine-rich (GAA-GAGAAGC) and pyrimidine-rich (GCTTCTCTTC or GCUUCUCUUC) DNA and RNA strands have been measured by UV melting, electrophoresis, circular dichroism, and NMR spectroscopy. The free energies of stabilization (-delta G) were in the order rR.rY > rR.dY > dR.dY > dR.rY. The two DNA.RNA hybrid duplexes showed conformational properties intermediate between those of DNA.DNA and RNA.RNA duplexes and also different from one another. Differences between 1H chemical shifts of the DNA strands in the two hybrid duplexes and those of the DNA duplex were larger than analogous shift differences for the RNA protons, and the differences were larger for the purine than the pyrimidines in the DNA strands. Detailed analysis of the nucleotide conformations using both NOE and scalar coupling data showed that the sugar conformations of the ribonucleotides are all near C3'-endo. The deoxyribonucleotides were in the "S" domain, i.e., near C2'-endo in the DNA duplex, and C1'-exo to C2'-endo in the two hybrids. However, the deoxyriboses in the two hybrids appear more flexible than in the DNA duplex, with the fraction in the "N" (C3'-endo) state increasing in the order dR.dY < dR.rY < rR.dY. Globally, the pure DNA duplex was B-form and the pure RNA duplex A form. The two DNA.RNA hybrids were neither A nor B, but closer globally to the A than the B form. The less stable dR.rY duplex has a significantly different conformation from rR.dY both at the local nucleotide level and globally.

Synthesis and properties of DNA-PNA chimeric oligomers

Adenine, thymine and cytosine PNA monomers have been prepared using 3-amino-1,2-propanediol as a starting material. The benzoyl group was used to protect the exocyclic amines of the heterocyclic bases of A and C PNA monomers and the backbone primary amine was protected with the monomethoxytrityl group. The thymine and cytosine PNA monomers were used in conjunction with standard DNA synthesis monomers to produce chimeric PNA DNA (PDC) oligomers. Ultraviolet melting studies confirmed that these oligomers form stable hybrids with complementary DNA strands and that mismatches in the DNA but more so in the PNA sections lead to duplex destabilisation.

Solution structure of (rGCGGACGC)2 by two-dimensional NMR and the iterative relaxation matrix approach

The three-dimensional solution structure of the RNA self-complementary duplex [sequence: see text] was derived from two-dimensional NMR and the iterative relaxation matrix approach. Each GA mismatch forms two hydrogen bonds: A-NH6 to G-O6 and A-N1 to G-NH1 (imino). This is the first three-dimensional RNA structure with imino hydrogen-bonded tandem GA mismatches. This GA structure is totally different from the sheared tandem GA structure in [sequence: see text] which also has two hydrogen bonds: A-N7 to G-NH2 and A-NH6 to G-N3 [SantaLucia, J., Jr., & Turner, D. H. (1993) Biochemistry 32, 12612-12623]. In particular, the sheared and imino GA mismatches produce a narrowing and widening of the backbone, respectively. The results show that substitutions of Watson-Crick base pairs can have dramatic effects on the three-dimensional structures of adjacent non-Watson-Crick paired regions; i.e., the structure depends on sequence context. Thus compensating substitutions in site-directed mutagenesis experiments may not always restore biological activities.

Thermodynamics of base pairing

The study of the thermodynamics of base pairing is a rapidly expanding field due to its importance for taking advantage of opportunities presented by the growing database of nucleic acid sequences. Both natural and unnatural nucleic acids are being studied. The determination of nearest-neighbor parameters from optical melting data for RNA-DNA and DNA-DNA duplexes with Watson-Crick base pairs is particularly noteworthy. The foundations are also being laid for a theoretical understanding of these thermodynamics.

Structural features of a six-nucleotide RNA hairpin loop found in ribosomal RNA

The hairpin loop GUAAUA occurs frequently in ribosomal RNA. Optical melting studies show that r(GGCGUAAUAGCC) folds into a hairpin containing this loop. The structural features of the r(GGCGUAAUAGCC) hairpin have been determined by NMR and molecular modeling. NOEs from G4-H1' to A9-H2 and from A9-H2 to G10-H1' show that G4 and A9 form a sheared base pair with two hydrogen bonds: A-N7 to G-NH2 and A-NH6 to G-N3. One-dimensional NOE data show no NOEs between the imino protons of U5 and U8, but NOEs are observed between the U5-H1' and the U8-H6 and U8-H5, thus orienting the U8 imino proton away from U5. Thus U5 and U8 do not form an imino hydrogen-bonded U.U pair. The U5-H2' exhibits NOEs to both the A6-H8 and A7-H8, and the 3' phosphorus resonances of U5 and A6 are shifted downfield. This suggests that the helix turn is between the U5 and A6 nucleotides. The JH1'-H2 and JH3'-H4' coupling constants indicate that the loop is dynamic, particularly at 35 degrees C, well below the melting temperature of 63 degrees C. Structures were generated using 75 distance and 46 dihedral angle restraints. In these structures, the U5 base is stacked on the sheared base pair formed by G4 and A9 and can initiate a uridine turn similar to that observed in the anticodon loop of tRNA. The A6, A7, and U8 bases can stack on one another with their hydrogen-bonding surfaces exposed to the solvent, suggesting that they are available for tertiary interactions or protein recognition in rRNA. A range of loop structures are consistent with the data, however. The lack of formation of a U.U mismatch is consistent with a recent model that predicts the stability of hairpin loops of six nucleotides on the basis of the closing base pair and first mismatch in the loop [Serra, M. J., Axenson, T. J., & Turner, D. H. (1994) Biochemistry 33, 14289-14296].

Nonplanar DNA base pairs

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.

High pressure effects on conformation of homo- and heteroduplexes of nucleic acids

Four different chemically synthesized single stranded complementary oligonucleotides: DNA I, d(GCGCGCATATAT); RNA I, r(AUAUAUGCGCGC): RNA II, r(GGCCGGUUAAUU); and RNA III, r(AAUUAACCGGCC) were studied in order that the effects of high pressure on heteroduplex and homoduplex structures could be understood. The oligonucleotides were subjected to a high pressure at low and/or high salt buffer and analyzed by circular dichroism spectroscopy. In these conditions, both DNA-RNA and RNA-RNA duplexes with different purine-pyrimidine sequences change their conformation. The heteroduplex DNA I-RNA I with the complementary alternating purine-pyrimidine sequence, does not change its conformation of A type at high salt alone or at high salt and high pressure applied together. The homoduplex RNA II-RNA III with purine-purine-pyrimidine-pyrimidine sequence does not change strongly its. A-RNA conformation either. However, a structure of the homoduplex is affected by high pressure alone or with high salt as concluded from shifting the maximum of the CD spectrum to around 265 nm and inducing higher Cotton effect. These observations clearly suggest some conformational changes of the homoduplex. A single stranded oligonucleotide (RNA I) and oligodeoxynucleotide (DNA I) alone showed up a different conformation. The CD spectrum of RNA I is similar to that of A-RNA structure, out that of DNA I shows a very small Cotton effect and has not an ordered structure.

A periodic table of symmetric tandem mismatches in RNA

The stabilities and structures of a series of RNA octamers containing symmetric tandem mismatches were studied by UV melting and imino proton NMR. The free energy increments for tandem mismatch formation are found to depend upon both mismatch sequence and adjacent base pairs. The observed sequence dependence of tandem mismatch stability is UGGU > GUUG > GAAG > or = AGGA > UUUU > CAAC > or = CUUC approximately UCCU approximately CCCC approximately ACCA approximately AAAA, and the closing base pair dependence is 5'G3'C > 5'C3'G > 5'U3'A approximately 5'A3'U. These results differ from expectations based on models used in RNA folding algorithms and from the sequence dependence observed for folding of RNA hairpins. Imino proton NMR results indicate the sequence dependence is partially due to hydrogen bonding within mismatches.