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.     

Thermodynamic characterization of single mismatches found in naturally occurring RNA

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

Thermodynamic contributions of single internal rA·dA, rC·dC, rG·dG and rU·dT mismatches in RNA/DNA duplexes

The thermodynamic contributions of rA·dA, rC·dC, rG·dG and rU·dT single internal mismatches were measured for 54 RNA/DNA duplexes in a 1 M NaCl buffer using UV absorbance thermal denaturation. Thermodynamic parameters were obtained by fitting absorbance versus temperature profiles using the curve-fitting program Meltwin. The weighted average thermodynamic data were fit using singular value decomposition to determine the eight non-unique nearest-neighbor parameters for each internal mismatch. The new parameters predict the ΔG°₃₇, ΔH° and melting temperature (T_m) of duplexes containing these single mismatches within an average of 0.33 kcal/mol, 4.5 kcal/mol and 1.4°C, respectively. The general trend in decreasing stability for the single internal mismatches is rG·dG > rU·dT > rA·dA > rC·dC. The stability trend for the base pairs 5′ of the single internal mismatch is rG·dC > rC·dG > rA·dT > rU·dA. The stability trend for the base pairs 3′ of the single internal mismatch is rC·dG > rG·dC >> rA·dT > rU·dA. These nearest-neighbor values are now a part of a complete set of single internal mismatch thermodynamic parameters for RNA/DNA duplexes that are incorporated into the nucleic acid assay development software programs Visual oligonucleotide modeling platform (OMP) and ThermoBLAST.     

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.     

Thermodynamics of single mismatches in RNA duplexes

The thermodynamic properties and structures of single mismatches in short RNA duplexes were studied in optical melting and imino proton NMR experiments. The free energy increments at 37 degrees C measured for non-GU single mismatches range from -2.6 to 1.7 kcal/mol. These increments depend on the identity of the mismatch, adjacent base pairs, and the position in the helix. UU and AA mismatches are more stable close to a helix end, but GG mismatch stability is essentially unaffected by the position in the helix. Approximations are suggested for predicting stabilities of single mismatches in short RNA duplexes.     

Thermodynamics of RNA duplexes modified with unlocked nucleic acid nucleotides

Thermodynamics provides insights into the influence of modified nucleotide residues on stability of nucleic acids and is crucial for designing duplexes with given properties. In this article, we introduce detailed thermodynamic analysis of RNA duplexes modified with unlocked nucleic acid (UNA) nucleotide residues. We investigate UNA single substitutions as well as model mismatch and dangling end effects. UNA residues placed in a central position makes RNA duplex structure less favourable by 4.0–6.6 kcal/mol. Slight destabilization, by ∼0.5–1.5 kcal/mol, is observed for 5′- or 3′-terminal UNA residues. Furthermore, thermodynamic effects caused by UNA residues are extremely additive with ΔG°₃₇ conformity up to 98%. Direct mismatches involving UNA residues decrease the thermodynamic stability less than unmodified mismatches in RNA duplexes. Additionally, the presence of UNA residues adjacent to unpaired RNA residues reduces mismatch discrimination. Thermodynamic analysis of UNA 5′- and 3′-dangling ends revealed that stacking interactions of UNA residues are always less favourable than that of RNA residues. Finally, circular dichroism spectra imply no changes in overall A-form structure of UNA–RNA/RNA duplexes relative to the unmodified RNA duplexes.     

The thermodynamics of 3'-terminal pyrene and guanosine for the design of isoenergetic 2'-O-methyl-RNA-LNA chimeric oligonucleotide probes of RNA structure

