Homopurine and homopyrimidine strands complementary in parallel orientation form an antiparallel duplex at neutral pH with A-C,G-T,and T-C mismatched base pairs

DNA sequences d-TGAGGAAAGAAGGT (a 14-mer) and d-CTCCTTTCTTCC (a 12-mer) are complementary in parallel orientation forming either Donahue (reverse Watson-Crick) base pairing at neutral pH or Hoogsteen base pairing at slightly acidic pH. The structure of the complex formed by dissolving the two strands in equimolar ratio in water has been investigated by nmr. At neutral pH, the system forms an ordered antiparallel duplex with five A : T and four G : C Watson-Crick base pairs and three mismatches, namely G-T, A-C, and T-C. The nuclear Overhauser effect cross-peak pattern suggests an overall B-DNA conformation with major structural perturbations near the mismatches. The duplex has a low melting point and dissociates directly into single strands with a broad melting profile. The hydrogen-bonding schemes in the mismatched base pairs have been investigated. It has been shown earlier that in acidic pH, the system prefers a triple-stranded structure with two pyrimidine strands and one purine strand. One of the pyrimidine strands has protonated cytosines, forms Hoogsteen base pairing, and is aligned parallel to the purine strand; the other has nonprotonated cytosines and has base-pairing scheme similar to the one discussed in this paper. The parallel duplex is therefore less stable than either the antiparallel duplex or the triplex, in spite of its perfect complementarity. © 1997 John Wiley & Sons, Inc. Biopoly 41: 773–784, 1997     

Stability of non-Watson-Crick G-A/A-G base pair in synthetic DNA and RNA oligonucleotides

A non-Watson-Crick G-A/A-G base pair is found in SECIS (selenocysteine-insertion sequence) element in the 3'-untranslated region of Se-protein mRNAs and in the functional site of the hammerhead ribozyme. We studied the stability of G-A/A-G base pair (bold) in 17mer GT(U)GACGGAAACCGGAAC synthetic DNA and RNA oligonucleotides by thermal melting experiments and gel electrophoresis. The measured Tm value of DNA oligonucleotide having G-A/A-G pair showed an intermediate value (58 degrees C) between that of Watson-Crick G-C/C-G base pair (75 degrees C) and that of G-G/A-A of non-base-pair (40 degrees C). Similar thermal melting patterns were obtained with RNA oligonucleotides. This result indicates that the secondary structure of oligonucleotide having G-A/A-G base pair is looser than that of the G-C type Watson-Crick base pair. In the comparison between RNA and DNA having G-A/A-G base pair, the Tm value of the RNA oligonucleotide was 11 degrees C lower than that of DNA, indicating that DNA has a more rigid structure than RNA. The stained pattern of oligonucleotide on polyacrylamide gel clarified that the mobility of the DNA oligonucleotide G-A/A-G base pair changed according to the urea concentration from the rigid state (near the mobility of G-C/C-G oligonucleotide) in the absence of urea to the random state (near the mobility of G-G/A-A oligonucleotide) in 7 M urea. However, the RNA oligonucleotide with G-A/A-G pair moved at an intermediate mobility between that of oligonucleotide with G-C/C-G and of the oligonucleotide with G-G/A-A, and the mobility pattern did not depend on urea concentration. Thus, DNA and RNA oligonucleotides with the G-A/A-G base pair showed a pattern indicating an intermediate structure between the rigid Watson-Crick base pair and the random structure of non-base pair. RNA with G-A/A-G base pair has the intermediate structure not influenced by urea concentration. Finally, this study indicated that the intermediate rigidity imparted by Non-Watson-Crick base pair in SECIS element plays an important role in the selenocysteine expression by UGA codon.     

Specific oligonucleotide invasion into the end of DNA duplex

