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.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

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.     

Thermal stability of RNA hairpins containing a four-membered loop and a bulge nucleotide

Fourteen RNA hairpins containing a four-membered loop and a bulge nucleotide were synthesized and their thermal stabilities determined. The combined contribution of a four-membered loop and bulge A to the free energy of a hairpin is calculated to be 9.3 kcal/mol at 37 degrees C and successfully predicts the stability of an independent RNA hairpin. The introduction of a bulge nucleotide to the helical stem of an RNA hairpin destabilizes the molecule in a sequence-dependent manner. The individual thermodynamic contributions of a four-membered loop and bulge A, G, and U residues to the stability of an RNA hairpin loop are presented.     

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.     

Characterization of RNA hairpin loop stability.

Fifteen RNA hairpins that share the same stem sequence and have homopolymer loops of A, C and U residues which vary in length from three to nine nucleotides were synthesized and their thermal stabilities determined. Tm varies as a function of loop size but is almost independent of loop composition. Loops of four or five nucleotides are found to be the most stable loop size. This is consistent with the observation that four-membered loops are the most prevalent loop size in 16S-like RNAs. The contribution of each loop to hairpin stability was calculated by subtracting the known contribution of the helical stem. These data should be useful for predicting the stability of other hairpins.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Thermodynamics of single mismatches in RNA duplexes

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

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

Influence of nearest neighbor sequence on the stability of base pair mismatches in long DNA; determination by temperature-gradient gel electrophoresis.

Temperature-gradient gel electrophoresis (TGGE) was employed to determine the thermal stabilities of 48 DNA fragments that differ by single base pair mismatches. The approach provides a rapid way for studying how specific base mismatches effect the stability of a long DNA fragment. Homologous 373 bp DNA fragments differing by single base pair substitutions in their first melting domain were employed. Heteroduplexes were formed by melting and reannealing pairs of DNAs, one of which was 32P-labeled on its 5'-end. Product DNAs were separated based on their thermal stability by parallel and perpendicular temperature-gradient gel electrophoresis. The order of stability was determined for all common base pairs and mismatched bases in four different nearest neighbor environments; d(GXT).d(AYC), d(GXG).d(CYC), d(CXA).d(TYG), and d(TXT).d(AYA) with X,Y = A, T, C, or G. DNA fragments containing a single mismatch were destabilized by 1 to 5 degrees C with respect to homologous DNAs with complete Watson-Crick base pairing. Both the bases at the mismatch site and neighboring stacking interactions influence the destabilization caused by a mismatch. G.T, G.G and G.A mismatches were always among the most stable mismatches for all nearest neighbor environments examined. Purine.purine mismatches were generally more stable than pyrimidine.pyrimidine mispairs. Our results are in very good agreement with data where available from solution studies of short DNA oligomers.     

Thermodynamic stability of RNA structures formed by CNG trinucleotide repeats. Implication for prediction of RNA structure

Trinucleotide repeat expansion diseases (TREDs) are correlated with elongation of CNG DNA and RNA repeats to pathological level. This paper shows, for the first time, complete data concerning thermodynamic stabilities of RNA with CNG trinucleotide repeats. Our studies include the stability of oligoribonucleotides composed of two to seven of CAG, CCG, CGG, and CUG repeats. The thermodynamic parameters of helix propagation correlated with the presence of multiple N-N mismatches within CNG RNA duplexes were also determined. Moreover, the total stability of CNG RNA hairpins, as well as the contribution of trinucleotide repeats placed only in the stem or loop regions, was evaluated. The improved thermodynamic parameters allow to predict much more accurately the thermodynamic stabilities and structures of CNG RNAs.     

Influence of sequence context and length on the structure and stability of triplet repeat DNA oligomers

