Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Deriving the Baseline Response at Equilibrium

Indirect Response Models · Tutorial 111 of 598

Deriving the Baseline Response at Equilibrium

Learn how an untreated equilibrium baseline constrains the ratio of response production to response loss, and how to check whether that interpretation is justified.

Advanced 10 min read

What You'll Learn

  • Derive the equilibrium baseline from the untreated response state equation.
  • Check the conditions needed for a unique, stable baseline.
  • Interpret the units of the production-to-loss rate ratio.
  • Explain why one baseline value cannot identify both rates separately.
  • Recognize when observations may not represent an untreated equilibrium.
  • Compare mechanistic baseline estimation with a direct empirical baseline.

Why the Baseline Is a Rate Balance

In an indirect response model, the response is maintained by production and reduced by loss. Even when no drug is present, these processes continue. The untreated baseline is therefore not a response with no underlying activity; it is the value at which the ongoing production and loss balance. This balance gives baseline observations a mechanistic interpretation, while also limiting what those observations can tell us about the separate rates.

The earlier material on state equations for indirect response models establishes the production-and-loss description. Here we derive its untreated equilibrium and examine what that equation implies for interpreting baseline data. The derivation is conditional: it applies when production and loss follow the stated model and the untreated response has reached equilibrium.

Definition: Let \(R(t)\) denote the response at time \(t\). In the untreated production-and-loss model, \(k_{\mathrm{in}}\) is the constant response production rate and \(k_{\mathrm{out}}\) is the first-order response loss rate constant. The untreated equilibrium baseline is \(R_0\), the constant response value at which the rate of change is zero. Under this model, \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\), provided \(k_{\mathrm{out}}>0\).

The corresponding untreated state equation is

$$ \frac{dR(t)}{dt}=k_{\mathrm{in}}-k_{\mathrm{out}}R(t). $$

The production term \(k_{\mathrm{in}}\) adds response per unit time. The loss term \(k_{\mathrm{out}}R(t)\) removes response at a rate proportional to the amount currently present. Thus \(k_{\mathrm{out}}\) has units of inverse time, while \(k_{\mathrm{in}}\) has units of response per time. Their ratio has units of response, as a baseline response must.

Deriving the Equilibrium Baseline

An equilibrium is a state that does not change over time. For the untreated response, this requires \(dR(t)/dt=0\). Setting the state equation's right-hand side to zero gives the rate-balance equation:

$$ 0=k_{\mathrm{in}}-k_{\mathrm{out}}R_0. $$

If \(k_{\mathrm{out}}>0\), rearranging yields the unique equilibrium:

$$ R_0=\frac{k_{\mathrm{in}}}{k_{\mathrm{out}}}. $$

This result does not say that either rate is zero at baseline. Instead, at equilibrium the production rate equals the loss rate:

$$ k_{\mathrm{in}}=k_{\mathrm{out}}R_0. $$

For example, if \(R_0\) is 10 response units and \(k_{\mathrm{out}}\) is \(0.2\ \mathrm{h}^{-1}\), the implied production rate is \(2\) response units per hour. Both processes are active, but their net contribution to the response is zero.

The equilibrium is stable under the model when \(k_{\mathrm{out}}>0\). To see why, define the deviation from equilibrium as \(D(t)=R(t)-R_0\). Subtract the equilibrium balance \(k_{\mathrm{in}}=k_{\mathrm{out}}R_0\) from the state equation:

$$ \frac{dD(t)}{dt} = k_{\mathrm{in}}-k_{\mathrm{out}}R(t) -\bigl(k_{\mathrm{in}}-k_{\mathrm{out}}R_0\bigr) = -k_{\mathrm{out}}\bigl(R(t)-R_0\bigr) = -k_{\mathrm{out}}D(t). $$

The solution is \(D(t)=D(0)e^{-k_{\mathrm{out}}t}\). Because \(e^{-k_{\mathrm{out}}t}\) decreases toward zero when \(k_{\mathrm{out}}>0\), deviations diminish and the response approaches \(R_0\). If \(k_{\mathrm{out}}=0\), a positive constant production rate has no finite equilibrium; if \(k_{\mathrm{out}}<0\), deviations do not decay. These cases do not support the stable baseline interpretation used here.

