Basal Steady States and Turnover Half-Lives

Why Basal Steady States Matter

In a gene-expression model, a basal state is the constant level of mRNA and protein maintained when synthesis and loss rates remain unchanged. Even in this simple setting, an observed concentration is not determined by synthesis alone: it reflects the balance between production and turnover. This balance is useful in quantitative systems pharmacology because a measured abundance can be translated into a synthesis requirement only when the relevant loss rate is known or estimated.

The earlier material on deriving coupled mRNA–protein ODEs supplies the model structure. Here, the focus is on the equilibrium concentrations, the relation between first-order loss and half-life, and what can be inferred from those quantities. The calculations assume constant basal rates; they do not describe how the system responds when a rate changes.

Model and assumptions: Let $m(t)$ be mRNA concentration and $P(t)$ be protein concentration. In the basal linear expression model, $$ \frac{dm}{dt}=s_m-k_m m,\qquad \frac{dP}{dt}=k_{\mathrm{tr}}m-k_P P. $$ Here $s_m$ is the constant mRNA synthesis rate, $k_m$ is the first-order mRNA degradation constant, $k_{\mathrm{tr}}$ is the translation coefficient, and $k_P$ is the first-order protein degradation constant. Assume the rates are constant and nonnegative, $k_m>0$ and $k_P>0$, and there are no other sources, sinks, feedbacks, or capacity limits in the model. If concentrations are used, $s_m$ has units of mRNA concentration per time, $k_m$ and $k_P$ have units of inverse time, and $k_{\mathrm{tr}}m$ has units of protein concentration per time.

Deriving the Basal Concentrations

At a steady state, neither state variable changes with time. Set both derivatives in the model equal to zero. The mRNA balance is $0=s_m-k_m m^*$, where $m^*$ denotes the basal steady-state mRNA concentration. Solving gives $m^*=s_m/k_m$. For protein, the steady-state balance is $0=k_{\mathrm{tr}}m^*-k_P P^*$. Substituting the mRNA result gives the basal protein concentration.

m^*=\frac{s_m}{k_m},\qquad P^*=\frac{k_{\mathrm{tr}}m^*}{k_P} =\frac{k_{\mathrm{tr}}s_m}{k_m k_P}.

These expressions have a direct flux interpretation. At basal equilibrium, mRNA synthesis equals mRNA degradation, $s_m=k_m m^*$. Likewise, protein synthesis by translation equals protein degradation, $k_{\mathrm{tr}}m^*=k_P P^*$. The synthesis rates here are fluxes, not concentrations: for example, $s_m$ is measured in concentration per time, whereas $m^*$ is measured in concentration.

With positive synthesis and positive degradation constants, both steady-state concentrations are positive and uniquely determined by the rates. The equilibrium is also stable for this model: the matrix of the linear system is triangular, with diagonal entries $-k_m$ and $-k_P$. Its eigenvalues are therefore $-k_m$ and $-k_P$, both negative. This establishes that small deviations from the equilibrium decay under the constant-rate model; it does not make the rates identifiable from concentrations alone.

First-Order Degradation and Turnover Half-Life

For a species undergoing first-order degradation alone, its concentration $C$ satisfies $dC/dt=-kC$, where $k>0$ is the degradation constant. The solution is $C(t)=C(0)e^{-kt}$. The half-life $t_{1/2}$ is the elapsed time at which $C(t_{1/2})=C(0)/2$. Substituting that condition and taking logarithms gives $e^{-kt_{1/2}}=1/2$, so $t_{1/2}=\ln(2)/k$.

Definition: Under first-order degradation, the turnover half-life is the time required for the amount of a species to fall to one-half its starting amount when no synthesis replenishes it. The degradation constant and half-life are related by $$ t_{1/2}=\frac{\ln(2)}{k},\qquad k=\frac{\ln(2)}{t_{1/2}}. $$ The half-life has units of time; the degradation constant has units of inverse time.

The same conversion applies separately to mRNA and protein: $k_m=\ln(2)/t_{1/2,m}$ and $k_P=\ln(2)/t_{1/2,P}$. A shorter half-life corresponds to a larger degradation constant. These relations concern a first-order loss process; they should not be applied as if a concentration’s observed decline were caused by degradation alone when synthesis or other processes also affect that decline.

Worked Example: From Half-Lives to Basal Synthesis Rates

Suppose a hypothetical cell has a basal mRNA concentration of $2.00$ nM and a basal protein concentration of $50.0$ nM. The mRNA and protein degradation half-lives are $4.00$ h and $10.0$ h, respectively. Calculate the degradation constants, the synthesis fluxes required at equilibrium, and the translation coefficient.

