18 Coaxial Stacking
Coaxial stacking is the stacking of two base pairs at the terminii of adjacent helices. This stacking aligns the two helices along a common axis. In unimolecular secondary structures, coaxial stacking occurs in multibranch and exterior loops.
This set of nearest neighbor parameters allows for two types of coaxial stacking, flush stacking of helices that are directly adjacent (no intervening unpaired nucleotides) and mismatch-mediated coaxial stacking in which a single mistmatch occurs between the stacked helices. Mismatch-mediated coaxial stacking of helices is only allowed when there is exactly one unpaired nucleotide between the helices that can form a non-canonical pair with a nucleotide on the other side of one of the two helices.
18.1 Flush Coaxial Stacking
In flush coaxial stacking, the free energy change of the coaxial stack is approximated using the helical nearest neighbor parameter Watson-Crick-Franklin or GU as though there was no break in the backbone.
18.2 Mismatch-Mediated Coaxial Stacking
In the case of mismatch-mediated coaxial stacking, there are two adjacent stacks. The stack of the mismatch on the adjacent helix, where there is no break in the backbone, is approximated using the sequence-dependent terminal mismatch parameter. The second stack is the stack of the mismatch on the second helix, where the backbone is not continuous. This stack is approximated using a sequence-independent term of –2.1 kcal/mol for folding free energy change.
18.3 Examples
Flush coaxial stacking
ΔG°37 coaxial stack = ΔG°37(Watson-Crick-Franklin Stack)
ΔG°37 coaxial stack = ΔG°37(AU pair followed by CG pair)
ΔG°37 coaxial stack = –2.24 kcal/mol
Note that at this interface, the terminal AU pair penalty would still apply when calculating the helix stability.
Mismatch-mediated coaxial stacking
ΔG°37 coaxial stack = ΔG°37(Continuous Backbone Stack) + ΔG°37(Discontinuous Backbone Stack)
ΔG°37 coaxial stack = ΔG°37(AU pair followed by GG mismatch) + ΔG°37(Discontinuous Backbone Stack)
ΔG°37 coaxial stack = –0.8 kcal/mol – 2.1 kcal/mol
ΔG°37 coaxial stack = –2.9 kcal/mol
Note that at this interface, the terminal AU pair penalty would still apply when calculating the helix stability.
18.4 Tables
A table summarizing the parameters is available in html format.
18.5 References
The coaxial stacking nearest neighbor parameters for free energy change were reported in:
Mathews, D.H., Sabina, J., Zuker, M. and Turner, D.H. (1999) Expanded sequence dependence of thermodynamic parameters provides improved prediction of RNA secondary structure. J. Mol. Biol., 288, 911-940.
The experimental data for the fit of the parameters were taken from:
- Walter, A.E., Turner, D.H., Kim, J., Lyttle, M.H., Müller, P., Mathews, D.H. and Zuker, M. (1994) Coaxial stacking of helixes enhances binding of oligoribonucleotides and improves predictions of RNA folding. Proc. Natl. Acad. Sci. USA., 91, 9218-9222.
- Kim, J., Walter, A.E. and Turner, D.H. (1996) Thermodynamics of coaxially stacked helices with GA and CC mismatches. Biochemistry, 35, 13753-13761.