Parallel analysis of oligodeoxyribonucleotide (oligonucleotide) interactions. I. Analysis of factors influencing oligonucleotide duplex formation.

A novel method for the analysis of oligonucleotide-oligonucleotide interactions is described. Oligonucleotides of different sequence are synthesised in situ as stripes on the surface of a glass slide (see accompanying paper). Multiple hybridizations are then carried out on each oligonucleotide simultaneously to determine the dependence of oligonucleotide duplex formation on duplex length, base composition, hybridisation solvent and sequence complexity.     

Stability of 3' double nucleotide overhangs that model the 3' ends of siRNA

Thermodynamic parameters are reported for duplex formation in 1 M NaCl for 16 RNA sequences, each containing a core tetramer duplex, GGCC, and a 3' overhang consisting of two bases. The results indicate additional double-helical stability is conferred by the double 3' terminal overhang relative to the single 3' terminal overhang. A nearest-neighbor analysis of the data indicates that the free energy contribution at 37 degrees C of the second base in the double 3' terminal overhang varies from 0 to 0.7 kcal/mol. The second base in the 3' double overhang can contribute nearly the same stability to a duplex as a base pair or a 3' dangling overhang. Stability contribution of a dangling base, two nucleotides removed from the 3' end of a duplex, is dependent upon both the identity of the base as well as that of the dangling base that it neighbors. A second dangling base only increases the stability of the duplex when it is neighboring a 3' purine dangling nucleotide. Furthermore, a second dangling pyrimidine provides a greater contribution to duplex stability than a purine. A nearest-neighbor model was developed to predict the influence of 3' double overhang on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against six sequences with different core duplexes.     

Thermodynamics of triple helix formation: spectrophotometric studies on the d(A)10.2d(T)10 and d(C+3T4C+3).d(G3A4G3).d(C3T4C3) triple helices.

We have stabilized the d(A)10.2d(T)10 and d(C+LT4C+3).d(G3A4G3).d(C3T4C3) triple helices with either NaCl or MgCl2 at pH 5.5. UV mixing curves demonstrate a 1:2 stoichiometry of purine to pyrimidine strands under the appropriate conditions of pH and ionic strength. Circular dichroic titrations suggest a possible sequence-independent spectral signature for triplex formation. Thermal denaturation profiles indicate the initial loss of the third strand followed by dissociation of the underlying duplex with increasing temperature. Depending on the base sequence and ionic conditions, the binding affinity of the third strand for the duplex at 25 degrees C is two to five orders of magnitude lower than that of the two strands forming the duplex. Thermodynamic parameters for triplex formation were determined for both sequences in the presence of 50 mM MgCl2 and/or 2.0 M NaCl. Hoogsteen base pairs are 0.22-0.64 kcal/mole less stable than Watson-Crick base pairs, depending on ionic conditions and base composition. C+.G and T.A Hoogsteen base pairs appear to have similar stability in the presence of Mg2+ ions at low pH.     

DNA sequence recognition by the antitumor antibiotic CC-1065 and analogs of CC-1065

The DNA base pair preferences of the antitumor antibiotic CC-1065 and two analogs of CC-1065 were studied by following the rate of covalent bond formation (N-3 adenine adduct) with DNA oligomers containing the 5'NNTTA* and 5'NNAAA* sequences (N = nucleotide, A* = alkylated adenine). The rate of adduct formation of CC-1065 is greatly affected by DNA base changes at the fourth and fifth positions of the bonding site for the 5'NNAAA sequences, but not the 5'NNTTA sequences. However, an analog of CC-1065 containing the same alkylating moiety as CC-1065, but not the third fused ring system or additional methylene and oxygen substituents, shows similar rates of adduct formation for all sequences. A second analog of CC-1065 containing three fused ring systems, but not the methylene and oxygen substituents of CC-1065, shows rates of adduct formation with the same sequence dependence as CC-1065, but does not distinguish between the sequences to the degree shown by CC-1065. Adduct formation of CC-1065, but not the analogs, competes with a reversibly bound species. Thymine bases to the 3' side of a potentially reactive adenine or a cytosine base at the fifth position from the bonding adenine create reversible binding sites which decrease the rate of adduct formation of CC-1065. The sequence 5'GCGAATT binds CC-1065 only reversibly. This sequence can compete for CC-1065 with covalent bonding sequences if the sites are located in different oligomers, or if the sites are located (overlapped or not overlapped) in the same oligomer. The results of these competitive binding experiments suggest that the transfer of CC-1065 from the reversible binding site to the covalent bonding site with both sites located on a single DNA duplex, not overlapped, occurs through an equilibrium of CC-1065 in solution, not by migration of CC-1065 in the minor groove.     

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)     