Genetic expansion diseases have been linked to the properties of triplet repeat DNA sequences during replication. The most common triplet repeats associated with such diseases are CAG, CCG, CGG, and CTG. It has been suggested that gene expansion occurs as a result of hairpin formation of long stretches of these sequences on the leading daughter strand synthesized during DNA replication [Gellibolian, R., Bacolla, A., and Wells, R. D. (1997) J. Biol. Chem. 272, 16793-7]. To test the biophysical basis for this model, oligonucleotides of general sequence (CNG)(n), where N = A, C, G, or T and n = 4, 5, 10, 15, or 25, were synthesized and characterized by circular dichroism (CD) spectropolarimetry, optical melting studies, and differential scanning calorimetry (DSC). The goal of these studies was to evaluate the influence of sequence context and oligomer length on their secondary structures and stabilities. The results indicate that all single oligomers, even those as short as 12 nucleotides, form stable hairpin structures at 25 degrees C. Such hairpins are characterized by the presence of N:N mismatched base pairs sandwiched between G:C base pairs in the stems and loops of three to four unpaired bases. Thermodynamic analysis of these structures reveals that their stabilities are influenced by both the sequence of the particular oligomer and its length. Specifically, the stability order of CGG > CTG > CAG > CCG was observed. In addition, longer oligomers were found to be more stable than shorter oligomers of the same sequence. However, a stability plateau above 45 nucleotides suggests that the length dependence reaches a maximum value where the stability of the G:C base pairs can no longer compensate the instability of the N:N mismatches in the stems of the hairpins. The results are discussed in terms of the above model proposed for gene expansion.     

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

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

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

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

RNA hairpin loop stability depends on closing base pair.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXAUAAUAYCC, where X and Y are CG, GC, AU, UA, GU, or UG. A nearest neighbor analysis of the data indicates the free energy change for loop formation at 37 degrees C, delta degrees Gl,37, averages 3.4 kcal/mol for hairpin loops closed with C.G, G.C, and G.U pairs. In contrast, delta G degree l,37 averages 4.6 kcal/mol for loops closed with A.U, U.A, or U.G pairs. Thus the stability of an RNA hairpin depends on the closing base pair. The hairpin with a GA mismatch that is formed by GGCGUAAUAGCC is more stable than the corresponding hairpin with an AA mismatch. Thus hairpin stability also depends on loop sequence. These effects are not included in current algorithms for prediction of RNA 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.     

Mutational analysis of the RNA pseudoknot involved in efficient ribosomal frameshifting in simian retrovirus-1.

Mutational effects on frameshifting efficiency of the RNA pseudoknot involved in ribosomal frameshifting in simian retrovirus-1 (SRV-1) have been investigated. The primary sequence and the proposed secondary structure of the SRV-1 pseudoknot are similar to those of other efficient frameshifting pseudoknots in mouse mammary tumor virus (MMTV) and feline immunodeficiency virus (FIV), where an unpaired adenine nucleotide intercalates between stem 1 and stem 2. In SRV-1 pseudoknot, the adenine nucleotide in between stem 1 and stem 2 has a potential to form an A*U base pair with the last uridine nucleotide in the loop 2, resulting in a continuous A-form helix with coaxially stacked stem 1 and stem 2. To test whether this A*U base pairing and coaxial stacking of stem 1 and stem 2 is absolutely required for efficient frameshifting in SRV-1, a series of mutants changing this potential A.U base pair to either G.C base pair or A.A, A.G, A.C, G.A, G.G mismatch is generated, and their frameshifting efficiencies are investigated in vitro using rabbit reticulocyte lysate translation assay. The frameshifting abilities of these mutant pseudoknots are similar to that of the wild-type pseudoknot, suggesting that the A*U base pair in between stem 1 and stem 2 is not necessary to promote efficient frameshifting in SRV-1. These results reveal that coaxial stacking of stem 1 and stem 2 with a Watson-Crick A.U base pair in between two stems is not a required structural feature of the pseudoknot for promoting efficient frameshifting in SRV-1. Our mutational data suggest that SRV-1 pseudoknot adopts similar structural features common to other efficient frameshifting pseudoknots as observed in MMTV and FIV.     

Stability and mismatch discrimination of locked nucleic acid-DNA duplexes

