Thermodynamic analysis of 5' and 3' single- and 3' double-nucleotide overhangs neighboring wobble terminal base pairs

Thermodynamic parameters are reported for duplex formation of 40 self-complementary RNA duplexes containing wobble terminal base pairs with all possible 3′ single and double-nucleotide overhangs, mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest neighbor analysis, the addition of a single 3′ dangling nucleotide increases the stability of duplex formation up to 1 kcal/mol in a sequence-dependent manner. The addition of a second dangling nucleotide increases the stability of duplexes closed with wobble base pairs in an idiosyncratic manner. The results allow for the development of a nearest neighbor model, which improves the predication of free energy and melting temperature for duplexes closed by wobble base pairs with 3′ single or double-nucleotide overhangs. Phylogenetic analysis of naturally occurring miRNAs was performed. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent on the orientation of the GU closing base pair rather than the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for the 5′ single terminal overhangs adjacent to wobble closing base pairs are also presented.     

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.     

HEN1 recognizes 21-24 nt small RNA duplexes and deposits a methyl group onto the 2' OH of the 3' terminal nucleotide

microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants bear a methyl group on the ribose of the 3′ terminal nucleotide. We showed previously that the methylation of miRNAs and siRNAs requires the protein HEN1 in vivo and that purified HEN1 protein methylates miRNA/miRNA* duplexes in vitro. In this study, we show that HEN1 methylates both miRNA/miRNA* and siRNA/siRNA* duplexes in vitro with a preference for 21–24 nt RNA duplexes with 2 nt overhangs. We also demonstrate that HEN1 deposits the methyl group on to the 2′ OH of the 3′ terminal nucleotide. Among various modifications that can occur on the ribose of the terminal nucleotide, such as 2′-deoxy, 3′-deoxy, 2′-O-methyl and 3′-O-methyl, only 2′-O-methyl on a small RNA inhibits the activity of yeast poly(A) polymerase (PAP). These findings indicate that HEN1 specifically methylates miRNAs and siRNAs and implicate the importance of the 2′-O-methyl group in the biology of RNA silencing.     

Thermodynamics of RNA duplexes modified with unlocked nucleic acid nucleotides

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

Nucleotide sequence of miRNA precursor contributes to cleavage site selection by Dicer

The ribonuclease Dicer excises mature miRNAs from a diverse group of precursors (pre-miRNAs), most of which contain various secondary structure motifs in their hairpin stem. In this study, we analyzed Dicer cleavage in hairpin substrates deprived of such motifs. We searched for the factors other than the secondary structure, which may influence the length diversity and heterogeneity of miRNAs. We found that the nucleotide sequence at the Dicer cleavage site influences both of these miRNA characteristics. With regard to cleavage mechanism, we demonstrate that the Dicer RNase IIIA domain that cleaves within the 3′ arm of the pre-miRNA is more sensitive to the nucleotide sequence of its substrate than is the RNase IIIB domain. The RNase IIIA domain avoids releasing miRNAs with G nucleotide and prefers to generate miRNAs with a U nucleotide at the 5′ end. We also propose that the sequence restrictions at the Dicer cleavage site might be the factor that contributes to the generation of miRNA duplexes with 3′ overhangs of atypical lengths. This finding implies that the two RNase III domains forming the single processing center of Dicer may exhibit some degree of flexibility, which allows for the formation of these non-standard 3′ overhangs.     

Contributions of dangling end stacking and terminal base-pair formation to the stabilities of XGGCCp, XCCGGp, XGGCCYp, and XCCGGYp helixes

The role of stacking in terminal base-pair formation was studied by comparison of the stability increments for dangling ends to those for fully formed base pairs. Thermodynamic parameters were measured spectrophotometrically for helix formation of the hexanucleotides AGGCCUp, UGGCCAp, CGGCCGp, GCCGGCp, and UCCGGAp and for the corresponding pentanucleotides containing a 5'-dangling end on the GGCCp or CCGGp core helix. In 1 M NaCl at 1 X 10(-4) M strands, a 5'-dangling nucleotide in this series increases the duplex melting temperature (Tm) only 0-4 degrees C, about the same as adding a 5'-phosphate. In contrast, a 3'-dangling nucleotide increases the Tm at 1 X 10(-4) M strands 7-23 degrees C, depending on the sequence [Freier, S. M., Burger, B. J., Alkema, D., Neilson, T., & Turner, D. H. (1983) Biochemistry 22, 6198-6206]. These results are consistent with stacking patterns observed in A-form RNA. The stability increments from terminal A.U, C.G, or U.A base pairs on GGCC or a terminal U.A pair on CCGG are nearly equal to the sums of the stability increments from the corresponding dangling ends. This suggests stacking plays a large role in nucleic acid stability. The stability increment from the terminal base pairs in GCCGGCp, however, is about 5 times the sum of the corresponding dangling ends, suggesting hydrogen bonding can also make important contributions.     