Worked Example: Inferring Production From a Baseline

Suppose a hypothetical biomarker has an untreated baseline of \(12.0\ \mathrm{ng/mL}\), and an independently specified loss rate constant of \(0.30\ \mathrm{h}^{-1}\). Under the constant production, first-order loss model, calculate the production rate required to maintain that baseline.

At equilibrium, \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\). Rearranging gives \(k_{\mathrm{in}}=k_{\mathrm{out}}R_0\). Substituting the values:

$$ k_{\mathrm{in}} =(0.30\ \mathrm{h}^{-1})(12.0\ \mathrm{ng/mL}) =3.60\ \mathrm{ng/mL/h}. $$

Thus, the model requires production of \(3.60\ \mathrm{ng/mL}\) per hour to balance loss at a baseline of \(12.0\ \mathrm{ng/mL}\). The units agree: inverse hours multiplied by response units gives response units per hour.

What Baseline Observations Constrain

The equilibrium equation directly constrains a ratio, not both rates separately. Rearranging it gives

$$ \frac{k_{\mathrm{in}}}{k_{\mathrm{out}}}=R_0. $$

If a baseline is known exactly, any positive pair of rates with this ratio gives the same equilibrium value. Multiplying both rates by the same positive constant changes the speed at which deviations return toward baseline, but leaves the ratio unchanged. A single equilibrium observation therefore cannot, by itself, determine the two rates separately.

Worked Example: Different Rates, the Same Baseline

Consider a hypothetical untreated response with a baseline of \(80\) response units. Compare two possible parameter pairs under the same model.

For the first pair, let \(k_{\mathrm{out}}=0.080\ \mathrm{h}^{-1}\). The production rate that gives the specified baseline is

$$ k_{\mathrm{in}}=R_0k_{\mathrm{out}} =(80\ \mathrm{units})(0.080\ \mathrm{h}^{-1}) =6.40\ \mathrm{units/h}. $$

For the second pair, let \(k_{\mathrm{out}}=0.200\ \mathrm{h}^{-1}\). Then

$$ k_{\mathrm{in}}=(80\ \mathrm{units})(0.200\ \mathrm{h}^{-1}) =16.0\ \mathrm{units/h}. $$

Both pairs give \(k_{\mathrm{in}}/k_{\mathrm{out}}=80\) response units. The baseline alone cannot distinguish these pairs. Information about how the response changes over time, or other independent information about a rate, would be needed to distinguish their dynamics.

With measured rather than exact baselines, the inference is less direct. Repeated untreated observations can help estimate a typical baseline, but their spread may reflect measurement error, within-person fluctuation, or genuine change over time. A sample average is not automatically the equilibrium value: it is a useful estimate of \(R_0\) only to the extent that the observations represent the same stable untreated condition and the measurement process is adequately characterized.

Worked Example: Checking a Proposed Equilibrium

In a hypothetical case, the response is \(10.0\) units, production is \(2.40\ \mathrm{units/h}\), and the loss rate constant is \(0.20\ \mathrm{h}^{-1}\). Determine whether this response is at equilibrium under the model, and calculate the implied equilibrium baseline.

At the observed response, the loss rate is

$$ k_{\mathrm{out}}R =(0.20\ \mathrm{h}^{-1})(10.0\ \mathrm{units}) =2.00\ \mathrm{units/h}. $$

The net rate of change is therefore

$$ \frac{dR}{dt} =2.40-2.00 =0.40\ \mathrm{units/h}. $$

Because this rate is not zero, \(10.0\) units is not the equilibrium value under the specified rates. The model-implied baseline is

$$ R_0=\frac{2.40\ \mathrm{units/h}}{0.20\ \mathrm{h}^{-1}} =12.0\ \mathrm{units}. $$

The response is below that equilibrium, and the positive net rate means it is increasing at that instant. This calculation checks the proposed equilibrium condition rather than treating any untreated observation as a baseline by definition.

Assumptions and Practical Checks

