Distinguishing Explicit Precursor Pools From Transit Chains

Why Distinguish the Two Structures?

A delayed biomarker response can arise because drug action changes a real biochemical pool that must form or convert before biomarker appears. It can also be represented empirically by a sequence of transit stages that smooths and delays an input signal. Both structures can produce gradual responses, so a plausible-looking curve alone does not establish which mechanism is correct. The distinction matters when interpreting parameters or predicting the result of perturbing a particular biochemical step.

The earlier material on precursor-product models and transit compartments establishes the basic mass-balance and chain equations. Here the central question is comparative: what biological commitments does each structure make, and which observations could reveal those commitments?

Model Definitions and Assumptions

Consider an irreversible precursor-to-product model with constant formation rate \(R\), precursor loss rate constant \(k_p\), conversion rate constant \(k_c\), and biomarker loss rate constant \(k_b\):

Amounts are measured in units, \(R\) in units per hour, and each \(k\) in \(\mathrm{h}^{-1}\). The assumptions are first-order loss and conversion, a well-mixed precursor pool, and no return from \(B\) to \(P\). At a constant baseline, the steady-state amounts are \(P_0=R/(k_p+k_c)\) and \(B_0=k_cP_0/k_b\). Thus \(k_p\) and \(k_c\) are distinct routes out of \(P\), not interchangeable labels for one rate.

For comparison, an empirical two-stage transit delay can be written with dimensionless signals \(z_1(t)\) and \(z_2(t)\):

Here \(u(t)\) is a normalized drug-modulated input, \(S_0\) is baseline biomarker production in units per hour, and \(k_{\mathrm{tr}}\) is the shared transit rate. At baseline \(u=z_1=z_2=1\), giving \(B_0=S_0/k_b\). The chain assumes sequential, first-order stages with the specified common rate; unless supported by independent biology or measurements, its intermediate signals are delay states rather than identified precursor amounts.

What the Equations Imply

For a constant input, the explicit pool approaches its new equilibrium with rate \(k_p+k_c\). Its mean residence time under the stated first-order assumptions is \(1/(k_p+k_c)\); the fraction of precursor leaving by conversion is \(k_c/(k_p+k_c)\). A transit chain instead distributes delay across stages. For \(m\) equal-rate stages, its mean transit time is \(m/k_{\mathrm{tr}}\), while the spread of transit times depends on both \(m\) and \(k_{\mathrm{tr}}\). Matching a mean delay does not make the two response shapes identical.

A short derivation shows how a formation step propagates through the explicit pool. Suppose \(R\) changes at time zero to a new constant \(R_1\), while the rate constants remain fixed. Let \(\lambda=k_p+k_c\). Solving the precursor equation gives \(P(t)=P_1+[P(0)-P_1]e^{-\lambda t}\), where \(P_1=R_1/\lambda\). Substituting this into the biomarker equation yields a constant term plus exponentials with rates \(\lambda\) and \(k_b\). When \(\lambda\ne k_b\):

The result follows by multiplying \(dB/dt+k_bB=k_cP(t)\) by the integrating factor \(e^{k_bt}\), integrating from zero to \(t\), and collecting terms. If \(\lambda=k_b\), the coincident exponential rates produce a term proportional to \(t e^{-k_bt}\), obtained by taking the limit. This response is a consequence of the specified two-pool mechanism, not a universal shape for delayed biomarker data.

Comparing Observable Predictions

The most useful distinctions often concern what changes immediately and what must change first. If formation \(R\) decreases in the explicit model, the precursor begins changing at once, but \(B\)'s instantaneous slope is initially unchanged: at the moment of the step, \(P\) and \(B\) still have their pre-step values, so \(dB/dt=k_cP_0-k_bB_0=0\). The biomarker changes after precursor availability changes. In a transit chain, an input step first changes \(z_1\), then \(z_2\), and only then changes the production term for \(B\). This is a sequential delay, but the transit states do not by themselves assert measurable precursor concentrations.

A selective change to the conversion rate makes a stronger mechanistic prediction. In the explicit model, changing \(k_c\) alters both the rate at which precursor is converted and the total rate at which precursor leaves. A transit chain whose input modulation acts only through \(u(t)\) has no corresponding conversion parameter unless one is deliberately added. Consequently, the two structures may agree for one type of input change yet predict different outcomes for another.

Diagnostics, Alternatives, and Limitations

Evidence favoring an explicit precursor interpretation is strongest when the precursor can be measured, its amount changes as predicted, and independent knowledge supports the proposed formation, conversion, and loss pathways. A selectively targeted perturbation can also be informative if its action on a particular rate is well justified. By contrast, response-shape information alone—such as a gradual onset or a long delay—usually supports a delay structure, not a unique biochemical identity for its stages. With only biomarker measurements under a limited set of inputs, different structures can produce similar curves.

The main alternative is not that one model is universally superior. A transit chain can be a useful parsimonious description when the delay is reproducible but the underlying stages are unknown or unmeasured. An explicit pool is preferable when its states and pathways have independent biological meaning and its additional predictions are relevant. Both models may be inadequate if rates vary over time, conversion is saturable, feedback is important, or the assumed input effect is wrong. A better visual fit alone does not validate a mechanistic interpretation.

• Warning sign: interpreting transit states as measured chemistry. A fitted stage amount has no biochemical units or direct interpretation unless the model defines and supports one.
• Warning sign: treating equal mean delays as equivalent models. Different stage structures can have different response shapes despite sharing a mean delay.
• Warning sign: assigning a drug effect to the wrong step. Formation inhibition and conversion inhibition can produce different early slopes and different precursor predictions.
• Diagnostic practice: compare predicted and observed trajectories for all measured states and perturbations, inspect whether early and late changes follow the assumed pathways, and consider whether conclusions depend on an unmeasured state.

A common mistake is to see a delayed biomarker decline and conclude that the transit stages must be literal precursor pools. That reasoning confuses a mathematical representation of delay with a biochemical claim. The correction is to ask what each model predicts for independently observed precursor amounts and for perturbations targeted to specific rates; without such information, describe the delay as an empirical structure rather than established mechanism.