RNA optical melting

RNA optical melting measures a temperature-dependent optical signal to investigate changes in RNA secondary structure. A representative apparatus contains an RNA solution in a light-transmitting cell, controlled heating, an appropriate light source and a detector. A spectrophotometer compares transmitted intensity $I$ with reference intensity $I_0$ to obtain absorbance, $A=-{\log }_{10}(I/I_0)$. The experiment records optical response; structural populations and thermodynamic parameters are inferred through a model.

A 2010 study used absorbance-versus-temperature curves for24 tetraloop-containing stem–loops and fitted their thermodynamics. Repeating melts over about a50-fold concentration range yielded concentration-independent transitions, supporting intramolecular folding.^[ref] The published study supplies an experimental use of the method, not universal conditions for every RNA.

For an original optical calculation, $I/I_0=0.50$ corresponds to $A\approx 0.3010$, while $I/I_0=0.25$ corresponds to $A\approx 0.6021$. Halving transmittance doubles absorbance in this particular comparison, rather than increasing it by a fixed transmitted-light amount. Absorbance is a logarithmic measurement.

Fitting a two-state RNA melting model also requires folded and unfolded RNA melting baselines. If $f_U(T)$ is the unfolded fraction, an ideal mixture givesA changing baseline can alter the apparent midpoint or slope even when the same population curve is used. Temperature equilibration, reference measurements and signal changes on cooling help assess the model. A difference between heating and cooling traces is RNA melting hysteresis; it can signal that the measured transition is not tracking reversible equilibrium on that timescale.

Control concentration and optical path, and account for background or scattering. Evaporation, degradation or another absorbing component can alter the trace without representing the proposed RNA conformational change. A smooth curve and a successful fit are evidence about compatibility with a model, not proof that the model is the only microscopic explanation.

RNA melting hysteresis is dependence of an RNA optical melting trace on thermal history: heating and cooling can give different absorbance at the same temperature. It warns that the observed trace may not follow reversible chemical equilibrium on the measurement timescale. Two different traces do not by themselves establish two equilibrium melting temperatures.

A primary experiment illustrates the importance of history. A 21-nucleotide RNA showed a low-temperature dimer-to-hairpin transition that did not reverse on cooling, whereas its higher-temperature RNA hairpin melting transition was reversible. Measurements used a Cary 100 spectrophotometer, a 1-cm quartz cell and 268-nm light, heating at $1.0{}^{\circ }\mathrm{C}/\min$ and cooling at $0.5$ or $0.2{}^{\circ }\mathrm{C}/\min$. An adjacent solvent cell held the temperature probe.^[ref] The authors interpreted the dimer as kinetically trapped under those conditions, rather than applying a reversible equilibrium fit to that transition.

For an original simple kinetic model, let $f$ be the unfolded fraction, $f_{\mathrm{eq}}(T)$ its equilibrium value and $\tau$ an RNA melting relaxation time:After a fixed-temperature step with constant $\tau$, the deviation undergoes exponential decay:With stipulated $f(0)=0$, $f_{\mathrm{eq}}=1$ and $\tau =2\min$, the fraction reaches only 0.632 after two minutes and 0.950 after six. During a continuing temperature ramp, this lag can put heating and cooling on different branches. This single-time model is illustrative; it is not a fit to the published dimer experiment.

Useful checks compare ramp speeds and include fixed-temperature holds and repeat cycles. Verify measurement calibration and exclude evaporation or chemical degradation before attributing a discrepancy to folding kinetics. Slower scans can reduce kinetic lag, but a smooth trace alone cannot establish reversibility.

An RNA melting relaxation time characterizes how quickly a perturbed RNA population approaches its new equilibrium. It belongs to a specified kinetic model, temperature, solution and measured signal. It is not an RNA melting temperature or, in general, the mean waiting time for an individual unfolding event.

For a reversible two-state RNA melting model, let $k_u$ be the first-order unfolding rate constant and $k_f$ the refolding rate constant. The unfolded fraction obeysAt fixed temperature,Both directions contribute to the return to equilibrium. For original stipulated values $k_u=20{\mathrm{s}}^{-1}$ and $k_f=80{\mathrm{s}}^{-1}$, $f_{U,\mathrm{eq}}=0.20$ and $\tau =10\mathrm{ms}$. The mean lifetime of an isolated folded state in this model is $1/k_u=50\mathrm{ms}$, which differs from the observed population exponential decay time.

A temperature jump shifts populations rapidly enough to reveal the subsequent motion. A primary study used 4-ns heating pulses, a 12-μm sample path and time-resolved infrared spectroscopy to follow melting and refolding of a UACG-containing RNA hairpin. Thin-cell cooling enabled direct observation of the return process. Its loop is UACG, not a UUCG tetraloop.^[ref]

