Journal of Nucleic AcidsVolume 2011 (2011), Article ID 513910, 9 pagesdoi:10.4061/2011/513910

Research Article

Influence of a Hairpin Loop on the Thermodynamic Stability of a DNA Oligomer

Jurij Lah, Mojca Seručnik, and Gorazd Vesnaver

Faculty of Chemistry and Chemical Technology, University of Ljubljana, Askerceva 5, 1000 Ljubljana, Slovenia

Received 13 April 2011; Accepted 16 May 2011

Academic Editor: Luis A. Marky

Copyright © 2011 Jurij Lah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

DSC was used to evaluate the mechanism of the thermally induced unfolding of the single-stranded hairpin HP = 5′-CGGAATTCCGTCTCCGGAATTCCG-3′ and its core duplex D (5′-CGGAATTCCG-3′)₂. The DSC melting experiments performed at several salt concentrations were successfully described for HP and D in terms of a three-state transition model (intermediate state) S (unfolded single-stranded state) and two state transition model , respectively. Comparison of the model-based thermodynamic parameters obtained for each HP and D transition shows that in unfolding of HP only the transition is affected by the TCTC loop. This observation suggests that in the intermediate state its TCTC loop part exhibits significantly more flexible structure than in the folded state while its duplex part remains pretty much unchanged.

1. Introduction

Hairpin loops are a common form of nucleic acid secondary structure and are crucial for tertiary structure and function [1]. They are known to play a key role in a number of biological processes such as gene expressions, DNA recombination, and DNA transposition [2–4]. In RNA molecules hairpins act as nucleation sites for RNA folding into final conformations [5–7] and play a critical role in RNA-protein recognition and gene regulation [8, 9]. Furthermore, due to the specificity of probe/target hybridization determined as a match-versus-mismatch discrimination, hairpin DNA oligomer probes have become an important tool in modern biotechnology and diagnostics [10, 11]. The thermodynamics and kinetics of hairpin formation, hairpin binding to complementary nucleic acids, and hairpin-ligand associations have been studied extensively [12–21]. There is no doubt that studies of hairpin-to-coil transitions and hairpin-ligand binding affinity and specificity have greatly enhanced our understanding of structural features and function of the naturally occurring nucleic acids [22, 23]. However, despite extensive biophysical research on the systems involving hairpin structures that produced a number of high-quality explanations and evaluations on properties and behavior of nucleic acids containing hairpin formations, there are still many unresolved questions.

As pointed out by Marky et al. [24] the most suitable hairpin molecules for studying the thermodynamics of their conformational transitions and ligand binding are the single-stranded hairpin molecules. They form stable partially paired duplexes that tend to melt in simple monomolecular transitions. Furthermore, their conformational stability and ligand binding properties are easily compared with those of the corresponding core duplexes. In this way one can evaluate the contributions of the loops to the thermal stability of the hairpins. Despite the simple structure of single-stranded hairpins it is not clear whether their monomolecular folding/unfolding transitions occur in a two-state or multistate manner. Measurements of their thermally induced unfolding transitions followed by UV, CD, and/or fluorescence spectroscopy as a rule result in sigmoidal melting curves suggesting that they may be considered as two-state processes. The same conclusion can be reached also on the basis of DSC measurements performed on the same sample solutions in the older generation of less sensitive DSC instruments (e.g., Microcal MC-2) which resulted in single-peak DSC thermograms. Recent measurements of conformational transitions of DNA quadruplex structures have shown, however, that the sigmoidal shape of UV or CD melting curves may be misleading. Namely, the DSC measurements performed on samples for which sigmoidal UV and CD melting curves were observed using the DSC of the latest generation (CSC, Microcal) resulted in thermograms containing two or three well-distinguished peaks thus indicating that the observed DNA melting process occurs in a multi-state manner [25, 26]. Furthermore, recent T-jump experiments performed on small hairpin molecules have produced a direct evidence that their unfolding transitions involve intermediate structures and thus cannot be considered as two-state processes [19, 27–29].

In our DSC study of the unfolding mechanism and stability of the 5′-CGGAATTCCGTCTCCGGAATTCCG-3′ hairpin we performed the DSC melting experiments on the hairpin and its core duplex, (5′-CGGAATTCCG-3′)₂ (see Figure 1), at several salt concentrations using an extremely sensitive microcalorimeter (CSC). To see to what extent the TCTC loop affects the hairpin unfolding process we attempted to describe for each oligonucleotide the measured DSC thermograms in terms of the simplest possible unfolding model. We derived the corresponding model functions and by fitting them to the experimental data we tested the appropriateness of the suggested models and obtained for each transition the characteristic thermodynamic quantities of transition , , , and . By comparing these values determined for the hairpin and the core duplex we tried to estimate the contribution of the TCTC loop to the stability of the hairpin.