mRNA degradation constant. Convert the half-life to a rate: $$ k_m=\frac{\ln(2)}{4.00\ \mathrm{h}} =\frac{0.6931}{4.00\ \mathrm{h}} =0.1733\ \mathrm{h}^{-1}. $$

mRNA synthesis rate. At steady state, mRNA synthesis equals mRNA degradation: $$ s_m=k_m m^* =(0.1733\ \mathrm{h}^{-1})(2.00\ \mathrm{nM}) =0.3466\ \mathrm{nM/h}. $$ Thus, maintaining the specified basal mRNA concentration requires a synthesis flux of $0.3466$ nM/h under the model assumptions.

Protein degradation constant and synthesis flux. The protein half-life gives $$ k_P=\frac{\ln(2)}{10.0\ \mathrm{h}} =0.06931\ \mathrm{h}^{-1}. $$ At the protein steady state, translation must supply the protein removed by degradation: $$ k_{\mathrm{tr}}m^*=k_P P^* =(0.06931\ \mathrm{h}^{-1})(50.0\ \mathrm{nM}) =3.466\ \mathrm{nM/h}. $$

Translation coefficient. Divide the required protein synthesis flux by the mRNA concentration: $$ k_{\mathrm{tr}}=\frac{3.466\ \mathrm{nM/h}}{2.00\ \mathrm{nM}} =1.733\ \mathrm{h}^{-1}. $$ Because both concentrations are expressed in nM, their concentration units cancel here. The calculation is a balance implied by the assumed model, not evidence that the hypothetical values were measured in an experiment.

How Steady-State Measurements Constrain Synthesis

The equilibrium equations can be rearranged to estimate synthesis rates when the corresponding steady-state concentrations and degradation constants are available: $$ s_m=k_m m^*,\qquad k_{\mathrm{tr}}m^*=k_P P^*. $$ Thus, an mRNA concentration and its degradation constant determine the mRNA synthesis flux. A protein concentration and its degradation constant determine the protein synthesis flux. If the mRNA concentration is also known, that protein flux can be expressed as a translation coefficient.

Half-life measurements can supply the degradation constants through the first-order conversion. But a basal concentration by itself does not determine both a synthesis rate and a degradation rate. For mRNA, the same $m^*$ can result from many pairs $(s_m,k_m)$ satisfying $s_m=k_m m^*$. The protein balance has the same structure: many combinations of translation and degradation can maintain the same $P^*$. Independent turnover or synthesis information is needed to distinguish among those combinations.

Worked Example: Calculating Synthesis From Half-Lives

Consider another hypothetical basal state with $s_m=0.600$ nM/h, an mRNA half-life of $3.00$ h, $k_{\mathrm{tr}}=2.50\ \mathrm{h}^{-1}$, and a protein half-life of $12.0$ h. Calculate the steady-state mRNA and protein concentrations.

mRNA concentration. First convert the mRNA half-life: $$ k_m=\frac{0.6931}{3.00\ \mathrm{h}}=0.2310\ \mathrm{h}^{-1}. $$ Then use $m^*=s_m/k_m$: $$ m^*=\frac{0.600\ \mathrm{nM/h}}{0.2310\ \mathrm{h}^{-1}} =2.597\ \mathrm{nM}. $$

Protein synthesis flux. Translation produces protein at the rate $$ k_{\mathrm{tr}}m^* =(2.50\ \mathrm{h}^{-1})(2.597\ \mathrm{nM}) =6.493\ \mathrm{nM/h}. $$

Protein concentration. The protein degradation constant is $$ k_P=\frac{0.6931}{12.0\ \mathrm{h}}=0.05776\ \mathrm{h}^{-1}. $$ Consequently, $$ P^*=\frac{k_{\mathrm{tr}}m^*}{k_P} =\frac{6.493\ \mathrm{nM/h}}{0.05776\ \mathrm{h}^{-1}} =112.4\ \mathrm{nM}. $$ At this equilibrium, mRNA synthesis is $0.600$ nM/h and mRNA degradation is $(0.2310)(2.597)=0.600$ nM/h. Protein translation and degradation are both $6.493$ nM/h, as required by the two balance equations.

Worked Example: Why Abundance Alone Does Not Identify Turnover

Suppose hypothetical measurements give $m^*=1.50$ nM and $P^*=30.0$ nM. Show that these steady-state concentrations do not uniquely determine all the rates in the model by constructing two rate sets that produce the same concentrations.