An RNA relaxation instrument-response correction may be needed if heating, detector response or cooling occur on comparable timescales. A changing sample temperature changes the equilibrium target during the measurement, invalidating a constant-temperature fit across the whole trace. Multiple kinetic states can also produce several decay components. One fitted exponential should therefore be identified as a model-dependent observed relaxation, rather than automatically assigned to a unique molecular step.

An RNA relaxation instrument-response correction accounts for the way an experiment's heating, detection and recording alter an underlying molecular transient. A measured RNA melting relaxation time can otherwise combine molecular kinetics with instrumental behavior. The response must be characterized for the actual apparatus and measurement conditions.

For an illustrative linear, time-invariant detector, let the normalized instrument response function be $h(t)={\tau }_I^{-1}{e}^{-t/{\tau }_I}$ for $t\ge 0$. Suppose the ideal molecular signal is $s(t)=A{e}^{-t/{\tau }_R}$ after excitation, and both vanish before it. Their convolution gives, for ${\tau }_R\ne {\tau }_I$,The recorded trace rises before decaying even though the assumed molecular transient begins at its maximum. Fitting the entire trace as one unblurred exponential would misdescribe the model.

An original example sets ${\tau }_R=10\mu \mathrm{s}$ and ${\tau }_I=5\mu \mathrm{s}$. The measured signal peaks at $10\ln 2\approx 6.93\mu \mathrm{s}$ with height $A/2$. Its normalized temporal centroid is ${\tau }_R+{\tau }_I=15\mu \mathrm{s}$, rather than 10. These different time summaries cannot be substituted for one another.

One useful approach predicts a molecular transient, convolves it with the characterized instrument response function, and fits that forward prediction to the data. Directly undoing a smoothing response can amplify measurement noise. Parameters that produce nearly indistinguishable recorded curves may remain difficult to determine separately.

Sample cooling is an additional physical process. A primary temperature jump study used a thin sample cell to accelerate cooling enough to observe RNA refolding, with separate infrared probe pulses following the transient.^[ref] If temperature changes during a trace, the molecular equilibrium and rate constants may change too; detector convolution alone does not account for that changing input. The two-exponential example is a stipulated teaching model, not a fit to that paper.

An RNA melting baseline is the predicted absorbance of a specified folded or unfolded RNA population as a function of temperature, before mixing those populations in an RNA optical melting model. It differs from a solvent blank: a folded RNA sample still absorbs light, and that response can change with temperature even without the modeled transition.

For a two-state RNA melting model with unfolded fraction $f_U$, writeThe folded and unfolded RNA melting baselines might be modeled as $A_F=a_F+b_{F}T$ and $A_U=a_U+b_{U}T$. Their continuation into the transition region is RNA baseline extrapolation. Solving the optical mixture givesFor original stipulated readings $A=0.55$, $A_F=0.40$ and $A_U=0.70$, the inferred fraction is 0.50. Changing only the folded RNA melting baseline to 0.45 gives 0.40; changing only the unfolded RNA melting baseline to 0.80 gives 0.375. The same measured absorbance can therefore imply different populations under different baseline models.

Sensitivity is explicit. With $D=A_U-A_F$,A small baseline separation amplifies measurement uncertainty. Correlated fitted baselines require covariance terms in uncertainty propagation, rather than independent error bars by assumption.^[ref]

Optical-melting studies infer thermodynamics by fitting measured curves, rather than directly counting unfolded molecules.^[ref] Inspect the pre-transition and post-transition regions and test the chosen RNA baseline extrapolation. An inferred fraction outside $[0,1]$ signals noise or incompatibility with the assumed optical model; it is not a negative molecular population.

RNA baseline extrapolation continues a fitted folded-state or unfolded-state RNA melting baseline from its observed temperature range into the transition region. It is necessary when estimating populations from RNA optical melting, because both hypothetical pure-state signals are needed at temperatures where the measured sample contains a mixture.

For a linear folded baseline, choose a reference $T_0$ and writeThe reference merely changes parameterization. An original fit with $T_0=300\mathrm{K}$, $a=0.400$ and $b=0.001{\mathrm{K}}^{-1}$ predicts 0.420 at $320\mathrm{K}$ and 0.450 at $350\mathrm{K}$. The first prediction might interpolate within measured data; the second extrapolates if all folded-state observations were below $320\mathrm{K}$.

The inferred population depends on that continuation. At $350\mathrm{K}$, suppose measured absorbance is 0.600 and the unfolded RNA melting baseline is 0.700. The linear folded-baseline prediction givesIf the actual pure folded-state response there were 0.490, the same optical-mixture equation would instead give $f_U\approx 0.524$. A good fit in a limited low-temperature region does not rule out curvature beyond it.