Figure 1: Schematic presentation of the model oligonucleotides: hairpin (HP), duplex (D).

2. Materials and Methods

2.1. Materials

Self-complementary oligonucleotide 5′-CGGAATTCCG-3′ and oligonucleotide 5′-CGGAATTCCGTCTCCGGAATTCCG-3′ that in solution at room temperature form a duplex (D) and a single-stranded hairpin structure (HP), respectively, were purchased HPLC pure from Invitrogen Co., Germany and used without any further purification. Their concentrations in buffer solution (10 mM phosphate buffer and 1 mM Na₂EDTA adjusted to pH = 7.0)) in the presence of 100 mM NaCl were determined at 25°C spectrophotometrically in the Cary Bio 100 UV-spectrophotometer. The molar extinction coefficients were determined using the nearest neighbor data of Cantor et al. [30] for single-stranded DNA at 25°C and the absorbance at 260 nm of thermally unfolded oligonucleotide extrapolated back to 25°C ( M⁻¹cm⁻¹,  M⁻¹cm⁻¹). The phosphate buffer solutions used in all experiment contained 0, 0.1, 0.3, or 1.0 M NaCl.

Differential scanning calorimetry (DSC). Thermally induced unfolding of duplex (D) and hairpin (H) in buffer solutions with different added NaCl concentrations was followed between 5 and 95°C in a Nano-II DSC calorimeter (CSC; UT) at the heating rate of 1°C/min and essentially the same results were obtained from several test-experiments performed at the heating rate of 0.25°C/min. The thermally induced unfolding of both oligonucleotides was monitored in terms of versus thermograms in which the differences between the partial molar heat capacity of the measured oligonucleotide (raw signals corrected for the solvent contributions) and the partial molar heat capacities of the corresponding folded states extrapolated from low temperatures over the whole measured temperature interval, , are normalized for the duplex or hairpin concentration. The total enthalpy of unfolding, ,was obtained from the measured thermograms as the area under the versus curve.

2.2. Analysis of the DSC Thermograms

The thermally induced conformational transitions can be experimentally followed in a model-independent way only by DSC. At relatively low concentrations used in DSC experiments the measured solute-normalized heat capacity of the sample solution, , with the subtracted baseline may be equalized with the oligonucleotide partial molar heat capacity, . Thus, the overall heat effect that accompanies the measured conformational transition from its initial folded state at the temperature to its final unfolded state at can be expressed as Since the enthalpy is the state function, the enthalpy change may be expressed also as where and are the partial molar heat capacities of the folded and unfolded DNA conformation, respectively, is the enthalpy of unfolding at which can be any temperature between and . By choosing where is the melting temperature at which a half of oligonucleotide molecules undergo unfolding transition (2) transforms into where is a model-independent enthalpy of transition at that can be easily determined by the appropriate integration of the experimental and curves as presented in (3).

According to the DSC thermograms of the measured hairpin (HP), its thermally induced unfolding involves at least two conformational transitions (Figure 2). Thus, the simplest suggested model to describe the observed thermal behavior would consist of two consecutive monomolecular transitions: HP (hairpin) I (intermediate state) S (unfolded single-stranded state). The enthalpy, H, of a solution containing an HP sample characterized by the suggested thermal unfolding can be expressed at given and as where and are the corresponding equilibrium constants, the quantities in brackets are the equilibrium molar concentrations of , and , is the number of moles of solvent, and is the number of moles of solute (oligonucleotide) that can be further expressed as: in (6) represents the number of moles of the oligonucleotide in the folded hairpin state, is the number of moles in the intermediate state, is the number of moles in the unfolded single-stranded state and , and are the corresponding partial molar enthalpies of the solvent, solute and folded, intermediate and unfolded oligonucleotide, respectively. By defining the molar fraction,, of the solute species, , as one obtains from (5) that Finally, by introducing into (7) and taking the temperature derivative of the modified (7) one obtains the model function for the measured DSC signal, , expressed as in which at any temperature is the measured (with subtracted baseline), is the partial molar heat capacity of HP extrapolated from low-temperature region over the entire measured temperature interval, (enthalpy of the hairpin to intermediate state transition), (enthalpy of the intermediate state to unfolded single stranded state transition), and .

