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.     

A DNA repair process in Escherichia coli corrects U:G and T:G mismatches to C:G at sites of cytosine methylation

Escherichia coli contains a base mismatch correction system called VSP repair that is known to correct T:G mismatches to C:G when they occur in certain sequence contexts. The preferred sequence context for this process is the site for methylation by the E. coli DNA cytosine methylase (Dcm). For this reason, VSP repair is thought to counteract potential mutagenic effects of deamination of 5-methylcytosine to thymine. We have developed a genetic reversion assay that quantitates the frequency of C to T mutations at Dcm sites and the removal of such mutations by DNA repair processes. Using this assay, we have studied the repair of U: G mismatches in DNA to C: G and have found that VSP repair is capable of correcting these mismatches. Although VSP repair substantially affects the reversion frequency, it may not be as efficient at correcting U: G mismatches as the uracil DNA glycosylase-mediated repair process.     

A DNA repair process in Escherichia coli corrects U:G and T:G mismatches to C:G at sites of cytosine methylation

Escherichia coli contains a base mismatch correction system called VSP repair that is known to correct T:G mismatches to C:G when they occur in certain sequence contexts. The preferred sequence context for this process is the site for methylation by the E. coli DNA cytosine methylase (Dcm). For this reason, VSP repair is thought to counteract potential mutagenic effects of deamination of 5-methylcytosine to thymine. We have developed a genetic reversion assay that quantitates the frequency of C to T mutations at Dcm sites and the removal of such mutations by DNA repair processes. Using this assay, we have studied the repair of U: G mismatches in DNA to C: G and have found that VSP repair is capable of correcting these mismatches. Although VSP repair substantially affects the reversion frequency, it may not be as efficient at correcting U: G mismatches as the uracil DNA glycosylase-mediated repair process.     

The thermal stability of DNA fragments with tandem mismatches at a d(CXYG).d(CY'X'G) site.

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.     

An intramolecular quadruplex of (GGA)(4) triplet repeat DNA with a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad, and its dimeric interaction.

The structure of d(GGAGGAGGAGGA) containing four tandem repeats of a GGA triplet sequence has been determined under physiological K(+) conditions. d(GGAGGAGGAGGA) folds into an intramolecular quadruplex composed of a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad. Four G-G segments of d(GGAGGAGGAGGA) are aligned parallel with each other due to six successive turns of the main chain at each of the GGA and GAGG segments. Two quadruplexes form a dimer stabilized through a stacking interaction between the heptads of the two quadruplexes. Comparison of the structure of d(GGAGGAGGAGGA) with the reported structure of d(GGAGGAN) (N=G or T) containing two tandem repeats of the GGA triplet revealed that although the two structures resemble each other to some extent, the extension of the repeats of the GGA triplet leads to distinct structural differences: intramolecular quadruplex for 12-mer versus intermolecular quadruplex for 7-mer; heptad versus hexad in the quadruplex; and three sheared G:A base-pairs versus two sheared G:A base-pairs plus one A:A base-pair per quadruplex. It was also suggested that d(GGAGGAGGAGGA) forms a similar quadruplex under low salt concentration conditions. This is in contrast to the case of d(GGAGGAN) (N=G or T), which forms a duplex under low salt concentration conditions. On the basis of these results, the structure of naturally occurring GGA triplet repeat DNA is discussed.     

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

Very stable mismatch duplexes: structural and thermodynamic studies on tandem G.A mismatches in DNA.

We have used ultraviolet melting techniques to compare the stability of several DNA duplexes containing tandem G.A mismatches to similar duplexes containing tandem A.G, I.A, and T.A base pairs. We have found that tandem G.A mismatches in 5'-Y-G-A-R-3' duplexes are more stable than their I.A counterparts and that they are sometimes more stable than tandem 5'-Y-T-A-R-3' sequences. This is not the case for tandem G.A mismatches in other base stacking environments, and it suggests that tandem G.A mismatches in 5'-Y-G-A-R-3' sequences have a unique configuration. In contrast to tandem 5'-G-A-3' mismatches, tandem 5'-A-G-3' mismatches were found to be unstable in all cases examined.     