Rate set A. Choose $k_m=0.200\ \mathrm{h}^{-1}$, $k_{\mathrm{tr}}=0.800\ \mathrm{h}^{-1}$, and $k_P=0.0400\ \mathrm{h}^{-1}$. To maintain $m^*=1.50$ nM, mRNA synthesis must be $$ s_m=k_m m^*=(0.200)(1.50)=0.300\ \mathrm{nM/h}. $$ The protein synthesis and degradation fluxes are both $$ k_{\mathrm{tr}}m^*=(0.800)(1.50)=1.20\ \mathrm{nM/h}, \qquad k_P P^*=(0.0400)(30.0)=1.20\ \mathrm{nM/h}. $$

Rate set B. Instead choose $k_m=0.100\ \mathrm{h}^{-1}$, $k_{\mathrm{tr}}=0.400\ \mathrm{h}^{-1}$, and $k_P=0.0200\ \mathrm{h}^{-1}$. The required mRNA synthesis is now $(0.100)(1.50)=0.150$ nM/h. The protein synthesis and degradation fluxes are both $(0.400)(1.50)=0.600$ nM/h and $(0.0200)(30.0)=0.600$ nM/h. This second set also maintains $m^*=1.50$ nM and $P^*=30.0$ nM.

The two distinct sets of rates give identical basal concentrations. Therefore, the concentrations alone cannot identify the individual synthesis and degradation constants. Additional measurements, such as a turnover half-life, would constrain the possible rates.

Assumptions, Diagnostics, and Alternatives

The algebra is simple, but its biological interpretation depends on the model assumptions. The synthesis rates must be constant over the period represented by the basal state; degradation must be first-order over the relevant concentration range; and the modeled species must have no unrepresented source or loss that materially affects the balance. If, for example, an effective loss includes both degradation and dilution, the relevant constant is the total first-order loss constant, not degradation alone. A saturable degradation process does not generally obey $k_P P^*$ with a fixed $k_P$.

Several checks help identify inconsistent calculations or a poor model fit:

Check units. A degradation constant must have inverse-time units, and multiplying it by a concentration must give a concentration-per-time flux.
Check both flux balances. At the proposed basal state, verify $s_m=k_m m^*$ and $k_{\mathrm{tr}}m^*=k_P P^*$ separately.
Check positivity and plausibility. Negative rates or concentrations are inconsistent with the stated model. Implausibly large inferred synthesis may point to a unit error, an incorrect half-life, or missing model processes.
Check the first-order assumption. A half-life treated as constant should not vary substantially with the starting concentration under the conditions represented by the model.

A useful alternative to solving the equations by substitution is to draw a flux balance for each species: list its production flux and loss flux, then set them equal at equilibrium. This gives the same result and makes missing terms easier to notice. A more detailed model may include additional sources, losses, or nonlinear rates, but then the simple ratios $s_m/k_m$ and $k_{\mathrm{tr}}m^*/k_P$ need not apply. The linear model is valuable as a transparent baseline, not a guarantee that every expression system follows first-order turnover.

Common pitfall: “The basal mRNA concentration is high, so mRNA synthesis must be high.” This does not follow from abundance alone. Since $m^*=s_m/k_m$, a high concentration can reflect a high synthesis rate, a slow degradation rate, or both. To infer synthesis from abundance, the relevant degradation constant must also be known or assumed.

Check Your Understanding

Answer each question using the linear basal model and show the balance or conversion used. Solutions follow the questions.

A species has a first-order degradation half-life of $6.00$ h. What is its degradation constant?
If $m^*=3.00$ nM and $k_m=0.250\ \mathrm{h}^{-1}$, what mRNA synthesis rate maintains that steady state?
A protein has $P^*=20.0$ nM and $k_P=0.0500\ \mathrm{h}^{-1}$. What protein synthesis flux is required at equilibrium?
Can a measured basal mRNA concentration by itself determine both $s_m$ and $k_m$? Explain briefly.

Solutions. (1) $k=\ln(2)/(6.00\ \mathrm{h})=0.1155\ \mathrm{h}^{-1}$. (2) $s_m=k_m m^*=(0.250\ \mathrm{h}^{-1})(3.00\ \mathrm{nM})=0.750\ \mathrm{nM/h}$. (3) The required protein synthesis flux equals degradation: $k_P P^*=(0.0500\ \mathrm{h}^{-1})(20.0\ \mathrm{nM})=1.00\ \mathrm{nM/h}$. (4) No. The steady-state equation gives $s_m=k_m m^*$, so many synthesis–degradation pairs yield the same measured concentration unless an additional rate or turnover measurement is available.