Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Stimulation of Response Loss

Indirect Response Models · Tutorial 116 of 598

Stimulation of Response Loss

Learn how increasing the loss rate, with production held constant, drives an indirect response downward and how to calculate its initial slope and steady state.

Advanced 10 min read

What You'll Learn

  • Write the stimulated-loss model while keeping the production rate constant.
  • Derive the baseline and drug-on steady-state responses.
  • Determine the initial response direction from the starting response and stimulated loss rate.
  • Calculate the transient response under a constant drug effect, including the equilibrium initial condition.
  • Distinguish stimulated loss from changes in production and identify assumptions that can fail.

How Increased Loss Changes an Indirect Response

A drug can lower a response without directly subtracting from it. If the response is continually produced and removed, increasing the removal rate can lower the amount that accumulates. This is the defining mechanism considered here: drug stimulates response loss while the production rate remains constant. Examples of response variables might include a biomarker concentration or a physiological quantity, but the model itself does not establish which biological process the drug acts on.

The key comparison is between the response already present and the rate at which the drug-modified system removes it. At baseline equilibrium, production and loss balance. After drug increases the loss rate, the existing response is initially removed faster than it is replaced, so the response begins to decline. If the drug effect is sustained, the response approaches a lower equilibrium.

Definition: Let \(R(t)\) be the response, \(k_{in}\) the constant production rate, and \(k_{out}\) the baseline first-order loss-rate constant. Let \(C(t)\) be drug concentration and \(S(C)\) a dimensionless, nonnegative stimulation function. The stimulated-loss model is
$$ \frac{dR}{dt}=k_{in}-k_{out}\bigl[1+S(C(t))\bigr]R(t). $$
The drug leaves production at \(k_{in}\) and multiplies the loss-rate constant by \(1+S(C(t))\). Thus \(S=0\) means no stimulation, while \(S>0\) means increased loss.

The units must be compatible: if \(R\) is measured in response units, then \(k_{in}\) has units of response units per time, and \(k_{out}\) has units of inverse time. The stimulation \(S\) is dimensionless. The term \(k_{out}[1+S]R\) therefore has units of response units per time, matching the production term.

Baseline and Drug-On Equilibria

Before drug exposure, suppose the response is at its baseline equilibrium and there is no drug stimulation. Setting the derivative to zero gives \(k_{in}=k_{out}R_0\), where \(R_0\) denotes the baseline response. Therefore, \(R_0=k_{in}/k_{out}\). This balance is the starting point for predicting what stimulation of loss will do.

Now suppose the drug effect is constant, \(S(C(t))=S\), with \(S>0\). A steady state \(R_{ss}\) occurs when the derivative is zero. Production must then equal stimulated loss:

$$ 0=k_{in}-k_{out}(1+S)R_{ss}, \qquad R_{ss}=\frac{k_{in}}{k_{out}(1+S)} =\frac{R_0}{1+S}. $$

Because \(1+S>1\), the drug-on steady state is below baseline. The fractional decrease from baseline is \(1-R_{ss}/R_0=S/(1+S)\). This is a steady-state comparison, not a statement that the response falls immediately to \(R_{ss}\): the response changes over time because the response pool has to turn over.

Worked Example: Predicting the Drug-On Steady State

Suppose a hypothetical response has \(k_{in}=12\) response units per hour and \(k_{out}=0.20\) per hour. At baseline, the response is at equilibrium. A sustained drug effect produces \(S=0.50\). Find baseline and drug-on steady-state responses.

First calculate baseline:

$$ R_0=\frac{k_{in}}{k_{out}} =\frac{12\ \text{response units/hour}}{0.20\ \text{hour}^{-1}} =60\ \text{response units}. $$

With stimulation, the loss-rate constant is \(0.20(1+0.50)=0.30\) per hour. Hence

$$ R_{ss}=\frac{12\ \text{response units/hour}}{0.30\ \text{hour}^{-1}} =40\ \text{response units}. $$

The sustained stimulation therefore lowers the equilibrium by \(60-40=20\) response units, or \(20/60=0.3333\), a 33.3% decrease. The units cancel correctly in each rate ratio, leaving response units.

Initial Direction and the Transient Response