To facilitate design of short isoenergetic hybridization probes for RNA, we report the influence of adding 5'- or 3'-terminal 2'-O-methylguanosine (GM), LNA-guanosine (GL), or 3'-terminal pyrene pseudo-nucleotide (PPN) on the thermodynamic stability of 2'-O-methyl-RNA/RNA (2'-O-Me-RNA/RNA) duplexes with sequences 5'CMGMGMCMAM/3'AAXGCCGUXAA, where X is A, C, G, or U. A 3'-terminal GM or GL added to the 2'-O-Me-RNA strand to form a G-A, G-G or G-U mismatch enhances thermodynamic stability (DeltaDeltaG degrees 37) of the 2'-O-Me-RNA/RNA duplexes on average by 0.7 and 1.5 kcal/mol, respectively. A 3'-terminal GM or GL in a GM-C or GL-C pair stabilizes the 2'-O-Me-RNA/RNA duplex by 2.6 and 3.4 kcal/mol, respectively. A 5'-terminal GM or GL in a G-A or G-G mismatch provided less stabilization in comparison with a 3'-terminal G-A or G-G mismatch, but more stabilization in a G-C or G-U pair. In contrast to guanosine derivatives, pyrene residue (P) as PPN at the 3'-terminal position enhances thermodynamic stability of the 2'-O-Me-RNA/RNA duplexes on average by 2.3 +/- 0.1 kcal/mol, relatively independent of the type of ribonucleotide placed in the opposite strand. The thermodynamic data can be applied to design 2'-O-Me-RNA/RNA duplexes with enhanced thermodynamic stability that is also sequence independent. This is useful for design of hybridization probes to interrogate RNA structure and/or expression by microarray and other methods.     

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).     

Effects of GA mismatches on the structure and thermodynamics of RNA internal loops

UV melting, CD and NMR studies indicate rGCGAGCG and rGCAGGCG from unusually stable duplexes of type a and b. The observed delta G degree 37's in 1 M NaCl are -6.7 and -6.3 kcal/mol, respectively. For the related duplex, c, delta G degree 37 is -4.2 kcal/mol. The predicted delta G degree 37 from nearest-neighbor parameters (formula; see text) for all three duplexes is -4.7 kcal/mol (Freier, S.M., Kierzek, R., Jaeger, J.A., Sugimoto, N., Caruthers, M.H., Neilson, T., & Turner, D.H. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373-9377). The results suggest a special interaction in the duplexes containing GA mismatches. Presumably, this is hydrogen bonding between G and A. While the thermodynamics for (rGCGAGCG)2 and (rGCAGGCG)2 are similar, CD and the imino region of the proton NMR spectra indicate their structures are different. In particular, (rGCAGGCG)2 exhibits a CD spectrum typical of A-form geometry with a weak negative band at 280 nm. In contrast, the CD spectrum for (rGCGAGCG)2 has an intense positive band at 285 nm. The NMR spectrum of (rGCAGGCG)2 has a resonance corresponding to a hydrogen-bonded GA mismatch, while for (rGCGAGCG)2 no hydrogen-bonded imino proton is observed for the mismatch. The glycosidic torsion angles of the bases in the GA mismatches of (rGCAGGCG)2 and (rCGCAGGCG)2 are anti. Duplexes of type d, where X is A, G, or U, are more stable than e, and the stability differences are similar to those (formula; see text) observed for f versus g. Thus, 3'-dangling ends in this system make contributions to duplex stability that are similar to contributions observed with fully paired duplexes.     

Thermodynamics of internal C.T mismatches in DNA.

Thermodynamics of 23 oligonucleotides with internal single C.T mismatches were obtained by measuring UV absorbance as a function of temperature. Results from these 23 duplexes were combined with three measurements from the literature to derive nearest-neighbor thermodynamic parameters for seven linearly independent trimer sequences with internal C.T mismatches. The data show that the nearest-neighbor model is adequate for predicting thermodynamics of oligonucleotides with internal C.T with average deviations for Delta G degrees37, Delta H degrees, Delta S degrees and T m of 6.4%, 9.9%, 10.6%, and 1.9 degreesC respectively. C.T mismatches destabilize the duplex in all sequence contexts. The thermodynamic contribution of C. T mismatches to duplex stability varies weakly depending on the orientation of the mismatch and its context and ranges from +1.02 kcal/mol for GCG/CTC and CCG/GTC to +1.95 kcal/mol for TCC/ATG.     