Evidence for a DNA triplex in a recombination-like motif: I. Recognition of Watson-Crick base pairs by natural bases in a high-stability triplex

Data are presented on a triplex type with two parallel homologous strands for which triplex formation is almost as strong as duplex formation at least for some sequences and even at pH 7 and 0.2 M NaCl. The evidence mainly rests upon comparing thermodynamic properties of similar systems. A paperclip oligonucleotide d(A12C4T12C4A12) with two linkers C4 obviously can form a triplex with parallel back-folded adenine strand regions, because the single melting transition of this complex splits in two transitions by introducing mismatches only in the third strand region. Respectively, a hairpin duplex d(A12C4T12) and a single strand d(A12) form a triplex as a 1:1 complex in which the second adenine strand is parallel oriented to the homologous one in the Watson-Crick paired duplex. In this system the melting temperature T(m) of the triplex is practically the same as that of the duplex d(A12)-d(T12), at least within a complex concentration range of 0.2-4.0 microM. The melting behaviour of complexes between triplex stabilizing ligand BePI and the system hairpin duplex plus single strand supports the triplex model. Non-denaturing gel electrophoresis suggests the existence of a triplex for a system in which five of the twelve A-T*A base triads are substituted by C-G*C base triads. The recognition between any substituted Watson-Crick base pair (X-Y) in the hairpin duplex d(A4XA7C4T7YT4) and the correspondingly replaced base (Z) in the third strand d(A4ZA7) is mutually selective. All triplexes with matching base substitutions (Z = X) have nearly the same stability (T(m) values from 29 to 33.5 degrees C), whereas triplexes with non-matching substitutions (Z not equal X) show a clearly reduced stability (T(m) values from 15 to 22 degrees C) at 2microM equimolar oligonucleotide concentration. Most nucleic acid triple helices hitherto known are limited to homopurine-homopyrimidine sequences in the target duplex. A stable triplex formation is demonstrated for inhomogeneous sequences tolerating at least 50% pyrimidine content in the homologous strands. On the basis of the surprisingly similar thermodynamic parameters for duplex and triplex, and of the fact that this triplex type seems to be more stable than many other natural DNA triplexes known, and on the basis of semiempirical and molecule mechanical calculations, we postulate bridging interactions of the third strand with the two other strands in the triplex according to the recombination motif. This triplex, denoted by us 'recombination-like form', tolerates heterogeneous base sequences.     

Effect of base-pair sequence on the conformations and thermally induced transitions in oligodeoxyribonucleotides containing only AT base pairs

Tm curves, CD spectra, and kinetics results of the self-complementary DNA dodecamers d(A6T6), d(A3T3A3T3), d(A2T2A2T2A2T2), d(ATATATATATAT), and d(T6A6) demonstrate that the thermal transitions of these oligomers at low salt concentration involve a hairpin intermediate. At high salt concentrations (greater than 0.1 M Na+) only a duplex to denatured-strand transition appears to occur. The temperature and salt-concentration regions of the transitions are very sequence dependent. Alternating-type AT sequences have a lower duplex stability and a greater tendency to form hairpins than sequences containing more nonalternating AT base pairs. Of the two nonalternating sequences, d(T6A6) is significantly less stable than d(A6T6). Both oligomers have CD curves that are very similar to the unusual CD spectrum of poly(dA).poly(dT). The Raman spectra of these two oligomers are also quite similar, but at low temperature, small intensity differences in two backbone modes and three nucleoside vibrations are obtained. The hairpin to duplex transition for the AT dodecamers was examined by salt-jump kinetics measurements. The transition is faster than transitions for palindromic-sequence oligomers containing terminal GC base pairs. Stopped-flow kinetics studies indicate that the transition is second order and has a relatively low activation energy. The reaction rate increases with increasing ionic strength. These results are consistent with a three-step mechanism for the hairpin to duplex reaction: (i) fraying of the hairpin oligomers' terminal base pairs, (ii) a rate-determining bimolecular step involving formation of a cruciform-type intermediate from two hairpin oligomers with open terminal base pairs, and (iii) base-pair migration and formation in the intermediate to give the duplex.     

