Salt dependence of nucleic acid hairpin stability - PubMed (original) (raw)

Salt dependence of nucleic acid hairpin stability

Zhi-Jie Tan et al. Biophys J. 2008 Jul.

Abstract

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na(+) and Mg(2+) in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na(+) and Mg(2+) concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.

PubMed Disclaimer

Figures

FIGURE 1

FIGURE 1

Illustrations of the structural model for single-stranded nucleotide chain used in the work. (A) The virtual-bond model for nucleotide chain (–85). (B) The coarse-grained chain model where P atoms are represented by a series of spheres with charge −e (electronic charge) at the centers (79,89). (C) A three-dimensional chain conformation for the coarse-grained chain model where coordinates of P atoms are produced from the virtual bond model on diamond lattice (–85). (D) The RNA and DNA hairpin sequences used in our calculations. The thermodynamic parameters at 1 M NaCl, and experimental references for these hairpins are listed in Table 2. The values i and j are the numbers of nucleotides of U's and T's in RH6 and DH1 loops, respectively.

FIGURE 2

FIGURE 2

The chain conformational ensemble produced from the virtual bond model on diamond lattice (–85). (A) The number Ω(N, x) of all conformations with end-to-end distance x; (B) the number ω(N, x) of lattice nodes visited by the flexible chain end with end-to-end distance x; (C) the number of chain conformations with x averaged to one direction formula image and (D) the number Ωcoil of chain coil states, which includes all the conformations produced by the self-avoiding walk on a diamond lattice.

FIGURE 3

FIGURE 3

The calculated loop formation free energy Δ_G_ as functions of end-to-end distance x for different chain-lengths N for the electrostatics-free case (calculated through Eq. 6), i.e., by assuming chain backbone is neutral. The dotted line denotes the DNA hairpin loop free energy at 1 M NaCl for a 13-nt chain (loop size = 11-nt) (36).

FIGURE 4

FIGURE 4

The partition function z(N, x) averaged for one chain conformation as a function of end-to-end distance x for different Na+ (A_–_C) and Mg2+ (D_–_F) concentrations. The chain lengths are 5-nt (A and D), 9-nt (B and E), and 13-nt (C and F). The averaging is over Ω0 (lines) and Ω0/2 (symbols) randomly selected chain conformations except for the fully stretched conformation. The agreement between lines and symbols show that the calculations are rather stable for different (large) numbers of conformations sampled. The Ω0 values used for calculations are listed in Table 1. For the comparisons, we also show z(N, x) calculated from the PB theory (dotted lines).

FIGURE 5

FIGURE 5

The three-dimensional loop free-energy landscapes Δ_G_/k_B_T as functions of end-to-end distance x for different chain lengths N and different [Na+] (A_–_D) and [Mg2+] (E_–_H). For the comparisons, we also show the loop free energy for different [Mg2+] computed from PB theory (I_–_L). The three-dimensional plots are produced from the empirical formulas for Na+ (Eq. 9), Mg2+ (Eq. 12), and for Mg2+ with the PB treatment (

Data S1

, Eq. S1). The red and blue colors represent the low and high loop free energies, respectively. In

Data S1

, Figs. S11–S13, we shows the comparisons between the calculated results and the respective empirical formulas.

FIGURE 6

FIGURE 6

The comparisons between the calculated loop free energy Δ_G_ for two typical ionic conditions: 0.01 M MgCl2 and 1 M NaCl. The chain lengths are N = 13, 9, and 5 nucleotides, respectively.

FIGURE 7

FIGURE 7

The calculated loop free energies Δ_G_ as functions of loop size (N – 2) for different [Na+] and [Mg2+]. (A and D) Hairpin loop; (B and E) bulge loop; and (C and F) internal loop. The symbols are experimental data for RNA (♦) (27) and DNA (□) (36) at 1 M NaCl, and the lines are calculated from Eqs. 15–17, respectively. From the bottom to top, ion concentrations are: (A_–_C) [Na+] = 1 M, 0.1 M, 0.01 M, and 0.001 M; (D_–_F) [Mg2+] = 0.1 M, 0.01 M, 0.001 M, and 0.0001 M, respectively.

FIGURE 8

FIGURE 8

The RNA hairpin folding free energy Δ_G_°37 (A) and melting temperature _T_m (B) as functions of [Na+] for three sequences RH1, RH2, and RH3. (Symbols) Experimental data: ⋄, RH3 in Na+ solution (20); ▴, RH3 in Na+ solution (95); ♦, RH2 in Na+ solution (49); and □, RH1 in Na+ solution (49). (Solid lines) Predictions with salt-dependent loop free energy (Δ_G_[Na+]); (dotted lines) predictions with invariable loop free energy at 1 M NaCl (Δ_G_[1 M Na+]). Hairpins RH1, RH2, and RH3 are shown in Fig. 1 and Table 2.

FIGURE 9

FIGURE 9

The RNA hairpin folding free energy Δ_G_°25 (A_–_C) and melting temperature _T_m (D_–_F) as functions of [Mg2+] for different [Na+]. Three sequences used are RH3, RH4, and RH5, which are shown in Fig. 1 and Table 2. (Symbols) Experimental data: (A) ♦, RH3 in 0.054 M [Na+] (20); (E) ♦, RH4 in 0.1 M [K+] (96); (F) ♦, RH5 in 0.1 M [K+] (97). (Solid lines) Predictions with salt (Na+/Mg2+) dependent loop free energy Δ_G_[Na+/Mg2+]; (dotted lines) predictions with invariable loop free energy at 1 M NaCl (Δ_G_[1 M Na+]). In panels A_–_C, we calculate Δ_G_°25 instead of Δ_G_°37 because of the available experimental data at 25°C. In the calculations, we use the same salt empirical formula for K+ as that for Na+, since Na+ and K+ have the similar electrostatic properties (36). For the comparisons, we also show the predictions with the loop free energy from the PB empirical formulas (

Data S1

, Eqs. S1 and S3) (dashed lines).

FIGURE 10

FIGURE 10

The melting temperatures _T_m of DNA (DH1, A and B) and RNA (RH6, C) hairpins as a function of loop size for different ionic conditions: 0.1 M Na+ (A), 2.5 mM and 33 mM Mg2+ (B), and 2.5 mM Mg2+ (C). (Lines) Predictions from our theory; (symbols) experimental data (93).

Similar articles

Cited by

References

    1. Bloomfield, V. A., D. M. Crothers, and I. Tinoco, Jr. 2000. Nucleic Acids: Structure, Properties and Functions. University Science Books, Sausalito, CA.
    1. Tinoco, I., and C. Bustamante. 1999. How RNA folds. J. Mol. Biol. 293:271–281. - PubMed
    1. Anderson, C. F., and T. M. Record. 1995. Salt-nucleic acid interactions. Annu. Rev. Phys. Chem. 46:657–700. - PubMed
    1. Bloomfield, V. A. 1997. DNA condensation by multivalent cations. Biopolymers. 44:269–282. - PubMed
    1. Brion, P., and E. Westhof. 1997. Hierarchy and dynamics of RNA folding. Annu. Rev. Biophys. Biomol. Struct. 26:113–137. - PubMed

Publication types

MeSH terms

Substances

Grants and funding

LinkOut - more resources