Phenomenon of the interaction of a double-stranded DNA fragment with an oligonucleotide complementary to the end of the duplex strand was demonstrated to occur via formation of three-stranded DNA structure with an oligonucleotide invasion. It was shown that oligonucleotides complementary to the duplex ends inhibit Holliday junction formation in solutions of homologous linear DNA fragments. This effect depends on the oligonucleotide concentration, sequence and their complementarity to the duplex ends. Formation of three-stranded complexes was demonstrated using radiolabeled oligonucleotides by agarose gel-electrophoresis followed by autoradiography. Analysis of three-stranded DNA structures by chemical cleavage of non-canonical base pairs revealed that oligonucleotide invades into duplex ends via a sequential displacement mechanism and that the level of the invasion may vary considerably.     

Mismatch formation in solution and on DNA microarrays: how modified nucleosides can overcome shortcomings of imperfect hybridization caused by oligonucleotide composition and base pairing

The DNA microarray technology is a well-established and widely used technology although it has several drawbacks. The accurate molecular recognition of the canonical nucleobases of probe and target is the basis for reliable results obtained from microarray hybridization experiments. However, the great flexibility of base pairs within the DNA molecule allows the formation of various secondary structures incorporating Watson-Crick base pairs as well as non-canonical base pair motifs, thus becoming a source of inaccuracy and inconsistence. The first part of this report provides an overview of unusual base pair motifs formed during molecular DNA interaction in solution highlighting selected secondary structures employing non-Watson-Crick base pairs. The same mispairing phenomena obtained in solution are expected to occur for immobilized probe molecules as well as for target oligonucleotides employed in microarray hybridization experiments the effect of base pairing and oligonucleotide composition on hybridization is considered. The incorporation of nucleoside derivatives as close shape mimics of the four canonical nucleosides into the probe and target oligonucleotides is discussed as a chemical tool to resolve unwanted mispairing. The second part focuses non-Watson-Crick base pairing during hybridization performed on microarrays. This is exemplified for the unusual stable dG.dA base pair.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

An innate twist between Crick’s wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.     

Analysis of RNA motifs

RNA motifs are directed and ordered stacked arrays of non-Watson-Crick base pairs forming distinctive foldings of the phosphodiester backbones of the interacting RNA strands. They correspond to the 'loops' - hairpin, internal and junction - that intersperse the Watson-Crick two-dimensional helices as seen in two-dimensional representations of RNA structure. RNA motifs mediate the specific interactions that induce the compact folding of complex RNAs. RNA motifs also constitute specific protein or ligand binding sites. A given motif is characterized by all the sequences that fold into essentially identical three-dimensional structures with the same ordered array of isosteric non-Watson-Crick base pairs. It is therefore crucial, when analyzing a three-dimensional RNA structure in order to identify and compare motifs, to first classify its non-Watson-Crick base pairs geometrically.     

Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

Mispairs with Watson‐Crick base‐pair geometry observed in ternary complexes of an RB69 DNA polymerase variant

Recent structures of DNA polymerase complexes with dGMPCPP/dT and dCTP/dA mispairs at the insertion site have shown that they adopt Watson‐Crick geometry in the presence of Mn²⁺ indicating that the tautomeric or ionization state of the base has changed. To see whether the tautomeric or ionization state of base‐pair could be affected by its microenvironment, we determined 10 structures of an RB69 DNA polymerase quadruple mutant with dG/dT or dT/dG mispairs at position n‐1 to n‐5 of the Primer/Template duplex. Different shapes of the mispairs, including Watson‐Crick geometry, have been observed, strongly suggesting that the local environment of base‐pairs plays an important role in their tautomeric or ionization states.     

Structure of the DNA interstrand cross-link of 4,5',8-trimethylpsoralen.