RNA baseline prediction covariance is important even within the chosen linear model. With $x=T-T_0$,For stipulated independent fit parameters $u_a=0.002$ and $u_b=0.0002{\mathrm{K}}^{-1}$, the prediction's standard uncertainty rises from 0.002 at $x=0$ to approximately 0.0102 at $x=50\mathrm{K}$. This is uncertainty propagation for estimated parameters; it does not include unmodeled curvature.^[ref]

Optical-melting experiments fit signal curves to infer thermodynamics.^[ref] Report which regions constrain each RNA melting baseline, inspect residuals and compare plausible baseline forms. Confidence in the fitted parameters and confidence in the functional form are different parts of the inference.

RNA baseline prediction covariance quantifies how errors in fitted RNA melting baselines link their predicted values at different temperatures. Reusing the same fit parameters creates shared uncertainty; predictions are not independent merely because they are evaluated at different temperatures.

For a linear RNA baseline extrapolation $A_F(T)=a+b(T-T_0)$, put $x_i=T_i-T_0$. ThenThis relation follows from the covariance of two linear combinations. The diagonal $C_{ii}$ gives a prediction variance; the off-diagonal entries determine how predictions move together.^[ref]

An original example assumes independent parameters with $u_a=0.002$ and $u_b=0.0002{\mathrm{K}}^{-1}$. At temperatures $20\mathrm{K}$ below and above $T_0$, both prediction variances are $20\times 10^{-6}$, while their covariance is $-12\times 10^{-6}$. The correlation coefficient is therefore $-0.6$. A slope error moves the predictions in opposite directions around the reference point, outweighing the common intercept error.

The difference between those two predictions cancels the intercept exactly:Its standard uncertainty is 0.008. Ignoring covariance would give $\sqrt{40\times 10^{-6}}\approx 0.00632$, which is too small for this model. Shared fitted inputs can produce negative as well as positive prediction correlations.

In vector form, a prediction covariance matrix is $C=X{\Sigma }_{\theta }{X}^{\mathsf{T}}$, with rows of $X$ equal to $(1,x_i)$ and parameter covariance ${\Sigma }_{\theta }$. A two-parameter fit supplies at most two independent shared-error directions, even if evaluated at a hundred temperatures. Additional independent readout errors contribute their own covariance.

When these predictions are used to infer populations in RNA optical melting, propagate their joint uncertainty rather than only separate error bars. This accounting does not test whether a linear baseline is chemically adequate; curvature or an incorrect pure-state model needs separate evidence.

A two-state RNA melting model treats an RNA transition as an equilibrium between folded F and unfolded U macrostates, with no substantially populated intermediate. Each macrostate can contain many microscopic conformations; two states does not mean exactly two atom-by-atom structures. The model relates populations to thermodynamics and then, separately, to an observed optical signal.

For unfolding free energy $\Delta G_u=G_U-G_F$,The population formula assumes an ideal unimolecular transition. With constant unfolding enthalpy and entropy,An RNA melting log-odds plot against $1/T$ is then a straight line. Nonzero changes in heat capacity, alternative states or incorrect population estimates can invalidate this simplified relation.

An original consistency test fits stipulated $f_U=0.10$ at $300\mathrm{K}$ and $f_U=0.50$ at $330\mathrm{K}$. These two points determineThe fitted model predicts approximately $f_U=0.862$ at $360\mathrm{K}$. A third observation tests that prediction; fitting two parameters to two points was not independent validation. Measurement uncertainty and baselines must be considered before attributing disagreement to a specific physical mechanism.

For an ideal optical mixture, $A(T)=[1-f_U(T)]A_F(T)+f_U(T)A_U(T)$. Thus the population model does not remove the need for folded and unfolded RNA melting baselines. Optical-melting studies infer thermodynamics by fitting measured curves under stated analysis assumptions.^[ref]

Reversibility, the absence of important intermediates and the unimolecular assumption are physical conditions to test. A two-strand association transition has a different concentration dependence even if its optical trace looks similarly sigmoidal. The model is a useful quantitative hypothesis, not a consequence of drawing one smooth curve through data.

An RNA melting log-odds plot graphs the natural logarithm of the unfolded-to-folded population ratio against reciprocal absolute temperature. For $0<f_U<1$, let $y=\ln [f_U/(1-f_U)]$ and $x=1/T$. An ideal unimolecular two-state RNA melting model with temperature-independent unfolding enthalpy and entropy predictsThe slope is $-\Delta H_u/R$; the intercept is $\Delta S_u/R$. The equation is a testable model relation, not a reason to treat every transformed optical curve as linear.

