Deriving Steady States in Turnover Models

Why Solve for Turnover Steady States?

Cell counts reflect a balance of flows, not just the rate at which cells are produced. A precursor population can be small because it matures rapidly, while the mature population can be large if its loss is slow. Solving for steady-state abundances makes these relationships explicit and helps distinguish changes in production from changes in maturation or survival.

Consider a precursor compartment $P$ that receives a constant input and produces mature cells $M$. Precursors can leave by maturing or by another first-order loss process; mature cells also have a first-order loss. The aim here is to find the abundances at which each compartment's inflows and outflows balance. This calculation identifies an equilibrium, but by itself does not establish how the system behaves over time or whether it approaches that equilibrium.

Model and assumptions: Let $I$ be the constant input rate into the precursor compartment, $k_m$ the first-order maturation rate constant, $\delta_P$ the other first-order precursor-loss rate constant, and $\delta_M$ the first-order mature-cell-loss rate constant. The compartment abundances are $P(t)$ and $M(t)$. Assume these rates are constant, cells move from $P$ to $M$ without branching, and there are no density-dependent effects, feedback, or other sources and sinks. The balances are $$ \frac{dP}{dt}=I-(k_m+\delta_P)P,\qquad \frac{dM}{dt}=k_mP-\delta_MM. $$ Here $I$ has units of cells per unit time, the rate constants have units of inverse time, and $P$ and $M$ have units of cells.

Maturation is an outflow from the precursor compartment and an equal inflow to the mature compartment. It is therefore included in the total precursor outflow, but it is not a loss from the two-compartment system as a whole. The term $\delta_P P$ is the precursor loss that does not enter $M$, and $\delta_M M$ is mature-cell loss.

Deriving the Steady-State Abundances

A steady state is a pair of constant abundances for which both time derivatives are zero. For positive $I$, $k_m+\delta_P$, and $\delta_M$, solve the precursor balance first and substitute its solution into the mature-cell balance.

Set the precursor derivative to zero.
At steady state, $0=I-(k_m+\delta_P)P_{\mathrm{ss}}$. Thus $P_{\mathrm{ss}}=I/(k_m+\delta_P)$, provided $k_m+\delta_P>0$.

Set the mature-cell derivative to zero.
The mature-cell balance gives $0=k_mP_{\mathrm{ss}}-\delta_MM_{\mathrm{ss}}$. If $\delta_M>0$, then $M_{\mathrm{ss}}=k_mP_{\mathrm{ss}}/\delta_M$.

Substitute and check the flows.
Substitution gives $M_{\mathrm{ss}}=k_mI/[\delta_M(k_m+\delta_P)]$. At this pair of abundances, input equals total precursor outflow, and mature-cell formation equals mature-cell loss.

Steady-state result: For $I\geq 0$, $k_m\geq 0$, $k_m+\delta_P>0$, and $\delta_M>0$, the finite steady-state solution is $$ P_{\mathrm{ss}}=\frac{I}{k_m+\delta_P},\qquad M_{\mathrm{ss}}=\frac{k_m I}{\delta_M(k_m+\delta_P)}. $$ When $I>0$, the steady-state maturation flux is $k_mP_{\mathrm{ss}}$, and it equals the mature-cell loss flux $\delta_MM_{\mathrm{ss}}$. If $k_m=0$, the solution has $M_{\mathrm{ss}}=0$: no cells mature into that compartment.

A useful ratio follows directly from the two steady-state balances. When $P_{\mathrm{ss}}>0$, $$ \frac{M_{\mathrm{ss}}}{P_{\mathrm{ss}}}=\frac{k_m}{\delta_M}. $$ Thus the mature-to-precursor abundance ratio is set by maturation relative to mature-cell loss. This ratio does not depend on the input or on the other precursor-loss rate, although those quantities do affect the absolute abundances.

Worked Examples

Worked Example: Calculating Both Compartments and Checking Flows

Suppose a hypothetical cell system has a constant precursor input of $I=120$ cells/day, maturation rate $k_m=0.20$ day$^{-1}$, other precursor-loss rate $\delta_P=0.10$ day$^{-1}$, and mature-cell-loss rate $\delta_M=0.30$ day$^{-1}$. Calculate both steady-state abundances and check the balances.

