First-Order mRNA Synthesis and Degradation

From mRNA Fluxes to a Balance Equation

A transcript abundance is determined by two opposing processes: new mRNA is synthesized, and existing mRNA is lost. A simple model can represent synthesis as a constant input and degradation as a loss proportional to the amount of mRNA present. This formulation is useful when the aim is to calculate a basal steady-state abundance from specified synthesis and degradation rates.

Let $M(t)$ denote the mRNA abundance at time $t$. It may be a number of transcripts per cell or a concentration, provided the units are used consistently. Let $k_{\mathrm{syn}}$ be the constant synthesis rate, expressed in units of abundance per unit time, and let $k_{\mathrm{deg}}$ be the first-order degradation rate constant, expressed in inverse time.

Definition: In the constant-synthesis, first-order-degradation model, the synthesis flux is $k_{\mathrm{syn}}$, and the degradation flux is $k_{\mathrm{deg}}M(t)$. The mRNA balance equation is $$ \frac{dM}{dt}=k_{\mathrm{syn}}-k_{\mathrm{deg}}M(t). $$ Here $k_{\mathrm{syn}}$ is constant and nonnegative, and $k_{\mathrm{deg}}$ is a nonnegative rate constant. When $M$ is measured as transcripts per cell, the synthesis flux has units of transcripts per cell per unit time, while $k_{\mathrm{deg}}$ has units of inverse time.

The equation states that the rate of change in mRNA abundance equals the rate of input minus the rate of loss. The first-order assumption means that doubling $M$, while keeping $k_{\mathrm{deg}}$ fixed, doubles the degradation flux. It does not mean that degradation removes a fixed number of transcripts per unit time; that would be a constant loss flux, not a first-order one.

This model is an application of the earlier material on state variables and fluxes: $M$ is the state variable, and the two terms on the right-hand side are fluxes with units matching $dM/dt$. In the biological setting considered here, “constant synthesis” means the transcriptional input is treated as unchanged over the modeled condition. It is not a claim that transcription is universally constant in a cell.

Steady State: Balance of the Two Fluxes

A steady state is an abundance at which the state variable no longer changes with time. For this model, set the derivative to zero and solve the resulting algebraic equation. If $k_{\mathrm{deg}}>0$, there is one steady-state abundance:

0=k_{\mathrm{syn}}-k_{\mathrm{deg}}M_{\mathrm{ss}} \qquad\Longrightarrow\qquad M_{\mathrm{ss}}=\frac{k_{\mathrm{syn}}}{k_{\mathrm{deg}}}.

At this abundance, synthesis and degradation fluxes are equal: $k_{\mathrm{syn}}=k_{\mathrm{deg}}M_{\mathrm{ss}}$. The equation also passes a dimensional check. Dividing abundance per unit time by inverse time gives abundance. For example, a synthesis rate in transcripts per cell per hour divided by a degradation rate constant in per hour gives transcripts per cell.

The steady-state calculation requires more than the algebraic balance. It assumes the synthesis and degradation rates are constant in the condition being represented, the degradation flux is adequately proportional to $M$, and the modeled mRNA pool can be treated as a single well-mixed state. The system must also have had sufficient time under those conditions to approach the steady state. If the abundance is still changing appreciably, the ratio $k_{\mathrm{syn}}/k_{\mathrm{deg}}$ is an equilibrium prediction, not necessarily the abundance at the observation time.

Why the Steady State Is Stable

For $k_{\mathrm{deg}}>0$, the steady state is not only a flux balance; small deviations from it decay under the model. To see this, define the deviation $\Delta(t)=M(t)-M_{\mathrm{ss}}$. Since $M_{\mathrm{ss}}=k_{\mathrm{syn}}/k_{\mathrm{deg}}$ is constant, subtract its balance equation from the mRNA balance equation:

\frac{d\Delta}{dt} =\frac{dM}{dt} =k_{\mathrm{syn}}-k_{\mathrm{deg}}M =-k_{\mathrm{deg}}\bigl(M-M_{\mathrm{ss}}\bigr) =-k_{\mathrm{deg}}\Delta.

