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.     

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

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.     

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.     

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

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.     

Thermodynamic stability of DNA tandem mismatches

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.     

Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A.A, C.C, G.G, and T.T mismatches

Thermodynamic measurements are reported for 51 DNA duplexes with A.A, C.C, G.G, and T.T single mismatches in all possible Watson-Crick contexts. These measurements were used to test the applicability of the nearest-neighbor model and to calculate the 16 unique nearest-neighbor parameters for the 4 single like with like base mismatches next to a Watson-Crick pair. The observed trend in stabilities of mismatches at 37 degrees C is G.G > T.T approximately A.A > C.C. The observed stability trend for the closing Watson-Crick pair on the 5' side of the mismatch is G.C >/= C.G >/= A.T >/= T.A. The mismatch contribution to duplex stability ranges from -2.22 kcal/mol for GGC.GGC to +2.66 kcal/mol for ACT.ACT. The mismatch nearest-neighbor parameters predict the measured thermodynamics with average deviations of DeltaG degrees 37 = 3.3%, DeltaH degrees = 7. 4%, DeltaS degrees = 8.1%, and TM = 1.1 degrees C. The imino proton region of 1-D NMR spectra shows that G.G and T.T mismatches form hydrogen-bonded structures that vary depending on the Watson-Crick context. The data reported here combined with our previous work provide for the first time a complete set of thermodynamic parameters for molecular recognition of DNA by DNA with or without single internal mismatches. The results are useful for primer design and understanding the mechanism of triplet repeat diseases.     

Thermodynamic characterization of tandem mismatches found in naturally occurring RNA

Although all sequence symmetric tandem mismatches and some sequence asymmetric tandem mismatches have been thermodynamically characterized and a model has been proposed to predict the stability of previously unmeasured sequence asymmetric tandem mismatches [Christiansen,M.E. and Znosko,B.M. (2008) Biochemistry, 47, 4329–4336], experimental thermodynamic data for frequently occurring tandem mismatches is lacking. Since experimental data is preferred over a predictive model, the thermodynamic parameters for 25 frequently occurring tandem mismatches were determined. These new experimental values, on average, are 1.0 kcal/mol different from the values predicted for these mismatches using the previous model. The data for the sequence asymmetric tandem mismatches reported here were then combined with the data for 72 sequence asymmetric tandem mismatches that were published previously, and the parameters used to predict the thermodynamics of previously unmeasured sequence asymmetric tandem mismatches were updated. The average absolute difference between the measured values and the values predicted using these updated parameters is 0.5 kcal/mol. This updated model improves the prediction for tandem mismatches that were predicted rather poorly by the previous model. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA duplexes containing tandem mismatches, and, furthermore, should allow for improved prediction of secondary structure from sequence.     

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.     

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.     

Thermodynamics of RNA internal loops with a guanosine-guanosine pair adjacent to another noncanonical pair

Thermodynamic parameters measured by optical melting are reported for formation of RNA duplexes containing tandem noncanonical pairs with at least one guanosine-guanosine (GG) pair. For selected sequences, imino proton NMR provides evidence that the desired duplex forms and that the structure of a GG pair adjacent to a noncanonical pair depends on context. A GG pair next to a different noncanonical pair is more stable than expected from measurements of adjacent GG pairs. This is likely due to an unfavorable stacking interaction between adjacent GG pairs, where areas of high negative charge probably overlap. The results suggest a model where tandem noncanonical pairs closed by two GC pairs are assigned the following free energy increments at 37 degrees C: 0.8 kcal/mol for adjacent GG pairs, 1.0 kcal/mol for GG next to UU, and -0.3 kcal/mol for all others. These values are adjusted by 0.65 kcal/mol for each closing AU pair.     

Recognition of G-U mismatches by tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III).

The coordination complex tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III) [Rh(DIP)3(3+)], which promotes RNA cleavage upon photoactivation, has been shown to target specifically guanine-uracil (G-U) mismatches in double-helical regions of folded RNAs. Photoactivated cleavage by Rh(DIP)3(3+) has been examined on a series of RNAs that contain G-U mismatches, yeast tRNA(Phe) and yeast tRNA(Asp), as well as on 5S rRNAs from Xenopus oocytes and Escherichia coli. In addition, a "microhelix" was synthesized, which consists of seven base pairs of the acceptor stem of yeast tRNA(Phe) connected by a six-nucleotide loop and contains a mismatch involving residues G4 and U69. A U4.G69 variant of this sequence was also constructed, and cleavage by Rh(DIP)3(3+) was examined. In each of these cases, specific cleavage is observed at the residue which lies to the 3'-side of the wobble-paired U; some cleavage by the rhodium complex is also evident in several structured RNA loops. The remarkable site selectivity for G-U mismatches within double-helical regions is attributed to shape-selective binding by the rhodium complex. This binding furthermore depends upon the orientation of the G-U mismatch, which produces different stacking interactions between the G-U base pair with the Watson-Crick base pair following it on the 5'-side of U compared to the Watson-Crick pair preceding it on the 3'-side of U. Rh(DIP)3(3+) therefore serves as a unique probe of G-U mismatches and may be useful both as a model and in probing RNA-protein interactions as well as in identifying G-U mismatches within double-helical regions of folded RNAs.     

Base pairing in wheat germ ribosomal 5S RNA as measured by ultraviolet absorption, circular dichroism, and Fourier-transform infrared spectrometry

Ultraviolet absorption (UV) and circular dichroism (CD) spectra of wheat germ 5S RNA, when compared to tRNAPhe, indicate a largely base-paired and base-stacked helical structure, containing up to 36 base pairs. Fourier-transform infrared (FT-IR) spectra of tRNAPhe and wheat germ ribosomal 5S RNA have been acquired at 30 and 90 degrees C. From the difference of the FT-IR spectra between 90 and 30 degrees C, the number of base pairs in both RNAs was determined by modification of a previously published procedure [Burkey, K. O., Marshall, A. G., & Alben, J. O. (1983) Biochemistry 22, 4223-4229]. The base-pair composition and total base-pair number from FT-IR data are now consistent for the first time with optical (UV, CD, Raman) and NMR results for ribosomal 5S RNA. Without added Mg2+, tRNAPhe gave 18 +/- 2 base pairs [7 A-U and 11 G-C], in good agreement with the number of secondary base pairs from X-ray crystallography [8 A-U, 12 G-C, and 1 G-U]. Within the 10% precision of the FT-IR method, wheat germ 5S RNA exhibits essentially the same number of base pairs [14 A-U, 17 G-C, and 5 G-U; for a total of 36] in the absence of Mg2+ as in the presence of Mg2+ [14 A-U, 18 G-C, and 3 G-U; for a total of 35], in agreement with the UV hyperchromism estimate of G-C/(A-U + G-C) = 0.58.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

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.     

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.     

Free energy contributions of G.U and other terminal mismatches to helix stability

Thermodynamic parameters of helix formation were measured spectroscopically for seven hexaribonucleotides containing a GC tetramer core and G.U or other terminal mismatches. The free energies of helix formation are compared with those for the tetramer core alone and with those for the hexamer with six Watson-Crick base pairs. In 1 M NaCl, at 37 degrees C, the free energy of a terminal G.U mismatch is about equal to that of the corresponding A.U pair. Although other terminal mismatches studied add between -1.0 and -1.6 kcal/mol to delta G0 37 for helix formation, all are less stable than the corresponding Watson-Crick pairs. Comparisons of the stability increments for terminal G.U mismatches and G.C pairs suggest when stacking is weak the additional hydrogen bond in the G.C pair adds roughly -1 kcal/mol to the favorable free energy of duplex formation.