For an original worked example, stipulate unfolded fractions 0.20 at $300\mathrm{K}$ and 0.50 at $320\mathrm{K}$. Their log-odds are $-\ln 4$ and 0, givingThe fitted enthalpy is approximately $55.3\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$, and the model predicts $f_U\approx 0.773$ at $340\mathrm{K}$. A third measurement tests this prediction; two chosen points always determine a line.

The transformation also changes the noise. First-order uncertainty propagation gives an RNA melting log-odds uncertaintyFor the same small standard uncertainty $u_f=0.001$, this gives $u_y\approx 0.004$ at $f_U=0.5$, but approximately 0.101 at $f_U=0.01$. Thus equal uncertainty in population fractions produces heteroscedasticity in the transformed data; unweighted linear regression can overemphasize unreliable endpoint values. Exactly 0 or 1 cannot be transformed, and replacing them with arbitrary small offsets manufactures finite values.^[ref]

Fractions inferred from absorbance also inherit the fitted RNA melting baselines. Optical-melting thermodynamics depends on this modeling step.^[ref] Curvature may motivate testing baseline errors, additional states or changes in heat capacity; the plot alone does not identify which explanation is correct.

RNA melting log-odds uncertainty quantifies uncertainty in $y=\ln [f_U/(1-f_U)]$, the transformed population used in an RNA melting log-odds plot. It depends on the unfolded fraction $f_U$ as well as its uncertainty. The transformation becomes particularly sensitive near 0 and 1.

For a sufficiently concentrated uncertainty distribution within $0<f_U<1$, first-order uncertainty propagation givesWith an original stipulated $u_f=0.010$, the standard uncertainty is 0.040 at $f_U=0.50$ but approximately 0.510 at $f_U=0.02$. Equal population uncertainties do not justify equal statistical treatment of transformed observations.^[ref]

Nonlinearity also creates asymmetry. At $f_U=0.02$, $y\approx -3.892$. Changing the fraction to 0.01 gives $y\approx -4.595$, while changing it to 0.03 gives $y\approx -3.476$: equal changes in fraction produce unequal transformed changes. These are deterministic perturbations, not asserted confidence limits. If the uncertainty materially reaches 0 or 1, the linear approximation needs reconsideration; taking the logarithm of those endpoints is undefined.

The underlying RNA melting baselines can induce shared errors across temperatures. For an optical mixture with $A_F<A<A_U$,soThese derivatives propagate the joint uncertainty of the measured absorbance and fitted baselines. Reusing baseline parameters for many observations creates RNA log-odds covariance; merely assigning separate error bars does not remove it.

A fit can model the original absorbance observations and their RNA melting baselines, or handle the transformed data with an appropriate joint error model. Optical curves are fitted quantities in primary melting analyses, rather than direct counts of unfolded molecules.^[ref] Precision in a fitted straight line cannot establish the correctness of the two-state RNA melting model used to interpret it.

RNA log-odds covariance describes correlated errors between transformed population estimates at different temperatures in an RNA melting log-odds plot. Shared RNA melting baselines can generate it even when individual detector readings are independent.

Let $f_i$ be unfolded fractions and $y_i=\ln [f_i/(1-f_i)]$. First-order propagation givesIn matrix form, $C_y\approx J{C}_f{J}^{\mathsf{T}}$ with diagonal $J_{ii}=1/[f_i(1-f_i)]$. This is the multivariate form of RNA melting log-odds uncertainty; the derivatives are evaluated at the estimated fractions.^[ref]

A concrete shared-input model uses $y_i=\ln (A_i-A_{F,i})-\ln (A_{U,i}-A_i)$. If all folded-state predictions contain the same uncertain offset $a$, thenThis component has one shared error direction. It is not equivalent to adding independent noise with variance $g_i^2{u}_a^2$ to every observation. Other baseline parameters contribute through RNA baseline prediction covariance.

For an original simplified example, nine observations each have independent log-odds standard uncertainty0.03 and the same shared-offset contribution 0.01. The standard uncertainty of their mean isIncorrectly treating the shared term as independent would give $\sqrt{(0.03^2+0.01^2)/9}\approx 0.0105$. Repetition reduces the independent component while leaving this common offset contribution intact.

A model fit can use the joint covariance matrix in generalized least squares, with an objective of the form $r^{\mathsf{T}}{C}_y^{-1}r$ when that inverse exists. Alternatively, fit the underlying absorbance and baseline parameters jointly. Shared-error components alone may be singular, and near population endpoints the logarithmic linearization may fail. A precise straight-line fit does not independently validate the two-state RNA melting model or its baseline assumptions.