Locked nucleic acids (LNA; symbols of bases, +A, +C, +G, and +T) are introduced into chemically synthesized oligonucleotides to increase duplex stability and specificity. To understand these effects, we have determined thermodynamic parameters of consecutive LNA nucleotides. We present guidelines for the design of LNA oligonucleotides and introduce free online software that predicts the stability of any LNA duplex oligomer. Thermodynamic analysis shows that the single strand-duplex transition is characterized by a favorable enthalpic change and by an unfavorable loss of entropy. A single LNA modification confines the local conformation of nucleotides, causing a smaller, less unfavorable entropic loss when the single strand is restricted to the rigid duplex structure. Additional LNAs adjacent to the initial modification appear to enhance stacking and H-bonding interactions because they increase the enthalpic contributions to duplex stabilization. New nearest-neighbor parameters correctly forecast the positive and negative effects of LNAs on mismatch discrimination. Specificity is enhanced in a majority of sequences and is dependent on mismatch type and adjacent base pairs; the largest discriminatory boost occurs for the central +C·C mismatch within the +T+C+C sequence and the +A·G mismatch within the +T+A+G sequence. LNAs do not affect specificity in some sequences and even impair it for many +G·T and +C·A mismatches. The level of mismatch discrimination decreases the most for the central +G·T mismatch within the +G+G+C sequence and the +C·A mismatch within the +G+C+G sequence. We hypothesize that these discrimination changes are not unique features of LNAs but originate from the shift of the duplex conformation from B-form to A-form.     

Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Crystal structure, stability and in vitro RNAi activity of oligoribonucleotides containing the ribo-difluorotoluyl nucleotide: insights into substrate requirements by the human RISC Ago2 enzyme

Short interfering RNA (siRNA) duplexes are currently being evaluated as antisense agents for gene silencing. Chemical modification of siRNAs is widely expected to be required for therapeutic applications in order to improve delivery, biostability and pharmacokinetic properties. Beyond potential improvements in the efficacy of oligoribonucleotides, chemical modification may also provide insight into the mechanism of mRNA downregulation mediated by the RNA–protein effector complexes (RNA-induced silencing complex or RISC). We have studied the in vitro activity in HeLa cells of siRNA duplexes against firefly luciferase with substitutions in the guide strand of U for the apolar ribo-2,4-difluorotoluyl nucleotide (rF) [Xia, J. et al. (2006) ACS Chem. Biol., 1, 176–183] as well as of C for rF. Whereas an internal rF:A pair adjacent to the Ago2 (‘slicer’ enzyme) cleavage site did not affect silencing relative to the native siRNA duplex, the rF:G pair and other mismatches such as A:G or A:A were not tolerated. The crystal structure at atomic resolution determined for an RNA dodecamer duplex with rF opposite G manifests only minor deviations between the geometries of rF:G and the native U:G wobble pair. This is in contrast to the previously found, significant deviations between the geometries of rF:A and U:A pairs. Comparison between the structures of the RNA duplex containing rF:G and a new structure of an RNA with A:G mismatches with the structures of standard Watson–Crick pairs in canonical duplex RNA leads to the conclusion that local widening of the duplex formed by the siRNA guide strand and the targeted region of mRNA is the most likely reason for the intolerance of human Ago2 (hAgo2), the RISC endonuclease, toward internal mismatch pairs involving native or chemically modified RNA. Contrary to the influence of shape, the thermodynamic stabilities of siRNA duplexes with single rF:A, A:A, G:A or C:A (instead of U:A) or rF:G pairs (instead of C:G) show no obvious correlation with their activities. However, incorporation of three rF:A pairs into an siRNA duplex leads to loss of activity. Our structural and stability data also shed light on the role of organic fluorine as a hydrogen bond acceptor. Accordingly, UV melting (T_M) data, osmotic stress measurements, X-ray crystallography at atomic resolution and the results of semi-empirical calculations are all consistent with the existence of weak hydrogen bonds between fluorine and the H-N1(G) amino group in rF:G pairs of the investigated RNA dodecamers.