DNA mismatch correction by Very Short Patch repair may have altered the abundance of oligonucleotides in the E. coli genome.

A base mismatch correction process in E. coli K-12 called Very Short Patch (VSP) repair corrects T:G mismatches to C:G when found in certain sequence contexts. Two of the substrate mismatches (5'-CTWGG/3'-GGW'CC; W = A or T) occur in the context of cytosine methylation in DNA and reduce the mutagenic effects of 5-methylcytosine deamination to thymine. However, VSP repair is also known to repair T:G mismatches that are not expected to arise from 5-methylcytosine deamination (example--CTAG/GGT-C). In these cases, if the original base pair were a T:A, VSP repair would cause a T to C transition. We have carried out Markov chain analysis of an E. coli sequence database to determine if repair at the latter class of sites has altered the abundance of the relevant tetranucleotides. The results are consistent with the prediction that VSP repair would tend to deplete the genome of the 'T' containing sequences (example--CTAG), while enriching it for the corresponding 'C' containing sequences (CCAG). Further, they provide an explanation for the known scarcity of CTAG containing restriction enzyme sites among the genomes of enteric bacteria and identify VSP repair as a force in shaping the sequence composition of bacterial genomes.     

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.     

Nearest-neighbor parameters for G.U mismatches: [formula; see text] is destabilizing in the contexts [formula; see text] and [formula; see text] but stabilizing in [formula; see text].

Thermodynamic parameters derived from optical melting studies are reported for duplex formation by a series of oligoribonucleotides containing G.U mismatches. The results are used to determine nearest-neighbor parameters for helix propagation by G.U mismatches. Surprisingly, the [formula; see text] nearest-neighbor free energy increment in unfavorable in the contexts [formula; see text], and [formula; see text] but favorable in the context [formula; see text]. This is a non-nearest-neighbor effect. In contrast, the [formula; see text] free energy increment is favorable and independent of context. Circular dichroism and imino proton NMR spectra of several sequences do not reveal an obvious structural basis for this dichotomy. For example, all the G.U mismatches have two slowly exchanging imino protons. The imino resonances for the G.U mismatches in GGAGUUCC, GUCGUGAC, and CCUGUAGG, however, broaden at lower temperature than the imino resonances for the interior Watson-Crick base pairs. In contrast, the imino resonances for the G.U mismatches in GGAUGUCC remain sharp at high temperature. The improved parameters for G.U mismatches should improve predictions of RNA structure from sequence.     

Comparative Studies of the Thermodynamic Stabilities Between Sheared A:G and Watson-Crick A:U(T) Base Pairs in RNA and DNA

Abstract

Thermodynamic parameters for duplex formation were determined from CD melting curves for r(GGACGAGUCC)₂ and d(GGACGAGTCC)₂, both of which form two consecutive ‘sheared’ A:G base pairs at the center [Katahira et al. (1993) Nucleic Acids Res. 21, 5418–5424; Katahira et al., (1994) Nucleic Acids Res. 22, 2752–27591. The parameters were determined also for r(GGACUAGUCC)₂ and d(GGACTAGTCC)₂, where the A:G mismatches are replaced by Watson-Crick A:U(T) base pairs. Thermodynamic properties for duplex formation are compared between the sheared and the Watson-Crick base pairs, and between RNA and DNA. Difference in the thermodynamic stability is analyzed and discussed in terms of enthalpy and entropy changes. The characteristic features in CD spectra of RNA and DNA containing the sheared A:G base pairs are also reported.

     

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.     

Structure of the ribozyme substrate hairpin of Neurospora VS RNA: a close look at the cleavage site

The cleavage site of the Neurospora VS RNA ribozyme is located in a separate hairpin domain containing a hexanucleotide internal loop with an A-C mismatch and two adjacent G-A mismatches. The solution structure of the internal loop and helix la of the ribozyme substrate hairpin has been determined by nuclear magnetic resonance (NMR) spectroscopy. The 2 nt in the internal loop, flanking the cleavage site, a guanine and adenine, are involved in two sheared G.A base pairs similar to the magnesium ion-binding site of the hammerhead ribozyme. Adjacent to the tandem G.A base pairs, the adenine and cytidine, which are important for cleavage, form a noncanonical wobble A+-C base pair. The dynamic properties of the internal loop and details of the high-resolution structure support the view that the hairpin structure represents a ground state, which has to undergo a conformational change prior to cleavage. Results of chemical modification and mutagenesis data of the Neurospora VS RNA ribozyme can be explained in context with the present three-dimensional structure.     

Thermodynamics of RNA hairpins containing single internal mismatches.

Thermodynamic parameters and circular dichroism spectra are presented for RNA hairpins containing single internal mismatches in the stem regions. Three different sequence contexts for the G*U mismatch and two contexts for C*A, G*A, U*U, A*C and U*G mismatches were examined and compared with Watson-Crick base-pair stabilities. The RNA hairpins employed were a microhelix and tetraloop representing the Escherichia coli tRNAAlaacceptor stem and sequence variants that have been altered at the naturally occurring G*U mismatch site. UV melting studies were carried out under different conditions to evaluate the effects of sodium ion concentration and pH on the stability of mismatch-containing hairpins. Our main findings are that single internal mismatches exhibit a range of effects on hairpin stability. In these studies, the size and sequence of the loop and stem are shown to influence the overall stability of the RNA, and have a minor effect on the relative mismatch stabilities. The relationship of these results to RNA-ligand interactions involving mismatch base-pairs is discussed.     

Substrate specificity and sequence-dependent activity of the Saccharomyces cerevisiae 3-methyladenine DNA glycosylase (Mag)

DNA glycosylases initiate base excision repair by first binding, then excising aberrant DNA bases. Saccharomyces cerevisiae encodes a 3-methyladenine (3MeA) DNA glycosylase, Mag, that recognizes 3MeA and various other DNA lesions including 1,N6-ethenoadenine (epsilon A), hypoxanthine (Hx) and abasic (AP) sites. In the present study, we explore the relative substrate specificity of Mag for these lesions and in addition, show that Mag also recognizes cisplatin cross-linked adducts, but does not catalyze their excision. Through competition binding and activity studies, we show that in the context of a random DNA sequence Mag binds epsilon A and AP-sites the most tightly, followed by the cross-linked 1,2-d(ApG) cisplatin adduct. While epsilon A binding and excision by Mag was robust in this sequence context, binding and excision of Hx was extremely poor. We further studied the recognition of epsilon A and Hx by Mag, when these lesions are present at different positions within A:T and G:C tracts. Overall, epsilon A was slightly less well excised from each position within the A:T and G:C tracts compared to excision from the random sequence, whereas Hx excision was greatly increased in these sequence contexts (by up to 7-fold) compared to the random sequence. However, given most sequence contexts, Mag had a clear preference for epsilon A relative to Hx, except in the TTXTT (X=epsilon A or Hx) sequence context from which Mag removed both lesions with almost equal efficiency. We discuss how DNA sequence context affects base excision by various 3MeA DNA glycosylases.     

Solution structure of Cobalt(III)hexammine complexed to the GAAA tetraloop, and metal-ion binding to G.A mismatches

The solution structure of a 22 nt RNA hairpin and its complex with Co(NH(3))(6)(3+) bound to the GAAA tetraloop has been determined by NMR spectroscopy. Co(NH(3))(6)(3+) has a similar geometry to Mg(H(2)O)(6)(2+) and can be used as a probe for binding sites of completely solvated magnesium ions. The hairpin contains tandem G.A mismatches, similar to the P5abc region of a group I intron, and is closed by a GAAA tetraloop. The tandem G.A mismatches are imino hydrogen bonded in contrast with the sheared G.A mismatches found in a different context in the crystal structure of the P4-P6 domain. Chemical shift changes of the imino protons upon titration of the RNA hairpin with Mg(2+) and with Co(NH(3))(6)(3+) were used to identify ion-binding sites. Paramagnetic resonance broadening upon titration with Mn(2+) was also used. The titration curves gave dissociation binding constants for the magnesium ions in the millimolar range, similar to the binding in the major groove of RNA at tandem G.U base-pairs. Although the largest chemical shift change occurred at an imino proton of one of the G.A base-pairs, no nuclear Overhauser enhancement cross-peaks between the cobalt ligand and neighboring RNA protons were seen, presumably due to the high mobility of the Co(NH(3))(6)(3+) at this site. Nuclear Overhauser enhancement cross-peaks between Co(NH(3))(6)(3+) and the GAAA tetraloop were observed, which allowed the determination of the structure of the tetraloop binding site. The Co(NH(3))(6)(3+) is bound in the major groove of the GAAA tetraloop with hydrogen bonds to guanine base N7 and to phosphate oxygen atoms of the tetraloop.     

Experimental demonstration of T:(G:G:G:G):T hexad and T:A:A:T tetrad alignments within a DNA quadruplex stem

A template-based approach was used to design unprecedented architectural motifs into a known DNA framework. The structure formed by the sequence d(GCGGTTGGAT) in 0.1 M Na(+) solution has been determined using molecular dynamics simulations constrained by distance and dihedral restraints derived from NMR experiments. The molecular topology has been previously observed for the sequence d(GCGGTGGAT) (Webba da Silva, M. (2003) Biochemistry 42, 14356-65). Insertion of a single thymine into the double chain reversal formed by the segment GGTGG results in the unprecedented experimental demonstration of a T:(G:G:G:G):T hexad. The bi-stranded hexad results from the pairing alignment of two G(T-G) triads. Each triad results from recognition of the sheared edge of a guanine by the Watson-Crick edge of a thymine of the segment GGTTGG. The alignment is stabilized by base-stacking of the thymine to the sugar pucker of the preceding thymine. The latter is involved in formation of the T:A:A:T tetrad alignment by forming a hydrogen bond with the free amino proton of a Watson-Crick aligned A:A mispair. We have thus established that residues in double chain reversal loops linking juxtaposed tetrads of a quadruplex stem may facilitate formation of yet unknown hydrogen bond alignments. By employing a systematic approach analysis of sequence motifs appearing in double chain reversals, bridging tetrad layers should allow for the prediction of topologies and architectural motifs appearing in biologically relevant genomic regions.     

Formation of sheared G:A base pairs in an RNA duplex modelled after ribozymes, as revealed by NMR.

The thermal stability and structure of an RNA duplex, r(GGACGAGUCC)2, the base sequence of which was modelled after both a hammerhead ribozyme and a lead ribozyme, were studied by CD and NMR. We previously demonstrated that the corresponding DNA duplex, d(GGACGAGTCC)2, formed unique 'sheared' G:A base pairs, 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' (Nucleic Acids Res. (1993) 21, 5418-5424). CD melting profiles showed that the RNA duplex is thermally more stable than the corresponding DNA duplex. NMR studies revealed that sheared G:A base pairs are formed in the RNA duplex, too, although the overall structure of the RNA is the A form, which differs from the B form taken on by the corresponding DNA. A model building study confirmed that sheared G:A base pairs can be accommodated in the double helical structure of the A form. A difference between the RNA and DNA duplexes in the stacking interaction involving G:A mismatch bases is also suggested. The demonstration that sheared G:A base pairs can be formed not only in DNA but also in RNA suggests that this base pairing plays an important role regarding the RNA structure.