The total first-order outflow rate from the precursor compartment is $0.20+0.10=0.30$ day$^{-1}$. Therefore, $$ P_{\mathrm{ss}}=\frac{120\ \text{cells/day}}{0.30\ \text{day}^{-1}}=400\ \text{cells}. $$ The mature-cell abundance is $$ M_{\mathrm{ss}}=\frac{(0.20\ \text{day}^{-1})(400\ \text{cells})}{0.30\ \text{day}^{-1}} =266.7\ \text{cells}. $$

The maturation flux is $0.20\times400=80$ cells/day, and the other precursor-loss flux is $0.10\times400=40$ cells/day. Their sum is $120$ cells/day, equal to the input. Mature-cell loss is $0.30\times266.7=80.01$ cells/day, which is $80$ cells/day up to rounding. The two balances therefore hold.

Worked Example: Increasing the Maturation Rate

Use the same hypothetical parameters as in the previous calculation, except that $k_m$ increases from $0.20$ to $0.40$ day$^{-1}$, with $I=120$ cells/day, $\delta_P=0.10$ day$^{-1}$, and $\delta_M=0.30$ day$^{-1}$ unchanged. Find the new steady state and compare it with the original one.

The precursor abundance is $$ P_{\mathrm{ss}}=\frac{120}{0.40+0.10} =\frac{120}{0.50}=240\ \text{cells}. $$ The mature-cell abundance is $$ M_{\mathrm{ss}}=\frac{(0.40)(240)}{0.30} =\frac{96}{0.30}=320\ \text{cells}. $$

The precursor pool falls from $400$ to $240$ cells because a larger fraction of its outflow is maturation. The mature-cell pool rises from $266.7$ to $320$ cells: the maturation flux is now $0.40\times240=96$ cells/day, balanced by $0.30\times320=96$ cells/day of mature-cell loss. The increase in maturation does not simply “increase every compartment”; it reduces the precursor abundance while increasing the mature abundance in this model.

Worked Example: Reducing Other Precursor Loss

Consider a different hypothetical system with $I=90$ cells/day, $k_m=0.15$ day$^{-1}$, $\delta_P=0.15$ day$^{-1}$, and $\delta_M=0.25$ day$^{-1}$. Calculate its steady state, then determine the effect of reducing $\delta_P$ to $0.05$ day$^{-1}$ while leaving the other parameters unchanged.

Initially, $$ P_{\mathrm{ss}}=\frac{90}{0.15+0.15}=300\ \text{cells}, \qquad M_{\mathrm{ss}}=\frac{(0.15)(300)}{0.25}=180\ \text{cells}. $$ After the reduction in precursor loss, $$ P_{\mathrm{ss}}=\frac{90}{0.15+0.05}=450\ \text{cells}, \qquad M_{\mathrm{ss}}=\frac{(0.15)(450)}{0.25}=270\ \text{cells}. $$

The maturation flux rises from $45$ to $67.5$ cells/day as more precursors remain available to mature. In the new steady state, the other precursor-loss flux is $0.05\times450=22.5$ cells/day; together with maturation, it equals the $90$-cell/day input. Mature-cell loss is $0.25\times270=67.5$ cells/day, matching maturation. In this example, reducing the non-maturation precursor loss increases both compartment abundances.

How Rate Changes Shift Abundances

The steady-state formulas give precise directions of change while the model assumptions and other parameter values are held fixed. The precursor abundance increases with input and decreases when either $k_m$ or $\delta_P$ increases. Its denominator, $k_m+\delta_P$, is the total rate at which precursors leave $P$.

The mature-cell abundance increases with input and decreases with mature-cell loss. The effect of maturation is more nuanced: increasing $k_m$ removes precursors more quickly, but also sends a larger share of the available precursor flow into $M$. For positive $\delta_P$, the net steady-state effect on mature-cell abundance is positive. For example, differentiating the expression for $M_{\mathrm{ss}}$ with respect to $k_m$, while holding the other parameters fixed, gives $$ \frac{\partial M_{\mathrm{ss}}}{\partial k_m} = \frac{I\delta_P}{\delta_M(k_m+\delta_P)^2}. $$ This derivative is positive when $I>0$, $\delta_P>0$, and $\delta_M>0$. If $\delta_P=0$, it is zero: with no other precursor loss, the steady-state maturation flux equals the input regardless of the maturation rate, so $M_{\mathrm{ss}}=I/\delta_M$.

