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.     

Solution conformation of a deoxynucleotide containing tandem G.A mismatched base pairs and 3'-overhanging ends in d(GTGAACTT)2.

We have used 31P and 1H NMR spectroscopy and circular dichroism to define the solution conformation of d(GTGAACTT)2 which contains tandem G.A mismatched base pairs and 3'-overhanging TT ends. Measurements of coupling constants and NOE intensities show that the sugar puckers of the nucleotides are predominantly in the south domain (i.e., near C2'-endo) and that the glycosidic torsion angles are anti. The sequential NOE intensities indicate the presence of a right-handed helix. Analysis of the 31P and 1H NMR spectra of the duplex shows that the tandem mismatch forms a block in which there are unusual backbone torsion angles (i.e., in the BII state), within an otherwise B-like structure. The chemical shift of the N1H of the mismatched guanosine and NOEs between the mismatched base pairs and their nearest neighbors are inconsistent with the imino pairing present in single A.G mismatches or in the X-ray structure of a tandem mismatch [Privé, G. G., et al. (1987) Science 238, 498-503] but the data are consistent with the amino pairing found by Li et al. (1991) [Li, Y., et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The strong base-base stacking both within the tandem G.A block and between the G.A mismatches and their other nearest neighbors offsets the intrinsic destabilizing effects of the mismatch. Further, the 3'-TT overhangs stack onto the ends of the helix and stabilize the duplex against fraying, which accounts for the observed increase in the melting temperature compared with the flush-ended duplex.     

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.     

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.     

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.     

Molecular recognition in purine-rich internal loops: thermodynamic,structural, and dynamic consequences of purine for adenine substitutions in 5'(rGGCAAGCCU)2

The contribution of amino groups to the thermodynamics, structure, and dynamics of tandem A.A mismatches is investigated by substitution of purine (P) for adenine (A) within the RNA duplex, 5'(rGGCAAGCCU)(2), to give 5'(rGGCPAGCCU)(2), 5'(rGGCAPGCCU)(2), and 5'(rGGCPPGCCU)(2). The 5'(rGGCAAGCCU)(2) duplex has sheared A(anti).A(anti) (A.A trans Hoogsteen/Sugar-edge) pairs in which the A5 amino group is involved in hydrogen bonds but the A4 amino group is not [Znosko, B. M., Burkard, M. E., Schroeder, S. J., Krugh, T. R., and Turner, D. H. (2002) Biochemistry 41, 14969-14977]. In comparison to 5'(rGGCAAGCCU)(2), replacing the amino group of A4 with a hydrogen stabilizes the duplex by 1.3 kcal/mol, replacement of the A5 amino group destabilizes the duplex by 0.6 kcal/mol, and replacement of both A4 and A5 amino groups destabilizes the duplex by 0.8 kcal/mol. In NMR structures, the P.A noncanonical pairs of the 5'(rGGCPAGCCU)(2) duplex have a sheared anti-anti structure (P.A trans Hoogsteen/Sugar-edge) with P4.A5 interstrand hydrogen bonding and A5 bases that interstrand stack, similar to the structure of 5'(rGGCAAGCCU)(2). In contrast, the A.P pairs of the 5'(rGGCAPGCCU)(2) duplex have a face-to-face conformation (A.P trans Watson-Crick/Watson-Crick) with intrastrand stacking resembling typical A-form geometry. Although the P5 bases in 5'(rGGCPPGCCU)(2) are involved in an interstrand stack, the loop region is largely undefined. The results illustrate that both hydrogen-bonded and non-hydrogen-bonded amino groups play important roles in determining the thermodynamic, structural, and dynamic characteristics of purine rich internal loops.     

Purine-purine mismatches in RNA helices: evidence for protonated G.A pairs and next-nearest neighbor effects.

Thermodynamic parameters are presented for 12 different RNA duplexes containing A.A, A.G, G.A and G.G mismatches flanked by C-G base pairs. UV melting studies were conducted under three different buffer conditions in order to evaluate the effects of salt concentration and pH on the stability of each mismatch-containing duplex. The main findings are: (i) the mismatches have a wide range of effects on duplex stability, decreasing delta G degrees 37 of denaturation by approximately 0-7 kcal/mol; (ii) the nearest-neighbor assumption commonly used to calculate helix stability breaks down for G.A mismatches; and (iii) G.A mismatches separated by 2 bp form a protonated structure.     

Crystal structure of an RNA helix recognized by a zinc-finger protein: an 18-bp duplex at 1.6 A resolution

The crystal structure of the 19-mer RNA, 5'-GAAUGCCUGCGAGCAUCCC-3' has been determined from X-ray diffraction data to 1.6 A resolution by the multiwavelength anomalous diffraction method from crystals containing a brominated uridine. In the crystal, this RNA forms an 18-mer self-complementary double helix with the 19th nucleotide flipped out of the helix. This helix contains most of the target stem recognized by the bacteriophage Mu Com protein (control of mom), which activates translation of an unusual DNA modification enzyme, Mom. The 19-mer duplex, which contains one A.C mismatch and one A.C/G.U tandem wobble pair, was shown to bind to the Com protein by native gel electrophoresis shift assay. Comparison of the geometries and base stacking properties between Watson-Crick base pairs and the mismatches in the crystal structure suggest that both hydrogen bonding and base stacking are important for stabilizing these mismatched base pairs, and that the unusual geometry adopted by the A.C mismatch may reveal a unique structural motif required for the function of Com.     

Thermodynamic characterization of the complete set of sequence symmetric tandem mismatches in RNA and an improved model for predicting the free energy contribution of sequence asymmetric tandem mismatches

Because of the availability of an abundance of RNA sequence information, the ability to rapidly and accurately predict the secondary structure of RNA from sequence is becoming increasingly important. A common method for predicting RNA secondary structure from sequence is free energy minimization. Therefore, accurate free energy contributions for every RNA secondary structure motif are necessary for accurate secondary structure predictions. Tandem mismatches are prevalent in naturally occurring sequences and are biologically important. A common method for predicting the stability of a sequence asymmetric tandem mismatch relies on the stabilities of the two corresponding sequence symmetric tandem mismatches [Mathews, D. H., Sabina, J., Zuker, M., and Turner, D. H. (1999) J. Mol. Biol. 288, 911-940]. To improve the prediction of sequence asymmetric tandem mismatches, the experimental thermodynamic parameters for the 22 previously unmeasured sequence symmetric tandem mismatches are reported. These new data, however, do not improve prediction of the free energy contributions of sequence asymmetric tandem mismatches. Therefore, a new model, independent of sequence symmetric tandem mismatch free energies, is proposed. This model consists of two penalties to account for destabilizing tandem mismatches, two bonuses to account for stabilizing tandem mismatches, and two penalties to account for A-U and G-U adjacent base pairs. This model improves the prediction of asymmetric tandem mismatch free energy contributions and is likely to improve the prediction of RNA secondary structure from sequence.     

Structural features and hydration of d(C-G-C-G-A-A-T-T-A-G-C-G); a double helix containing two G.A mispairs

Single crystal X-ray diffraction techniques have been used to characterise the molecular structure of the title compound to 2.5A resolution. The structure consists of ten standard Watson-Crick base pairs and two G.A mismatched base pairs. The purine-purine mismatches have guanine in the usual anti orientation with respect to the sugar and adenine in syn orientation. There are two hydrogen bonds formed between the mismatch bases, N-1 and O-6 of guanine with N-7 and N-6 of adenine respectively. The bulky purine-purine mismatches are accommodated with minor perturbation of the sugar-phosphate backbone. There is a slight improvement in base pair overlap at the mismatch sites. Details of the backbone conformation, base stacking interactions and hydration are presented and compared with those of the parent compound d(C-G-C-G-A-A-T-T-C-G-C-G).     

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     

Thermodynamic characterization of naturally occurring RNA single mismatches with G-U nearest neighbors

Due to their prevalence and roles in biological systems, single mismatches adjacent to G-U pairs are important RNA structural elements. Since there are only limited experimental values for the stability of single mismatches adjacent to G-U pairs, current algorithms using free energy minimization to predict RNA secondary structure from sequence assign predicted thermodynamic values to these types of single mismatches. Here, thermodynamic data are reported for frequently occurring single mismatches adjacent to at least one G-U pair. This experimental data can be used in place of predicted thermodynamic values in algorithms that predict secondary structure from sequence using free energy minimization. When predicting the thermodynamic contributions of previously unmeasured single mismatches, most algorithms apply the same thermodynamic penalty for an A-U pair adjacent to a single mismatch and a G-U pair adjacent to a single mismatch. A recent study, however, suggests that the penalty for a G-U pair adjacent to a tandem mismatch should be 1.2 +/- 0.1 kcal/mol, and the penalty for an A-U pair adjacent to a tandem mismatch should be 0.5 +/- 0.2 kcal/mol [Christiansen, M. E. and Znosko, B. M. (2008) Biochemistry 47, 4329-4336]. Therefore, the data reported here are combined with the existing thermodynamic dataset of single mismatches, and nearest neighbor parameters are derived for an A-U pair adjacent to a single mismatch (1.1 +/- 0.1 kcal/mol) and a G-U pair adjacent to a single mismatch (1.4 +/- 0.1 kcal/mol).     

Thermodynamics-structure relationship of single mismatches in RNA/DNA duplexes

Thermodynamics of 66 RNA/DNA duplexes containing single mismatches were measured by UV melting methods. Stability enhancements for rG. dT mismatches were the largest of all mismatches examined here, while rU.dG mismatches were not as stable. The methyl group on C5 of thymine enhanced the stability by 0.12 approximately 0.53 kcal mol(-)(1) depending on the identity of adjacent Watson-Crick base pairs, whereas the 2'-hydroxyl group in ribouridine stabilized the duplex by approximately 0.6 kcal mol(-)(1) regardless of the adjacent base pairs. Stabilities induced by the methyl group in thymine, the 2'-hydroxyl group of ribouridine, and an nucleotide exchange at rG.dT and rU.dG mismatches were found to be independent of each other. The order for the mismatch stabilities is rG.dT > rU. dG approximately rG.dG > rA.dG approximately rG.dA approximately rA. dC > rA.dA approximately rU.dT approximately rU.dC > rC.dA approximately rC.dT, although the identity of the adjacent base pairs slightly altered the order. The pH dependence stability and structural changes were suggested for the rA.dG but not for rG.dA mismatches. Comparisons of trinucleotide stabilities for G.T and G.U pairs in RNA, DNA, and RNA/DNA duplexes indicate that stable RNA/DNA mismatches exhibit a stability similar to RNA mismatches while unstable RNA/DNA mismatches show a stability similar to that of DNA mismatches. These results would be useful for the design of antisense oligonucleotides.     

Structure, dynamics, and thermodynamics of mismatched DNA oligonucleotide duplexes d(CCCAGGG)2 and d(CCCTGGG)2

The structures and hydrogen exchange properties of the mismatched DNA oligonucleotide duplexes d(CCCAGGG)2 and d(CCCTGGG)2 have been studied by high-resolution nuclear magnetic resonance. Both the adenine-adenine and thymine-thymine mismatches are intercalated in the duplexes. The structures of these self-complementary duplexes are symmetric, with the two strands in equivalent positions. The evidence indicates that these mismatches are not stably hydrogen bonded. The mismatched bases in both duplexes are in the anti conformation. The mismatched thymine nucleotide in d(CCCTGGG)2 is intercalated in the duplex with very little distortion of the bases or sugar-phosphate backbone. In contrast, the bases of the adenine-adenine mismatch in d(CCCAGGG)2 must tilt and push apart to reduce the overlap of the amino groups. The thermodynamic data show that the T-T mismatch is less destabilizing than the A-A mismatch when flanked by C-G base pairs in this sequence, in contrast to their approximately equal stabilities when flanked by A-T base pairs in the sequence d(CAAAXAAAG.CTTTYTTTG) where X and Y = A, C, G, and T [Aboul-ela, F., Koh, D., & Tinoco, I., Jr. (1985) Nucleic Acids Res. 13, 4811]. Although the mechanism cannot be determined conclusively from the limited data obtained, exchange of the imino protons with solvent in these destabilized heteroduplexes appears to occur by a cooperative mechanism in which half the helix dissociates.     

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.     

C5-(1-propynyl)-2'-deoxy-pyrimidines enhance mismatch penalties of DNA:RNA duplex formation

UV melting experiments show that C5-(1-propynyl)ation of seven pyrimidines to give a fully propynylated oligodeoxynucleotide (PrODN) heptamer increases the thermodynamic stability of six Watson-Crick paired DNA:RNA duplexes by 8.2 kcal/mol, on average, at 37 degrees C. About 2.5 kcal/mol of this enhancement is due to long-range cooperativity between the propynylated pyrimidines, Y(p)'s. On average, penalties for dU(p):rG, dC(p):rA, dU(p):rC, and dC(p):rC mismatches are enhanced by 2.9 kcal/mol in PrODN:RNA duplexes over those in unmodified duplexes. This results in penalties as large as 10 kcal/mol for a single mismatch. Removing a single propyne two base pairs away from a mismatch in a PrODN:RNA duplex eliminates the enhancement in specificity. Evidently, enhanced specificity is directly linked to long-range cooperativity between Y(p)'s. In most cases, the enhanced specificity is larger for internal than for terminal mismatches. PrODN:RNA duplexes are destabilized by full phosphorothioate backbone substitution to give S-PrODN:RNA duplexes. The S-PrODN:RNA duplexes retain enhanced mismatch penalties, however. These results provide insight for utilizing long-range cooperativity and enhanced specificity to improve nucleic acid based probe and drug design.     

Nearest neighbor parameters for inosine x uridine pairs in RNA duplexes

An enzyme family known as adenosine deaminases that act on RNA (ADARs) catalyzes adenosine deamination in RNA. ADARs act on RNA that is largely double-stranded and convert adenosine to inosine, resulting, in many cases, in an I x U pair. Thermodynamic parameters derived from optical melting studies are reported for a series of 14 oligoribonucleotides containing single I x U pairs adjacent to Watson-Crick pairs. In order to determine unique linearly independent nearest neighbor parameters for I x U pairs, four duplexes containing 3'-terminal I x U pairs and four duplexes containing 5'-terminal I x U pairs have also been thermodynamically characterized. This data was combined with previously published data of seven duplexes containing internal, terminal, or tandem I x U pairs from Strobel et al. [Strobel, S. A., Cech, T. R., Usman, N., and Beigelman, L. (1994) Biochemistry 33, 13824-13838] and Serra et al. [Serra, M. J., Smolter, P. E., and Westhof, E. (2004) Nucleic Acids Res. 32, 1824-1828]. On average, a duplex with an internal I x U pair is 2.3 kcal/mol less stable than the same duplex with an A-U pair, however, a duplex with a terminal I x U pair is 0.8 kcal/mol more stable than the same duplex with an A-U pair. Although isosteric with a G-U pair, on average, a duplex with an internal I x U pair is 1.9 kcal/mol less stable than the same duplex with a G-U pair, however, a duplex with a terminal I x U pair is 0.9 kcal/mol more stable than the same duplex with a G-U pair. Duplexes with tandem I x U pairs are on average 5.9 and 3.8 kcal/mol less stable than the same duplex with tandem A-U or tandem G-U pairs, respectively. Using the combined thermodynamic data and a complete linear least-squares fitting routine, nearest neighbor parameters for all nearest neighbor combinations of I x U pairs and an additional parameter for terminal I x U pairs have been derived.     

Structural Features and Hydration of d(C-G-C-G-A-A-T-T-A-G-C-G); a Double Helix Containing Two G.A Mispairs

Abstract

Single crystal X-ray diffraction techniques have been used to characterise the molecular structure of the title compound to 2.5Å resolution. The structure consists of ten standard Watson-Crick base pairs and two G.A mismatched base pairs. The purine-purine mismatches have guanine in the usual anti orientation with respect to the sugar and adenine in syn orientation. There are two hydrogen bonds formed between the mismatch bases, N-l and 0–6 of guanine with N-7 and N-6 of adenine respectively. The bulky purine-purine mismatches are accommodated with minor perturbation of the sugar-phosphate backbone. There is a slight improvement in base pair overlap at the mismatch sites. Details of the backbone conformation, base stacking interactions and hydration are presented and compared with those of the parent compound d(C-G-C-G-A-A-T-T-C-G-C-G).     

Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.     

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)