(a)

(b)

(c)

(d)

Figure 2: DSC thermograms and their model analysis: hairpin (HP) unfolding characterized in terms of a three-state model HP I S (a) and the corresponding fractions of species (b); duplex (D) unfolding characterized in terms of a two-state model D 2S (c) and the corresponding fractions of species (d). In panels (a) and (c) symbols represent experimental data points while lines of the same color correspond to the best-fit model functions ((14) and (18)).

can be obtained experimentally simply by subtracting the hairpin versus curve extrapolated from low- region over the entire measured temperature interval from the corresponding measured versus curve. The model-based , however, can be calculated from the right-hand side term of (8). According to the suggested model (4) the total solute molar concentration, , and the fractions of the solute species present in the solution can be expressed as and , and . Since one obtains from (4) that For the description of the DSC experiment with the model function (8) one needs also the temperature derivatives of and . By using for each transition, , the van’t Hoff relation one obtains Assuming that for each transition, , the corresponding does not depend on the standard free energy of that transition, , can be obtained at any from the integrated form of the Gibbs-Helmholtz relation as where is the temperature at which the values sof species participating in transition are the same. The corresponding equilibrium constant, , is related to as and for the suggested mechanism of the hairpin unfolding (4) it can be easily seen that for each suggested monomolecular transition . Finally, according to the DSC experiments performed at different oligonucleotide concentrations the values appear to be concentration independent thus indicating that one may assume for each transition that and . Using these assumptions and (8)–(13) one can express the model function (14) only in terms of parameters , , and , characteristic for each of the suggested transitions. Their “best fit” values are obtained by fitting the model function (14) to the experimental versus curves. Furthermore, since for each transition, , the corresponding and quantities can be expressed as the “best fit” parameters , and can be used also to obtain the and values at any .

In contrast to HP unfolding, the measured thermally induced duplex (D) to single strand (S) transition appears to be a simpler, all-or-none process that can be described in terms of the total oligonucleotide concentration, , the concentrations of the duplex form and the single strands , the fraction of duplex molecules that undergo the unfolding transition at a given temperature, , and the equilibrium constant interrelated as A similar, though much simpler derivation of the model function than the one presented for unfolding of the hairpin structure (14) leads for the suggested D 2S transition to where is the measured of the sample solution with subtracted baseline, is the heat capacity of the duplex form extrapolated from the low- region over the whole measured temperature interval, and . From the suggested model (16) and (17) it follows that and The corresponding expressions for    and are the same as those shown for each transition in the suggested hairpin unfolding mechanism (12), (13) and (15). Similarly, in deriving (18) the and are assumed to be independent on the oligonucleotide concentration and thus equal to , and .

Inspection of (18) and (19) shows that the model function (18) is expressed in terms of adjustable parameters, , and that can be determined by fitting the model function to the experimental versus curve and further used to determine the , and values at any .

To obtain a set of the “best fit” adjustable parameters , and describing the hairpin and duplex thermal unfolding at each of the added salt concentrations the iterative nonlinear Levenberg-Marquardt regression procedure was used [31]. Furthermore, assuming that for the observed transitions the accompanying quantities do not depend on the added NaCl concentration their values may be determined also from the slopes of the versus curves constructed from the “best fit” and parameters determined at different added salt concentrations [32]. These data can be also used to construct the corresponding versus plots from which the amount of the ions released upon thermal unfolding of the hairpin or duplex structure can be estimated (see discussion, (20)).

3. Results and Discussion

According to the measured DSC thermograms presented in Figure 2 the thermally induced unfolding of the hairpin HP consists of at least two conformational transitions while the one observed for the duplex D occurs in a simpler “all or none” manner. In addition, UV melting experiments (not shown) resulted for HP in biphasic melting curves that exhibit transitions independent on HP concentration (monomolecular transitions) while for D monophasic melting curves dependent on D concentration (nonmonomolecular transition) were observed. Moreover, excellent repeatability of the consecutive measured heating and cooling versus curves and the observed independence of the measured DSC peaks on the applied heating rate (several test experiments) clearly shows that the studied thermal unfolding events may be considered as reversible processes. Model analysis of the measured thermograms shows that the hairpin thermal unfolding can be well described in terms of a three state model involving (hairpin) (intermediate state) (unfolded single-stranded structure) transitions and the corresponding model function (14) characterized for each of the suggested transitions with the corresponding “best fit” adjustable parameters , , and (Table 1). However, analysis of the applied fitting procedure indicates that the parameter is highly correlated to some other adjustable parameters. Thus, to obtain safe estimate of another method of its determination has to be used. Assuming that it does not depend on the simple salt concentration was estimated as a slope of the versus plot (Figure 3(a)) constructed from the “best fit” parameters determined at different NaCl concentrations (Table 1). This method of determining was justified by a good agreement between the values for other transitions obtained by the described fitting procedure and the values determined as the slopes of the corresponding versus plots (Table 1). By using the parameters presented in Table 1 one can calculate for the duplex and hairpin solutions the relative populations of the model-predicted species in the measured temperature interval and at all added salt concentrations (Figure 2). Evidently, the thermal stability of the folded state of the measured duplex and the hairpin is substantially enhanced by increasing the added salt concentration. At low salt concentrations, however, a small fraction of the hairpin molecules undergoes transition into the intermediate state already at physiological temperatures.