Molecular and Functional Characterization of RecD,a Novel Member of the SF1 Family of Helicases,from Mycobacterium tuberculosis

The annotated whole-genome sequence of Mycobacterium tuberculosis revealed the presence of a putative recD gene; however, the biochemical characteristics of its encoded protein product (MtRecD) remain largely unknown. Here, we show that MtRecD exists in solution as a stable homodimer. Protein-DNA binding assays revealed that MtRecD binds efficiently to single-stranded DNA and linear duplexes containing 5′ overhangs relative to the 3′ overhangs but not to blunt-ended duplex. Furthermore, MtRecD bound more robustly to a variety of Y-shaped DNA structures having ≥18-nucleotide overhangs but not to a similar substrate containing 5-nucleotide overhangs. MtRecD formed more salt-tolerant complexes with Y-shaped structures compared with linear duplex having 3′ overhangs. The intrinsic ATPase activity of MtRecD was stimulated by single-stranded DNA. Site-specific mutagenesis of Lys-179 in motif I abolished the ATPase activity of MtRecD. Interestingly, although MtRecD-catalyzed unwinding showed a markedly higher preference for duplex substrates with 5′ overhangs, it could also catalyze significant unwinding of substrates containing 3′ overhangs. These results support the notion that MtRecD is a bipolar helicase with strong 5′ → 3′ and weak 3′ → 5′ unwinding activities. The extent of unwinding of Y-shaped DNA structures was ∼3-fold lower compared with duplexes with 5′ overhangs. Notably, direct interaction between MtRecD and its cognate RecA led to inhibition of DNA strand exchange promoted by RecA. Altogether, these studies provide the first detailed characterization of MtRecD and present important insights into the type of DNA structure the enzyme is likely to act upon during the processes of DNA repair or homologous recombination.     

DNA strand transfer catalyzed by vaccinia topoisomerase: ligation of DNAs containing a 3' mononucleotide overhang

The specificity of vaccinia topoisomerase for transesterification to DNA at the sequence 5′-CCCTT and its versatility in strand transfer have illuminated the recombinogenic properties of type IB topoisomerases and spawned topoisomerase-based strategies for DNA cloning. Here we characterize a pathway of topoisomerase-mediated DNA ligation in which enzyme bound covalently to a CCCTT end with an unpaired +1T nucleotide rapidly and efficiently joins the CCCTT strand to a duplex DNA containing a 3′ A overhang. The joining reaction occurs with high efficiency, albeit slowly, to duplex DNAs containing 3′ G, T or C overhangs. Strand transfer can be restricted to the correctly paired 3′ A overhang by including 0.5 M NaCl in the ligation reaction mixture. The effects of base mismatches and increased ionic strength on the rates of 3′ overhang ligation provide a quantitative picture of the relative contributions of +1 T:A base pairing and electrostatic interactions downstream of the scissile phosphate to the productive binding of an unlinked acceptor DNA to the active site. The results clarify the biochemistry underlying topoisomerase-cloning of PCR products with non-templated 3′ overhangs.     

Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.     

The influence of locked nucleic acid residues on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes

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.     

Propynyl groups in duplex DNA: stability of base pairs incorporating 7-substituted 8-aza-7-deazapurines or 5-substituted pyrimidines

Oligonucleotides incorporating the 7-propynyl derivatives of 8-aza-7-deaza-2′-deoxyguanosine (3b) and 8-aza-7-deaza-2′-deoxyadenosine (4b) were synthesized and their duplex stability was compared with those containing the 5-propynyl derivatives of 2′-deoxycytidine (1) and 2′-deoxyuridine (2). For this purpose phosphoramidites of the 8-aza- 7-deazapurine (pyrazolo[3,4-d]pyrimidine) nucleosides were prepared and employed in solid-phase synthesis. All propynyl nucleosides exert a positive effect on the DNA duplex stability because of the increased polarizability of the nucleobase and the hydrophobic character of the propynyl group. The propynyl residues introduced into the 7-position of the 8-aza-7-deazapurines are generally more stabilizing than those at the 5-position of the pyrimidine bases. The duplex stabilization of the propynyl derivative 4b was higher than for the bromo nucleoside 4c. The extraordinary stability of duplexes containing the 7-propynyl derivative of 8-aza-7- deazapurin-2,6-diamine (5b) is attributed to the formation of a third hydrogen bond, which is apparently not present in the base pair of the purin-2,6-diamine 2′-deoxyribonucleoside with dT.     