Azole substituted oligonucleotides promote antiparallel triplex formation at non-homopurine duplex targets.

The ability of certain azole substituted oligodeoxy-ribonucleotides to promote antiparallel triple helix formation with duplex targets having CG or TA interruptions in the otherwise homopurine sequence was examined. 2'-Deoxyribonucleosides of the azoles, which include pyrazole, imidazole, 1,2,4-triazole and 1,2,3,4-tetrazole were synthesized using the stereo-specific sodium salt glycosylation procedure. These nucleosides were successfully incorporated using solid-support, phosphoramidite chemistry, into oligonucleotides designed to interact with the non-homopurine duplex targets. The interaction of these modified oligonucleotides with all four possible base pairs was evaluated and compared to similar data for a series of natural oligonucleotides. The oligonucleotides containing simple azoles enhanced the triplex forming ability considerably at non-homopurine targets. Binding of these modified oligonucleotides to duplex targets containing TA inversion sites was particularly noteworthy, and compare favorably to unmodified oligonucleotides for binding to duplex targets containing CG as well as TA base pairs. The selectivity exhibited by certain azoles is suggestive of base pair specific interactions. Thus, the azoles evaluated during this study show considerable promise for efforts to develop generalized triplex formation at non-homopurine duplex sequences.     

Recognition of 5-methyl-CG and CG base pairs in duplex DNA with high stability using antiparallel-type triplex-forming oligonucleotides with 2-guanidinoethyl-2′-deoxynebularine

The formation of triplex DNA is a site-specific recognition method that directly targets duplex DNA. However, triplex DNA formation is generally formed for the GC and AT base pairs of duplex DNA, and there are no natural nucleotides that recognize the CG and TA base pairs, or even the 5-methyl-CG (^5mCG) base pair. Moreover, duplex DNA, including ^5mCG base pairs, epigenetically regulates gene expression in vivo, and thus targeting strategies are of biological importance. Therefore, the development of triplex-forming oligonucleotides (TFOs) with artificial nucleosides that selectively recognize these base pairs with high affinity is needed. We recently reported that 2′-deoxy-2-aminonebularine derivatives exhibited the ability to recognize ^5mCG and CG base pairs in triplex formation; however, this ability was dependent on sequences. Therefore, we designed and synthesized new nucleoside derivatives based on the 2′-deoxy-nebularine (dN) skeleton to shorten the linker length connecting to the hydrogen-bonding unit in formation of the antiparallel motif triplex. We successfully demonstrated that TFOs with 2-guanidinoethyl-2′-deoxynebularine (guanidino-dN) recognized ^5mCG and CG base pairs with very high affinity in all four DNA sequences with different adjacent nucleobases of guanidino-dN as well as in the promoter sequences of human genes containing ^5mCG base pairs with a high DNA methylation frequency.     

Temperature dependence of thermodynamic properties for DNA/DNA and RNA/DNA duplex formation.

A clear difference in the enthalpy changes derived from spectroscopic and calorimetric measurements has recently been shown. The exact interpretation of this deviation varied from study to study, but it was generally attributed to the non-two-state transition and heat capacity change. Although the temperature-dependent thermodynamics of the duplex formation was often implied, systemic and extensive studies have been lacking in universally assigning the appropriate thermodynamic parameter sets. In the present study, the 24 DNA/DNA and 41 RNA/DNA oligonucleotide duplexes, designed to avoid the formation of hairpin or slipped duplex structures and to limit the base pair length less than 12 bp, were selected to evaluate the heat capacity changes and temperature-dependent thermodynamic properties of duplex formation. Direct comparison reveals that the temperature-independent thermodynamic parameters could provide a reasonable approximation only when the temperature of interest has a small deviation from the mean melting temperature over the experimental range. The heat capacity changes depend on the base composition and sequences and are generally limited in the range of -160 to approximately -40 cal.mol-1.K-1 per base pair. In contrast to the enthalpy and entropy changes, the free energy change and melting temperature are relatively insensitive to the heat capacity change. Finally, the 16 NN-model free energy parameters and one helix initiation at physiological temperature were extracted from the temperature-dependent thermodynamic data of the 41 RNA/DNA hybrids.     