Effect of base stacking on the relative thermodynamic stability of oligonucleotide complexes: a spectroscopic study

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be -1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Determining the origin of the stabilization of DNA by 5-aminopropynylation of pyrimidines

DNA duplexes are stabilized by aminopropynyl modification of pyrimidines at the 5 position. A combination of thermodynamic analyses as a function of ionic strength, NMR, and molecular modeling has been applied to determine the origin of the stabilization. UV melting studies of a dodecamer bearing one, two, or three nonadjacent modified dU and dC and of a single dU(8) in the Dickerson-Drew dodecamer revealed that the modifications are essentially additive in terms of T(m), DeltaG, and DeltaH, and there is little difference between dU and dC. The free energy change was parsed into electrostatic and nonelectrostatic components, which showed a significant contribution from charge interactions at physiological ionic strength but also a nonelectrostatic contribution that arises in part from hydration. NMR spectroscopy of the modified Dickerson-Drew dodecamer revealed that the conformation of the duplexes is not significantly altered by the modifications, though (31)P NMR shows that the positive charge may affect ionic interactions with the oxygen atoms of the neighboring phosphates. The modified duplex showed significant hydration in both major and minor grooves. The single strands were also analyzed by NMR, which showed evidence of significant stacking interactions in the modified oligonucleotide. Parsing the energy contribution has shown that electrostatics and hydration can produce substantial increases in thermodynamic stability without significant changes in the conformation of the duplex state. These considerations have significance for the design of oligonucleotides used for hybridization.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Effect of Base Stacking on the Relative Thermodynamic Stability of Oligonucleotide Complexes: A Spectroscopic Study

Abstract

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be ?1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Mismatch recognition and uracil excision provide complementary paths to both Ig switching and the A/T-focused phase of somatic mutation

AID-mediated deamination of dC residues within the immunoglobulin locus generates dU:dG lesions whose resolution leads to class-switch recombination and somatic hypermutation. The dU:dG pair is a mismatch and comprises a base foreign to DNA and is, thus, recognized by proteins from both base excision (uracil-DNA glycosylase, UNG) and mismatch recognition (MSH2/MSH6) pathways. Strikingly, while antibody diversification is perturbed by single deficiency in either UNG or MSH2, combined UNG/MSH2 deficiency leads to a total ablation both of switch recombination and of IgV hypermutation at dA:dT pairs. The initiating dU:dG lesions appear not to be recognized and are simply replicated over. The results indicate that the major pathway for switch recombination occurs through uracil excision with mismatch recognition of dU:dG providing a backup; the second phase of hypermutation (essentially introducing mutations solely at dA:dT pairs) is triggered by mismatch recognition of the dU:dG lesion with uracil excision providing a backup.     

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.     

Studies on nucleic acid interactions. I. Stabilities of mini-duplexes (dG2A4XA4G2-dC2T4YT4C2) and self-complementary d(GGGAAXYTTCCC) containing deoxyinosine and other mismatched bases.

The thermal stability of DNA duplexes containing deoxyinosine in a pairing position in turn with each of the four major deoxynucleotides has been investigated by measuring ultraviolet-absorbance at different temperatures. d(G2A4 X A4G2) and d(C2T4YT4C2) were prepared by the solid-phase phosphotriester method. When X is deoxyinosine, the Tm values of the duplexes are in the order Y = dC greater than dA greater than dG greater than dT greater than dU. The Tm of other duplexes containing dG, dA and dT at X were also measured. Self-complementary duplexes d(GGGAAINTTCCC) showed the same order of stability with N being dC, dA, dG and dT. Thermal stabilities of duplexes containing dG instead of dI were compared with other matched and mismatched duplexes. The Tm values of sequence isomers containing purine-pyrimidine combinations were compared. Self-complementary duplexes containing G-C and A-T in the central positions showed Tm values ca. 10 degrees higher than those containing C-G and T-A in the same positions. Thermodynamic parameters and circular dichroism spectra of these oligonucleotides were compared.     

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.