If $\Delta$ is positive, its derivative is negative, so the abundance moves downward toward the steady state. If $\Delta$ is negative, its derivative is positive, so the abundance moves upward. More formally, the deviation has the form $\Delta(t)=\Delta(0)e^{-k_{\mathrm{deg}}t}$, which tends to zero as time increases when $k_{\mathrm{deg}}>0$. This establishes stability for the stated model; it does not establish that the biological assumptions are appropriate.

Worked Example: Calculating Transcript Abundance per Cell

Suppose a hypothetical cell type has a constant mRNA synthesis rate of 120 transcripts per cell per hour and a first-order degradation rate constant of $0.40\ \mathrm{h}^{-1}$. Calculate the model’s steady-state abundance.

The units are compatible: the synthesis rate is abundance per hour, and the degradation constant is per hour. The steady-state formula is $M_{\mathrm{ss}}=k_{\mathrm{syn}}/k_{\mathrm{deg}}$. Substituting the given values,

M_{\mathrm{ss}} =\frac{120\ \text{transcripts per cell per hour}} {0.40\ \mathrm{h}^{-1}} =300\ \text{transcripts per cell}.

The predicted degradation flux at this abundance is $0.40\ \mathrm{h}^{-1}\times300$ transcripts per cell, or 120 transcripts per cell per hour. This equals the synthesis flux, as required at steady state. The answer is a model prediction under the specified constant-rate assumptions, not evidence that a measured cell population has exactly this abundance.

Worked Example: Using Concentration Units

For a hypothetical mRNA measured as a cellular concentration, let the constant synthesis rate be $0.18\ \mathrm{nM\,h}^{-1}$ and the first-order degradation rate constant be $0.060\ \mathrm{h}^{-1}$. Find the steady-state concentration and verify the balance.

Because $k_{\mathrm{syn}}$ is concentration per hour and $k_{\mathrm{deg}}$ is inverse hours, their ratio has concentration units:

M_{\mathrm{ss}} =\frac{0.18\ \mathrm{nM\,h}^{-1}}{0.060\ \mathrm{h}^{-1}} =3.0\ \mathrm{nM}.

At $3.0\ \mathrm{nM}$, the degradation flux is $0.060\ \mathrm{h}^{-1}\times3.0\ \mathrm{nM} =0.18\ \mathrm{nM\,h}^{-1}$, equal to the synthesis rate. The calculation is valid only if the concentration-based synthesis and loss terms refer to the same modeled compartment and use compatible time units.

Comparing Conditions Without Losing the Rate Information

The steady-state formula makes clear that abundance depends on the ratio of synthesis to degradation, not on either rate alone. For two conditions that each satisfy the model assumptions, the ratio of their predicted steady states is

\frac{M_{\mathrm{ss},2}}{M_{\mathrm{ss},1}} = \frac{k_{\mathrm{syn},2}}{k_{\mathrm{syn},1}} \frac{k_{\mathrm{deg},1}}{k_{\mathrm{deg},2}}.

Thus an abundance can rise because synthesis is higher, degradation is lower, or both. An observed steady-state abundance by itself generally identifies only the ratio $k_{\mathrm{syn}}/k_{\mathrm{deg}}$, not the two rates separately. Independent information about at least one rate is needed to calculate the other from that abundance.

Worked Example: Comparing Two Hypothetical Conditions

In a hypothetical comparison, condition 1 has $k_{\mathrm{syn},1}=90$ transcripts per cell per hour and $k_{\mathrm{deg},1}=0.30\ \mathrm{h}^{-1}$. Condition 2 has $k_{\mathrm{syn},2}=54$ transcripts per cell per hour and $k_{\mathrm{deg},2}=0.15\ \mathrm{h}^{-1}$. Assuming both conditions meet the model assumptions and have reached steady state, compare their predicted abundances.

For condition 1,

M_{\mathrm{ss},1} =\frac{90}{0.30} =300\ \text{transcripts per cell}.

For condition 2,

M_{\mathrm{ss},2} =\frac{54}{0.15} =360\ \text{transcripts per cell}.

The ratio is $360/300=1.20$, so condition 2 has a predicted steady-state abundance 20% higher than condition 1. This occurs even though its synthesis rate is 40% lower: its degradation rate constant is 50% lower. The same conclusion follows from the rate ratio: $\left(54/90\right)\left(0.30/0.15\right)=0.60\times2=1.20$. This model-based comparison does not, by itself, establish which biological mechanism caused the difference.