The initial direction follows directly from the differential equation. At any starting value \(R(0)\), under a constant effect \(S\), the initial slope is \(k_{in}-k_{out}(1+S)R(0)\). It is negative if the starting response exceeds the drug-on steady state, positive if it is below that steady state, and zero if it equals that steady state.

In the common case that dosing begins when the system is at baseline, substitute \(R(0)=R_0=k_{in}/k_{out}\). The initial slope becomes \(-k_{out}S R_0=-k_{in}S\), which is negative whenever \(S>0\). This result gives the direction at onset; it does not imply a constant downward slope throughout the response trajectory.

Key result: With constant production and constant stimulation \(S>0\), a response starting at baseline initially declines at rate \(k_{in}S\). Its drug-on equilibrium is \(R_0/(1+S)\). More generally, the initial direction depends on whether \(R(0)\) is above or below that drug-on equilibrium.

For a sustained constant effect, the transient can also be obtained from the first-order differential equation. Write it as \(dR/dt=-k_{out}(1+S)(R-R_{ss})\). When \(R(0)\ne R_{ss}\), separating variables and integrating from time zero to \(t\) gives

$$ \int_{R(0)}^{R(t)}\frac{dR}{R-R_{ss}} =-k_{out}(1+S)\int_0^t du, $$

so that \(\ln|R(t)-R_{ss}|-\ln|R(0)-R_{ss}|=-k_{out}(1+S)t\). Exponentiating and retaining the initial difference's sign yields

$$ R(t)=R_{ss}+\bigl[R(0)-R_{ss}\bigr]e^{-k_{out}(1+S)t}. $$

The derivation by separation divides by \(R-R_{ss}\), so it requires \(R(0)\ne R_{ss}\). If instead \(R(0)=R_{ss}\), the differential equation gives \(dR/dt=0\), and the response remains constant at \(R_{ss}\). The displayed transient formula also gives this constant solution when \(R(0)=R_{ss}\), so it holds in both cases. The exponential factor decays to zero because \(k_{out}(1+S)>0\); consequently, the response approaches \(R_{ss}\) from its starting value.

Worked Example: Calculating the Initial Slope and a Later Response

Use \(k_{in}=12\) response units per hour, \(k_{out}=0.20\) per hour, and \(S=0.50\), with the response initially at baseline. Find the initial slope and the model-predicted response after 3 hours of sustained stimulation.

The baseline is \(R_0=12/0.20=60\) response units, and the drug-on steady state is \(R_{ss}=40\) response units. The initial slope at baseline is

$$ \left.\frac{dR}{dt}\right|_{t=0} =k_{in}-k_{out}(1+S)R_0 =12-(0.20)(1.50)(60) =-6.00\ \text{response units/hour}. $$

The negative slope confirms an initial decline. The effective rate constant is \(0.20(1.50)=0.30\) per hour. Using the transient expression,

$$ R(3)=40+(60-40)e^{-(0.30)(3)} =40+20e^{-0.90} \approx 48.13\ \text{response units}. $$

After 3 hours, the response is predicted to be about 48.13 response units: below its baseline value of 60, but still above its eventual drug-on equilibrium of 40. The response has not reached steady state at that time.

Worked Example: Direction From a Nonbaseline Starting Value

A hypothetical system has \(k_{in}=10\) response units per day, \(k_{out}=0.25\) per day, and constant stimulation \(S=1.0\). The response is 15 units when the drug effect becomes constant. Determine the drug-on steady state and the initial direction.

The stimulated loss-rate constant is \(0.25(1+1.0)=0.50\) per day, so

$$ R_{ss}=\frac{10\ \text{response units/day}}{0.50\ \text{day}^{-1}} =20\ \text{response units}. $$

Since \(R(0)=15\) is below \(R_{ss}=20\), the initial slope is positive. Direct substitution verifies this:

$$ \left.\frac{dR}{dt}\right|_{t=0} =10-(0.50)(15) =2.50\ \text{response units/day}. $$

Thus stimulated loss does not guarantee a decline from every possible starting value. Here the drug-on equilibrium exceeds the starting response, so the response initially rises toward it. The decrease relative to baseline applies when the system starts at baseline and stimulation is positive.