Crystal structures of coaxially stacked kissing complexes of the HIV-1 RNA dimerization initiation site

We describe the crystal structures of the RNA dimerization initiation sites (DIS) of HIV-1 subtypes A and B. Both molecules adopt a hairpin conformation, with loop sequences consisting of 272-AGGUGCACA-280 and 272-AAGCGCGCA-280, respectively. This includes a six-base self-complementary stretch (underlined) that allows homodimerization through the formation of a loop-loop, or 'kissing-loop', complex. The DISs for the two sequences have identical conformations, and the two interacting hairpins show a perfect coaxial alignment. The conserved purines, A272 and R273, are stacked in a bulged-out conformation and symmetrically join the upward and downward strands of each hairpin by crossing the helix major groove. There is good agreement between these structures and previous results from chemical probing in solution, as well as with differences in magnesium dependence for dimerization. The overall shape of the kissing-loop complex is very similar to that of the previously determined subtype A DIS duplex form. Unexpectedly, the purine R273 is the only base seen at a different position and is responsible for the difference in topology between the two forms. We propose that the transition from kissing-loop duplex could occur by a recombination mechanism based on symmetrical chain cleavage at R273 of each hairpin and subsequent cross-religation, rather than by base-pair melting and subsequent reannealing.     

Solution conformation of an RNA hairpin loop.

The hairpin conformation adopted by the RNA sequence 5'GCGAUUUCUGACCGCC3' has been studied by one- and two-dimensional NMR spectroscopy. Exchangeable imino spectra in 60 mM Na+ indicate that the hairpin has a stem of six base pairs (indicated by boldface type) and a loop of three nucleotides. NOESY spectra of nonexchangeable protons confirm the formation of the stem region. The duplex has an A-conformation and contains an A.C apposition; a G.U base pair closes the loop region. The stem nucleotides have C3'-endo sugar conformations, as expected of an A-form duplex, whereas the three loop nucleotides adopt C2'-endo sugar puckers. Stacking within the loop, C8 upon the sugar of U7, stabilizes the structure. The pH dependence of both the exchangeable and nonexchangeable NMR spectra is consistent with the formation of an A+.C base pair, protonated at the N1 position of adenine. The stability of the hairpin was probed by using absorbance melting curves. The hairpin structure with the A+.C base pair is about +2 kcal/mol less stable in free energy at 37 degrees C than the hairpin formed with an A.U pair replacing the A+.C pair.     

Long range molecular wire behaviour in a metal complex of DNA

M-DNA is a complex of metal ions such as Zn(2+) with duplex DNA. Previous results showed that the fluorescence of a donor fluorophore was quenched when an acceptor fluorophore was placed at the opposite end of a short M-DNA duplex. In order to investigate further the molecular wire behaviour of M-DNA, 30-mer duplexes were constructed with fluorescein as donor and rhodamine, pyrene and the cyanine dyes, Cy5 and Cy5.5 as acceptors. Good quenching was observed in all cases even though the efficiency of resonance energy transfer was calculated to be < 5%. The distance dependence of quenching was investigated by preparing doubly-labelled duplexes ranging in length from 20 to 1,000 base pairs. Upon formation of M-DNA significant quenching of the fluorescence of the donor fluorophore was observed in duplexes up to 500 base pairs in length. The amount of quenching decreased with increasing length of the duplexes with a shallow distance dependence. The results are consistent with an electron transfer mechanism in which the electron hops between metal centers. This process can occur efficiently over long distances.     

Proton NMR study and conformational analysis of d(CGT), d(TCG) and d(CGTCG) in aqueous solution. The effect of a dangling thymidine and of a thymidine mismatch on DNA mini-duplexes

