Why Similar Response Curves Can Imply Different Parameters
In an indirect response model, a drug changes a production or loss rate, while the measured response changes over time through turnover. Two parameters can therefore influence the same observed trajectory in different ways: the drug-effect amplitude sets how far the response is driven from baseline, and the turnover rate sets how quickly it moves. If observations cover only a narrow part of the trajectory, a larger amplitude paired with slower turnover can resemble a smaller amplitude paired with faster turnover.
This is a parameter-confounding problem, not evidence that the parameters are biologically interchangeable. The distinction matters when interpreting an estimated Emax or predicting how quickly a response develops or recovers. The examples below use a stimulatory effect on production; an inhibitory effect has the same amplitude-versus-rate issue, with the response moving below baseline instead of above it.
Model and Assumptions
Let \(R(t)\) be the response, measured in response units, and let \(k_{\mathrm{in}}\) and \(k_{\mathrm{out}}\) be the baseline production and loss-rate constants. Under baseline conditions, the response is at steady state. During exposure, assume drug concentration \(C\) is constant and the fractional stimulation of production follows an Emax function. For the purpose of isolating the confounding, take \(C\) and \(EC_{50}\) as known, so the concentration-dependent factor \(f(C)\) is known.
This model assumes a constant baseline state before dosing, a known and constant concentration during the interval being analyzed, and a constant effect amplitude while exposure is maintained. It also assumes that the measured response follows the stated turnover equation, without an additional delay or a separate drug effect on loss. These assumptions matter: if concentration changes appreciably, the input is not constant; if the drug also changes loss, the model’s interpretation of \(E_{\max}\) may be wrong.
Where the Confounding Comes From
Suppose exposure begins at time zero and \(y(0)=1\). Solving the constant-exposure equation gives
The long-term response is \(1+q\), so \(q\) controls the plateau. The rate constant controls how quickly that plateau is approached. At early times, however, the trajectory is governed mainly by their product. The initial slope is
For small \(k_{\mathrm{out}}t\), the expansion \(1-e^{-k_{\mathrm{out}}t}\approx k_{\mathrm{out}}t\) gives \(y(t)-1\approx qk_{\mathrm{out}}t\). Thus an early slope alone identifies a combination of amplitude and turnover, not each parameter separately.
The formula also shows the model’s limits. For an increasing response under this model, \(0<\Delta<q\) at a finite positive time. If an observed change is at or above the proposed plateau change, the proposed parameters and assumptions cannot explain that observation exactly. In practice, measurement error can produce small apparent violations, but systematic disagreement calls for checking the exposure profile, model structure, and data.
Worked Example: One Observation, Two Plausible Parameter Pairs
Suppose two hypothetical parameter sets are considered, with concentration such that the known factor is \(f(C)=0.5\). Set A has \(E_{\max}=1.0\) and \(k_{\mathrm{out}}=0.200\ \mathrm{h}^{-1}\), so \(q=0.500\). At two hours, its predicted response change is
Now keep \(t=2\ \mathrm{h}\) and \(k_{\mathrm{out}}=0.100\ \mathrm{h}^{-1}\). To give the same response change, set
Both parameter sets predict \(y(2\ \mathrm{h})\approx1.1648\). They do not predict the same later response. At eight hours, Set A gives \(y\approx1+0.500(1-e^{-1.6})=1.3991\), whereas Set B gives \(y\approx1+0.9094(1-e^{-0.8})=1.5008\). A lone two-hour observation cannot distinguish these pairs; later observations can be informative because their predicted trajectories separate. All values here are hypothetical and the response is normalized to baseline.
What Baseline and Post-Dose Data Contribute
Stable baseline measurements establish the response level \(R_0\), help quantify its measurement variability, and permit normalization. They do not, by themselves, identify \(k_{\mathrm{out}}\) or \(k_{\mathrm{in}}\): the steady-state relationship is \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\), so many production-loss pairs have the same baseline. If \(k_{\mathrm{out}}\) is established from informative time-course data, then \(k_{\mathrm{in}}=R_0k_{\mathrm{out}}\) follows under the model. A precise baseline is useful, but it is not a substitute for observing response dynamics.
Post-dose data can separate the parameters when they capture distinct features of the trajectory. A maintained exposure with observations approaching a plateau informs \(q\); observations during the transition inform \(k_{\mathrm{out}}\). Alternatively, after the drug effect has ended and the response returns toward baseline, the recovery rate informs \(k_{\mathrm{out}}\), while the amount of response displacement at the start of recovery depends on both parameters.
Worked Example: Combining a Plateau and a Transient Observation
Consider hypothetical normalized-response data under constant exposure: the plateau is \(1.400\), and at \(t=4.00\ \mathrm{h}\), \(y=1.2203\). Under the model, the plateau gives \(q=1.400-1=0.400\). The fraction of the plateau change reached by four hours is
The step-response equation therefore gives
The plateau supplies amplitude information; the fraction of that plateau change reached by a specified time supplies rate information. For example, if the baseline were \(100\ \mathrm{ng/mL}\), a normalized plateau of \(1.400\) would correspond to \(140\ \mathrm{ng/mL}\). Once the turnover rate is supported by the time course, the model implies \(k_{\mathrm{in}}=(100\ \mathrm{ng/mL})(0.2000\ \mathrm{h}^{-1})=20.00\ \mathrm{ng/mL/h}\). Baseline alone would not establish that rate.
Worked Example: Using Recovery to Estimate Turnover
Suppose hypothetical exposure is maintained for six hours, then the drug effect ends. With \(q=0.400\) and \(k_{\mathrm{out}}=0.200\ \mathrm{h}^{-1}\), the response change at the end of exposure is
After the effect ends, assume the system returns to baseline with no residual drug effect. For elapsed recovery time \(\tau\), the excess response is \(y(6+\tau)-1=[y(6)-1]e^{-k_{\mathrm{out}}\tau}\). Suppose the measured hypothetical excesses are \(0.1875\) at two hours after the end of exposure and \(0.1256\) at four hours. Their ratio removes the unknown starting excess:
The recovery data identify the decay rate under the stated assumptions without requiring the response at the end of dosing to equal the plateau. The recovery level at a particular time still depends on the displacement at the start of recovery, which carries information about the drug effect and the prior dosing duration.
Design Choices, Diagnostics, and Alternatives
The most useful design depends on what is feasible and what the data actually cover. A long, stable exposure can reveal the plateau, but may be impractical or unsafe. A post-exposure recovery phase can inform turnover, provided drug effect has ended and there is enough measurable response displacement. Multiple observations during onset and recovery provide more shape information than one early observation. When the concentration is changing, the constant-exposure expression is not appropriate; interpretation should use the time-varying exposure in the model rather than treating it as a step.
| Data feature | Most direct information | Important limitation |
|---|---|---|
| Stable baseline level | \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\) | Does not separate production from loss |
| Early onset observations | Approximately \(qk_{\mathrm{out}}\), through the initial slope | Amplitude and rate can trade off |
| Exposure plateau | Effect amplitude \(q\) | Requires sufficiently sustained exposure and a credible plateau |
| Onset plus recovery time course | Trajectory shape and turnover rate | Depends on the assumed mechanism and drug-effect offset |
Warning signs include highly different fitted parameter pairs that predict almost the same observed values, large uncertainty in both parameters, or estimates that change substantially when early or late observations are omitted. A fitted curve may appear adequate over the observed interval yet imply very different plateaus or recovery times. Check whether the data include observations that inform both amplitude and rate, whether exposure is correctly represented, and whether the baseline is plausibly stable. A model that systematically misses the observed plateau or recovery should not be rescued simply by interpreting its parameter estimates more strongly.
The main alternative to relying on a full onset-and-recovery profile is to obtain one parameter from an independent source or a separate informative experiment. For example, a turnover rate supported by recovery data can reduce the amplitude-rate trade-off when analyzing onset data. Conversely, a well-supported plateau can anchor amplitude while the transition informs rate. Neither strategy removes the need to check that the model assumptions apply. If the observation window contains only early responses and no independent information, reporting the combined slope information may be more defensible than treating both parameters as precisely determined.
Common Mistake
The same caution applies in the opposite direction: a large early slope does not by itself prove a large Emax. The concentration-dependent factor must also be accounted for, and the turnover rate contributes to the slope. Parameter interpretation should follow from the entire informative trajectory and the model assumptions, not from one visual feature.
Check Your Understanding
Answer each question using the model assumptions stated above. Short solutions follow the questions.
- Why can a single early response slope fail to distinguish \(q\) from \(k_{\mathrm{out}}\)?
- A normalized plateau is \(1.30\). What is \(q\) under the stimulatory model?
- If the baseline response is \(80\ \mathrm{mg/L}\), does that value alone determine \(k_{\mathrm{out}}\)? Explain.
- What post-exposure pattern can inform \(k_{\mathrm{out}}\) without knowing the recovery starting displacement exactly?
Solutions. (1) The initial slope is \(qk_{\mathrm{out}}\), so different values of the two parameters can have the same product. (2) \(q=1.30-1=0.300\). (3) No. Baseline provides \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}=80\ \mathrm{mg/L}\), a ratio; it does not give either rate separately. (4) Two or more positive excess-response measurements during recovery can be compared: their ratio cancels the initial excess, and the log ratio divided by the elapsed-time difference estimates \(k_{\mathrm{out}}\), if the drug effect has ended and recovery follows the assumed turnover model.