Sequential Precursor Maturation With Irreversible Conversion

Why Model Sequential Precursor Maturation?

A biomarker may not be produced directly from its initial source. Instead, material can pass through several precursor pools, with each pool either converting to the next or being lost by another route. This sequence matters: the terminal biomarker’s steady-state concentration depends not only on its own clearance, but also on the fraction of input that survives each upstream pool and reaches it.

The earlier material on precursor-product models from biochemical mass balance provides the basis for writing pool-specific input and output terms. Here, the focus is a sequence of pools with irreversible, forward-only conversion, and the calculation of steady-state concentrations from the rate constants.

All \(P_i\) and \(B\) are concentrations in the same units. The input \(R\) is a concentration per unit time, and \(k_i\), \(\lambda_i\), and \(k_B\) have units of inverse time. The model assumes constant input, constant rate constants, well-mixed pools, and one-for-one conversion in concentration units. It also assumes no feedback between pools and no saturation of conversion or loss. Rate constants are nonnegative, and every precursor pool has a positive total removal rate \(k_i+\lambda_i\); \(k_B\) must also be positive for a finite biomarker steady state.

Mass Balance and the Steady-State Result

For the first pool, input is \(R\), while total output is the sum of conversion and competing loss. Each later pool receives the preceding pool’s conversion flux. The terminal biomarker receives conversion from the last precursor and is cleared at rate \(k_B\).

At steady state, each derivative is zero. Setting the first equation to zero gives \(P_1^*=R/(k_1+\lambda_1)\). Setting each subsequent precursor equation to zero gives its concentration as its incoming conversion flux divided by its total removal rate. The biomarker concentration is its formation flux divided by its clearance rate.

Repeated substitution gives a closed form. An empty product, as occurs for \(P_1^*\), is defined to equal 1.

The equations show why the distinction between conversion and loss is important. For pool \(i\), \(k_i+\lambda_i\) determines how quickly material leaves that pool, while the fraction \(k_i/(k_i+\lambda_i)\) is the fraction of its outflow that continues along the modeled sequence rather than being lost. Thus, the terminal biomarker formation flux is the input multiplied by the successive conversion fractions.

This flux expression is a useful interpretation and calculation check. It is not itself the biomarker concentration: dividing the formation flux by \(k_B\) is still necessary to obtain \(B^*\).

Worked Steady-State Calculations

Derivation, Checks, and Model Diagnostics

The recursion follows directly from setting each mass-balance equation to zero. For \(P_1\), \(0=R-(k_1+\lambda_1)P_1^*\), so \(P_1^*=R/(k_1+\lambda_1)\). If \(P_{i-1}^*\) is known, then \(0=k_{i-1}P_{i-1}^*-(k_i+\lambda_i)P_i^*\), which rearranges to \(P_i^*=k_{i-1}P_{i-1}^*/(k_i+\lambda_i)\). Substituting this expression repeatedly produces the product formula. Finally, \(0=k_mP_m^*-k_BB^*\) yields \(B^*=k_mP_m^*/k_B\). Positive total removal rates make each division valid and determine one finite steady state for a given nonnegative input.

Several checks are useful when implementing or interpreting the model:

• Check units. Each concentration is a rate divided by a rate, or a concentration multiplied by a rate and divided by a rate. The result must have concentration units.
• Check direction and indexing. \(k_i\) is the conversion rate out of pool \(i\), so \(k_iP_i\) is that pool’s forward conversion flux. The terminal conversion is \(k_mP_m\), not a rate indexed beyond the last pool.
• Check fluxes at steady state. For each precursor, incoming flux equals conversion plus competing loss. For the biomarker, formation equals clearance.
• Check nonnegative concentrations. Nonnegative input and rates, together with positive total removal rates, cannot produce a negative steady-state concentration. A negative value signals an algebraic, coding, or parameter error.
• Check the whole-system balance. With one-for-one conversion and no other inputs or outputs, \(R\) must equal the sum of all competing precursor losses and biomarker clearance at steady state.

A structural diagnostic is whether the assumed path is plausible: if measured or mechanistic evidence supports return flux, branching inputs, changing synthesis, or saturable rates, this forward-only linear model omits relevant processes. A numerical diagnostic is whether the reported “steady state” satisfies every equation when the computed concentrations are substituted back in. Small rounding differences are expected; a material mismatch is not.

How This Model Compares With Alternatives

A single-pool turnover model has one production input and one loss process. It can summarize a terminal biomarker’s net turnover, but it cannot separately represent intermediate precursor concentrations or the loss that competes with conversion at each maturation step. The sequential model adds those explicit pools, at the cost of requiring more rate constants and, where data are available, more information to estimate them.

The defining contrast with a reversible precursor-product exchange model is pathway direction: here, conversion is irreversible and there is no backward flux. If backward exchange is biologically important, the one-way equations do not describe the mass balance. Likewise, this tutorial’s steady-state expressions rely on linear first-order rates; nonlinear or time-varying systems generally require their own equations and may need numerical solution rather than these formulas.

These expressions calculate steady states, not the time needed to reach them or the full response after a changing input. They also do not establish that the rate constants can be uniquely estimated from a particular data set. A mathematically defined steady state and an empirically well-determined set of parameters are separate matters.

Common Mistakes and Limitations

Another frequent error is to report the last conversion flux as the biomarker concentration. A flux has units of concentration per time; the concentration requires division by \(k_B\). Indexing errors can cause a similar problem: in an \(m\)-pool system, \(k_m\) converts the final precursor pool to the biomarker. A rate constant for a nonexistent next precursor should not be introduced.

The model’s simplicity is useful when its assumptions are defensible, but it is a limitation when they are not. Constant rate constants exclude saturation and time-dependent regulation. A single forward path excludes branches and return pathways. One-for-one conversion excludes changes in yield or units between pools. If these features matter, the mass-balance equations must be modified before steady-state calculations are made.