The ratio \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\) is valid when the following conditions are reasonable over the baseline period:

  • The response follows constant production and first-order loss, so its rate is \(k_{\mathrm{in}}-k_{\mathrm{out}}R(t)\).
  • The production and loss rates are stable over the period being interpreted; important time-varying influences are absent or accounted for.
  • The observations represent an untreated condition, without a relevant prior intervention still affecting the response.
  • The response is at, or sufficiently close to, equilibrium that its mean rate of change is approximately zero.
  • \(k_{\mathrm{out}}>0\), giving a finite and stable equilibrium for nonnegative production.

These are modeling assumptions, not facts established merely because a value is labelled “baseline.” Useful checks include plotting serial untreated observations against time, checking for persistent upward or downward trends, comparing measurements across relevant times or conditions, and examining whether replicate measurements vary more than expected from the observation process. A marked trend, systematic cycling, or abrupt shift warns that a single constant \(R_0\) may not describe the data.

A baseline taken soon after a prior intervention may also fail to represent the untreated equilibrium. In that situation, the response can still be moving even if no intervention is being administered at the time of measurement. The relevant question is whether the system has returned to the stable untreated condition, not simply whether treatment is currently absent.

Mechanistic Baseline Versus Direct Estimation

There are two useful but different approaches. A direct empirical approach estimates baseline from untreated response measurements, for example by summarizing stable repeated observations. A mechanistic approach uses \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\) to relate that baseline to production and loss. The direct approach can provide a clear data-based summary without requiring the rate model to be correct. The mechanistic approach makes the rate balance explicit and can connect baseline information to a dynamic model, but depends on the stated assumptions.

The approaches need not compete. An empirical baseline estimate can be compared with the ratio implied by proposed rates. Disagreement may indicate parameter inconsistency, measurement or sampling issues, failure to reach equilibrium, or a misspecified production-and-loss structure. Agreement is a useful consistency check, but does not prove that the two rates are separately correct: many rate pairs can have the same ratio.

The equation also has limits. If production varies with time, loss is not proportional to response, or other processes contribute to the state equation, the equilibrium must be derived from that expanded model's zero-rate condition. The simple ratio should not be carried over automatically. Even within the simple model, uncertain baseline observations constrain the ratio with uncertainty rather than fixing it exactly.

Common Mistake

Wrong reasoning: “The baseline is \(R_0\), so the production rate must equal \(R_0\), and the loss rate must be zero.”

Correction: \(R_0\) is a response amount, while \(k_{\mathrm{in}}\) is a production rate and \(k_{\mathrm{out}}\) is a loss rate constant. Their units differ. At equilibrium, the production rate equals the loss rate evaluated at the baseline: \(k_{\mathrm{in}}=k_{\mathrm{out}}R_0\). Neither process must be zero.

Key takeaway: Under constant production, first-order loss, and untreated equilibrium, the baseline is \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\). Baseline data constrain this rate ratio; they do not identify both rates separately.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A response has an untreated equilibrium of \(15\) units and \(k_{\mathrm{out}}=0.10\ \mathrm{h}^{-1}\). What production rate is required?
  2. Why can a single exact equilibrium baseline not identify \(k_{\mathrm{in}}\) and \(k_{\mathrm{out}}\) separately?
  3. Under the model, \(k_{\mathrm{in}}=4.0\ \mathrm{units/h}\), \(k_{\mathrm{out}}=0.50\ \mathrm{h}^{-1}\), and \(R=6.0\) units. Is the response at equilibrium? What is its instantaneous rate of change?
  4. Name two features of baseline data that would make the equilibrium interpretation questionable.

Solutions. (1) \(k_{\mathrm{in}}=R_0k_{\mathrm{out}}=(15\ \mathrm{units})(0.10\ \mathrm{h}^{-1})=1.50\ \mathrm{units/h}\). (2) The observation fixes only the ratio \(k_{\mathrm{in}}/k_{\mathrm{out}}=R_0\); multiplying both rates by the same positive factor preserves that ratio. (3) The net rate is \(4.0-(0.50)(6.0)=1.0\ \mathrm{units/h}\), so the response is not at equilibrium and is increasing at that instant. (4) Examples include a persistent time trend, systematic cycling, an abrupt shift, or measurements taken before return to the untreated stable condition.

Clinical Biostats