The Turner Energy Model is the standard nearest-neighbor empirical framework that assigns experimentally derived thermodynamic parameters—enthalpy and entropy changes—to RNA base pair stacks, loops, and other structural motifs. By summing these individual contributions, the model calculates the total free energy change of a given RNA secondary structure, enabling the prediction of the most thermodynamically stable fold via free energy minimization algorithms.
Glossary
Turner Energy Model

What is the Turner Energy Model?
The standard empirical framework for predicting RNA secondary structure stability from sequence using nearest-neighbor thermodynamic parameters.
Developed through decades of optical melting experiments on short synthetic oligoribonucleotides, the model parameterizes the stability of Watson-Crick and wobble base pairs based on their sequence context. It forms the computational backbone of widely used software suites like RNAstructure and ViennaRNA, providing the energy rules that drive both Minimum Free Energy (MFE) predictions and partition function calculations for ensemble analysis.
Core Thermodynamic Parameters
The Turner model decomposes RNA folding energetics into a set of empirically derived, sequence-dependent parameters. Each structural motif—from stacked base pairs to hairpin loops—contributes an additive free energy term, enabling the calculation of total folding stability.
Nearest-Neighbor Base Pair Stacks
The foundational energy unit of the Turner model. Rather than treating base pairs independently, the model assigns a unique free energy change to each adjacent stacking interaction (e.g., 5'-AU-3' stacked on 5'-UA-3'). These parameters, derived from optical melting experiments on short duplexes, capture the sequence-dependent enthalpic and entropic contributions of base stacking, hydrogen bonding, and solvation effects. The total helical stability is the sum of all nearest-neighbor increments.
Loop Free Energy Penalties
Unpaired regions impose a thermodynamic cost that opposes folding. The Turner model categorizes and parameterizes these motifs distinctly:
- Hairpin loops: Destabilizing penalty that scales with loop size, with a minimum for tetraloops (e.g., GNRA, UNCG) which are exceptionally stable.
- Internal loops & bulges: Penalized based on size, asymmetry, and sequence-dependent closure energetics.
- Multi-branch loops: Modeled using a linear approximation with an intercept penalty plus a per-helix branching term.
Dangling End Contributions
Unpaired nucleotides at the termini of helices can stack coaxially on the adjacent base pair, providing additional stabilization. The Turner model includes parameters for 5' and 3' dangling ends, where the identity of the dangling nucleotide and the closing base pair determine the magnitude of the favorable free energy increment. This effect is critical for accurately predicting the stability of short helices and multi-loop junctions where coaxial stacking mimics continuous helices.
Mismatch and Tandem Mismatch Energetics
Non-canonical base pairs within a helix are not simply destabilizing errors; they have specific thermodynamic signatures. The Turner model includes parameters for single mismatches (e.g., G·U wobble pairs) and tandem mismatches (two consecutive non-canonical pairs). Tandem mismatches can be surprisingly stabilizing due to cross-strand stacking and specific hydrogen bonding networks, making their accurate parameterization essential for predicting biologically relevant structures like internal loops.
Ion and Temperature Corrections
The standard Turner parameters are defined at 1 M NaCl and 37°C. To predict folding under physiological or non-standard conditions, correction terms are applied:
- Salt correction: A term derived from polyelectrolyte theory adjusts the entropic penalty of folding based on monovalent and divalent cation concentration, reflecting the reduced charge repulsion as counterions condense on the backbone.
- Temperature extrapolation: Standard enthalpy and entropy values allow free energy calculation at any temperature using the Gibbs-Helmholtz equation, enabling prediction of melting temperatures (Tm).
Coaxial Stacking in Multi-Loops
When two helices meet within a multi-branch loop, their terminal base pairs can stack coaxially to form a pseudo-continuous helix. This interaction is not explicitly parameterized in the original linear multi-loop model but is a critical determinant of 3D architecture. Modern implementations incorporate coaxial stacking free energy increments—often approximated by the nearest-neighbor parameter for the stacked interface—to favor junction topologies observed in phylogenetic structures and crystallography.
Frequently Asked Questions
Clarifying the empirical rules and computational logic behind the nearest-neighbor model that powers RNA secondary structure prediction.
The Turner Energy Model is the standard empirical nearest-neighbor model that assigns experimentally derived thermodynamic parameters to RNA base pair stacks and loops, forming the basis for minimum free energy (MFE) and partition function calculations. Rather than treating individual base pairs in isolation, the model operates on the principle that the stability of a base pair depends on the identity and orientation of its immediate neighbors. The total free energy change of an RNA secondary structure is calculated as the sum of incremental free energy terms for each structural motif: stacking interactions between adjacent base pairs, hairpin loops, internal loops, bulge loops, and multi-branch loops. These parameters, measured through optical melting experiments on model oligonucleotides, are compiled in the Turner Rules database and are the default energy function in software packages like RNAstructure, ViennaRNA, and mfold.
Enabling Efficiency, Speed & Accuracy
Intelligent Analysis, Decision & Execution
We build AI systems for teams that need search across company data, workflow automation across tools, or AI features inside products and internal software.
Talk to Us
Search across company data
Give teams answers from docs, tickets, runbooks, and product data with sources and permissions.
Useful when people spend too long searching or get different answers from different systems.

