Why Baseline Concentration Does Not Determine Baseline Input
A biomarker concentration at baseline is often treated as if it directly determines the rate at which the biomarker is produced. It does not. At steady state, the observed concentration reflects a balance between endogenous input and loss: a high concentration can result from high input, slow loss, or both. If baseline is estimated rather than fixed, the model may fit the observed concentrations well while leaving the input and turnover rate poorly separated.
The practical question is therefore not just whether a fitted curve passes through the data. It is whether the observations contain enough distinct information to estimate the baseline concentration, endogenous input, and turnover rate separately. Sensitivity calculations and profile plots help reveal when estimates depend strongly on assumptions or sparse observations.
The units provide an immediate check: multiplying concentration by an inverse-time rate gives concentration per time, the units of input. The equilibrium identity also makes clear that one baseline concentration supplies one relationship between two rates. Without further information, it cannot uniquely determine both.
What the Data Can Identify
For constant input and first-order loss, the turnover equation is
At equilibrium, the derivative is zero, giving \(R_{\mathrm{base}}=k_{\mathrm{in}}/k_{\mathrm{out}}\). If the only observations are steady-state baseline concentrations, the data may estimate this ratio, but do not separately identify the numerator and denominator. A post-perturbation time course can provide information about \(k_{\mathrm{out}}\) through its rate of change, but that information depends on the perturbation being adequately specified and on the sampling times capturing the response.
Worked Example: Two Turnover Rates With the Same Baseline
Suppose a hypothetical biomarker has an observed baseline concentration of \(100\) nmol/L. Consider two candidate turnover rates: \(k_{\mathrm{out}}=0.10\) day\(^{-1}\) and \(k_{\mathrm{out}}=0.20\) day\(^{-1}\). For each candidate, calculate the input required to maintain the same baseline.
For \(k_{\mathrm{out}}=0.10\) day\(^{-1}\),
For \(k_{\mathrm{out}}=0.20\) day\(^{-1}\),
Both parameter pairs predict the same steady-state concentration, \(100\) nmol/L. The second pair has twice the input and twice the loss rate. Baseline data alone cannot choose between them; observations that capture a known perturbation or another source of information are needed.
Profiling Baseline Instead of Fixing It
Consider a hypothetical washout in which endogenous input is completely suppressed at time zero and loss remains first-order. If the concentration at the start of washout is the unknown \(R_{\mathrm{base}}\), then
This form assumes the suppression is complete and immediate, \(k_{\mathrm{out}}\) is constant during observation, and the measured concentrations correspond to the stated times. For observations \(y_i\) at times \(t_i\), suppose the residuals are analyzed by unweighted least squares, with independent errors of common variance. At any fixed candidate value of \(k_{\mathrm{out}}\), the least-squares estimate of \(R_{\mathrm{base}}\) is available in closed form.
This is called profiling: hold one parameter fixed, optimize the other parameter(s), and record the best achievable fit. A profile plot displays the resulting objective against the fixed parameter. A shallow or extended low region means that a range of turnover rates can be paired with different baseline estimates and still fit nearly as well. A sharp minimum indicates stronger information about the rate, though it does not by itself guarantee that the model assumptions are correct.
Worked Example: Profiling the Baseline in a Hypothetical Washout
Suppose hypothetical concentrations are \(100\), \(85\), and \(75\) nmol/L at times \(0\), \(1\), and \(2\) days, respectively. To see how baseline is estimated conditional on turnover, fix \(k_{\mathrm{out}}=0.20\) day\(^{-1}\). The exponential factors are \(1\), \(e^{-0.20}=0.81873\), and \(e^{-0.40}=0.67032\). Substitution in the profiling formula gives
The fitted predictions at days \(0\), \(1\), and \(2\) are therefore
The residuals, observed minus predicted, are approximately \(-3.727606\), \(0.075019\), and \(5.469306\) nmol/L. Thus the sum of squared errors is
This calculation is conditional on the selected turnover rate; it does not establish that \(0.20\) day\(^{-1}\) is uniquely supported. A profile plot would repeat the calculation over a range of rates, re-estimating baseline at every rate. If the curve is broad, reporting only the best-fitting rate and a standard error can conceal the trade-off.
Sensitivity Shows the Direction of the Trade-Off
A profile summarizes fit across parameter values. Local sensitivity gives a complementary view: it measures how much a prediction changes when a parameter changes slightly. Suppose a known fraction \(I\) of endogenous production is inhibited from time zero, with the baseline steady state \(R_{\mathrm{base}}\) before inhibition. The equation becomes \(dR/dt=k_{\mathrm{in}}(1-I)-k_{\mathrm{out}}R\). Substituting \(k_{\mathrm{in}}=k_{\mathrm{out}}R_{\mathrm{base}}\) and solving from \(R(0)=R_{\mathrm{base}}\) yields
This expression assumes a constant known inhibition fraction, unchanged first-order loss, and a pre-perturbation steady state. It shows that baseline scales the response, while turnover controls how quickly the response approaches its new level. At limited sampling times, changes in the two parameters can compensate, especially when the data do not clearly resolve the response’s curvature.
Worked Example: Comparing Baseline and Turnover Sensitivity
Take a hypothetical response with \(R_{\mathrm{base}}=100\) nmol/L, \(I=0.50\), \(k_{\mathrm{out}}=0.20\) day\(^{-1}\), and \(t=2\) days. Since \(e^{-0.40}=0.67032\), the response factor is
The local sensitivity to baseline is \(\partial R/\partial R_{\mathrm{base}}=0.83516\): a small increase of \(1\) nmol/L in baseline changes this prediction by about \(0.83516\) nmol/L. The sensitivity to turnover is
A local increase of \(0.01\) day\(^{-1}\) in turnover would therefore lower this prediction by about \(0.67032\) nmol/L, using a first-order approximation. The two sensitivities have different units, so their magnitudes should not be compared without considering parameter scales and uncertainty. Their signs and sizes nevertheless help explain how baseline and turnover adjustments can offset one another in a fitted response.
Reading Profiles and Diagnosing Weak Separation
A useful diagnostic is to plot the profiled objective against \(k_{\mathrm{out}}\), with the corresponding \(\widehat R_{\mathrm{base}}(k_{\mathrm{out}})\) shown on a second panel or as a companion curve. If the objective remains nearly unchanged while the estimated baseline moves substantially, the data weakly distinguish the parameter combinations. For a likelihood-based fit, the analogous display is the profile likelihood; uncertainty limits are obtained using an appropriate likelihood-based criterion and its assumptions, not by treating every low point as equally plausible.
Check whether the apparent profile shape is stable to reasonable choices of observation error model, starting values, and parameter bounds. A minimum located at a bound may indicate that the data favor values outside the explored range, or that the model is weakly identified. A narrow minimum can still be misleading if a major assumption—such as complete input suppression—is false. Inspect observed and predicted time courses and residuals as well as the profile.
The main alternatives answer somewhat different questions. Fixing baseline from an external, reliable measurement can reduce the trade-off, but treats that value as known and may understate uncertainty if it is noisy. Jointly estimating baseline and turnover uses the study data but must expose their dependence rather than hide it. Adding early and late observations, or measuring baseline repeatedly, can improve information when those observations target response scale and response rate, respectively. No fitting method can recover separate parameters if the experimental design contains no information to distinguish them.
Once \(R_{\mathrm{base}}\) and \(k_{\mathrm{out}}\) are estimated, the derived input \(k_{\mathrm{in}}=k_{\mathrm{out}}R_{\mathrm{base}}\) inherits uncertainty from both and their covariance. Treating the two estimates as independent can misrepresent uncertainty in input. Likewise, a precisely fitted baseline does not imply that the endogenous input is precise. The equilibrium relationship is a model-based conversion, not an additional independent observation.
Common Mistake and Practical Limits
A common mistake is to say, “The baseline concentration is measured, so baseline input is known; the remaining parameter is turnover.” That conclusion is only justified if turnover is known independently, because \(k_{\mathrm{in}}=k_{\mathrm{out}}R_{\mathrm{base}}\). The correction is to state what the data directly identify, profile the uncertain parameters, and report the resulting dependence or uncertainty.
Profiles are diagnostics, not repairs. They depend on the model, the error structure, and the parameters allowed to vary. A one-dimensional profile can miss confounding with other estimated parameters if those parameters were held fixed rather than re-optimized. Sparse sampling, uncertain perturbation magnitude, measurement error near baseline, or departures from steady state can all weaken the interpretation. In particular, a slow baseline change is not addressed by a constant-baseline model; the data and assumptions should be assessed before interpreting a turnover estimate mechanistically.
Check Your Understanding
Answer each question, then compare with the brief solutions.
- If \(R_{\mathrm{base}}=80\) nmol/L and \(k_{\mathrm{out}}=0.15\) day\(^{-1}\), what baseline input is implied?
- Why can a steady-state concentration alone not identify both \(k_{\mathrm{in}}\) and \(k_{\mathrm{out}}\)?
- In a profiled least-squares plot, what does a broad, nearly flat region indicate?
- For \(R(t)=R_{\mathrm{base}}e^{-k_{\mathrm{out}}t}\), which parameter is re-estimated at each fixed turnover rate when profiling?
- Name one assumption of the complete-suppression washout model whose violation could distort a turnover estimate.
Solutions. (1) \(k_{\mathrm{in}}=(80)(0.15)=12\) nmol/L/day. (2) The observation supplies only \(R_{\mathrm{base}}=k_{\mathrm{in}}/k_{\mathrm{out}}\), so multiple input–loss pairs give the same concentration. (3) A range of turnover values, paired with optimized baseline values, fits nearly as well; turnover is weakly separated by the data. (4) \(R_{\mathrm{base}}\) is re-estimated conditionally at each fixed \(k_{\mathrm{out}}\). (5) Examples include incomplete or gradual input suppression, a changing loss rate, or failure to begin at steady state.