Assumptions, Diagnostics, and Model Comparisons

The direction and equilibrium conclusions require specific assumptions. Production is constant and unaffected by drug; loss is first-order in \(R\); the drug increases the loss-rate constant through a nonnegative \(S\); and a drug-on steady state is relevant only if the effect is held constant long enough for the response to approach it. The transient calculation additionally assumes a constant stimulation over the interval. If concentration changes, \(S(C(t))\) changes too, and the constant-effect exponential expression cannot simply be used with one fixed \(S\).

The observed response may not follow the model even when it initially falls. Useful warning signs include a response that continues to drift before dosing, a drug-on plateau inconsistent with \(R_0/(1+S)\), or systematic residual patterns suggesting that production also changes or loss is not first-order. A measured concentration-effect delay can also make the observed response lag the assumed effect; the concentration and effect history should not be treated as interchangeable without justification.

A useful conceptual comparison is with stimulated production. In that mechanism, drug raises the input while the loss process remains unchanged, so the drug-on equilibrium moves above baseline. In stimulated loss, the input is unchanged and the drug-on equilibrium moves below baseline. Both mechanisms can produce delayed changes because the response pool turns over. A falling response alone does not prove that the drug stimulates loss: reduced production or another mechanism could also lower the response. Discriminating mechanisms requires the model to be evaluated against the exposure and response data, not inferred from direction alone.

The model is intentionally simple. It does not specify the shape of \(S(C)\), establish a biological mechanism, or guarantee that response measurements represent a single well-mixed pool. If production varies over time, the baseline is not at equilibrium, the response has multiple turnover processes, or the effect is delayed, the simple predictions may be inadequate. In those situations, a more appropriate model should represent the relevant changing input, delay, or response structure rather than treating the mismatch as evidence for a different value of \(k_{out}\).

Common Mistakes

  • Changing production as well as loss. The stimulated-loss model changes \(k_{out}\) by the factor \(1+S\); \(k_{in}\) remains fixed. Altering both rates describes a different model.
  • Assuming the response always falls. A positive stimulation makes the drug-on equilibrium lower than baseline, but a response starting below that equilibrium initially rises. Check the starting value against \(R_{ss}\).
  • Equating initial slope with long-term change. The initial slope is a rate at onset. The eventual response is an equilibrium reached over time, not the response after an arbitrarily short interval.
  • Dividing by zero in the transient derivation. Separation using \(1/(R-R_{ss})\) applies only when \(R(0)\ne R_{ss}\). When they are equal, the solution is constant; the final formula still covers this case.
  • Reading the model as proof of mechanism. A decrease is compatible with stimulated loss but does not establish it. Alternative explanations must be assessed against the data and design.
Key takeaway: With production held constant, stimulation of loss raises the first-order loss-rate constant. From baseline, any positive stimulation produces a negative initial slope and a lower sustained equilibrium, \(R_0/(1+S)\). The initial direction from another starting point depends on its position relative to that equilibrium.

Check Your Understanding

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

  1. If \(k_{in}=8\) response units per hour and \(k_{out}=0.40\) per hour, what is the baseline equilibrium?
  2. For the same rates and \(S=0.25\), calculate the drug-on equilibrium and state whether it is above or below baseline.
  3. With \(k_{in}=8\), \(k_{out}=0.40\), \(S=0.25\), and \(R(0)=15\) response units, determine the initial slope and direction.
  4. Why does the separation-of-variables derivation need a separate argument when \(R(0)=R_{ss}\)? What is the response in that case?
  5. Name two assumptions that would make the constant-effect transient expression inappropriate.

Solutions. (1) \(R_0=k_{in}/k_{out}=8/0.40=20\) response units. (2) \(R_{ss}=8/[0.40(1.25)]=16\) response units, below baseline. (3) The initial slope is \(8-[0.40(1.25)(15)]=8-7.50=0.50\) response units per hour, so the response initially rises toward 16. (4) Separation divides by \(R-R_{ss}\), which is zero at the initial value. Substituting \(R=R_{ss}\) into the differential equation gives zero slope, so the response stays constant at \(R_{ss}\). (5) Examples include time-varying stimulation during the interval, drug effects on production, or loss that is not first-order in the response.

Clinical Biostats