4,5',8-Trimethylpsoralen (TMP) cross-links a 5' TpA or a 5' ApT site by photoreacting with one thymine moiety in each DNA strand. We are interested in whether psoralen interstrand cross-links all share one structure or whether there are significant differences. In this paper, we employed a rapid method for probing the structure of the cross-link by making a series of TMP cross-linked duplexes containing specific base-pair mismatches. The relative stability provided by a base pair can be correlated with neighboring base pairs by comparing the extents of gel retardation when base-pair mismatches happen in each position. From our studies, we infer that with respect to the furan-side strand, the 5'T.A base pair of the two T.A base pairs in the TpA site is not hydrogen bonded. Immediately on each side of the cross-linked TpA site is a highly stabilized base pair. Next, a region of decreased stability occurs in each arm of a cross-linked duplex and these base pairs of least stability are located farther away from the cross-linked thymines as the lengths of the arms of the cross-linked helix increase. Finally, even in 7 M urea at 49 degrees C the cross-linked helix is hydrogen bonded at both ends of a duplex of 22 base pairs. We propose that the structures of interstrand cross-links in DNA vary appreciably with the DNA sequence, the length of the DNA duplex, and the structures of the DNA cross-linking agents.     

Temperature dependence of thermodynamic properties for DNA/DNA and RNA/DNA duplex formation.

A clear difference in the enthalpy changes derived from spectroscopic and calorimetric measurements has recently been shown. The exact interpretation of this deviation varied from study to study, but it was generally attributed to the non-two-state transition and heat capacity change. Although the temperature-dependent thermodynamics of the duplex formation was often implied, systemic and extensive studies have been lacking in universally assigning the appropriate thermodynamic parameter sets. In the present study, the 24 DNA/DNA and 41 RNA/DNA oligonucleotide duplexes, designed to avoid the formation of hairpin or slipped duplex structures and to limit the base pair length less than 12 bp, were selected to evaluate the heat capacity changes and temperature-dependent thermodynamic properties of duplex formation. Direct comparison reveals that the temperature-independent thermodynamic parameters could provide a reasonable approximation only when the temperature of interest has a small deviation from the mean melting temperature over the experimental range. The heat capacity changes depend on the base composition and sequences and are generally limited in the range of -160 to approximately -40 cal.mol-1.K-1 per base pair. In contrast to the enthalpy and entropy changes, the free energy change and melting temperature are relatively insensitive to the heat capacity change. Finally, the 16 NN-model free energy parameters and one helix initiation at physiological temperature were extracted from the temperature-dependent thermodynamic data of the 41 RNA/DNA hybrids.     

Application of synthetic oligonucleotides to the diagnosis of human genetic diseases

Under appropriate conditions synthetic oligonucleotide hybridization probes display essentially absolute hybridization specificity. That is, every nucleotide must form a Watson-Crick base pair in order that the probe forms a stable duplex. All of the non-Watson-Crick base pairs, including G-T, have a destabilizing effect. Thus, it is possible to choose stringent conditions of hybridization such that, while a perfectly matched duplex between an oligonucleotide and complementary DNA will form, duplexes mismatched at one or more position will not. Mutations in a single base in the DNA sequence of a gene can and do result in genetic diseases. The hybridization of oligonucleotides to the region of DNA containing these base changes would be affected by the mutations and thus, oligonucleotide hybridization provides a means of detecting single base changes. In an attempt to develop a non-radioactive method for the detection of human genetic diseases, we have prepared biotinylated-oligonucleotides by an enzymatic method. An oligonucleotide probe (23-mer) containing a single biotinylated deoxyuridine residue at the 3'-terminus was prepared by a primer extention reaction using E. coli DNA polymerase I (Klenow fragment). The probe could be specifically and tightly bound with Avidin D in 1 M NaCl. It could be hybridized to a plasmid DNA containing a perfectly matched complementary sequence, but not to a DNA containing 5 non-consecutive non-complementary bases. The hybridized biotinylated probe could be visualized by Avidin D and biotinylated alkaline phosphatase, even when 1.8 ng of the plasmid DNA (0.5 fmol) was used. A general approach to the enzymatic synthesis of oligonucleotides containing a single biotinylated deoxyuridine at the 3' end is described.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

Effect of LNA- and OMeN-modified oligonucleotide probes on the stability and discrimination of mismatched base pairs of duplexes

