Concentration dependence of RNA melting 2026-11-05
The concentration dependence of RNA melting is the change of an RNA melting temperature or population curve with total strand molar concentration, while other solution conditions are held fixed. It helps distinguish intramolecular RNA hairpin melting from association of separate strands, provided the system reaches chemical equilibrium.
For an ideal unimolecular folded-to-unfolded transition, , the equilibrium constant is . Multiplying both populations by the same amount leaves their ratio unchanged: the two-state RNA melting model has no total-concentration term. By contrast, consider self-complementary RNA duplex melting, with association . Let count strands and let denote the reference molar concentration of the standard state. The dimensionless association equilibrium constant isAt the midpoint defined by half the strands belonging to duplexes, and , giving . This factor follows from strand counting; a different reaction or midpoint convention requires its own derivation.
An original ideal calculation assumes constant association enthalpy and a midpoint at . Increasing strand molar concentration tenfold givesor approximately . Here is the molar gas constant. Greater strand availability favors association, so more heating is needed to reach the same strand-based midpoint.
Experimentally, one records RNA optical melting curves at several strand molar concentrations, with matched solvent conditions. A tetraloop study found concentration-independent transitions over approximately a 50-fold range, supporting intramolecular stem–loop formation.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, Optical melting experiments and thermodynamics.
This test is conditional on reversibility: a kinetically trapped association can fail to show the expected equilibrium concentration response. Justin P. Sheehy 2026-11-05
Justin P. Sheehy is first author of a2010 study measuring the thermodynamic stability of common RNA tetraloops. These four-nucleotide loops cap RNA hairpins, and their contribution to stability depends on sequence and structural context.[ref]
The study used RNA optical melting: measuring absorbance as temperature changed, then fitting the transitions to obtain thermodynamic parameters. Its Results distinguish33 stem–loop sequences from24 unique tetraloops, because some loops were tested in multiple stems. Repeated melting curves over an approximately50-fold concentration range showed single transitions and concentration independence, supporting an intramolecular hairpin interpretation.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Thermodynamic parameters,Table2; Optical melting experiments and thermodynamics
An original two-state RNA melting model makes the measured quantities interpretable. If folding has temperature-independent enthalpy change and entropy change , then its folding Gibbs free energy isThe RNA melting temperature of this stipulated unimolecular model satisfies , givingabout60.2 degrees Celsius. At310 kelvin the model instead gives , favoring folding. These invented parameters illustrate the calculation; they are not reported values for one of Sheehy's sequences.
The experiment links a temperature-dependent optical signal with a model of RNA stability. An absorbance transition alone is not an atomic structure, and changing the stem can change the inferred loop contribution. The author identity here follows the primary scientific byline; it does not assume a career or portrait from a same-name directory entry.
Noncanonical base pair 2026-11-05
A noncanonical base pair joins two nucleobases in an arrangement other than the standard Watson–Crick base pairing geometry. For RNA, the canonical Watson–Crick combinations are adenine–uracil and guanine–cytosine. Other combinations and contact geometries can contribute to an organized, functional RNA secondary structure; noncanonical does not mean an experimental mistake.
RNA base-pair geometry describes more than the two base letters. Which parts of the nucleobases contact each other, their relative orientation and their hydrogen bonds distinguish arrangements that a letter-only sequence description cannot specify. For a concrete structural example, solution NMR spectroscopy studies of GCAA and GAAA GNRA tetraloops identified a guanine–adenine base pair between the first and last loop residues, alongside other stabilizing contacts.[ref]
Structural features that give rise to the unusual stability of RNA hairpins containing GNRA loops by H. A. Heus and others. Abstract, loop G–A pairing and solution structures.
An original counting example makes the distinction useful. A three-pair stem capped by a GAAA RNA tetraloop contains ten nucleotides: six in the stem and four in the loop. A contact between the loop's first guanine and last adenine can add a noncanonical base pair without adding any nucleotides or turning the four-residue loop into a different-length sequence. Sequence length, stem-pair count and the complete three-dimensional contact inventory answer different questions.
Conventions should be stated. Some RNA analyses include the common guanine–uracil wobble arrangement among their canonical pairs even though it differs from Watson–Crick geometry. The tetraloop survey, for example, explicitly adopts that broader convention.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, Compiling and searching a database for RNA tetraloops.
Comparing reported counts therefore requires checking what each author counted as canonical. RNA baseline extrapolation 2026-11-05
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 and writeThe reference merely changes parameterization. An original fit with , and predicts 0.420 at and 0.450 at . The first prediction might interpolate within measured data; the second extrapolates if all folded-state observations were below .
The inferred population depends on that continuation. At , 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 . 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 ,For stipulated independent fit parameters and , the prediction's standard uncertainty rises from 0.002 at to approximately 0.0102 at . This is uncertainty propagation for estimated parameters; it does not include unmodeled curvature.[ref]
Evaluation of measurement data — Guide to the expression of uncertainty in measurement by JCGM. §§5.1–5.2, linear sensitivities and correlated inputs.
Optical-melting experiments fit signal curves to infer thermodynamics.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, individual-curve thermodynamic fitting.
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 hairpin 2026-11-05
An RNA hairpin is a stem-loop formed when two regions of one RNA strand make base pairs, leaving an intervening unpaired loop. It is a local RNA secondary structure motif: covalent continuity runs through the loop while pairing brings separated parts of the sequence together. The motif does not require a second RNA molecule and need not describe the fold of the whole RNA.
For an original sequence illustration, consider
5'-GGCCUUCGGGCC-3'. Pair positions1–12,2–11,3–10 and4–9. The four pairs are G–C, G–C, C–G and C–G, leaving positions5–8, UUCG, as a four-nucleotide loop. Thus four pairs use eight nucleotides, and the loop supplies four more, giving twelve in total. The corresponding illustrative dot-bracket notation is ((((....)))). This proposed pairing is a structural example, not an experimental demonstration that an isolated sequence occupies only that conformation.A four-residue loop is an RNA tetraloop. Both the paired stem and the loop sequence influence stability. In a 2010 experiment, researchers optically melted24 stem–loop sequences containing common tetraloops and determined their thermodynamic parameters. Measuring absorbance as temperature changes connects sequence differences to an experimental folding transition rather than inferring stability from the number of pairs alone.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Abstract and Materials and Methods: optical melting and thermodynamic analysis of24 sequences.
A two-state model makes RNA hairpin melting quantitative. If only folded F and unfolded U states are allowed, set , with and the molar gas constant. The folded fraction is . A stipulated gives and folded fraction ; equal free energies give . These model fractions are not the paper's measurements. Alternative folds, non-two-state transitions and assay conditions can require a richer interpretation.
RNA hairpin melting 2026-11-05
RNA hairpin melting is the temperature-dependent loss of the folded stem-loop structure of an RNA hairpin. Here melting describes a molecular conformational transition in solution, rather than conversion of a solid sample into a liquid. Changing temperature shifts the relative populations of folded and unfolded RNA.
In optical measurements, a spectrophotometer tracks absorbance while the RNA solution is heated under controlled conditions. Changes in base stacking can alter ultraviolet absorbance, an effect associated with RNA hyperchromicity. In a 2010 tetraloop study, researchers recorded absorbance-versus-temperature curves, repeated measurements over an approximately 50-fold RNA concentration range, and observed concentration-independent transitions. That behavior supported intramolecular hairpin formation rather than association of separate strands.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, Optical melting experiments and thermodynamics; repeated curves and concentration dependence.
An explicit two-state model lets F and U denote folded and unfolded RNA. Define unfolding free energy and take temperature-independent unfolding parameters and :Here is absolute temperature and the molar gas constant. The RNA melting temperature satisfies , hence and in this model.
For original stipulated values and , the predicted unfolded fractions at , and are approximately , and . They are model calculations, not measurements from the cited study.
Raw absorbance is not automatically the unfolded fraction. Even in a two-state interpretation, folded and unfolded baselines can vary with temperature: . Baseline estimates, reversibility and alternative conformations affect the inference. A smooth single transition is useful evidence, but it does not independently prove that every RNA molecule follows one microscopic pathway.
RNA hyperchromicity 2026-11-05
RNA hyperchromicity is an increase in an RNA sample's ultraviolet absorption associated with a change in its molecular environment or conformation. The general hyperchromic effect concerns increased spectral-band intensity, rather than a shift of the band to another wavelength.[ref] In RNA melting measurements, altered base stacking can contribute to the change in absorption; signal interpretation still requires a defined wavelength and sample conditions.
The Beer-Lambert law writes , where is absorbance, a molar absorption coefficient, path length and concentration. A change in can therefore change absorbance even when the amount of RNA is unchanged. Conversely, changed concentration, path length or other optical contributions can mimic a structural signal.
For an original comparison at one stated temperature, suppose extrapolated folded and unfolded absorbances are and . Defining fractional hyperchromicity relative to the folded state gives , or50%. Using the unfolded state as denominator would instead give . A reported percentage therefore needs its convention.
In a linear two-population optical model,provided . With those baselines, gives . But with temperature-appropriate baselines and , the same absorbance gives . An RNA melting baseline is consequently part of the population estimate, not decorative plotting.
A 2010 tetraloop study measured absorbance-versus-temperature curves to infer folding thermodynamics.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, optical melting and curve fitting.
Such measurements connect an optical response to a model of RNA structure. An absorption increase alone does not count broken base pairs, establish a unique unfolded structure or distinguish every possible chemical change. RNA melting baseline 2026-11-05
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 , writeThe folded and unfolded RNA melting baselines might be modeled as and . Their continuation into the transition region is RNA baseline extrapolation. Solving the optical mixture givesFor original stipulated readings , and , 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 ,A small baseline separation amplifies measurement uncertainty. Correlated fitted baselines require covariance terms in uncertainty propagation, rather than independent error bars by assumption.[ref]
Evaluation of measurement data — Guide to the expression of uncertainty in measurement by JCGM. §§5.1–5.2, sensitivities and correlated inputs.
Optical-melting studies infer thermodynamics by fitting measured curves, rather than directly counting unfolded molecules.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, individual-curve fitting.
Inspect the pre-transition and post-transition regions and test the chosen RNA baseline extrapolation. An inferred fraction outside signals noise or incompatibility with the assumed optical model; it is not a negative molecular population. RNA melting log-odds plot 2026-11-05
An RNA melting log-odds plot graphs the natural logarithm of the unfolded-to-folded population ratio against reciprocal absolute temperature. For , let and . An ideal unimolecular two-state RNA melting model with temperature-independent unfolding enthalpy and entropy predictsThe slope is ; the intercept is . 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 and 0.50 at . Their log-odds are and 0, givingThe fitted enthalpy is approximately , and the model predicts at . 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 , this gives at , but approximately 0.101 at . 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]
Evaluation of measurement data — Guide to the expression of uncertainty in measurement by JCGM. §5, first-order sensitivity-coefficient method; derivatives here are for the stated transformation.
Fractions inferred from absorbance also inherit the fitted RNA melting baselines. Optical-melting thermodynamics depends on this modeling step.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, optical curves and thermodynamic fitting.
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 2026-11-05
RNA melting log-odds uncertainty quantifies uncertainty in , the transformed population used in an RNA melting log-odds plot. It depends on the unfolded fraction as well as its uncertainty. The transformation becomes particularly sensitive near 0 and 1.
For a sufficiently concentrated uncertainty distribution within , first-order uncertainty propagation givesWith an original stipulated , the standard uncertainty is 0.040 at but approximately 0.510 at . Equal population uncertainties do not justify equal statistical treatment of transformed observations.[ref]
Evaluation of measurement data — Guide to the expression of uncertainty in measurement by JCGM. §5, sensitivity-coefficient method; derivative here follows from the stated logarithmic transformation.
Nonlinearity also creates asymmetry. At , . Changing the fraction to 0.01 gives , while changing it to 0.03 gives : 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 ,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]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, optical curves and thermodynamic fitting.
Precision in a fitted straight line cannot establish the correctness of the two-state RNA melting model used to interpret it. RNA melting temperature 2026-11-05
An RNA melting temperature, , characterizes the midpoint of a specified RNA conformational transition under stated conditions. In a reversible, unimolecular two-state RNA melting model, it is the temperature at which folded and unfolded populations are equal. A midpoint assigned to an optical trace requires a model connecting absorbance to those populations.
Let unfolding free energy be . If unfolding enthalpy and entropy are treated as temperature independent,The formula assumes a finite relevant transition and uses absolute temperature. With original stipulated values and , . Holding enthalpy fixed but changing entropy to gives approximately . A melting temperature reflects both energetic and entropic contributions.
Molecularity matters. An ideal intramolecular hairpin transition need not change its midpoint with total strand concentration. A transition involving association of separate strands generally does: the equilibrium relation includes their concentrations. The concentration dependence of RNA melting can therefore help test which structural model is appropriate. In the 2010 tetraloop measurements, concentration-independent curves across a broad concentration range supported stem–loop formation.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, repeated melts and concentration dependence.
A raw halfway absorbance is not automatically this thermodynamic midpoint. If the folded and unfolded RNA melting baselines at one temperature are and , an absorbance of implies unfolded fraction , not . Intermediates, temperature-dependent parameters or RNA melting hysteresis can further complicate midpoint assignment. Report the solution conditions and definition of the measured transition rather than treating as a context-free constant of an RNA sequence.
RNA optical melting 2026-11-05
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 with reference intensity to obtain absorbance, . 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]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Abstract and Materials and Methods, optical-melting experiments.
The published study supplies an experimental use of the method, not universal conditions for every RNA.For an original optical calculation, corresponds to , while corresponds to . 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 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 tetraloop 2026-11-05
An RNA tetraloop is the four-nucleotide loop capping the paired stem of an RNA hairpin. The four count excludes the two nucleotides of the closing base pair. A tetraloop is therefore a local RNA secondary structure motif, not necessarily an entire four-nucleotide RNA molecule.
For an original sequence illustration,
5'-GGCGAAAGCC-3' can be drawn with pairs 1–10, 2–9 and 3–8. The central GAAA at positions 4–7 forms a four-residue loop; three stem pairs plus four loop residues give nucleotides. Whether a particular sequence actually favors that conformation requires experimental or structural evidence.Comparative analysis of ribosomal RNA identified recurrent motif families, including GNRA tetraloops and UNCG tetraloops, as well as the less common CUUG motif. The study also found relationships between loop sequence and its closing pair.[ref]
Architecture of ribosomal RNA: constraints on the sequence of "tetra-loops" by Carl R. Woese and others. Abstract, recurring motifs and closing-base-pair preferences.
In GNRA, N means any of the four ordinary RNA bases and R means a purine, A or G; GAAA is one example. UUCG belongs to UNCG, with its N position occupied by U.A simple sequence count separates a motif pattern from a frequency claim. Four unrestricted positions give possible sequences. GNRA specifies the first and last bases and permits sequences; UNCG permits 4. These counts assume the usual A, C, G, U alphabet. They do not say that the sequences occur equally often or have identical stability.
The stem and loop environment matter together. Optical-melting measurements on 24 selected stem–loops quantified sequence-dependent thermodynamics rather than treating every four-residue loop as one energetic object.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Abstract and sequence-design/thermodynamic-analysis sections.
A frequently observed motif and a particularly stable isolated hairpin are different scientific observations; both can inform how RNA structures are understood. Self-complementary RNA duplex melting 2026-11-05
Self-complementary RNA duplex melting dissociates a duplex formed by two copies of the same RNA strand into separate strands. Self-complementarity means that the strand's reverse complement can match its own nucleotide sequence. The reaction has molecularity two in the association direction, unlike intramolecular RNA hairpin melting.
Use the ideal association model . Let total strand molar concentration be and paired-strand fraction be . With reference molar concentration , the dimensionless equilibrium constant isSubstituting and givesThe final expression is the physical solution for a self-complementary RNA duplex population; the algebraically equivalent subtractive form loses numerical precision at small .
At , half the strands are paired. For an original example, this means and . Duplexes are only one-third of the independently diffusing molecular objects, . Half the strands is not half the objects. At fixed temperature and , increasing tenfold makes and .
Heating changes and thus the population. At the strand-based midpoint, , which explains the concentration dependence of RNA melting. A fit should state whether its constant describes association or dissociation and whether concentration counts strands or duplexes.
Real sequences can also form competing hairpins. A primary tetraloop study explicitly compared predicted hairpin and duplex equilibria when designing its optical-melting sequences.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, Design of sequences for optical melting studies.
The equations above assume only separated strands and duplexes; they do not automatically cover competing folds or RNA melting hysteresis. Self-complementary RNA duplex population 2026-11-05
A self-complementary RNA duplex population describes how identical RNA strands are distributed between free strands and duplexes under specified conditions. A fraction of strands in duplexes differs from a fraction of molecular objects that are duplexes. The distinction matters when comparing chemical mass balances with observed signals.
For the ideal equilibrium , let and . With dimensionless association constant and reference molar concentration , define . The strand balance givesThis is the population model underlying ideal self-complementary RNA duplex melting. It assumes separated strands and duplexes are the only relevant species and that chemical equilibrium is reached.
The concentration regimes are illuminating. For , ; for , . An original calculation gives at and at . Even a large association parameter leaves some free strands in this equilibrium model.
An original finite-count illustration has 100 strands,80 of which belong to 40 duplexes, with 20 free strands left. The paired-strand fraction is 0.8, but duplexes are of the 60 molecular objects. A detector responding equally to each object would not give the same signal as a detector responding equally per strand. RNA optical melting needs a model of the folded and unfolded optical responses to infer populations.
Differentiating the equilibrium relation givesAt , a small 10% relative change in changes by approximately 0.0167. Population is therefore not generally proportional to concentration.
Competing RNA hairpins require an expanded species balance. A primary tetraloop study evaluated hairpin and duplex formation explicitly when designing sequences for melting experiments.[ref]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, Design of sequences for optical melting studies.
Ignoring a populated competitor can produce a precise solution to the wrong balance equation. Two-state RNA melting model 2026-11-05
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 ,The population formula assumes an ideal unimolecular transition. With constant unfolding enthalpy and entropy,An RNA melting log-odds plot against 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 at and at . These two points determineThe fitted model predicts approximately at . 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, . 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]
Thermodynamic characterization of naturally occurring RNA tetraloops by Justin P. Sheehy and others. Materials and Methods, optical curves and individual-curve thermodynamic fitting.
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.