Table 1: Thermodynamic parametersa obtained from fitting the model functions ((14) and (18)) to the duplex (D) and hairpin (HP) DSC thermograms presented in Figure 2.

Transition

 b
 c

HP → I
59.2
37

0.34
0.9

I → S
75.0
74

0.37
0.25
1.5

HP → S

117d
111

0.59
2.4

D → 2S
56.5
71
71
0.25
0.30
1.7

Error
±0.2
±2
±2
±0.05
±0.05
±0.2

aUnits: , , ; unless stated otherwise the values are those obtained at ; bobtained from fitting the model function; cobtained as the slope of versus curves (Figure 3(a)); d the total enthalpy of the of the transition was calculated as where and are the melting temperatures of the and transitions and and are the corresponding enthalpies of transition.

(a)

(b)

Figure 3: Estimation of heat capacity changes and number of released Na⁺ ions: (a) was determined for each transition as the slope of the versus plot constructed from the model based and values determined at different salt concentrations; (b) the corresponding versus plots from which the values were determined according to (20).

A standard way of testing the quality of a suggested model is to compare the enthalpy of unfolding determined at a given temperature directly by an appropriate integration of the experimental versus curve (, see (3)) with the corresponding model-based value calculated at the same temperature using the reported “best fit” adjustable parameters. As shown in Table 1 a good agreement was obtained which clearly supports the appropriateness of the suggested unfolding model.

It is well known that DNA unfolding is accompanied by release of counterions. The number of the released ions, , upon each HP and D transition, expressed per mole of oligonucleotide, may be estimated from [33] in which is the melting temperature at a given concentration, , is the corresponding enthalpy of transition at . The values presented in Table 1 were determined from the slopes of the versus . plots (Figure 3(b)).

At any in the measured range of physiological temperatures the difference between the given property characterizing the total unfolding of the hairpin and the duplex (for ex. ) reflects the contribution of the TCTC loop to that property relative to the core duplex (Figure 4, Table 2). Thus, the observed indicates a favorable energy contribution of the TCTC loop to the stability of the hairpin that results, very likely, from the increased number of stacking interactions (in the first place from those occurring at the core duplex-loop connections) [34]. The corresponding is consistent with the highly positive indicating that the unfavorable entropy contribution of the TCTC loop to the hairpin stability arises largely from a substantial disruption of the loop structure accompanying the unfolding of the hairpin. The observed and show, however, that the loop contributions to and may be, to a certain extent, determined also by hydration [35] and electrostatic interactions. The suggests that within the folded hairpin conformation, not only the core duplex but also the TCTC loop are less exposed to water than in the unfolded state. In addition, the observed may be ascribed to a decrease in the surface charge density accompanying the unfolding of the oligonucleotide which is significantly more pronounced in the case of hairpin unfolding. Comparison of the values for , , , , and quantities determined for the and transitions of the hairpin with the corresponding values determined for the transition of the duplex shows that the values for the and transitions are very close (Figure 4, Table 2). Evidently, one may speculate that in the hairpin structure only the transition is influenced by the TCTC loop. In other words, it seems that the observed transition reflects mainly the changes in the TCTC conformation. Thus, the intermediate state may be considered as a state in which the core duplex remains more or less unchanged while the TCTC loop occurs as a more flexible structure characterized by the additional stacking interactions and the freedom of the neighboring water molecules and ions similar to the one in the unfolded state.

Table 2: Difference thermodynamic stability parametersa at 25°C exhibiting the influence of the TCTC loop to the unfolding features of the hairpin forming oligonucleotide (HP).

Transition

(HP → S)-(D → 2S)
2.4
26.6
24.3
0.29
0.7