Proton NMR studies of d(CGT), d(TCG) and d(CGTCG) were carried out at 300 and 500 MHz. The temperature and concentration dependence of the chemical shifts of various resonances indicates duplex formation only in the cases of d(TCG) and d(CGTCG). It is concluded that d(TCG) forms a mini-duplex stabilized by a 5'-dangling thymine base. Thermodynamic parameters of the duplex-to-coil equilibrium of the d(TCG) duplex are: delta H0 = -22.3 kcal/mol and delta S0 = -70 cal/mol. K, which correspond to approximately 40% duplex formation at 0 degrees C in a 2 mM nucleotide solution. Comparison of these data with thermodynamic parameters given earlier [Borer, P.N., Dengler, B., Tinoco, I. and Uhlenbeck, O.C. (1974) J. Mol. Biol. 86, 843-853] leads to the conclusion that the dangling base stabilization observed here is approximately equivalent to the stabilization caused by one or two additional A . T base pairs. The chemical shift behaviour of various resonances in d(CGTCG) indicates duplex formation without looping out of the thymine bases. The T X T mismatch does not seem to disturb the helical structure to a large extent. Analysis of the vicinal proton-proton coupling constants of the three compounds yielded geometrical data for the sugar rings. The data are interpreted in terms of N and S pseudorotational ranges. It is shown that a distinct conformation-transmission effect is exerted by the guanosine residues in a 5'----3' direction.     

d-CpCpGpG and d-GpGpCpC self-complementary duplexes: Nmr studies of the helix-coil transition

We have monitored the helix-coil transition of the self-complementary d-CpCpGpG and d-GpGpCpC sequences (20mM strand concentration) at the base pairs, sugar rings, and backbone phosphates by 360-MHz proton and 145.7-MHz phosphorus nmr spectroscopy in 0.1M phosphate solution between 5 and 95°C. The guanine 1-imino Watson-Crick hydrogen-bonded protons, characteristic of the duplex state, are observed below 10°C, with solvent exchange occurring by transient opening of the tetranucleotide duplexes. The cytosine 4-amino Watson-Crick hydrogen-bonded protons resonate 1.5 ppm downfield from the exposed protons at the same position in the tetranucleotide duplexes, with slow exchange indicative of restricted rotation about the C-N bond below 15°C. The guanine 2-amino exchangeable protons in the tetranucleotide sequence exhibit very broad resonances at low temperatures and narrow average resonances above 20°C, corresponding to intermediate and fast rotation about the C-N bond, respectively. Solvent exchange is slower at the amino protons compared to the imino protons since the latter broaden out above 10°C. The well-resolved nonexchangeable base proton chemical shifts exhibit helix-coil transition midpoints between 37 and 42°C. The transition midpoints and the temperature dependence of the chemical shifts at low temperatures were utilized to differentiate between resonances located at the terminal and internal base pairs while the H-5 and H-6 doublets of individual cytosines were related by spin decoupling studies. For each tetranucleotide duplex, the cytosine H-5 resonances exhibit the largest chemical shift change associated with the helix-coil transition, a result predicted from calculations based on nearest-neighbor atomic diamagnetic anisotropy and ring current contributions for a B-DNA duplex. There is reasonable agreement between experimental and calculated chemical shift changes for the helix-coil transition at the internal base pairs but the experimental shifts exceed the calculated values at the terminal base pairs due to end-to-end aggregation at low temperatures. Since the guanine H-8 resonances of the CpCpGpG and d-CpCpGpG sequences exhibit upfield shifts of 0.6–0.8 and <0.1 ppm, respectively, on duplex formation, these RNA and DNA tetranucleotides with the same sequence must adopt different base-pair overlap geometries. The large chemical shift changes associated with duplex formation at the sugar H-1′ triplets are not detected at the other sugar protons and emphasize the contribution of the attached base at the 1′ position. The coupling sum between the H-1′ and the H-2′ and H-2″ protons equals 15–17 Hz at all four sugar rings for the d-CpCpGpG and d-GpGpCpC duplexes (25°C), consistent with a C-3′ exo sugar ring pucker for the deoxytetranucleotides in solution. The temperature dependent phosphate chemical shifts monitor changes in the ω,ω′ angles about the O-P backbone bonds, in contrast to the base-pair proton chemical shifts, which monitor stacking interactions.     

Thermodynamic comparison of the base pairs formed by the carcinogenic lesion O6-methylguanine with reference both to Watson-Crick pairs and to mismatched pairs