Automate internal workflows
Use AI to route work, draft outputs, trigger actions, and keep approvals and logs in place.
Useful when repetitive work moves across multiple tools and teams.

Add AI to products and internal tools
Build assistants, guided actions, or decision support into the software your team or customers already use.
Useful when AI needs to be part of the product, not a separate tool.
Related Terms
The Turner Energy Model is the empirical backbone of RNA secondary structure prediction. These interconnected concepts define how thermodynamic parameters are calculated, applied, and extended in modern algorithms.
Minimum Free Energy (MFE)
The thermodynamic principle that predicts the single most stable RNA secondary structure by minimizing the sum of empirically derived loop and stacking energy parameters. The Turner model provides the nearest-neighbor parameters that make MFE calculations possible.
- Uses dynamic programming (Zuker algorithm) to efficiently search the space of possible structures
- The structure with the lowest ΔG° is considered the biologically relevant fold
- MFE prediction accuracy depends directly on the quality of Turner parameters
Partition Function
A statistical mechanics calculation that sums the Boltzmann-weighted free energies of all possible RNA secondary structures to derive base pairing probabilities and ensemble properties. The Turner energy parameters serve as the Hamiltonian for this calculation.
- Computes the equilibrium ensemble rather than a single structure
- Yields per-base pairing probabilities used in centroid structure prediction
- Enables calculation of ensemble diversity and Shannon entropy of the structural landscape
Nearest-Neighbor Model
The fundamental assumption that the stability of an RNA base pair depends on the identity and orientation of adjacent base pairs, not just the pair itself. The Turner model quantifies this with experimentally measured parameters.
- Each stacking interaction between adjacent base pairs has a distinct ΔG° value
- Parameters are derived from optical melting experiments on short model duplexes
- Accounts for sequence-dependent effects like the 5'-GC-3' vs 5'-CG-3' asymmetry
Dot-Bracket Notation
A standard string representation of RNA secondary structure where matching parentheses denote canonical base pairs and dots represent unpaired nucleotides. This format is the direct output of Turner-model-based folding algorithms.
- Example:
(((...)))represents a simple stem-loop - Serves as the training target for deep learning models predicting secondary structure
- Extended notation uses brackets
[]and braces{}to annotate pseudoknots
Loop Free Energy Parameters
The Turner model assigns distinct thermodynamic penalties to different loop types based on size, sequence, and closure geometry. These parameters are critical for accurately predicting multi-branch junctions.
- Hairpin loops: Sequence-dependent penalties for tetraloops like GNRA and UNCG
- Internal loops: Asymmetric and symmetric loops with size-dependent initiation costs
- Multi-branch loops: Linear approximation using junction initiation penalty and per-helix bonus
- Bulge loops: Single-nucleotide bulges have specific stacking context parameters
SHAPE-Directed Folding
An extension of Turner-model folding that incorporates experimental reactivity data as a pseudo-free energy term to constrain structure prediction. SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) measures nucleotide flexibility.
- Reactivity values are converted to ΔG° penalties for base-paired nucleotides that appear flexible
- Dramatically improves prediction accuracy, especially for long RNAs
- Bridges the gap between purely thermodynamic prediction and experimental observation

About the author
Prasad Kumkar
CEO & MD, Inference Systems
Prasad Kumkar is the CEO & MD of Inference Systems and writes about AI systems architecture, LLM infrastructure, model serving, evaluation, and production deployment. Over 5+ years, he has worked across computer vision models, L5 autonomous vehicle systems, and LLM research, with a focus on taking complex AI ideas into real-world engineering systems.
His work and writing cover AI systems, large language models, AI agents, multimodal systems, autonomous systems, inference optimization, RAG, evaluation, and production AI engineering.
Partnered with leading AI, data, and software stack.
How We Work
Custom AI workflows for your Business
One-fit-all AI don't work for modern businesses. At Inferensys, we aim to understand your business & custom requirements; which we use to define most efficient agentic workflows, the data, and the tools for your business.
01
Review the use case
We understand the task, the users, and where AI can actually help.
Read more02
Pick the right approach
We define what needs search, automation, or product integration.
Read more03
Build the first useful version
We implement the part that proves the value first.
Read more04
Improve from there
We add the checks and visibility needed to keep it useful.
Read moreThe first call is a practical review of your use case and the right next step.
Talk to Us