Selective microRNA uridylation by Zcchc6 (TUT7) and Zcchc11 (TUT4)

Recent small RNA sequencing data has uncovered 3′ end modification of mature microRNAs (miRNAs). This non-templated nucleotide addition can impact miRNA gene regulatory networks through the control of miRNA stability or by interfering with the repression of target mRNAs. The miRNA modifying enzymes responsible for this regulation remain largely uncharacterized. Here we describe the ability for two related terminal uridyl transferases (TUTases), Zcchc6 (TUT7) and Zcchc11 (TUT4), to 3′ mono-uridylate a specific subset of miRNAs involved in cell differentiation and Homeobox (Hox) gene control. Zcchc6/11 selectively uridylates these miRNAs in vitro, and we biochemically define a bipartite sequence motif that is necessary and sufficient to confer Zcchc6/11 catalyzed uridylation. Depletion of these TUTases in cultured cells causes the selective loss of 3′ mono-uridylation of many of the same miRNAs. Upon TUTase-dependent loss of uridylation, we observe a concomitant increase in non-templated 3′ mono-adenylation. Furthermore, TUTase inhibition in Zebrafish embryos causes developmental defects and aberrant Hox expression. Our results uncover the molecular basis for selective miRNA mono-uridylation by Zcchc6/11, highlight the precise control of different 3′ miRNA modifications in cells and have implications for miRNA and Hox gene regulation during development.     

Flap Endonuclease Activity of Gene 6 Exonuclease of Bacteriophage T7

Flap endonucleases remove flap structures generated during DNA replication. Gene 6 protein of bacteriophage T7 is a 5′–3′-exonuclease specific for dsDNA. Here we show that gene 6 protein also possesses a structure-specific endonuclease activity similar to known flap endonucleases. The flap endonuclease activity is less active relative to its exonuclease activity. The major cleavage by the endonuclease activity occurs at a position one nucleotide into the duplex region adjacent to a dsDNA-ssDNA junction. The efficiency of cleavage of the flap decreases with increasing length of the 5′-overhang. A 3′-single-stranded tail arising from the same end of the duplex as the 5′-tail inhibits gene 6 protein flap endonuclease activity. The released flap is not degraded further, but the exonuclease activity then proceeds to hydrolyze the 5′-terminal strand of the duplex. T7 gene 2.5 single-stranded DNA-binding protein stimulates the exonuclease and also the endonuclease activity. This stimulation is attributed to a specific interaction between the two proteins because Escherichia coli single-stranded DNA binding protein does not produce this stimulatory effect. The ability of gene 6 protein to remove 5′-terminal overhangs as well as to remove nucleotides from the 5′-termini enables it to effectively process the 5′-termini of Okazaki fragments before they are ligated.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.     

Evidence That the DNA Endonuclease ARTEMIS also Has Intrinsic 5′-Exonuclease Activity

ARTEMIS is a member of the metallo-β-lactamase protein family. ARTEMIS has endonuclease activity at DNA hairpins and at 5′- and 3′-DNA overhangs of duplex DNA, and this endonucleolytic activity is dependent upon DNA-PKcs. There has been uncertainty about whether ARTEMIS also has 5′-exonuclease activity on single-stranded DNA and 5′-overhangs, because this 5′-exonuclease is not dependent upon DNA-PKcs. Here, we show that the 5′-exonuclease and the endonuclease activities co-purify. Second, we show that a point mutant of ARTEMIS at a putative active site residue (H115A) markedly reduces both the endonuclease activity and the 5′-exonuclease activity. Third, divalent cation effects on the 5′-exonuclease and the endonuclease parallel one another. Fourth, both the endonuclease activity and 5′-exonuclease activity of ARTEMIS can be blocked in parallel by small molecule inhibitors, which do not block unrelated nucleases. We conclude that the 5′-exonuclease is intrinsic to ARTEMIS, making it relevant to the role of ARTEMIS in nonhomologous DNA end joining.     

Small RNA Sequence Analysis of Adenovirus VA RNA-Derived MiRNAs Reveals an Unexpected Serotype-Specific Difference in Structure and Abundance