A set of 10 non-self-complementary nonadeoxyribonucleoside octaphosphates, d(GGTTXTTGG) and d(CCAAYAACC), where X and Y are A, C, G, T, or O6MeG, has been synthesized by a large-scale, automated, phosphoramidite procedure. Purification was effected by reversed-phase HPLC, and the base composition was verified by analytical HPLC after enzymatic degradation to the constituent deoxynucleosides. This set of molecules was designed to allow evaluation of the nearest-neighbor dependence of each base pair. The thermal stability, expressed as Tmax, of each duplex containing one of the O6MeG base pairs, a Watson-Crick pair, or one of the mismatches possible with this set of molecules was determined over a concentration range of 5.7-200 microM. From these data the delta H degree, delta S degree, and delta G degree of each combination were calculated. In general, the relative thermal stabilities observed for the O6-methylguanine combinations confirm our previous findings that the most stable base pair is formed with cytosine rather than thymine and that all O6MeG pairs are much weaker than Watson-Crick base pairs [Kuzmich, S., Marky, L. A., & Jones, R. A. (1983) Nucleic Acids Res. 11, 3393-3404; Gaffney, B. L., Marky, L. A., & Jones, R. A. (1984) Biochemistry 23, 5686-5691]. Moreover, the nine combinations containing O6-methylguanine are all of similar thermal stability, cover a much smaller range in Tmax than do the mismatches, and show little sequence dependence.     

Base- and sequence-dependent binding of aristololactam beta-D-glucoside to deoxyribonucleic acid

The dependence on base-pair composition and sequence specificity of the (aristololactam beta-D-glucoside)-DNA interaction was examined by spectrophotometric, spectrofluorometric, spectropolarimetric, thermal melting, thermodynamic, and viscometric studies. Binding of this alkaloid to various natural and synthetic DNAs was dependent upon the base composition and sequences of DNA. The binding parameters obtained from spectrophotometric analysis, according to an excluded-site model, indicated a relatively high affinity of the alkaloid binding to GC-rich DNA and alternating GC polymer. This affinity was further evidenced by the quenching of fluorescence intensity, decrease in quantum yield, and perturbations in circular dichroic spectrum. The alkaloid stabilized all DNAs against thermal denaturation. The temperature dependence of the binding constants was used to estimate the thermodynamic parameters involved in the complex formation of the alkaloid with various DNAs. The negative enthalpy and entropy change increased with increasing GC content of DNA and also compensated one another to produce a relatively small Gibbs free energy change. Viscometric studies showed that in the strong binding region the increase of contour length of DNA depended strongly on its base composition and sequence of bases, being larger for GC-rich DNA and alternating GC polymer. On the basis of these observations, it is concluded that the alkaloid binds to DNA by a mechanism of intercalation and exhibits considerable specificity toward alternating GC polymer.     

Comprehensive thermodynamic analysis of 3' double-nucleotide overhangs neighboring Watson-Crick terminal base pairs

Thermodynamic parameters are reported for duplex formation of 48 self-complementary RNA duplexes containing Watson–Crick terminal base pairs (GC, AU and UA) with all 16 possible 3′ double-nucleotide overhangs; mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest-neighbor analysis, the addition of a second dangling nucleotide to a single 3′ dangling nucleotide increases stability of duplex formation up to 0.8 kcal/mol in a sequence dependent manner. Results from this study in conjunction with data from a previous study [A. S. O'Toole, S. Miller and M. J. Serra (2005) RNA, 11, 512.] allows for the development of a refined nearest-neighbor model to predict the influence of 3′ double-nucleotide overhangs on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against five oligomers with various core duplex sequences. Phylogenetic analysis of naturally occurring miRNAs was performed to support our results. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent upon the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for 3′ single terminal overhangs adjacent to a UA pair are also presented.     

Transition characteristics and thermodynamic analysis of DNA duplex formation: a quantitative consideration for the extent of duplex association

Transition characteristics and thermodynamic properties of the single-stranded self-transition and the double-stranded association were investigated and analyzed for 9-, 15- and 21-bp non-self-complementary DNA sequences. The multiple transition processes for the single-stranded self-transition and the double-stranded association were further put forth. The experimental results confirmed that the double-stranded association transition was generally imperfect and the thermodynamic properties of the single-stranded self-transition would exert an influence on a duplex formation. Combining ultraviolet melting experiments in various molar ratios, the extent of duplex association was estimated for three double-stranded DNAs. In our experimental range, the extent of duplex association decreases with increasing the number of base pairs in DNA sequences, which suggest that the short oligonucleotides may proceed in a two-state transition while the long oligonucleotides may not. When the extent of duplex association was considered, the true transition enthalpies of a duplex formation derived from UV and differential scanning calorimetry measurements were in good agreement.