Likewise, with $k_m>0$, increasing $\delta_P$ decreases $M_{\mathrm{ss}}$, because more precursors are lost before they mature. Increasing $\delta_M$ reduces mature-cell abundance in inverse proportion: doubling $\delta_M$ halves $M_{\mathrm{ss}}$ if the other parameters remain fixed. These are conditional comparisons within the specified model, not general claims about every biological system.

Assumptions, Diagnostics, and Alternatives

The calculation requires a constant input and constant first-order rate constants. It also assumes the maturation flux is exactly $k_mP$, that this flux enters $M$ without delay or branching, and that mature-cell loss is exactly $\delta_MM$. These assumptions make the balances linear. If rates depend on abundance, if input varies with time, or if cells enter through additional routes, the stated steady-state formulas need not apply.

A practical diagnostic is to substitute the computed abundances back into both original equations. Check that $I-(k_m+\delta_P)P_{\mathrm{ss}}=0$ and $k_mP_{\mathrm{ss}}-\delta_MM_{\mathrm{ss}}=0$, allowing only rounding differences. Also check dimensions: dividing cells/time by a rate in $1/\text{time}$ gives cells, and multiplying a rate by an abundance gives cells/time. Negative abundances or inconsistent units usually indicate an algebraic or parameter-definition error.

A finite unique solution requires positive total precursor outflow, $k_m+\delta_P>0$, and positive mature-cell loss, $\delta_M>0$. If $\delta_M=0$ while $I>0$ and $k_m>0$, mature cells have a continuing positive inflow and no loss, so a finite mature-cell steady state is impossible. Such edge cases should be checked rather than hidden by division by zero.

Directly setting each derivative to zero is usually the clearest method for this two-compartment chain. An equivalent alternative is to write the balances as a linear system and solve the simultaneous equations; this can be convenient for larger linear compartment systems. Numerical integration can illustrate a trajectory under the model, but it is not needed to calculate the equilibrium and should not replace checking the algebraic balances. Whether trajectories approach the equilibrium is a separate dynamical question.

Common Pitfall and Why the Result Matters

A common mistake is to treat maturation as if it were a loss from the whole system, or to assume that a faster maturation rate must increase both compartments. Maturation is a loss from $P$, but it is simultaneously an input to $M$. The worked comparison shows why the precursor pool can shrink while the mature pool grows. Keeping each flow attached to its correct compartment prevents this sign and interpretation error.

Steady-state calculations are useful in QSP because they translate a proposed mechanism into measurable abundance predictions. The equations show, for example, that the same mature-cell abundance might result from greater production, faster maturation, or slower mature-cell loss. An abundance alone may therefore not identify which mechanism changed; additional measurements or experiments may be needed. The formulas describe the equilibrium implied by the model and parameters, not proof that a biological system is at equilibrium.

Check Your Understanding

Try each question, then compare with the solutions below.

For $I=60$ cells/hour, $k_m=0.10$ hour$^{-1}$, $\delta_P=0.20$ hour$^{-1}$, and $\delta_M=0.50$ hour$^{-1}$, calculate $P_{\mathrm{ss}}$ and $M_{\mathrm{ss}}$.
At steady state, which flux balances the input into $P$? State the balance equation in words or symbols.
If $\delta_M$ doubles and all other rates and input stay fixed, what happens to $M_{\mathrm{ss}}$?
What happens to the mature-cell abundance as $k_m$ increases when $\delta_P=0$, with $I$ and $\delta_M$ positive and fixed?

Solutions. (1) $P_{\mathrm{ss}}=60/(0.10+0.20)=200$ cells. Then $M_{\mathrm{ss}}=(0.10)(200)/0.50=40$ cells. (2) Input balances both precursor outflows: $I=k_mP_{\mathrm{ss}}+\delta_PP_{\mathrm{ss}}$. (3) Since $M_{\mathrm{ss}}$ is inversely proportional to $\delta_M$, doubling it halves the mature-cell abundance. (4) It does not change: when $\delta_P=0$, $M_{\mathrm{ss}}=I/\delta_M$, independent of $k_m$ for positive $k_m$.