Human adenoviruses (HAds) encode for one or two highly abundant virus-associated RNAs, designated VA RNAI and VA RNAII, which fold into stable hairpin structures resembling miRNA precursors. Here we show that the terminal stem of the VA RNAs originating from Ad4, Ad5, Ad11 and Ad37, all undergo Dicer dependent processing into virus-specific miRNAs (so-called mivaRNAs). We further show that the mivaRNA duplex is subjected to a highly asymmetric RISC loading with the 3′-strand from all VA RNAs being the favored strand, except for the Ad37 VA RNAII, where the 5′-mivaRNAII strand was preferentially assembled into RISC. Although the mivaRNA seed sequences are not fully conserved between the HAds a bioinformatics prediction approach suggests that a large fraction of the VA RNAII-, but not the VA RNAI-derived mivaRNAs still are able to target the same cellular genes. Using small RNA deep sequencing we demonstrate that the Dicer processing event in the terminal stem of the VA RNAs is not unique and generates 3′-mivaRNAs with a slight variation of the position of the 5′ terminal nucleotide in the RISC loaded guide strand. Also, we show that all analyzed VA RNAs, except Ad37 VA RNAI and Ad5 VA RNAII, utilize an alternative upstream A start site in addition to the classical +1 G start site. Further, the 5′-mivaRNAs with an A start appears to be preferentially incorporated into RISC. Although the majority of mivaRNA research has been done using Ad5 as the model system our analysis demonstrates that the mivaRNAs expressed in Ad11- and Ad37-infected cells are the most abundant mivaRNAs associated with Ago2-containing RISC. Collectively, our results show an unexpected variability in Dicer processing of the VA RNAs and a serotype-specific loading of mivaRNAs into Ago2-based RISC.     

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T₂-biotin·dT-T₂ loop. The third base was a biotinylated uracil (U_B) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3′ dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3′ end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5′-FITC, or radiolabeled with [γ-³³P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45°C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin–target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.     

Thermodynamic stability of base pairs between 2-hydroxyadenine and incoming nucleotides as a determinant of nucleotide incorporation specificity during replication

We investigated the thermodynamic stability of double-stranded DNAs with an oxidative DNA lesion, 2-hydroxyadenine (2-OH-Ade), in two different sequence contexts (5′-GA*C-3′ and 5′-TA*A-3′, A* represents 2-OH-Ade). When an A*–N pair (N, any nucleotide base) was located in the center of a duplex, the thermodynamic stabilities of the duplexes were similar for all the natural bases except A (N = T, C and G). On the other hand, for the duplexes with the A*–N pair at the end, which mimic the nucleotide incorporation step, the stabilities of the duplexes were dependent on their sequence. The order of stability is T > G > C >> A in the 5′-GA*C-3′ sequences and T > A > C > G in the 5′-TA*A-3′ sequences. Because T/G/C and T/A are nucleotides incorporated opposite to 2-OH-Ade in the 5′-GA*C-3′ and 5′-TA*A-3′ sequences, respectively, these results agree with the tendency of mutagenic misincorporation of the nucleotides opposite to 2-OH-Ade in vitro. Thus, the thermodynamic stability of the A*–N base pair may be an important factor for the mutation spectra of 2-OH-Ade.     

Stability of XGCGCp, GCGCYp, and XGCGCYp helixes: an empirical estimate of the energetics of hydrogen bonds in nucleic acids

The stabilizing effects of dangling ends and terminal base pairs on the core helix GCGC are reported. Enthalpy and entropy changes of helix formation were measured spectrophotometrically for AGCGCU, UGCGCA, GGCGCCp, CGCGCGp, and the corresponding pentamers XGCGCp and GCGCYp containing the GCGC core plus a dangling end. Each 5' dangling end increases helix stability at 37 degrees C roughly 0.2 kcal/mol and each 3' end from 0.8 to 1.7 kcal/mol. The free energy increments for dangling ends on GCGC are similar to the corresponding increments reported for the GGCC core [Freier, S. M., Alkema, D., Sinclair, A., Neilson, T., & Turner, D. H. (1985) Biochemistry 24, 4533-4539], indicating a nearest-neighbor model is adequate for prediction of stabilization due to dangling ends. Nearest-neighbor parameters for prediction of the free energy effects of adding dangling ends and terminal base pairs next to G.C pairs are presented. Comparison of these free energy changes is used to partition the free energy of base pair formation into contributions of "stacking" and "pairing". If pairing contributions are due to hydrogen bonding, the results suggest stacking and hydrogen bonding make roughly comparable favorable contributions to the stability of a terminal base pair. The free energy increment associated with forming a hydrogen bond is estimated to be -1 kcal/mol of hydrogen bond.