Assumptions, Diagnostics, and Alternative Formulations

Before interpreting $M_{\mathrm{ss}}$, check whether the model matches the question and data. A useful sequence is to verify units, identify whether the rates are constant over the relevant condition, and assess whether the measurements are consistent with a stable abundance. These checks cannot prove the model, but they can reveal common mismatches.

Check the units. The two terms in the balance equation must have the same units as $dM/dt$. A synthesis rate expressed per hour cannot be divided directly by a degradation constant expressed per day without converting one time unit.
Check the proposed steady state. Repeated measurements under the same condition should be compatible with an abundance that has stopped changing, allowing for measurement variability. A clear trend over the observation period warns against treating the measured value as steady.
Check first-order loss. If degradation is not proportional to abundance across the range of interest, one constant $k_{\mathrm{deg}}$ may not describe the loss flux adequately.
Check what counts as the modeled pool. If important transport between compartments or distinct transcript pools affects the measured abundance, a single-state equation may omit relevant fluxes.

Alternative formulations should be chosen to address a specific failed assumption. A time-varying synthesis input is needed when transcription changes during the modeled period; a non-first-order loss term may be needed when loss is not proportional to abundance; and multiple state variables may be needed when distinct compartments or pools matter. These alternatives add assumptions and parameters, so additional complexity should be supported by the biological question and available data rather than added by default.

Common pitfall: “The mRNA abundance is 300 transcripts per cell, so the synthesis rate must be 300 transcripts per cell per hour.” This confuses abundance with flux. In this model, $M_{\mathrm{ss}}=k_{\mathrm{syn}}/k_{\mathrm{deg}}$; the same abundance could result from many synthesis and degradation rate pairs. For example, $300$ transcripts per cell is consistent with $k_{\mathrm{syn}}=120$ transcripts per cell per hour and $k_{\mathrm{deg}}=0.40\ \mathrm{h}^{-1}$, but not with the claim that abundance alone specifies either rate.

There is also an important boundary case. If $k_{\mathrm{deg}}=0$ and $k_{\mathrm{syn}}>0$, the equation gives $dM/dt=k_{\mathrm{syn}}$, so no finite steady state exists. If both rates are zero, every abundance has zero derivative, and the model does not select a unique steady state. The ratio formula therefore requires $k_{\mathrm{deg}}>0$.

The main limitation is that a correct algebraic balance is only as useful as the representation behind it. A single steady-state measurement cannot distinguish a high-synthesis/high-degradation system from a low-synthesis/low-degradation system when their ratios match. Nor does a steady-state calculation establish that the abundance has actually equilibrated. Treat $M_{\mathrm{ss}}$ as a conditional model prediction: it follows from the specified rates and assumptions, not from the word “mRNA” alone.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

A hypothetical cell has $k_{\mathrm{syn}}=75$ transcripts per cell per hour and $k_{\mathrm{deg}}=0.25\ \mathrm{h}^{-1}$. What steady-state abundance does the model predict?
A model uses a synthesis rate in $\mathrm{nM\,h}^{-1}$ and a degradation constant in $\mathrm{h}^{-1}$. What are the units of the predicted steady-state abundance?
If $k_{\mathrm{syn}}=0$ and $k_{\mathrm{deg}}>0$, what steady-state abundance follows? Explain using the balance equation.
Two conditions have identical steady-state abundances. Does that imply their synthesis and degradation rates are each identical? Explain.

Solutions. (1) $M_{\mathrm{ss}}=75/0.25=300$ transcripts per cell. The units reduce from transcripts per cell per hour divided by per hour to transcripts per cell. (2) The predicted abundance is in nM, because $(\mathrm{nM\,h}^{-1})/(\mathrm{h}^{-1})=\mathrm{nM}$. (3) Setting $dM/dt=0$ gives $0=0-k_{\mathrm{deg}}M_{\mathrm{ss}}$. Since $k_{\mathrm{deg}}>0$, $M_{\mathrm{ss}}=0$. (4) No. Equal abundance means the ratios $k_{\mathrm{syn}}/k_{\mathrm{deg}}$ are equal; different synthesis and degradation rates can have the same ratio.