(I → S)-(D → 2S)
−0.9
0.8
1.7
−0.05
−0.2

Error
±3
±3
±3
±0.07
±0.3

aUnits: (, (); for temperature dependence of , , and see Figure 4.

(a)

(b)

(c)

Figure 4: Thermodynamics of hairpin and duplex unfolding: standard Gibbs free energy (a), standard enthalpy (b), and the corresponding entropy contribution (c) presented for each model predicted hairpin (, ) and duplex () transitions as functions of temperature at [Na⁺] = 0.13M.

To the best of our knowledge this is the first time that a three-state unfolding of a simple hairpin structure, observed by DSC, has been reported and characterized thermodynamically. We believe that the main reason for this is that in most studies of thermal unfolding of hairpins too high starting temperatures have been chosen and therefore the low-temperature transitions have been overlooked.

Acknowledgments

This paper was supported by the Slovenian Research Agency through Grant no. P1-0201 and by the COST action MP0802.

References

P. Brion and E. Westhof, “Hierarchy and dynamics of RNA folding,” Annual Review of Biophysics and Biomolecular Structure, vol. 26, pp. 113–137, 1997.

D. B. Roth, J. P. Menetski, P. B. Nakajima, M. J. Bosma, and M. Gellert, “V(D)J recombination: broken DNA molecules with covalently sealed (hairpin) coding ends in scid mouse thymocytes,” Cell, vol. 70, no. 6, pp. 983–991, 1992.

A. Bhasin, I. Y. Goryshin, and W. S. Reznikoff, “Hairpin formation in Tn5 transposition,” Journal of Biological Chemistry, vol. 274, no. 52, pp. 37021–37029, 1999.

A. K. Kennedy, A. Guhathakurta, N. Kleckner, and D. B. Haniford, “Tn10 transposition via a DNA hairpin intermediate,” Cell, vol. 95, no. 1, pp. 125–134, 1998.

O. C. Uhlenbeck, “Nucleic-acid structure—tetraloops and RNA folding,” Nature, vol. 346, no. 6285, pp. 613–614, 1990.

C. Cheong, G. Varani, and I. Tinoco, “Solution structure of an unusually stable RNA heparin—GGAC(UUCG)GUCC,” Nature, vol. 346, no. 6285, pp. 680–682, 1990.

H. A. Heus and A. Pardi, “Structural features that give rise to the unusual stability of RNA hairpins containing GNRA loops,” Science, vol. 253, no. 5016, pp. 191–194, 1991.

D. E. Draper, “Themes in RNA-protein recognition,” Journal of Molecular Biology, vol. 293, no. 2, pp. 255–270, 1999.

P. Svoboda and A. Di Cara, “Hairpin RNA: a secondary structure of primary importance,” Cellular and Molecular Life Sciences, vol. 63, no. 7-8, pp. 901–918, 2006.

G. Bonnet, S. Tyagi, A. Libchaber, and F. R. Kramer, “Thermodynamic basis of the enhanced specificity of structured DNA probes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 11, pp. 6171–6176, 1999.

Y. Kim, D. Sohn, and W. Tan, “Molecular beacons in biomedical detection and clinical diagnosis,” International Journal of Clinical and Experimental Pathology, vol. 1, pp. 105–116, 2008.

I. E. Scheffler, E. L. Elson, and R. L. Baldwin, “Helix formation by d(TA) oligomers. II. Analysis of the helix-coil transitions of linear and circular oligomers,” Journal of Molecular Biology, vol. 48, no. 1, pp. 145–171, 1970.

J. Gralla and D. M. Crothers, “Free energy of imperfect nucleic acid helices. II. Small Hairpin Loops,” Journal of Molecular Biology, vol. 73, no. 4, pp. 497–511, 1973.

C. A. Haasnoot, C. W. Hilbers, G. A. van der Marel et al., “On loop folding in nucleic acid hairpin-type structures,” Journal of Biomolecular Structure &amp; Dynamics, vol. 3, no. 5, pp. 843–857, 1986.

M. M. Senior, R. A. Jones, and K. J. Breslauer, “Influence of loop residues on the relative stabilities of DNA hairpin structures,” Proceedings of the National Academy of Sciences of the United States of America, vol. 85, no. 17, pp. 6242–6246, 1988.

M. Zuker, “On findings all suboptimal foldings of an RNA molecule,” Science, vol. 244, no. 4900, pp. 48–52, 1989.

D. Rentzeperis, K. Alessi, and L. A. Marky, “Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding,” Nucleic Acids Research, vol. 21, no. 11, pp. 2683–2689, 1993.

M. W. Freyer, R. Buscaglia, B. Nguyen, W. David Wilson, and E. A. Lewis, “Binding of netropsin and 4,6-diamidino-2-phenylindole to an A2T2 DNA hairpin: a comparison of biophysical techniques,” Analytical Biochemistry, vol. 355, no. 2, pp. 259–266, 2006.

S. V. Kuznetsov, C. C. Ren, S. A. Woodson, and A. Ansari, “Loop dependence of the stability and dynamics of nucleic acid hairpins,” Nucleic Acids Research, vol. 36, no. 4, pp. 1098–1112, 2008.

J. Lah, I. Drobnak, M. Dolinar, and G. Vesnaver, “What drives the binding of minor groove-directed ligands to DNA hairpins?” Nucleic Acids Research, vol. 36, no. 3, pp. 897–904, 2008.

J. Lah, N. Carl, I. Drobnak, B. Šumiga, and G. Vesnaver, “Competition of some minor groove binders for a single DNA binding site,” Acta Chimica Slovenica, vol. 53, no. 3, pp. 284–291, 2006.

P. C. Bevilacqua and J. M. Blose, “Structures, kinetics, thermodynamics, and biological functions of RNA hairpins,” Annual Review of Physical Chemistry, vol. 59, pp. 79–103, 2008.

D. Bikard, C. Loot, Z. Baharoglu, and D. Mazel, “Folded DNA in action: hairpin formation and biological functions in prokaryotes,” Microbiology and Molecular Biology Reviews, vol. 74, no. 4, pp. 570–588, 2010.

L. A. Marky, D. Rentzeperis, and J. Ho, “Contribution of loops and nicks to the formation of DNA dumbbells: melting behavior and ligand binding,” Biochemistry, vol. 32, no. 10, pp. 2564–2572, 1993.

I. Prislan, J. Lah, and G. Vesnaver, “Diverse polymorphism of G-quadruplexes as a kinetic phenomenon,” Journal of the American Chemical Society, vol. 130, no. 43, pp. 14161–14169, 2008.

I. Prislan, J. Lah, M. Milanic, et al., “Kinetically governed polymorphism of d(G4T4G3) quadruplexes in solutions,” Nucleic Acids Research, vol. 39, pp. 1933–1942, 2011.

A. Ansari, S. V. Kuznetsov, and Y. Shen, “Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 14, pp. 7771–7776, 2001.

Y. Shen, S. V. Kuznetsov, and A. Ansari, “Loop dependence of the dynamics of DNA hairpins,” Journal of Physical Chemistry B, vol. 105, no. 48, pp. 12202–12211, 2001.

H. Ma, D. J. Proctor, E. Kierzek, R. Kierzek, P. C. Bevilacqua, and M. Gruebele, “Exploring the energy landscape of a small RNA hairpin,” Journal of the American Chemical Society, vol. 128, no. 5, pp. 1523–1530, 2006.

C. R. Cantor, M. M. Warshaw, and H. Shapiro, “Oligonucleotide interactions. 3. Circular dichroism studies of the conformation of deoxyoligonucleotides,” Biopolymers, vol. 9, no. 9, pp. 1059–1077, 1970.

W. H. Press, B. P. Flannery, S. A. Teukolsky, et al., Numerical Recipies, Cambridge University Press, Oxford, UK, 2007.

I. Drobnak, M. Seručnik, J. Lah, and G. Vesnaver, “Stability of a short DNA duplex as a function of temperature: the effect of ΔCp and added salt concentration,” Acta Chimica Slovenica, vol. 54, no. 3, pp. 445–451, 2007.

M. T. Record, C. F. Anderson, and T. M. Lohman, “Thermodynamic analysis of ion effects on the binding and conformational equilibria of proteins and nucleic acids: the roles of ion association or release, screening, and ion effects on water activity,” Quarterly Reviews of Biophysics, vol. 11, no. 2, pp. 103–178, 1978.

C. W. Hilbers, C. A. G. Haasnoot, S. H. de Bruin, J. J. M. Joordens, G. A. Van Der Marel, and J. H. Van Boom, “Hairpin formation in synthetic oligonucleotides,” Biochimie, vol. 67, no. 7-8, pp. 685–695, 1985.

P. M. Vallone and A. S. Benight, “Melting studies of short DNA hairpins containing the universal base 5-nitroindole,” Nucleic Acids Research, vol. 27, pp. 3589–3596, 1999.