Locked nucleic acid (LNA) and 2'-O-methyl nucleotide (OMeN) are the most extensively studied nucleotide analogues. Although both LNA and OMeN are characterized by the C3'-endo sugar pucker conformation, which is dominant in A-form DNA and RNA nucleotides, they demonstrate different binding behaviours. Previous studies have focused attention on their properties of duplex stabilities, hybridization kinetics and resistance against nuclease digestion; however, their ability to discriminate mismatched hybridizations has been explored much less. In this study, LNA- and OMeN-modified oligonucleotide probes have been prepared and their effects on the DNA duplex stability have been examined: LNA modifications can enhance the duplex stability, whereas OMeN modifications reduce the duplex stability. Next, we studied how the LNA:DNA and OMeN:DNA mismatches reduced the duplex stability. Melting temperature measurement showed that different LNA:DNA or OMeN:DNA mismatches indeed influence the duplex stability differently. LNA purines can discriminate LNA:DNA mismatches more effectively than LNA pyrimidines as well as DNA nucleotides. Furthermore, we designed five LNA- and five OMeN-modified oligonucleotide probes to simulate realistic situations where target-probe duplexes contain a complementary LNA:DNA or OMeN:DNA base pairs and a DNA:DNA mismatch simultaneously. The measured collective effect showed that the duplex stability was enhanced by the complementary LNA:DNA base pair but decreased by the DNA:DNA mismatch in a position-dependent manner regardless of the chemical identity and position of the complementary LNA:DNA base pair. On the other hand, the OMeN-modified probes also showed that the duplex stability was reduced by both the OMeN modification and the OMeN:DNA mismatch in a position-dependent manner.     

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

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

The solution structure of a locked nucleic acid (LNA) hybridized to DNA

LNA (Locked Nucleic Acids) is a novel oligonucleotide analogue containing a conformationally restricted nucleotide with a 2'-O, 4'-C-methylene bridge that induces unprecedented thermal affinities when mixed with complementary single stranded DNA and RNA. We have used two-dimensional 1H NMR spectroscopy obtained at 750 and 500 MHz to determine a high resolution solution structure of an LNA oligonucleotide hybridized to the complementary DNA strand. The determination of the structure was based on a complete relaxation matrix analysis of the NOESY cross peaks followed by restrained molecular dynamics calculations. Forty final structures were generated for the duplex from A-type and B-type dsDNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the forty structures of the complex was 0.32A. The structures were analysed by use of calculated helix parameters. This showed that the values for rise and buckle in the LNA duplex is markedly different from canonical B-DNA at the modification site. A value of twist similar to A-DNA is also observed at the modification site. The overall length of the helix which is 27.3 A. The average twist over the sequence are 35.9 degrees +/- 0.3 degrees. Consequently, the modification does not cause the helix to unwind. The bis-intercalation of the thiazole orange dye TOTO to the LNA duplex was also investigated by 1H NMR spectroscopy to sense the structural change from the unmodified oligonucleotide. We observed that the bis-intercalation of TOTO is much less favourable in the 5'-CT(L)AG-3' site than in the unmodified 5'-CTAG-3' site. This was related to the change in the base stacking of the LNA duplex compared to the unmodified duplex.     

Properties of circular dumbbell RNA/DNA chimeric oligonucleotides containing antisense phosphodiester oligonucleotides.

We have designed a new class of oligonucleotides, "dumbbell RNA/DNA chimeric phosphodiesters", containing two alkyl loop structures with RNA/DNA base pairs (sense (RNA) and antisense (DNA) in the double helical stem. The reaction of nicked (NDRDON) and circular (CDRDON) dumbbell RNA/DNA chimeric oligonucleotides with RNaseH gave the corresponding antisense phosphodiester oligonucleotide together with the sense RNA cleavage products. The liberated antisense phosphodiester oligodeoxynucleotide was bound to the target 35mer RNA, which gave 35mer RNA cleavage products by treatment with RNaseH. The circular dumbbell RNA/DNA chimeric oligonucleotide showed more nuclease resistance than the linear antisense phosphodiester oligodeoxynucleotide(anti-ODN) and the nicked dumbbell RNA/DNA chimeric oligonucleotide.