The influence of locked nucleic acid residues on the thermodynamic pro
http://www.100md.com
《核酸研究医学期刊》
1Department of Chemistry and Department of Pediatrics, University of Rochester RC Box 270216, Rochester, NY 14627-0216, USA 2Institute of Bioorganic Chemistry, Polish Academy of Sciences 60-714 Poznan, Noskowskiego 12/14, Poland 3Center for Molecular Pediatric Research, University of Rochester School of Medicine and Dentistry Rochester, New York, 14642, USA 4Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry Rochester, New York, 14642, USA
*To whom correspondence should be addressed. Tel: +1 585 275 3207; Fax: +1 585 276 0205; Email: turner@chem.rochester.edu
ABSTRACT
The influence of locked nucleic acid (LNA) residues on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes is reported. Optical melting studies indicate that LNA incorporated into an otherwise 2'-O-methyl RNA oligonucleotide usually, but not always, enhances the stabilities of complementary duplexes formed with RNA. Several trends are apparent, including: (i) a 3' terminal U LNA and 5' terminal LNAs are less stabilizing than interior and other 3' terminal LNAs; (ii) most of the stability enhancement is achieved when LNA nucleotides are separated by at least one 2'-O-methyl nucleotide; and (iii) the effects of LNA substitutions are approximately additive when the LNA nucleotides are separated by at least one 2'-O-methyl nucleotide. An equation is proposed to approximate the stabilities of complementary duplexes formed with RNA when at least one 2'-O-methyl nucleotide separates LNA nucleotides. The sequence dependence of 2'-O-methyl RNA/RNA duplexes appears to be similar to that of RNA/RNA duplexes, and preliminary nearest-neighbor free energy increments at 37°C are presented for 2'-O-methyl RNA/RNA duplexes. Internal mismatches with LNA nucleotides significantly destabilize duplexes with RNA.
INTRODUCTION
Understanding the thermodynamics of nucleic acid duplexes is important for many reasons. For example, such knowledge facilitates design of ribozymes (1), antisense and RNAi oligonucleotides (2–9), diagnostic probes including those employed on microarrays (10–23) and structures useful for nanotechnology (24–27). Many modified residues have been developed for such applications. Examples include propynylated bases (28–30), peptide nucleic acids (5,31–33), N3'–P5' phosphoramidates (34–38) and 2'-O-alkyl RNA (39–43). A modification that is particularly stabilizing in DNA and RNA duplexes (44–51) is a methyl bridge between the 2' oxygen and 4' carbon of ribose to form a ‘locked nucleic acid’ or LNA as shown in Figure 1. McTigue et al. (48) have shown that the enhanced stability due to a single LNA residue in a DNA duplex can be predicted from a nearest-neighbor model.
Figure 1 Covalent structures of sugars used.
Hybridization of oligonucleotides to RNA is important for applications, such as antisense therapeutics (4,8,21,46,52–54), diagnostics (32,33,42,55), profiling gene expression with microarrays (18–20,56), identifying bands by Northern blots of gels (57,58) and probing RNA structure (1,3,15,59–61). Oligonucleotides with 2'-O-alkyl modifications can be particularly useful for these applications because they are easily synthesized (39,43), chemically stable and bind relatively tightly to RNA (39–42). However, for many applications, it is desirable to modulate the binding affinity. For example, sequence independent duplex stabilities would benefit applications that involve multiplex detection, such as microarrays. Here, we show that introduction of LNA into 2'-O-methyl RNA oligonucleotides can increase stabilities of 2'-O-methyl RNA/RNA hybrid duplexes and that the enhancements in stability can usually be predicted with a simple model.
MATERIALS AND METHODS
General methods
High-performance liquid chromatography (HPLC) was performed on a Hewlett Packard series 1100 HPLC with a reverse-phase Supelco RP-18 column (4.6 x 250 mm). Mass spectra were obtained on an LC MS Hewlett Packard series 1100 MSD with API-ES detector or on an AMD 604/402. Thin-layer chromatography (TLC) was carried out on Merck 60 F254 TLC plates with the mixture 1-propanol/aqueous ammonia/water = 55:35:10 (v/v/v).
Synthesis and purification of oligonucleotides
Oligoribonucleotides were synthesized on an Applied Biosystems DNA/RNA synthesizer, using ?-cyanoethyl phosphoramidite chemistry (62). For synthesis of standard RNA oligonucleotides, the commercially available phosphoramidites with 2'-O-tertbutyldimethylsilyl groups were used (Glen Research). For synthesis of 2'-O-methyl RNA oligonucleotides, the 3'-O-phosphoramidites of 2'-O-methylnucleotides were used (Glen Research and Proligo). The 3'-O-phosphoramidites of LNA nucleotides were synthesized according to the published procedures with some minor modifications (44,47,63). The details of deprotection and purification of oligoribonucleotides were described previously (64).
UV melting
Oligonucleotides were melted in buffer containing 100 mM NaCl, 20 mM sodium cacodylate, 0.5 mM Na2EDTA, pH 7.0. The relatively low NaCl concentration kept melting temperatures in the reasonable range even when there were multiple LNA substitutions. Oligonucleotide single-strand concentrations were calculated from absorbencies above 80°C and single-strand extinction coefficients were approximated by a nearest-neighbor model (65,66). It was assumed that 2'-O-methyl RNA and RNA strands with identical sequences have identical extinction coefficients. Absorbancy versus temperature melting curves were measured at 260 nm with a heating rate of 1°C/min from 0 to 90°C on a Beckman DU 640 spectrophotometer with a water cooled thermoprogrammer. Melting curves were analyzed and thermodynamic parameters were calculated from a two-state model with the program MeltWin 3.5 (67). For almost all sequences, the H° derived from versus ln (CT/4) plots is within 15% of that derived from averaging the fits to individual melting curves, as expected if the two-state model is reasonable.
Parameter fitting
Free energy parameters for predicting stabilities of 2'-O-methyl RNA/RNA and 2'-O-methyl RNA–LNA/RNA duplexes with the Individual Nearest-Neighbor Hydrogen Bonding (INN-HB) model (64) were obtained by multiple linear regression with the program Analyse-it v.1.71 (Analyse-It Software, Ltd, Leeds, England; www.analyse-it.com) which expands Microsoft Excel. Analyse-It was also used to obtain parameters for enhancement of stabilities of 2'-O-methyl RNA/RNA duplexes by substitution of LNA nucleotides internally and/or at the 3' end when the LNAs are separated by at least one 2'-O-methyl nucleotide. Results from versus ln (CT/4) plots were used as the data for the calculations.
RESULTS
Figures 2 and 3 show typical data from optical melting curves, and Table 1 lists the thermodynamic parameters for the helix to coil transition with either no or one LNA nucleotide in the primarily 2'-O-methyl strand of a hybrid with a Watson–Crick complementary RNA strand.
Figure 2 Representative UV melting curves 5'AMCMUMALCMCMAM/5'UGGCAGU (1, 52 μM), 5'AMCMUMAMCMCMAM/5'UGGUAGU (2, 46 μM), 5'AMCMUMALCMCMAM/5'UGGUAGU (3, 51 μM) and 5'AMCLUMALCMCLAM/5'UGGUAGU (4, 51 μM). Note that sequence 1 has an ALC mismatch in the middle.
Figure 3 Reciprocal melting temperature versus ln (CT/4) for 5'AMCMUMALCMCMAM/5'UGGCAGU (solid squares), 5'AMCMUMAMCMCMAM/5'UGGUAGU (solid circles), 5'AMCMUMALCMCMAM/5'UGGUAGU (solid stars) and 5'AMCLUMALCMCLAM/5'UGGUAGU (solid triangles) in 100 mM NaCl, 20 mM sodium cacodylate, 0.5 mM Na2EDTA, pH 7.0. The thermodynamics of duplex formation were obtained by fitting the data points to the equation, = (R/H°) ln (CT/4) + S°/H°. Here, is the inverse melting temperature in kelvin, R is the gas constant, 1.987 cal K–1 mol–1, CT is the total oligonucleotide strand concentration, and both strands have the same concentration.
Table 1 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without a single LNA substitutiona
Single LNA substitutions at the 5' end of heptamer duplexes have little effect on stability
The effects of single LNA substitutions at the 5' end of the 2'-O-methyl strand were studied in duplexes of the form, 5'NCMUMAMCMCMAM/3'r(QGAUGGU), where superscript M denotes a 2'-O-methyl sugar, N is A, C, G, or U with a 2'-O-methyl or LNA sugar, r denotes ribose sugars, and Q is the Watson–Crick complement to N. As summarized in Table 1, 5' terminal LNA substitutions make duplex stability more favorable by 0.3–0.6 kcal/mol at 37°C with an average enhancement of 0.45 kcal/mol. Thus, 5' terminal LNA substitutions increase the binding constant for duplex formation by 2-fold at 37°C.
The effects of single LNA substitutions at the 3' ends of heptamer duplexes is idiosyncratic
The effects of single LNA substitutions at the 3' end of the 2'-O-methyl strand was studied in duplexes of the form, 5'AMCMUMAMCMCMN/3'r(UGAUGGQ) and in 5'AMCMUMAMCMGMUM/L/3'r(UGAUGCA) (Table 1). If N is A, C or G, then LNA substitutions have similar effects. On average, an LNA substitution makes duplex stability more favorable by 1.2 kcal/mol at 37°C. In the two sequences with a 3' terminal LNA U on the 2'-O-methyl strand, duplex stability is, however, affected little, averaging a destabilization of 0.08 kcal/mol at 37°C. In both cases, the terminal U is preceded by a GC pair, but both orientations of the GC pair give similar destabilization upon LNA substitution at the 3' terminal U.
Single LNA substitutions in the interior of AMCMUMAMCMCMAM enhance the stability of the duplex formed with its complementary RNA by 1.4 kcal/mol
The effect of interior position on the free energy increment for a single LNA substitution for a 2'-O-methyl RNA was studied for the duplex 5'AMCMUMAMCMCMAM/3'r(UGAUGGU). As summarized in Table 1, a single interior LNA substitution makes duplex stability more favorable by 1.2–1.7 kcal/mol at 37°C, with an average of 1.4 kcal/mol. This corresponds to roughly a 10-fold increase in binding constant. Thus, interior and 3' terminal LNA substitutions usually improve binding more than 5' terminal LNA substitutions.
The effects of single LNA substitutions as a function of sequence context in the middle of heptamer duplexes
The effects from single LNA substitutions in the middle of heptamer duplexes were studied as a function of the LNA base and the 3' adjacent 2'-O-methyl nucleotide in duplexes of the form, 5'AMCMUMNVMCMAM/3'r(UGAQWGU), where VM is 2'-O-methyl A, C, G or U, and W is the ribose Watson–Crick complement to V. The effects on the for duplex formation of substituting LNA for 2'-O-methyl at position N are summarized in Table 1. For 13 of 16 sequences, the LNA substitution makes duplex stability more favorable by 1.0–1.5 kcal/mol at 37°C, with an average enhancement of 1.3 kcal/mol. The enhancement for the other three sequences averages 2.1 kcal/mol at 37°C.
The dependence on the 5' nearest-neighbor nucleotide of effects from substituting UL for UM was studied in duplexes of the form, 5'AMCMVMUL/MCMCMAM/3'r(UGWAGGU) (Table 1). For three of the four sequences, the range for enhancement of duplex stability is 0.8–1.3 kcal/mol and the average is 1.0 kcal/mol. The exception is for the nearest neighbor 5'GMUL/3'r(CA), where LNA substitution destabilizes the duplex by 0.2 kcal/mol.
The results for single LNA substitutions in other contexts provide additional insight into the dependence of free energy increments on the 5' nearest neighbor 2'-O-methyl nucleotide. In particular, the Watson–Crick complementary duplexes with 5'AMCLUMAMCMCMAM and 5'AMCMUMCLUMCMAM strands have a 5'CLUM/3'r(GA) nearest neighbor that is preceded by AM and UM, respectively. The LNA substitution for 2'-O-methyl in these cases results in duplex formation being more stable by 1.29 and 1.48 kcal/mol, respectively, at 37°C. Similarly, the duplexes with 5'AMGLUMAMCMCMAM and 5'AMCMUMGLUMCMAM each have a 5'GLUM/3'r(CA) nearest neighbor that is preceded by AM and UM, respectively. In both cases, the LNA substitution enhances duplex stability by 1.14 kcal/mol at 37°C. Thus, for seven duplexes, the enhanced stability from an LNA substitution is relatively independent of the nearest-neighbor nucleotide 5' to the LNA. The one exception is for the nearest neighbor 5'GMUL/3'r(CA). Interestingly, this nearest-neighbor combination is also destabilized by LNA substitution at a 3' terminal U (Table 1). Evidently, an LNA substitution in the middle of a 2'-O-methyl strand usually affects heteroduplex stability with an RNA strand by about the same amount as an LNA substitution at a 3' terminus.
The effects of LNA substitutions are approximately additive when LNA nucleotides are spaced by at least one 2'-O-methyl nucleotide
Table 2 contains thermodynamic parameters measured for duplexes having more than one LNA substitution and Table 3 compares the stabilities at 37°C with those predicted from four simple models. The first model, labeled ‘additivity’, predicts the for duplex formation in the 5'ACUACCA/3'UGAUGGU series by adding the free energy increments measured for single LNA substitutions in the same context to the for duplex formation in the absence of LNA nucleotides. The second model predicts the (kcal/mol) for duplex formation with the following equation as deduced from fitting the data in Tables 1 and 2 excluding sequences with adjacent LNAs.
(1)
Here, (2'-O-MeRNA/RNA) is the free energy change at 37°C for duplex formation in the absence of any LNA nucleotides, n5'tL is the number of 5' terminal LNAs, niAL/UL and niGL/CL are the number of internal LNAs in AU and GC pairs, respectively, n3'tU and n3'tAL/CL/GL are the number of 3' terminal LNAs that are U or not U, respectively. Both methods that use experimental data for (2'-O-MeRNA/RNA) provide reasonable predictions that are within 1 kcal/mol of the measured value (Table 3). Two other methods that use nearest-neighbor models to approximate (2'-O-MeRNA/RNA) provide somewhat less accurate, but still reasonable predictions as described below. The duplex with the worst prediction, 5'GMULUMCLGMGL/3'CAAGCC has a 5'GMUL/3'CA nearest neighbor, consistent with this motif being unusually unstable by 1.2 kcal/mol. Thus, it is likely that the of Equation 1 should be made less favorable by 1.2 kcal/mol for every internal 5'GMUL/3'CA nearest neighbor in a duplex. Evidently, the effects of multiple LNA substitutions are approximately additive when the LNAs are spaced by at least 1 nt.
Table 2 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without multiple LNA substitutionsa
Table 3 Comparison of measured and predicted stabilities
The data may also be fit to a nearest-neighbor model containing 30 of the LNA enhancement parameters associated with duplexes of RNA strands bound to 2'-O-methyl RNA/LNA chimeras. These parameters are listed in Supplementary Material. The number of occurrences for each nearest neighbor is limited, however, so the values are only roughly determined.
Predictions for RNA/RNA duplexes at 1 M NaCl can be used to approximate stabilities of 2'-O-methyl RNA/RNA duplexes at 0.1 M NaCl
The stabilities of RNA/RNA duplexes at 37°C and 1 M NaCl are predicted well by an Independent Nearest-Neighbor Hydrogen Bonding (INN-HB) model (64). In this model, the stability of an RNA/RNA duplex is approximated by:
(2)
Here, is the free energy change for initiating a helix; each is the free energy increment of the jth type nearest neighbor (see Table 4) with nj occurrences in the sequence; mterm-AU is the number of terminal AU pairs; is the free energy increment per terminal AU pair; Gsym is 0.43 kcal/mol at 37°C for self-complementary duplexes and 0 for non-self-complementary duplexes.
Table 4 Preliminary nearest-neighbor free energy parameters (kcal/mol at 37°C) for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl compared with parameters for RNA/RNA duplexes in 1 M NaCl
Similar sequence dependent parameters may also be applicable to 2'-O-methyl RNA/RNA heteroduplexes because they are expected to have A-form conformations similar to those of RNA/RNA homoduplexes (68). This was tested by comparing the predicted stabilities of RNA/RNA duplexes in 1 M NaCl at 37°C with those measured for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl at 37°C. The predicted thermodynamics are listed in parentheses in Tables 1 and 2. On average at 37°C, the RNA/RNA duplexes in 1 M NaCl are 0.12 ± 0.01 kcal/mol of phosphate pairs more stable than the 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl. Presumably, much of this difference is due to a sequence independent effect of salt concentration, which would primarily affect the S° for duplex formation (22,69). Thus, a reasonable approximation for the first term on the right hand side of Equation 1 is:
(3)
Note that Gsym from the RNA/RNA calculation is subtracted because a 2'-O-methyl RNA/RNA duplex cannot be self-complementary because the backbones differ. For the duplexes studied here, the number of phosphate pairs is one less than the number of base pairs.
The effects of LNA substitutions are likely not very dependent on salt concentration. Thus, it is probable that in 1 M NaCl or in the presence of Mg2+ (70) that (2'-O-MeRNA/RNA) can be approximated by (RNA/RNA, 1 M NaCl).
Table 3 compares measured values for duplexes with more than one LNA to predictions from combining Equation 1–3. The measured values average –10.5 kcal/mol and the root-mean-square difference between measured and predicted values is 0.6 kcal/mol with the largest difference being 1.7 kcal/mol. Again, the sequence with the largest difference contains a 5'GMUL/3'CA nearest neighbor so the prediction would be improved if Equation 1 was corrected for the apparent instability of this motif.
The results for 2'-O-methyl RNA/RNA duplexes provide preliminary nearest-neighbor free energy increments for predicting stabilities of such duplexes
The comparison of predicted RNA/RNA stabilities with those measured for 2'-O-methyl RNA/RNA duplexes suggests that the INN-HB model will also be applicable to 2'-O-methyl RNA/RNA duplexes (71). The results in Tables 1 and 2 were combined with those for several other duplexes (E. Kierzek, R. Kierzek, and D.H. Turner, unpublished data) to give preliminary INN-HB parameters for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl (see Table 4). Three nearest neighbors are only represented once or twice in the database, and these parameters are in parentheses. The parameters for 2'-O-methyl RNA/RNA and RNA/RNA duplexes are similar, especially if the RNA/RNA Watson–Crick nearest-neighbor parameters are each made less favorable by 0.12 kcal/mol, which largely accounts for the difference in salt concentration as suggested above. Evidently, the first term on the right hand side of Equation 1 can also be approximated by:
(4)
Table 3 compares predictions from combining Equations 1 and 4 with measured values for duplexes with more than one LNA. The root-mean-square difference between measured and predicted values is 0.6 kcal/mol with the largest difference being the 1.7 kcal/mol associated with the duplex containing a 5'GMUL/3'CA nearest neighbor. Undoubtedly, this model can be expanded and refined by more measurements, but it appears sufficient to aid sequence design for many applications.
Complete LNA substitution is no more stabilizing than substitution at every other nucleotide starting at the second nucleotide from the 5' end
The effect of complete LNA substitution for a 2'-O-methyl RNA backbone was studied for the sequences 5'ALCLULALCLCLAL/3'r(UGAUGGU) and 5'GLCLULALCLULGL/3'r(CGAUGAC). As summarized in Table 2, the stabilities of these duplexes at 37°C are within experimental error of those measured for 5'AMCLUMALCMCLAM/3'r(UGAUGGU) and 5'GMCLUMALCMULGM/3'r(CGAUGAC), respectively. Evidently, the most effective use of LNA nucleotides is to space them every other nucleotide with the first LNA placed at the second nucleotide from the 5' end.
Internal mismatches make duplex formation less favorable
Table 5 contains thermodynamic parameters measured for the formation of duplexes containing single mismatches and the difference in stabilities relative to completely Watson–Crick complementary duplexes (Tables 1 and 2). All internal mismatches make duplex formation less favorable by at least 2 kcal/mol at 37°C corresponding to at least a 25-fold less favorable equilibrium constant for duplex formation. In general, terminal mismatches destabilize much less than internal mismatches. In fact, when the 3' terminal UL of 5'AMCMUMAMCMCMUL makes a GU pair, the duplex is stabilized by 0.14 kcal/mol at 37°C relative to a terminal AU pair.
Table 5 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without LNA substitutions: effects of mismatchesa
For four cases, the effect of a mismatch with an LNA nucleotide was compared with that for the equivalent 2'-O-methyl nucleotide. In each case, the mismatch penalty for the LNA was less than that for 2'-O-methyl RNA. However, for an AM-G mismatch flanked by LNAs in the context 5'ALCMULAMCLCMAL/3'r(UGAGGGU), the LNAs enhanced the mismatch penalty by 1 kcal/mol relative to a completely 2'-O-methyl RNA strand. Thus, oligonucleotides containing LNA may discriminate best against mismatches flanked by LNAs.
DISCUSSION
Oligonucleotide hybridization to RNA has many applications, ranging from quantifying gene expression (18–20,56) to designing therapeutics (4,8,21,46,52–54). LNA nucleotides have characteristics useful for these purposes. For example, LNA usually stabilizes duplexes (4,44,48,51) and is more resistant than RNA and DNA to nuclease digestion (4,6,51). The results presented here provide insights that are useful for designing 2'-O-methyl RNA/LNA chimeric oligonucleotides for various purposes. Some trends may be general for RNA A-form helixes and thus may also be relevant to other chimeras with nucleotides that favor A-form conformations.
The results suggest several principles for the design of 2'-O-methyl RNA/LNA chimeras for hybridization to RNA
The database in Tables 1 and 2 is too small to generate a complete model for the design of 2'-O-methyl RNA/LNA chimeras. Nevertheless, several trends are apparent: (i) LNA substitution at the 5' end has little effect on duplex stability. (ii) Except for a 3' terminal UL, interior and 3' terminal LNA substitutions have similar effects. (iii) Most of the stability enhancement is achieved when LNA nucleotides are separated by one 2'-O-methyl nucleotide. (iv) The effects of LNA substitutions are approximately additive when the LNA nucleotides are separated by at least one 2'-O-methyl nucleotide. Thus, in most cases, LNA nucleotides are used most effectively for duplex stabilization when separated by at least one 2'-O-methyl nucleotide and not placed on the 5' end of the chimera. The additivity observed suggests that a nearest-neighbor model (64,72–75) will be able to predict well the stabilities of duplexes formed between RNA and 2'-O-methyl RNA/LNA chimeras. While the database is not sufficient to determine all the parameters for such a model, the simpler model of Equation 1 provides reasonable approximations for stabilities at 37°C, though a correction should probably be applied for duplexes containing a 5'GMUL/3'CA nearest neighbor. It also appears that approximations for the stabilities of 2'-O-methyl RNA/RNA duplexes can be provided by nearest-neighbor models for 2'-O-methyl RNA/RNA duplexes (Equation 4) or even for RNA/RNA duplexes (Equations 2 and 3).
The magnitude and sequence dependence of the stabilization due to LNAs are surprising. Ribose and therefore probably 2'-O-methyl ribose sugars in single strands are typically found in roughly equal fractions in C2'-endo and C3'-endo conformations. If the methylene bridge of an LNA only locks the sugar into the C3'-endo conformation, then the expected stabilization due to preorganization would be: G° = –RT ln 2, which is –0.4 kcal/mol at 37°C (310.15 K). The stabilization observed for a 5' terminal LNA is roughly –0.4 kcal/mol, but the average stabilizations for internal LNAs and 3' terminal AL, CL and GL are more favorable at –1.3 and –1.2 kcal/mol, respectively. Moreover, if stabilization was only due to preorganization of an LNA sugar, then the effect would not saturate when alternate sugars are LNA. Evidently, the LNA substitution also affects the 5' neighboring base pair in a way that enhances the stabilization beyond that expected from preorganization of a single sugar. Interestingly, NMR structures of DNA/LNA chimeras bound to RNA show that only the DNA sugar 3' of the LNA is driven to a C3'-endo conformation for the sequence d(5'CTGATLATGC)/3'GACUAUACG, but all non-terminal DNA sugars are C3'-endo when all three Ts are LNAs (76).
Comparison with single LNA substitutions in DNA/DNA duplexes
McTigue et al. (48) measured the thermodynamic effects of single internal LNA substitutions in 100 DNA/DNA duplexes. The free energy increments at 37°C for LNA substitutions ranged from +0.83 to –1.90 kcal/mol with an average of –0.55 kcal/mol. This compares with a range from +0.18 to –2.17 kcal/mol and an average of –1.32 kcal/mol for the single internal LNA substitutions in Table 1. The comparision suggests that single LNA substitutions are on average more stabilizing to 2'-O-methyl RNA/RNA duplexes than to DNA/DNA duplexes. This may reflect the expectation that LNA substitutions do not have a large effect on the conformations of 2'-O-methyl RNA/RNA duplexes, but alter the conformations of DNA/DNA duplexes.
LNA substitutions should be useful for probing RNA with short 2'-O-methyl RNA oligonucleotides
RNA structure can be probed with short oligonucleotides on microarrays (3). To optimize such methods, it is necessary to have tight binding that is sequence independent and that discriminates against mismatches. It appears that LNA nucleotides can be used to achieve this. For example, free energy increments for 2'-O-methyl RNA/RNA nearest neighbors range from –0.7 to –3.5 kcal/mol, corresponding to 5'AMUM/3'UA and 5'GMCM/3'CG, respectively (Table 4). The average increment of –1.3 kcal/mol of internal and 3' terminal LNA can help compensate for such less favorable stability of AU relative to GC pairs. The stability enhancement from LNA can also allow the use of shorter oligonucleotides.
The potential disadvantage to LNA substitutions in 2'-O-methyl RNA oligonucleotides is that discrimination against mismatches containing an LNA may be less than with a complete 2'-O-methyl RNA backbone. This was clearly true for three of the four cases where such direct comparisons were made. Nevertheless, internal mismatches with LNA nucleotides are considerably destabilizing, averaging a penalty of 4.1 kcal/mol at 37°C (Table 5), which translates to almost a 1000-fold weaker binding due to a single mismatch. When LNAs flanked an AM-G mismatch, the mismatch penalty at 37°C was 4.4 kcal/mol compared with 3.3 kcal/mol in the absence of LNAs. Such an effect may reflect enhanced rigidity due to LNA, which thereby prevents a mismatch from adopting a favorable conformation. Thus, it may be advantageous to use LNAs to flank nucleotides likely to give small mismatch penalties.
SUPPLEMENTARY DATA
Supplementary Data is available at NAR Online.
ACKNOWLEDGEMENTS
This work was supported by Polish State Committee for Scientific Research (KBN) grant PZB-KBN 059/T09/14 to R.K., NIH grant 1R03 TW1068 to R.K. and D.H.T., and NIH grant GM 22939 to D.H.T. E.K. was a recipient of the fellowship for young scientist from Polish Science Foundation. A.C. is a recipient of a fellowship from the President of the Polish Academy of Sciences. Funding to pay the Open Access publication charges for this article was provided by NIH.
REFERENCES
Brown, T.S., Chadalavada, D.M., Bevilacqua, P.C. (2004) Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength J. Mol. Biol., 341, 695–712 .
Zamecnik, P.C. and Stephenson, M.L. (1978) Inhibition of Rous-sarcoma virus-replication and cell transformation by a specific oligodeoxynucleotide Proc. Natl Acad. Sci. USA, 75, 280–284 .
Ooms, M., Verhoef, K., Southern, E., Huthoff, H., Berkhout, B. (2004) Probing alternative foldings of the HIV-1 leader RNA by antisense oligonucleotide scanning arrays Nucleic Acids Res., 32, 819–827 .
Kurreck, J., Wyszko, E., Gillen, C., Erdmann, V.A. (2002) Design of antisense oligonucleotides stabilized by locked nucleic acids Nucleic Acids Res., 30, 1911–1918 .
Kushon, S.A., Jordan, J.P., Seifert, J.L., Nielsen, H., Nielsen, P.E., Armitage, B.A. (2001) Effect of secondary structure on the thermodynamics and kinetics of PNA hybridization to DNA hairpins J. Am. Chem. Soc., 123, 10805–10813 .
Elmen, J., Thonberg, H., Ljungberg, K., Frieden, M., Westergaard, M., Xu, Y.H., Wahren, B., Liang, Z.C., Urum, H., Koch, T., et al. (2005) Locked nucleic acid (LNA) mediated improvements in siRNA stability and functionality Nucleic Acids Res., 33, 439–447 .
Pancoska, P., Moravek, Z., Moll, U.M. (2004) Efficient RNA interference depends on global context of the target sequence: quantitative analysis of silencing efficiency using Eulerian graph representation of siRNA Nucleic Acids Res., 32, 1469–1479 .
Matveeva, O.V., Mathews, D.H., Tsodikov, A.D., Shabalina, S.A., Gesteland, R.F., Atkins, J.F., Freier, S.M. (2003) Thermodynamic criteria for high hit rate antisense oligonucleotide design Nucleic Acids Res., 31, 4989–4994 .
Schwarz, D.S., Hutvagner, G., Du, T., Xu, Z.S., Aronin, N., Zamore, P.D. (2003) Asymmetry in the assembly of the RNAi enzyme complex Cell, 115, 199–208 .
Bernard, P.S. and Wittwer, C.T. (2002) Real-time PCR technology for cancer diagnostics Clin. Chem., 48, 1178–1185 .
Bodrossy, L. and Sessitsch, A. (2004) Oligonucleotide microarrays in microbial diagnostics Curr. Opin. Microbiol., 7, 245–254 .
Heller, R.A., Schena, M., Chai, A., Shalon, D., Bedilion, T., Gilmore, J., Woolley, D.E., Davis, R.W. (1997) Discovery and analysis of inflammatory disease-related genes using cDNA microarrays Proc. Natl Acad. Sci. USA, 94, 2150–2155 .
Herr, A., Grutzmann, R., Matthaei, A., Artelt, J., Schrock, E., Rump, A., Pilarsky, C. (2005) High-resolution analysis of chromosomal imbalances using the Affymetrix 10K SNP genotyping chip Genomics, 85, 392–400 .
Lemieux, B., Aharoni, A., Schena, M. (1998) Overview of DNA chip technology Mol. Breed., 4, 277–289 .
Martell, M., Briones, C., de Vicente, A., Piron, M., Esteban, J.I., Esteban, R., Guardia, J., Gomez, J. (2004) Structural analysis of hepatitis C RNA genome using DNA microarrays Nucleic Acids Res., 32, e90 .
Long, W.H., Xiao, H.S., Gu, X.M., Zhang, Q.H., Yang, H.J., Zhao, G.P., Liu, J.H. (2004) A universal microarray for detection of SARS coronavirus J. Virol. Methods, 121, 57–63 .
Mei, R., Hubbell, E., Bekiranov, S., Mittmann, M., Christians, F.C., Shen, M.M., Lu, G., Fang, J., Liu, W.M., Ryder, T., et al. (2003) Probe selection for high-density oligonucleotide arrays Proc. Natl Acad. Sci. USA, 100, 11237–11242 .
Schena, M., Heller, R.A., Theriault, T.P., Konrad, K., Lachenmeier, E., Davis, R.W. (1998) Microarrays: biotechnology's discovery platform for functional genomics Trends Biotechnol., 16, 301–306 .
Schena, M., Shalon, D., Davis, R.W., Brown, P.O. (1995) Quantitative monitoring of gene-expression patterns with a complementary-DNA microarray Science, 270, 467–470 .
Schena, M., Shalon, D., Heller, R., Chai, A., Brown, P.O., Davis, R.W. (1996) Parallel human genome analysis: Microarray-based expression monitoring of 1000 genes Proc. Natl Acad. Sci. USA, 93, 10614–10619 .
Matveeva, O.V., Shabalina, S.A., Nemtsov, V.A., Tsodikov, A.D., Gesteland, R.F., Atkins, J.F. (2003) Thermodynamic calculations and statistical correlations for oligo-probes design Nucleic Acids Res., 31, 4211–4217 .
SantaLucia, J. and Hicks, D. (2004) The thermodynamics of DNA structural motifs Annu. Rev. Biophys. Biomol. Struct., 33, 415–440 .
Li, H.X. and Rothberg, L.J. (2004) Label-free colorimetric detection of specific sequences in genomic DNA amplified by the polymerase chain reaction J. Am. Chem. Soc., 126, 10958–10961 .
Chworos, A., Severcan, I., Koyfman, A.Y., Weinkam, P., Oroudjev, E., Hansma, H.G., Jaeger, L. (2004) Building programmable jigsaw puzzles with RNA Science, 306, 2068–2072 .
Pancoska, P., Moravek, Z., Moll, U.M. (2004) Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs Nucleic Acids Res., 32, 4630–4645 .
Wengel, J. (2004) Nucleic acid nanotechnology—towards Angstrom-scale engineering Org. Biomol. Chem., 2, 277–280 .
Seeman, N.C. (2005) From genes to machines: DNA nanomechanical devices Trends Biochem. Sci., 30, 119–125 .
Froehler, B.C., Wadwani, S., Terhorst, T.J., Gerrard, S.R. (1992) Oligodeoxynucleotides containing C-5 propyne analogs of 2'-deoxyuridine and 2'-deoxycytidine Tetrahedron Lett., 33, 5307–5310 .
Barnes, T.W. and Turner, D.H. (2001) C5-(1-Propynyl)-2'-deoxy-pyrimidines enhance mismatch penalties of DNA: RNA duplex formation Biochemistry, 40, 12738–12745 .
Barnes, T.W. and Turner, D.L.H. (2001) Long-range cooperativity due to C5-propynylation of oligopyrimidines enhances specific recognition by uridine of ribo-adenosine over ribo-guanosine J. Am. Chem. Soc., 123, 9186–9187 .
Nielsen, P.E. (2004) PNA technology Mol. Biotechnol., 26, 233–248 .
Pellestor, F. and Paulasova, P. (2004) The peptide nucleic acids (PNAs), powerful tools for molecular genetics and cytogenetics Eur. J. Hum. Genet., 12, 694–700 .
Pellestor, F. and Paulasova, P. (2004) The peptide nucleic acids, efficient tools for molecular diagnosis Int. J. Mol. Med., 13, 521–525 .
Gryaznov, S.M. (1999) Oligonucleotide N3'P5' phosphoramidates as potential therapeutic agents Biochim. Biophys. Acta, 1489, 131–140 .
Egli, M. and Gryaznov, S.M. (2000) Synthetic oligonucleotides as RNA mimetics: 2'-modified RNAs and N3'P5' phosphoramidates Cell. Mol. Life Sci., 57, 1440–1456 .
Gryaznov, S.M. and Winter, H. (1998) RNA mimetics: oligoribonucleotide N3'P5' phosphoramidates Nucleic Acids Res., 26, 4160–4167 .
Testa, S.M., Gryaznov, S.M., Turner, D.H. (1998) Antisense binding enhanced by tertiary interactions: Binding of phosphorothioate and N3'P5' phosphoramidate hexanucleotides to the catalytic core of a group I ribozyme from the mammalian pathogen Pneumocystis carinii Biochemistry, 37, 9379–9385 .
Testa, S.M., Gryaznov, S.M., Turner, D.H. (1999) In vitro suicide inhibition of self-splicing of a group I intron from Pneumocystis carinii by an N3'P5'phosphoramidate hexanucleotide Proc. Natl Acad. Sci. USA, 96, 2734–2739 .
Cramer, H. and Pfleiderer, W. (2000) Nucleotides LXIV: Synthesis, hybridization and enzymatic degradation studies of 2'-O-methyl-oligoribonucleotides and 2'-O-methyl/deoxy gapmers Nucleosides Nucleotides Nucleic Acids, 19, 1765–1777 .
Freier, S.M. and Altmann, K.H. (1997) The ups and downs of nucleic acid duplex stability: structure-stability studies on chemically-modified DNA:RNA duplexes Nucleic Acids Res., 25, 4429–4443 .
Lesnik, E.A. and Freier, S.M. (1998) What affects the effect of 2'-alkoxy modifications? 1. Stabilization effect of 2'-methoxy substitutions in uniformly modified DNA oligonucleotides Biochemistry, 37, 6991–6997 .
Majlessi, M., Nelson, N.C., Becker, M.M. (1998) Advantages of 2'-O-methyl oligoribonucleotide probes for detecting RNA targets Nucleic Acids Res., 26, 2224–2229 .
Sproat, B.S., Lamond, A.I., Beijer, B., Neuner, P., Ryder, U. (1989) Highly efficient chemical synthesis of 2'-O-methyloligoribonucleotides and tetrabiotinylated derivatives—novel probes that are resistant to degradation by RNA or DNA specific nucleases Nucleic Acids Res., 17, 3373–3386 .
Koshkin, A.A., Singh, S.K., Nielsen, P., Rajwanshi, V.K., Kumar, R., Meldgaard, M., Olsen, C.E., Wengel, J. (1998) LNA (locked nucleic acids): Synthesis of the adenine, cytosine, guanine, 5-methylcytosine, thymine and uracil bicyclonucleoside monomers, oligomerisation, and unprecedented nucleic acid recognition Tetrahedron, 54, 3607–3630 .
Jepsen, J.S., Sorensen, M.D., Wengel, J. (2004) Locked nucleic acid: a potent nucleic acid analog in therapeutics and biotechnology Oligonucleotides, 14, 130–146 .
Jepsen, J.S. and Wengel, J. (2004) LNA-antisense rivals siRNA for gene silencing Curr. Opin. Drug Discov. Dev., 7, 188–194 .
Koshkin, A.A., Fensholdt, J., Pfundheller, H.M., Lomholt, C. (2001) A simplified and efficient route to 2'-O, 4'-C-methylene-linked bicyclic ribonucleosides (locked nucleic acid) J. Org. Chem., 66, 8504–8512 .
McTigue, P.M., Peterson, R.J., Kahn, J.D. (2004) Sequence-dependent thermodynamic parameters; for locked nucleic acid (LNA)–DNA duplex formation Biochemistry, 43, 5388–5405 .
Singh, S.K., Nielsen, P., Koshkin, A.A., Wengel, J. (1998) LNA (locked nucleic acids): synthesis and high-affinity nucleic acid recognition Chem. Commun., 455–456 .
Frieden, M., Hansen, H.F., Koch, T. (2003) Nuclease stability of LNA oligonucleotides and LNA-DNA chimeras Nucleosides Nucleotides Nucleic Acids, 22, 1041–1043 .
Vester, B. and Wengel, J. (2004) LNA (locked nucleic acid): high-affinity targeting of complementary RNA and DNA Biochemistry, 43, 13233–13241 .
Morita, K., Hasegawa, C., Kaneko, M., Tsutsumi, S., Sone, J., Ishikawa, T., Imanishi, T., Koizumi, M. (2002) 2'-O,4'-C-ethylene-bridged nucleic acids (ENA): highly nuclease-resistant and thermodynamically stable oligonucleotides for antisense drug Bioorg. Med. Chem. Lett., 12, 73–76 .
Milligan, J.F., Matteucci, M.D., Martin, J.C. (1993) Current concepts in antisense drug design J. Med. Chem., 36, 1923–1937 .
Goodchild, J. (2004) Oligonucleotide therapeutics: 25 years agrowing Curr. Opin. Mol. Ther., 6, 120–128 .
Forrest, S., Cotton, R., Landegren, U., Southern, E. (1995) How to find all those mutations Nature Genet., 10, 375–376 .
Huber, M., Wei, T.F., Muller, U.R., Lefebvre, P.A., Marla, S.S., Bao, Y.P. (2004) Gold nanoparticle probe-based gene expression analysis with unamplified total human RNA Nucleic Acids Res., 32, e137 .
Chan, V.T.W., Fleming, K.A., McGee, J.O.D. (1988) Simultaneous extraction from clinical biopsies of high-molecular-weight DNA and RNA—comparative characterization by biotinylated and P-32-labeled probes on southern and northern blots Anal. Biochem., 168, 16–24 .
Oberbaumer, I. and Speth, C. (1990) Detection of RNA on northern blots by negative staining with aurintricarboxylic acid Anal. Biochem., 185, 77–79 .
Southern, E., Mir, K., Shchepinov, M. (1999) Molecular interactions on microarrays Nature Genet., 21, 5–9 .
Smith, R.M., Walton, C.M., Wu, C.H., Wu, G.Y. (2002) Secondary structure and hybridization accessibility of hepatitis C virus 3'-terminal sequences J. Virol., 76, 9563–9574 .
Uhlenbeck, O.C., Baller, J., Doty, P. (1970) Complementary oligonucleotide binding to the anticodon loop of fMet-transfer RNA Nature, 225, 508–510 .
McBride, L.J. and Caruthers, M.H. (1983) An investigation of several deoxynucleoside phosphoramidites useful for synthesizing deoxyoligonucleotides Tetrahedron Lett., 24, 245–248 .
Pedersen, D.S., Rosenbohm, C., Koch, T. (2002) Preparation of LNA phosphoramidites Synthesis, 802–808 .
Xia, T.B., SantaLucia, J., Burkard, M.E., Kierzek, R., Schroeder, S.J., Jiao, X.Q., Cox, C., Turner, D.H. (1998) Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson–Crick base pairs Biochemistry, 37, 14719–14735 .
Borer, P.N. (1975) Optical properties of nucleic acids, absorption and circular dichroism spectra In Fasman, G.D. (Ed.). Handbook of Biochemistry and Molecular Biology: Nucleic Acids, 3rd edn, Cleveland, OH CRC Press Vol. 1, pp. 589–595 .
Richards, E.G. (1975) Use of tables in calculation of absorption, optical rotatory dispersion and circular dichroism of polyribonucleotides In Fasman, G.D. (Ed.). Handbook of Biochemistry and Molecular Biology: Nucleic Acids, 3rd edn, Cleveland, OH CRC Press Vol. 1, pp. 596–603 .
McDowell, J.A. and Turner, D.H. (1996) Investigation of the structural basis for thermodynamic stabilities of tandem GU mismatches: Solution structure of (rGAGGUCUC)2 by two-dimensional NMR and simulated annealing Biochemistry, 35, 14077–14089 .
Adamiak, D.A., Rypniewski, W.R., Milecki, J., Adamiak, R.W. (2001) The 1.19 angstrom X-ray structure of 2'-O-Me(CGCGCG)2 duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove Nucleic Acids Res., 29, 4144–4153 .
Manning, G.S. (1978) Molecular theory of polyelectrolyte solutions with applications to electrostatic properties of polynucleotides Q. Rev. Biophys., 11, 179–246 .
Williams, A.P., Longfellow, C.E., Freier, S.M., Kierzek, R., Turner, D.H. (1989) Laser temperature-jump, spectroscopic, and thermodynamic study of salt effects on duplex formation by dGCATGC Biochemistry, 28, 4283–4291 .
SantaLucia, J. and Turner, D.H. (1997) Measuring the thermodynamics of RNA secondary structure formation Biopolymers, 44, 309–319 .
Tinoco, I., Borer, P.N., Dengler, B., Levine, M.D., Uhlenbeck, O.C., Crothers, D.M., Gralla, J. (1973) Improved estimation of secondary structure in ribonucleic-acids Nature New Biol., 246, 40–41 .
Borer, P.N., Dengler, B., Tinoco, I., Uhlenbeck, O.C. (1974) Stability of ribonucleic-acid double-stranded helices J. Mol. Biol., 86, 843–853 .
Gray, D.M. (1997) Derivation of nearest-neighbor properties from data on nucleic acid oligomers.1. Simple sets of independent sequences and the influence of absent nearest neighbors Biopolymers, 42, 783–793 .
Gray, D.M. (1997) Derivation of nearest-neighbor properties from data on nucleic acid oligomers.2. Thermodynamic parameters of DNA–RNA hybrids and DNA duplexes Biopolymers, 42, 795–810 .
Petersen, M., Bondensgaard, K., Wengel, J., Jacobsen, J.P. (2002) Locked nucleic acid (LNA) recognition of RNA: NMR solution structures of LNA: RNA hybrids J. Am. Chem. Soc., 124, 5974–5982 .(Elzbieta Kierzek1,2, Anna Ciesielska2, K)
*To whom correspondence should be addressed. Tel: +1 585 275 3207; Fax: +1 585 276 0205; Email: turner@chem.rochester.edu
ABSTRACT
The influence of locked nucleic acid (LNA) residues on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes is reported. Optical melting studies indicate that LNA incorporated into an otherwise 2'-O-methyl RNA oligonucleotide usually, but not always, enhances the stabilities of complementary duplexes formed with RNA. Several trends are apparent, including: (i) a 3' terminal U LNA and 5' terminal LNAs are less stabilizing than interior and other 3' terminal LNAs; (ii) most of the stability enhancement is achieved when LNA nucleotides are separated by at least one 2'-O-methyl nucleotide; and (iii) the effects of LNA substitutions are approximately additive when the LNA nucleotides are separated by at least one 2'-O-methyl nucleotide. An equation is proposed to approximate the stabilities of complementary duplexes formed with RNA when at least one 2'-O-methyl nucleotide separates LNA nucleotides. The sequence dependence of 2'-O-methyl RNA/RNA duplexes appears to be similar to that of RNA/RNA duplexes, and preliminary nearest-neighbor free energy increments at 37°C are presented for 2'-O-methyl RNA/RNA duplexes. Internal mismatches with LNA nucleotides significantly destabilize duplexes with RNA.
INTRODUCTION
Understanding the thermodynamics of nucleic acid duplexes is important for many reasons. For example, such knowledge facilitates design of ribozymes (1), antisense and RNAi oligonucleotides (2–9), diagnostic probes including those employed on microarrays (10–23) and structures useful for nanotechnology (24–27). Many modified residues have been developed for such applications. Examples include propynylated bases (28–30), peptide nucleic acids (5,31–33), N3'–P5' phosphoramidates (34–38) and 2'-O-alkyl RNA (39–43). A modification that is particularly stabilizing in DNA and RNA duplexes (44–51) is a methyl bridge between the 2' oxygen and 4' carbon of ribose to form a ‘locked nucleic acid’ or LNA as shown in Figure 1. McTigue et al. (48) have shown that the enhanced stability due to a single LNA residue in a DNA duplex can be predicted from a nearest-neighbor model.
Figure 1 Covalent structures of sugars used.
Hybridization of oligonucleotides to RNA is important for applications, such as antisense therapeutics (4,8,21,46,52–54), diagnostics (32,33,42,55), profiling gene expression with microarrays (18–20,56), identifying bands by Northern blots of gels (57,58) and probing RNA structure (1,3,15,59–61). Oligonucleotides with 2'-O-alkyl modifications can be particularly useful for these applications because they are easily synthesized (39,43), chemically stable and bind relatively tightly to RNA (39–42). However, for many applications, it is desirable to modulate the binding affinity. For example, sequence independent duplex stabilities would benefit applications that involve multiplex detection, such as microarrays. Here, we show that introduction of LNA into 2'-O-methyl RNA oligonucleotides can increase stabilities of 2'-O-methyl RNA/RNA hybrid duplexes and that the enhancements in stability can usually be predicted with a simple model.
MATERIALS AND METHODS
General methods
High-performance liquid chromatography (HPLC) was performed on a Hewlett Packard series 1100 HPLC with a reverse-phase Supelco RP-18 column (4.6 x 250 mm). Mass spectra were obtained on an LC MS Hewlett Packard series 1100 MSD with API-ES detector or on an AMD 604/402. Thin-layer chromatography (TLC) was carried out on Merck 60 F254 TLC plates with the mixture 1-propanol/aqueous ammonia/water = 55:35:10 (v/v/v).
Synthesis and purification of oligonucleotides
Oligoribonucleotides were synthesized on an Applied Biosystems DNA/RNA synthesizer, using ?-cyanoethyl phosphoramidite chemistry (62). For synthesis of standard RNA oligonucleotides, the commercially available phosphoramidites with 2'-O-tertbutyldimethylsilyl groups were used (Glen Research). For synthesis of 2'-O-methyl RNA oligonucleotides, the 3'-O-phosphoramidites of 2'-O-methylnucleotides were used (Glen Research and Proligo). The 3'-O-phosphoramidites of LNA nucleotides were synthesized according to the published procedures with some minor modifications (44,47,63). The details of deprotection and purification of oligoribonucleotides were described previously (64).
UV melting
Oligonucleotides were melted in buffer containing 100 mM NaCl, 20 mM sodium cacodylate, 0.5 mM Na2EDTA, pH 7.0. The relatively low NaCl concentration kept melting temperatures in the reasonable range even when there were multiple LNA substitutions. Oligonucleotide single-strand concentrations were calculated from absorbencies above 80°C and single-strand extinction coefficients were approximated by a nearest-neighbor model (65,66). It was assumed that 2'-O-methyl RNA and RNA strands with identical sequences have identical extinction coefficients. Absorbancy versus temperature melting curves were measured at 260 nm with a heating rate of 1°C/min from 0 to 90°C on a Beckman DU 640 spectrophotometer with a water cooled thermoprogrammer. Melting curves were analyzed and thermodynamic parameters were calculated from a two-state model with the program MeltWin 3.5 (67). For almost all sequences, the H° derived from versus ln (CT/4) plots is within 15% of that derived from averaging the fits to individual melting curves, as expected if the two-state model is reasonable.
Parameter fitting
Free energy parameters for predicting stabilities of 2'-O-methyl RNA/RNA and 2'-O-methyl RNA–LNA/RNA duplexes with the Individual Nearest-Neighbor Hydrogen Bonding (INN-HB) model (64) were obtained by multiple linear regression with the program Analyse-it v.1.71 (Analyse-It Software, Ltd, Leeds, England; www.analyse-it.com) which expands Microsoft Excel. Analyse-It was also used to obtain parameters for enhancement of stabilities of 2'-O-methyl RNA/RNA duplexes by substitution of LNA nucleotides internally and/or at the 3' end when the LNAs are separated by at least one 2'-O-methyl nucleotide. Results from versus ln (CT/4) plots were used as the data for the calculations.
RESULTS
Figures 2 and 3 show typical data from optical melting curves, and Table 1 lists the thermodynamic parameters for the helix to coil transition with either no or one LNA nucleotide in the primarily 2'-O-methyl strand of a hybrid with a Watson–Crick complementary RNA strand.
Figure 2 Representative UV melting curves 5'AMCMUMALCMCMAM/5'UGGCAGU (1, 52 μM), 5'AMCMUMAMCMCMAM/5'UGGUAGU (2, 46 μM), 5'AMCMUMALCMCMAM/5'UGGUAGU (3, 51 μM) and 5'AMCLUMALCMCLAM/5'UGGUAGU (4, 51 μM). Note that sequence 1 has an ALC mismatch in the middle.
Figure 3 Reciprocal melting temperature versus ln (CT/4) for 5'AMCMUMALCMCMAM/5'UGGCAGU (solid squares), 5'AMCMUMAMCMCMAM/5'UGGUAGU (solid circles), 5'AMCMUMALCMCMAM/5'UGGUAGU (solid stars) and 5'AMCLUMALCMCLAM/5'UGGUAGU (solid triangles) in 100 mM NaCl, 20 mM sodium cacodylate, 0.5 mM Na2EDTA, pH 7.0. The thermodynamics of duplex formation were obtained by fitting the data points to the equation, = (R/H°) ln (CT/4) + S°/H°. Here, is the inverse melting temperature in kelvin, R is the gas constant, 1.987 cal K–1 mol–1, CT is the total oligonucleotide strand concentration, and both strands have the same concentration.
Table 1 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without a single LNA substitutiona
Single LNA substitutions at the 5' end of heptamer duplexes have little effect on stability
The effects of single LNA substitutions at the 5' end of the 2'-O-methyl strand were studied in duplexes of the form, 5'NCMUMAMCMCMAM/3'r(QGAUGGU), where superscript M denotes a 2'-O-methyl sugar, N is A, C, G, or U with a 2'-O-methyl or LNA sugar, r denotes ribose sugars, and Q is the Watson–Crick complement to N. As summarized in Table 1, 5' terminal LNA substitutions make duplex stability more favorable by 0.3–0.6 kcal/mol at 37°C with an average enhancement of 0.45 kcal/mol. Thus, 5' terminal LNA substitutions increase the binding constant for duplex formation by 2-fold at 37°C.
The effects of single LNA substitutions at the 3' ends of heptamer duplexes is idiosyncratic
The effects of single LNA substitutions at the 3' end of the 2'-O-methyl strand was studied in duplexes of the form, 5'AMCMUMAMCMCMN/3'r(UGAUGGQ) and in 5'AMCMUMAMCMGMUM/L/3'r(UGAUGCA) (Table 1). If N is A, C or G, then LNA substitutions have similar effects. On average, an LNA substitution makes duplex stability more favorable by 1.2 kcal/mol at 37°C. In the two sequences with a 3' terminal LNA U on the 2'-O-methyl strand, duplex stability is, however, affected little, averaging a destabilization of 0.08 kcal/mol at 37°C. In both cases, the terminal U is preceded by a GC pair, but both orientations of the GC pair give similar destabilization upon LNA substitution at the 3' terminal U.
Single LNA substitutions in the interior of AMCMUMAMCMCMAM enhance the stability of the duplex formed with its complementary RNA by 1.4 kcal/mol
The effect of interior position on the free energy increment for a single LNA substitution for a 2'-O-methyl RNA was studied for the duplex 5'AMCMUMAMCMCMAM/3'r(UGAUGGU). As summarized in Table 1, a single interior LNA substitution makes duplex stability more favorable by 1.2–1.7 kcal/mol at 37°C, with an average of 1.4 kcal/mol. This corresponds to roughly a 10-fold increase in binding constant. Thus, interior and 3' terminal LNA substitutions usually improve binding more than 5' terminal LNA substitutions.
The effects of single LNA substitutions as a function of sequence context in the middle of heptamer duplexes
The effects from single LNA substitutions in the middle of heptamer duplexes were studied as a function of the LNA base and the 3' adjacent 2'-O-methyl nucleotide in duplexes of the form, 5'AMCMUMNVMCMAM/3'r(UGAQWGU), where VM is 2'-O-methyl A, C, G or U, and W is the ribose Watson–Crick complement to V. The effects on the for duplex formation of substituting LNA for 2'-O-methyl at position N are summarized in Table 1. For 13 of 16 sequences, the LNA substitution makes duplex stability more favorable by 1.0–1.5 kcal/mol at 37°C, with an average enhancement of 1.3 kcal/mol. The enhancement for the other three sequences averages 2.1 kcal/mol at 37°C.
The dependence on the 5' nearest-neighbor nucleotide of effects from substituting UL for UM was studied in duplexes of the form, 5'AMCMVMUL/MCMCMAM/3'r(UGWAGGU) (Table 1). For three of the four sequences, the range for enhancement of duplex stability is 0.8–1.3 kcal/mol and the average is 1.0 kcal/mol. The exception is for the nearest neighbor 5'GMUL/3'r(CA), where LNA substitution destabilizes the duplex by 0.2 kcal/mol.
The results for single LNA substitutions in other contexts provide additional insight into the dependence of free energy increments on the 5' nearest neighbor 2'-O-methyl nucleotide. In particular, the Watson–Crick complementary duplexes with 5'AMCLUMAMCMCMAM and 5'AMCMUMCLUMCMAM strands have a 5'CLUM/3'r(GA) nearest neighbor that is preceded by AM and UM, respectively. The LNA substitution for 2'-O-methyl in these cases results in duplex formation being more stable by 1.29 and 1.48 kcal/mol, respectively, at 37°C. Similarly, the duplexes with 5'AMGLUMAMCMCMAM and 5'AMCMUMGLUMCMAM each have a 5'GLUM/3'r(CA) nearest neighbor that is preceded by AM and UM, respectively. In both cases, the LNA substitution enhances duplex stability by 1.14 kcal/mol at 37°C. Thus, for seven duplexes, the enhanced stability from an LNA substitution is relatively independent of the nearest-neighbor nucleotide 5' to the LNA. The one exception is for the nearest neighbor 5'GMUL/3'r(CA). Interestingly, this nearest-neighbor combination is also destabilized by LNA substitution at a 3' terminal U (Table 1). Evidently, an LNA substitution in the middle of a 2'-O-methyl strand usually affects heteroduplex stability with an RNA strand by about the same amount as an LNA substitution at a 3' terminus.
The effects of LNA substitutions are approximately additive when LNA nucleotides are spaced by at least one 2'-O-methyl nucleotide
Table 2 contains thermodynamic parameters measured for duplexes having more than one LNA substitution and Table 3 compares the stabilities at 37°C with those predicted from four simple models. The first model, labeled ‘additivity’, predicts the for duplex formation in the 5'ACUACCA/3'UGAUGGU series by adding the free energy increments measured for single LNA substitutions in the same context to the for duplex formation in the absence of LNA nucleotides. The second model predicts the (kcal/mol) for duplex formation with the following equation as deduced from fitting the data in Tables 1 and 2 excluding sequences with adjacent LNAs.
(1)
Here, (2'-O-MeRNA/RNA) is the free energy change at 37°C for duplex formation in the absence of any LNA nucleotides, n5'tL is the number of 5' terminal LNAs, niAL/UL and niGL/CL are the number of internal LNAs in AU and GC pairs, respectively, n3'tU and n3'tAL/CL/GL are the number of 3' terminal LNAs that are U or not U, respectively. Both methods that use experimental data for (2'-O-MeRNA/RNA) provide reasonable predictions that are within 1 kcal/mol of the measured value (Table 3). Two other methods that use nearest-neighbor models to approximate (2'-O-MeRNA/RNA) provide somewhat less accurate, but still reasonable predictions as described below. The duplex with the worst prediction, 5'GMULUMCLGMGL/3'CAAGCC has a 5'GMUL/3'CA nearest neighbor, consistent with this motif being unusually unstable by 1.2 kcal/mol. Thus, it is likely that the of Equation 1 should be made less favorable by 1.2 kcal/mol for every internal 5'GMUL/3'CA nearest neighbor in a duplex. Evidently, the effects of multiple LNA substitutions are approximately additive when the LNAs are spaced by at least 1 nt.
Table 2 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without multiple LNA substitutionsa
Table 3 Comparison of measured and predicted stabilities
The data may also be fit to a nearest-neighbor model containing 30 of the LNA enhancement parameters associated with duplexes of RNA strands bound to 2'-O-methyl RNA/LNA chimeras. These parameters are listed in Supplementary Material. The number of occurrences for each nearest neighbor is limited, however, so the values are only roughly determined.
Predictions for RNA/RNA duplexes at 1 M NaCl can be used to approximate stabilities of 2'-O-methyl RNA/RNA duplexes at 0.1 M NaCl
The stabilities of RNA/RNA duplexes at 37°C and 1 M NaCl are predicted well by an Independent Nearest-Neighbor Hydrogen Bonding (INN-HB) model (64). In this model, the stability of an RNA/RNA duplex is approximated by:
(2)
Here, is the free energy change for initiating a helix; each is the free energy increment of the jth type nearest neighbor (see Table 4) with nj occurrences in the sequence; mterm-AU is the number of terminal AU pairs; is the free energy increment per terminal AU pair; Gsym is 0.43 kcal/mol at 37°C for self-complementary duplexes and 0 for non-self-complementary duplexes.
Table 4 Preliminary nearest-neighbor free energy parameters (kcal/mol at 37°C) for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl compared with parameters for RNA/RNA duplexes in 1 M NaCl
Similar sequence dependent parameters may also be applicable to 2'-O-methyl RNA/RNA heteroduplexes because they are expected to have A-form conformations similar to those of RNA/RNA homoduplexes (68). This was tested by comparing the predicted stabilities of RNA/RNA duplexes in 1 M NaCl at 37°C with those measured for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl at 37°C. The predicted thermodynamics are listed in parentheses in Tables 1 and 2. On average at 37°C, the RNA/RNA duplexes in 1 M NaCl are 0.12 ± 0.01 kcal/mol of phosphate pairs more stable than the 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl. Presumably, much of this difference is due to a sequence independent effect of salt concentration, which would primarily affect the S° for duplex formation (22,69). Thus, a reasonable approximation for the first term on the right hand side of Equation 1 is:
(3)
Note that Gsym from the RNA/RNA calculation is subtracted because a 2'-O-methyl RNA/RNA duplex cannot be self-complementary because the backbones differ. For the duplexes studied here, the number of phosphate pairs is one less than the number of base pairs.
The effects of LNA substitutions are likely not very dependent on salt concentration. Thus, it is probable that in 1 M NaCl or in the presence of Mg2+ (70) that (2'-O-MeRNA/RNA) can be approximated by (RNA/RNA, 1 M NaCl).
Table 3 compares measured values for duplexes with more than one LNA to predictions from combining Equation 1–3. The measured values average –10.5 kcal/mol and the root-mean-square difference between measured and predicted values is 0.6 kcal/mol with the largest difference being 1.7 kcal/mol. Again, the sequence with the largest difference contains a 5'GMUL/3'CA nearest neighbor so the prediction would be improved if Equation 1 was corrected for the apparent instability of this motif.
The results for 2'-O-methyl RNA/RNA duplexes provide preliminary nearest-neighbor free energy increments for predicting stabilities of such duplexes
The comparison of predicted RNA/RNA stabilities with those measured for 2'-O-methyl RNA/RNA duplexes suggests that the INN-HB model will also be applicable to 2'-O-methyl RNA/RNA duplexes (71). The results in Tables 1 and 2 were combined with those for several other duplexes (E. Kierzek, R. Kierzek, and D.H. Turner, unpublished data) to give preliminary INN-HB parameters for 2'-O-methyl RNA/RNA duplexes in 0.1 M NaCl (see Table 4). Three nearest neighbors are only represented once or twice in the database, and these parameters are in parentheses. The parameters for 2'-O-methyl RNA/RNA and RNA/RNA duplexes are similar, especially if the RNA/RNA Watson–Crick nearest-neighbor parameters are each made less favorable by 0.12 kcal/mol, which largely accounts for the difference in salt concentration as suggested above. Evidently, the first term on the right hand side of Equation 1 can also be approximated by:
(4)
Table 3 compares predictions from combining Equations 1 and 4 with measured values for duplexes with more than one LNA. The root-mean-square difference between measured and predicted values is 0.6 kcal/mol with the largest difference being the 1.7 kcal/mol associated with the duplex containing a 5'GMUL/3'CA nearest neighbor. Undoubtedly, this model can be expanded and refined by more measurements, but it appears sufficient to aid sequence design for many applications.
Complete LNA substitution is no more stabilizing than substitution at every other nucleotide starting at the second nucleotide from the 5' end
The effect of complete LNA substitution for a 2'-O-methyl RNA backbone was studied for the sequences 5'ALCLULALCLCLAL/3'r(UGAUGGU) and 5'GLCLULALCLULGL/3'r(CGAUGAC). As summarized in Table 2, the stabilities of these duplexes at 37°C are within experimental error of those measured for 5'AMCLUMALCMCLAM/3'r(UGAUGGU) and 5'GMCLUMALCMULGM/3'r(CGAUGAC), respectively. Evidently, the most effective use of LNA nucleotides is to space them every other nucleotide with the first LNA placed at the second nucleotide from the 5' end.
Internal mismatches make duplex formation less favorable
Table 5 contains thermodynamic parameters measured for the formation of duplexes containing single mismatches and the difference in stabilities relative to completely Watson–Crick complementary duplexes (Tables 1 and 2). All internal mismatches make duplex formation less favorable by at least 2 kcal/mol at 37°C corresponding to at least a 25-fold less favorable equilibrium constant for duplex formation. In general, terminal mismatches destabilize much less than internal mismatches. In fact, when the 3' terminal UL of 5'AMCMUMAMCMCMUL makes a GU pair, the duplex is stabilized by 0.14 kcal/mol at 37°C relative to a terminal AU pair.
Table 5 Thermodynamic parameters of duplex formation between RNA and 2'-O-methyl oligoribonucleotides with and without LNA substitutions: effects of mismatchesa
For four cases, the effect of a mismatch with an LNA nucleotide was compared with that for the equivalent 2'-O-methyl nucleotide. In each case, the mismatch penalty for the LNA was less than that for 2'-O-methyl RNA. However, for an AM-G mismatch flanked by LNAs in the context 5'ALCMULAMCLCMAL/3'r(UGAGGGU), the LNAs enhanced the mismatch penalty by 1 kcal/mol relative to a completely 2'-O-methyl RNA strand. Thus, oligonucleotides containing LNA may discriminate best against mismatches flanked by LNAs.
DISCUSSION
Oligonucleotide hybridization to RNA has many applications, ranging from quantifying gene expression (18–20,56) to designing therapeutics (4,8,21,46,52–54). LNA nucleotides have characteristics useful for these purposes. For example, LNA usually stabilizes duplexes (4,44,48,51) and is more resistant than RNA and DNA to nuclease digestion (4,6,51). The results presented here provide insights that are useful for designing 2'-O-methyl RNA/LNA chimeric oligonucleotides for various purposes. Some trends may be general for RNA A-form helixes and thus may also be relevant to other chimeras with nucleotides that favor A-form conformations.
The results suggest several principles for the design of 2'-O-methyl RNA/LNA chimeras for hybridization to RNA
The database in Tables 1 and 2 is too small to generate a complete model for the design of 2'-O-methyl RNA/LNA chimeras. Nevertheless, several trends are apparent: (i) LNA substitution at the 5' end has little effect on duplex stability. (ii) Except for a 3' terminal UL, interior and 3' terminal LNA substitutions have similar effects. (iii) Most of the stability enhancement is achieved when LNA nucleotides are separated by one 2'-O-methyl nucleotide. (iv) The effects of LNA substitutions are approximately additive when the LNA nucleotides are separated by at least one 2'-O-methyl nucleotide. Thus, in most cases, LNA nucleotides are used most effectively for duplex stabilization when separated by at least one 2'-O-methyl nucleotide and not placed on the 5' end of the chimera. The additivity observed suggests that a nearest-neighbor model (64,72–75) will be able to predict well the stabilities of duplexes formed between RNA and 2'-O-methyl RNA/LNA chimeras. While the database is not sufficient to determine all the parameters for such a model, the simpler model of Equation 1 provides reasonable approximations for stabilities at 37°C, though a correction should probably be applied for duplexes containing a 5'GMUL/3'CA nearest neighbor. It also appears that approximations for the stabilities of 2'-O-methyl RNA/RNA duplexes can be provided by nearest-neighbor models for 2'-O-methyl RNA/RNA duplexes (Equation 4) or even for RNA/RNA duplexes (Equations 2 and 3).
The magnitude and sequence dependence of the stabilization due to LNAs are surprising. Ribose and therefore probably 2'-O-methyl ribose sugars in single strands are typically found in roughly equal fractions in C2'-endo and C3'-endo conformations. If the methylene bridge of an LNA only locks the sugar into the C3'-endo conformation, then the expected stabilization due to preorganization would be: G° = –RT ln 2, which is –0.4 kcal/mol at 37°C (310.15 K). The stabilization observed for a 5' terminal LNA is roughly –0.4 kcal/mol, but the average stabilizations for internal LNAs and 3' terminal AL, CL and GL are more favorable at –1.3 and –1.2 kcal/mol, respectively. Moreover, if stabilization was only due to preorganization of an LNA sugar, then the effect would not saturate when alternate sugars are LNA. Evidently, the LNA substitution also affects the 5' neighboring base pair in a way that enhances the stabilization beyond that expected from preorganization of a single sugar. Interestingly, NMR structures of DNA/LNA chimeras bound to RNA show that only the DNA sugar 3' of the LNA is driven to a C3'-endo conformation for the sequence d(5'CTGATLATGC)/3'GACUAUACG, but all non-terminal DNA sugars are C3'-endo when all three Ts are LNAs (76).
Comparison with single LNA substitutions in DNA/DNA duplexes
McTigue et al. (48) measured the thermodynamic effects of single internal LNA substitutions in 100 DNA/DNA duplexes. The free energy increments at 37°C for LNA substitutions ranged from +0.83 to –1.90 kcal/mol with an average of –0.55 kcal/mol. This compares with a range from +0.18 to –2.17 kcal/mol and an average of –1.32 kcal/mol for the single internal LNA substitutions in Table 1. The comparision suggests that single LNA substitutions are on average more stabilizing to 2'-O-methyl RNA/RNA duplexes than to DNA/DNA duplexes. This may reflect the expectation that LNA substitutions do not have a large effect on the conformations of 2'-O-methyl RNA/RNA duplexes, but alter the conformations of DNA/DNA duplexes.
LNA substitutions should be useful for probing RNA with short 2'-O-methyl RNA oligonucleotides
RNA structure can be probed with short oligonucleotides on microarrays (3). To optimize such methods, it is necessary to have tight binding that is sequence independent and that discriminates against mismatches. It appears that LNA nucleotides can be used to achieve this. For example, free energy increments for 2'-O-methyl RNA/RNA nearest neighbors range from –0.7 to –3.5 kcal/mol, corresponding to 5'AMUM/3'UA and 5'GMCM/3'CG, respectively (Table 4). The average increment of –1.3 kcal/mol of internal and 3' terminal LNA can help compensate for such less favorable stability of AU relative to GC pairs. The stability enhancement from LNA can also allow the use of shorter oligonucleotides.
The potential disadvantage to LNA substitutions in 2'-O-methyl RNA oligonucleotides is that discrimination against mismatches containing an LNA may be less than with a complete 2'-O-methyl RNA backbone. This was clearly true for three of the four cases where such direct comparisons were made. Nevertheless, internal mismatches with LNA nucleotides are considerably destabilizing, averaging a penalty of 4.1 kcal/mol at 37°C (Table 5), which translates to almost a 1000-fold weaker binding due to a single mismatch. When LNAs flanked an AM-G mismatch, the mismatch penalty at 37°C was 4.4 kcal/mol compared with 3.3 kcal/mol in the absence of LNAs. Such an effect may reflect enhanced rigidity due to LNA, which thereby prevents a mismatch from adopting a favorable conformation. Thus, it may be advantageous to use LNAs to flank nucleotides likely to give small mismatch penalties.
SUPPLEMENTARY DATA
Supplementary Data is available at NAR Online.
ACKNOWLEDGEMENTS
This work was supported by Polish State Committee for Scientific Research (KBN) grant PZB-KBN 059/T09/14 to R.K., NIH grant 1R03 TW1068 to R.K. and D.H.T., and NIH grant GM 22939 to D.H.T. E.K. was a recipient of the fellowship for young scientist from Polish Science Foundation. A.C. is a recipient of a fellowship from the President of the Polish Academy of Sciences. Funding to pay the Open Access publication charges for this article was provided by NIH.
REFERENCES
Brown, T.S., Chadalavada, D.M., Bevilacqua, P.C. (2004) Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength J. Mol. Biol., 341, 695–712 .
Zamecnik, P.C. and Stephenson, M.L. (1978) Inhibition of Rous-sarcoma virus-replication and cell transformation by a specific oligodeoxynucleotide Proc. Natl Acad. Sci. USA, 75, 280–284 .
Ooms, M., Verhoef, K., Southern, E., Huthoff, H., Berkhout, B. (2004) Probing alternative foldings of the HIV-1 leader RNA by antisense oligonucleotide scanning arrays Nucleic Acids Res., 32, 819–827 .
Kurreck, J., Wyszko, E., Gillen, C., Erdmann, V.A. (2002) Design of antisense oligonucleotides stabilized by locked nucleic acids Nucleic Acids Res., 30, 1911–1918 .
Kushon, S.A., Jordan, J.P., Seifert, J.L., Nielsen, H., Nielsen, P.E., Armitage, B.A. (2001) Effect of secondary structure on the thermodynamics and kinetics of PNA hybridization to DNA hairpins J. Am. Chem. Soc., 123, 10805–10813 .
Elmen, J., Thonberg, H., Ljungberg, K., Frieden, M., Westergaard, M., Xu, Y.H., Wahren, B., Liang, Z.C., Urum, H., Koch, T., et al. (2005) Locked nucleic acid (LNA) mediated improvements in siRNA stability and functionality Nucleic Acids Res., 33, 439–447 .
Pancoska, P., Moravek, Z., Moll, U.M. (2004) Efficient RNA interference depends on global context of the target sequence: quantitative analysis of silencing efficiency using Eulerian graph representation of siRNA Nucleic Acids Res., 32, 1469–1479 .
Matveeva, O.V., Mathews, D.H., Tsodikov, A.D., Shabalina, S.A., Gesteland, R.F., Atkins, J.F., Freier, S.M. (2003) Thermodynamic criteria for high hit rate antisense oligonucleotide design Nucleic Acids Res., 31, 4989–4994 .
Schwarz, D.S., Hutvagner, G., Du, T., Xu, Z.S., Aronin, N., Zamore, P.D. (2003) Asymmetry in the assembly of the RNAi enzyme complex Cell, 115, 199–208 .
Bernard, P.S. and Wittwer, C.T. (2002) Real-time PCR technology for cancer diagnostics Clin. Chem., 48, 1178–1185 .
Bodrossy, L. and Sessitsch, A. (2004) Oligonucleotide microarrays in microbial diagnostics Curr. Opin. Microbiol., 7, 245–254 .
Heller, R.A., Schena, M., Chai, A., Shalon, D., Bedilion, T., Gilmore, J., Woolley, D.E., Davis, R.W. (1997) Discovery and analysis of inflammatory disease-related genes using cDNA microarrays Proc. Natl Acad. Sci. USA, 94, 2150–2155 .
Herr, A., Grutzmann, R., Matthaei, A., Artelt, J., Schrock, E., Rump, A., Pilarsky, C. (2005) High-resolution analysis of chromosomal imbalances using the Affymetrix 10K SNP genotyping chip Genomics, 85, 392–400 .
Lemieux, B., Aharoni, A., Schena, M. (1998) Overview of DNA chip technology Mol. Breed., 4, 277–289 .
Martell, M., Briones, C., de Vicente, A., Piron, M., Esteban, J.I., Esteban, R., Guardia, J., Gomez, J. (2004) Structural analysis of hepatitis C RNA genome using DNA microarrays Nucleic Acids Res., 32, e90 .
Long, W.H., Xiao, H.S., Gu, X.M., Zhang, Q.H., Yang, H.J., Zhao, G.P., Liu, J.H. (2004) A universal microarray for detection of SARS coronavirus J. Virol. Methods, 121, 57–63 .
Mei, R., Hubbell, E., Bekiranov, S., Mittmann, M., Christians, F.C., Shen, M.M., Lu, G., Fang, J., Liu, W.M., Ryder, T., et al. (2003) Probe selection for high-density oligonucleotide arrays Proc. Natl Acad. Sci. USA, 100, 11237–11242 .
Schena, M., Heller, R.A., Theriault, T.P., Konrad, K., Lachenmeier, E., Davis, R.W. (1998) Microarrays: biotechnology's discovery platform for functional genomics Trends Biotechnol., 16, 301–306 .
Schena, M., Shalon, D., Davis, R.W., Brown, P.O. (1995) Quantitative monitoring of gene-expression patterns with a complementary-DNA microarray Science, 270, 467–470 .
Schena, M., Shalon, D., Heller, R., Chai, A., Brown, P.O., Davis, R.W. (1996) Parallel human genome analysis: Microarray-based expression monitoring of 1000 genes Proc. Natl Acad. Sci. USA, 93, 10614–10619 .
Matveeva, O.V., Shabalina, S.A., Nemtsov, V.A., Tsodikov, A.D., Gesteland, R.F., Atkins, J.F. (2003) Thermodynamic calculations and statistical correlations for oligo-probes design Nucleic Acids Res., 31, 4211–4217 .
SantaLucia, J. and Hicks, D. (2004) The thermodynamics of DNA structural motifs Annu. Rev. Biophys. Biomol. Struct., 33, 415–440 .
Li, H.X. and Rothberg, L.J. (2004) Label-free colorimetric detection of specific sequences in genomic DNA amplified by the polymerase chain reaction J. Am. Chem. Soc., 126, 10958–10961 .
Chworos, A., Severcan, I., Koyfman, A.Y., Weinkam, P., Oroudjev, E., Hansma, H.G., Jaeger, L. (2004) Building programmable jigsaw puzzles with RNA Science, 306, 2068–2072 .
Pancoska, P., Moravek, Z., Moll, U.M. (2004) Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs Nucleic Acids Res., 32, 4630–4645 .
Wengel, J. (2004) Nucleic acid nanotechnology—towards Angstrom-scale engineering Org. Biomol. Chem., 2, 277–280 .
Seeman, N.C. (2005) From genes to machines: DNA nanomechanical devices Trends Biochem. Sci., 30, 119–125 .
Froehler, B.C., Wadwani, S., Terhorst, T.J., Gerrard, S.R. (1992) Oligodeoxynucleotides containing C-5 propyne analogs of 2'-deoxyuridine and 2'-deoxycytidine Tetrahedron Lett., 33, 5307–5310 .
Barnes, T.W. and Turner, D.H. (2001) C5-(1-Propynyl)-2'-deoxy-pyrimidines enhance mismatch penalties of DNA: RNA duplex formation Biochemistry, 40, 12738–12745 .
Barnes, T.W. and Turner, D.L.H. (2001) Long-range cooperativity due to C5-propynylation of oligopyrimidines enhances specific recognition by uridine of ribo-adenosine over ribo-guanosine J. Am. Chem. Soc., 123, 9186–9187 .
Nielsen, P.E. (2004) PNA technology Mol. Biotechnol., 26, 233–248 .
Pellestor, F. and Paulasova, P. (2004) The peptide nucleic acids (PNAs), powerful tools for molecular genetics and cytogenetics Eur. J. Hum. Genet., 12, 694–700 .
Pellestor, F. and Paulasova, P. (2004) The peptide nucleic acids, efficient tools for molecular diagnosis Int. J. Mol. Med., 13, 521–525 .
Gryaznov, S.M. (1999) Oligonucleotide N3'P5' phosphoramidates as potential therapeutic agents Biochim. Biophys. Acta, 1489, 131–140 .
Egli, M. and Gryaznov, S.M. (2000) Synthetic oligonucleotides as RNA mimetics: 2'-modified RNAs and N3'P5' phosphoramidates Cell. Mol. Life Sci., 57, 1440–1456 .
Gryaznov, S.M. and Winter, H. (1998) RNA mimetics: oligoribonucleotide N3'P5' phosphoramidates Nucleic Acids Res., 26, 4160–4167 .
Testa, S.M., Gryaznov, S.M., Turner, D.H. (1998) Antisense binding enhanced by tertiary interactions: Binding of phosphorothioate and N3'P5' phosphoramidate hexanucleotides to the catalytic core of a group I ribozyme from the mammalian pathogen Pneumocystis carinii Biochemistry, 37, 9379–9385 .
Testa, S.M., Gryaznov, S.M., Turner, D.H. (1999) In vitro suicide inhibition of self-splicing of a group I intron from Pneumocystis carinii by an N3'P5'phosphoramidate hexanucleotide Proc. Natl Acad. Sci. USA, 96, 2734–2739 .
Cramer, H. and Pfleiderer, W. (2000) Nucleotides LXIV: Synthesis, hybridization and enzymatic degradation studies of 2'-O-methyl-oligoribonucleotides and 2'-O-methyl/deoxy gapmers Nucleosides Nucleotides Nucleic Acids, 19, 1765–1777 .
Freier, S.M. and Altmann, K.H. (1997) The ups and downs of nucleic acid duplex stability: structure-stability studies on chemically-modified DNA:RNA duplexes Nucleic Acids Res., 25, 4429–4443 .
Lesnik, E.A. and Freier, S.M. (1998) What affects the effect of 2'-alkoxy modifications? 1. Stabilization effect of 2'-methoxy substitutions in uniformly modified DNA oligonucleotides Biochemistry, 37, 6991–6997 .
Majlessi, M., Nelson, N.C., Becker, M.M. (1998) Advantages of 2'-O-methyl oligoribonucleotide probes for detecting RNA targets Nucleic Acids Res., 26, 2224–2229 .
Sproat, B.S., Lamond, A.I., Beijer, B., Neuner, P., Ryder, U. (1989) Highly efficient chemical synthesis of 2'-O-methyloligoribonucleotides and tetrabiotinylated derivatives—novel probes that are resistant to degradation by RNA or DNA specific nucleases Nucleic Acids Res., 17, 3373–3386 .
Koshkin, A.A., Singh, S.K., Nielsen, P., Rajwanshi, V.K., Kumar, R., Meldgaard, M., Olsen, C.E., Wengel, J. (1998) LNA (locked nucleic acids): Synthesis of the adenine, cytosine, guanine, 5-methylcytosine, thymine and uracil bicyclonucleoside monomers, oligomerisation, and unprecedented nucleic acid recognition Tetrahedron, 54, 3607–3630 .
Jepsen, J.S., Sorensen, M.D., Wengel, J. (2004) Locked nucleic acid: a potent nucleic acid analog in therapeutics and biotechnology Oligonucleotides, 14, 130–146 .
Jepsen, J.S. and Wengel, J. (2004) LNA-antisense rivals siRNA for gene silencing Curr. Opin. Drug Discov. Dev., 7, 188–194 .
Koshkin, A.A., Fensholdt, J., Pfundheller, H.M., Lomholt, C. (2001) A simplified and efficient route to 2'-O, 4'-C-methylene-linked bicyclic ribonucleosides (locked nucleic acid) J. Org. Chem., 66, 8504–8512 .
McTigue, P.M., Peterson, R.J., Kahn, J.D. (2004) Sequence-dependent thermodynamic parameters; for locked nucleic acid (LNA)–DNA duplex formation Biochemistry, 43, 5388–5405 .
Singh, S.K., Nielsen, P., Koshkin, A.A., Wengel, J. (1998) LNA (locked nucleic acids): synthesis and high-affinity nucleic acid recognition Chem. Commun., 455–456 .
Frieden, M., Hansen, H.F., Koch, T. (2003) Nuclease stability of LNA oligonucleotides and LNA-DNA chimeras Nucleosides Nucleotides Nucleic Acids, 22, 1041–1043 .
Vester, B. and Wengel, J. (2004) LNA (locked nucleic acid): high-affinity targeting of complementary RNA and DNA Biochemistry, 43, 13233–13241 .
Morita, K., Hasegawa, C., Kaneko, M., Tsutsumi, S., Sone, J., Ishikawa, T., Imanishi, T., Koizumi, M. (2002) 2'-O,4'-C-ethylene-bridged nucleic acids (ENA): highly nuclease-resistant and thermodynamically stable oligonucleotides for antisense drug Bioorg. Med. Chem. Lett., 12, 73–76 .
Milligan, J.F., Matteucci, M.D., Martin, J.C. (1993) Current concepts in antisense drug design J. Med. Chem., 36, 1923–1937 .
Goodchild, J. (2004) Oligonucleotide therapeutics: 25 years agrowing Curr. Opin. Mol. Ther., 6, 120–128 .
Forrest, S., Cotton, R., Landegren, U., Southern, E. (1995) How to find all those mutations Nature Genet., 10, 375–376 .
Huber, M., Wei, T.F., Muller, U.R., Lefebvre, P.A., Marla, S.S., Bao, Y.P. (2004) Gold nanoparticle probe-based gene expression analysis with unamplified total human RNA Nucleic Acids Res., 32, e137 .
Chan, V.T.W., Fleming, K.A., McGee, J.O.D. (1988) Simultaneous extraction from clinical biopsies of high-molecular-weight DNA and RNA—comparative characterization by biotinylated and P-32-labeled probes on southern and northern blots Anal. Biochem., 168, 16–24 .
Oberbaumer, I. and Speth, C. (1990) Detection of RNA on northern blots by negative staining with aurintricarboxylic acid Anal. Biochem., 185, 77–79 .
Southern, E., Mir, K., Shchepinov, M. (1999) Molecular interactions on microarrays Nature Genet., 21, 5–9 .
Smith, R.M., Walton, C.M., Wu, C.H., Wu, G.Y. (2002) Secondary structure and hybridization accessibility of hepatitis C virus 3'-terminal sequences J. Virol., 76, 9563–9574 .
Uhlenbeck, O.C., Baller, J., Doty, P. (1970) Complementary oligonucleotide binding to the anticodon loop of fMet-transfer RNA Nature, 225, 508–510 .
McBride, L.J. and Caruthers, M.H. (1983) An investigation of several deoxynucleoside phosphoramidites useful for synthesizing deoxyoligonucleotides Tetrahedron Lett., 24, 245–248 .
Pedersen, D.S., Rosenbohm, C., Koch, T. (2002) Preparation of LNA phosphoramidites Synthesis, 802–808 .
Xia, T.B., SantaLucia, J., Burkard, M.E., Kierzek, R., Schroeder, S.J., Jiao, X.Q., Cox, C., Turner, D.H. (1998) Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson–Crick base pairs Biochemistry, 37, 14719–14735 .
Borer, P.N. (1975) Optical properties of nucleic acids, absorption and circular dichroism spectra In Fasman, G.D. (Ed.). Handbook of Biochemistry and Molecular Biology: Nucleic Acids, 3rd edn, Cleveland, OH CRC Press Vol. 1, pp. 589–595 .
Richards, E.G. (1975) Use of tables in calculation of absorption, optical rotatory dispersion and circular dichroism of polyribonucleotides In Fasman, G.D. (Ed.). Handbook of Biochemistry and Molecular Biology: Nucleic Acids, 3rd edn, Cleveland, OH CRC Press Vol. 1, pp. 596–603 .
McDowell, J.A. and Turner, D.H. (1996) Investigation of the structural basis for thermodynamic stabilities of tandem GU mismatches: Solution structure of (rGAGGUCUC)2 by two-dimensional NMR and simulated annealing Biochemistry, 35, 14077–14089 .
Adamiak, D.A., Rypniewski, W.R., Milecki, J., Adamiak, R.W. (2001) The 1.19 angstrom X-ray structure of 2'-O-Me(CGCGCG)2 duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove Nucleic Acids Res., 29, 4144–4153 .
Manning, G.S. (1978) Molecular theory of polyelectrolyte solutions with applications to electrostatic properties of polynucleotides Q. Rev. Biophys., 11, 179–246 .
Williams, A.P., Longfellow, C.E., Freier, S.M., Kierzek, R., Turner, D.H. (1989) Laser temperature-jump, spectroscopic, and thermodynamic study of salt effects on duplex formation by dGCATGC Biochemistry, 28, 4283–4291 .
SantaLucia, J. and Turner, D.H. (1997) Measuring the thermodynamics of RNA secondary structure formation Biopolymers, 44, 309–319 .
Tinoco, I., Borer, P.N., Dengler, B., Levine, M.D., Uhlenbeck, O.C., Crothers, D.M., Gralla, J. (1973) Improved estimation of secondary structure in ribonucleic-acids Nature New Biol., 246, 40–41 .
Borer, P.N., Dengler, B., Tinoco, I., Uhlenbeck, O.C. (1974) Stability of ribonucleic-acid double-stranded helices J. Mol. Biol., 86, 843–853 .
Gray, D.M. (1997) Derivation of nearest-neighbor properties from data on nucleic acid oligomers.1. Simple sets of independent sequences and the influence of absent nearest neighbors Biopolymers, 42, 783–793 .
Gray, D.M. (1997) Derivation of nearest-neighbor properties from data on nucleic acid oligomers.2. Thermodynamic parameters of DNA–RNA hybrids and DNA duplexes Biopolymers, 42, 795–810 .
Petersen, M., Bondensgaard, K., Wengel, J., Jacobsen, J.P. (2002) Locked nucleic acid (LNA) recognition of RNA: NMR solution structures of LNA: RNA hybrids J. Am. Chem. Soc., 124, 5974–5982 .(Elzbieta Kierzek1,2, Anna Ciesielska2, K)