Why classify the mechanism before interpreting the response?
An indirect response does not change instantaneously just because drug concentration changes. Instead, the drug alters the rate at which a measured response is produced or lost, and the response then changes as those rates become unbalanced. A rise in response, by itself, therefore does not identify the mechanism: it could follow increased production or decreased loss. Likewise, a fall could follow decreased production or increased loss.
The four canonical indirect effects are the combinations of two choices: the drug acts on production or loss, and it inhibits or stimulates that process. The drug-free equilibrium provides a common reference for predicting the initial response direction. These classifications describe the structure of the model; they do not, on their own, establish that a particular biological mechanism is true.
Baseline as the reference point
The drug-free baseline response, denoted \(R_0\), is the equilibrium value when no drug effect is present. At equilibrium the response is not changing, so its rate of change is zero:
Here \(k_{\mathrm{in}}\) has units of response per unit time and \(k_{\mathrm{out}}\) has units of inverse time, so their ratio has the units of response. The baseline equation is important for classification because it lets us evaluate the drug-induced rate of change at the moment treatment begins, when \(R=R_0\). If the rate is positive at that point, the response initially rises; if it is negative, the response initially falls.
For the equations below, let \(C\) denote drug concentration and let \(I(C)\) and \(S(C)\) denote dimensionless inhibition and stimulation effects. Assume these effects are zero when there is no drug effect. For ordinary fractional inhibition, \(0\leq I(C)<1\); this keeps the affected rate nonnegative. Assume \(S(C)\geq0\). These functions could be concentration-dependent in different ways; their exact shapes are not needed to classify the mechanisms.
The four canonical rate equations
Inhibition or stimulation of production changes the input rate while leaving the loss term unchanged. Inhibition or stimulation of loss changes the coefficient multiplying the response, while leaving production unchanged.
| Mechanism | Rate equation | Change at \(R=R_0\) | Initial direction |
|---|---|---|---|
| Inhibition of production | \(\dfrac{dR}{dt}=k_{\mathrm{in}}[1-I(C)]-k_{\mathrm{out}}R\) | \(-k_{\mathrm{in}}I(C)\) | Down |
| Inhibition of loss | \(\dfrac{dR}{dt}=k_{\mathrm{in}}-k_{\mathrm{out}}[1-I(C)]R\) | \(+k_{\mathrm{in}}I(C)\) | Up |
| Stimulation of production | \(\dfrac{dR}{dt}=k_{\mathrm{in}}[1+S(C)]-k_{\mathrm{out}}R\) | \(+k_{\mathrm{in}}S(C)\) | Up |
| Stimulation of loss | \(\dfrac{dR}{dt}=k_{\mathrm{in}}-k_{\mathrm{out}}[1+S(C)]R\) | \(-k_{\mathrm{in}}S(C)\) | Down |
The entries in the third column are obtained by substituting \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\) into each equation and simplifying. Thus, an effect on production and an effect on loss can produce the same initial direction: inhibition of loss and stimulation of production both initially raise the response, while inhibition of production and stimulation of loss both initially lower it. Direction alone cannot distinguish those pairs.
Worked Example: Classifying a Model From Its Equation
A hypothetical model is $$ \frac{dR}{dt}=k_{\mathrm{in}}-k_{\mathrm{out}}[1+S(C)]R. $$ Classify the drug effect and determine the direction of change from the drug-free baseline when \(S(C)>0\).
Step 1: Locate the modified rate. The drug-dependent factor multiplies \(k_{\mathrm{out}}R\), the loss term. Production remains \(k_{\mathrm{in}}\). The drug therefore acts on loss, not production.
Step 2: Determine whether the affected process is inhibited or stimulated. Since \(S(C)>0\), the multiplier \(1+S(C)\) exceeds 1. The drug increases the loss rate, so this is stimulation of response loss.
Step 3: Check the direction at baseline. At \(R=R_0\), use \(k_{\mathrm{in}}=k_{\mathrm{out}}R_0\): $$ \left.\frac{dR}{dt}\right|_{R=R_0} =k_{\mathrm{in}}-k_{\mathrm{out}}[1+S(C)]R_0 =k_{\mathrm{in}}-k_{\mathrm{in}}[1+S(C)] =-k_{\mathrm{in}}S(C)<0. $$ The response initially falls. The sign agrees with the interpretation that stimulated loss removes response faster than it is produced at the drug-free baseline.
Justifying the direction predictions
The prediction can be established with one substitution for each model. At the drug-free baseline, \(k_{\mathrm{out}}R_0=k_{\mathrm{in}}\). For production inhibition, for example,
Because \(k_{\mathrm{in}}>0\) and \(I(C)\geq0\), this rate is nonpositive, and it is strictly negative when the inhibition is active. The same substitution gives the other signs:
These signs describe the instantaneous direction at baseline when the drug effect is active. With constant concentration and constant effect, each model has a stable drug-on equilibrium toward which the response moves from baseline. For a time-varying concentration, the effect may change while the response is still adjusting, so a response trajectory need not continue in one direction at every later time. The baseline sign is a local prediction, not a guarantee of a globally monotone trajectory under any exposure history.
Worked Example: Predicting Direction With Hypothetical Values
Suppose a hypothetical response has \(k_{\mathrm{in}}=12\) response units per hour and \(k_{\mathrm{out}}=0.30\) per hour. A drug effect at the current concentration has magnitude \(0.25\), interpreted as \(I(C)=0.25\) for an inhibition or \(S(C)=0.25\) for a stimulation. Find the baseline and calculate the initial rate at baseline for each canonical mechanism.
First calculate the drug-free baseline: $$ R_0=\frac{k_{\mathrm{in}}}{k_{\mathrm{out}}} =\frac{12\ \text{response units/hour}}{0.30\ \text{per hour}} =40\ \text{response units}. $$
For inhibition of production, the initial rate is $$ -k_{\mathrm{in}}I(C)=-(12)(0.25)=-3.00\ \text{response units/hour}. $$ The response initially falls.
For inhibition of loss, the initial rate is $$ +k_{\mathrm{in}}I(C)=(12)(0.25)=+3.00\ \text{response units/hour}. $$ The response initially rises.
For stimulation of production, the initial rate is $$ +k_{\mathrm{in}}S(C)=(12)(0.25)=+3.00\ \text{response units/hour}. $$ The response initially rises.
For stimulation of loss, the initial rate is $$ -k_{\mathrm{in}}S(C)=-(12)(0.25)=-3.00\ \text{response units/hour}. $$ The response initially falls. These are rates of change at the instant the active effect is applied to a response at baseline, not predicted changes in response units after one hour. The equations would be needed to calculate the response at a later time.
Assumptions, diagnostics, and limits of classification
The simple sign predictions rely on a particular model setup. The response is at the drug-free equilibrium immediately before the effect begins; \(k_{\mathrm{in}}\) and \(k_{\mathrm{out}}\) are positive; and the drug changes only the specified production or loss term. The drug effect has the stated sign and is active. If the response is not at baseline when treatment begins, its initial direction also depends on how far it is from the current balance of production and loss.
A useful diagnostic is to write the equation as production minus loss before assigning a mechanism. Then check which term contains the drug-dependent factor. For example, an equation that reduces the input is production inhibition even if the observed response happens to rise over some interval because the response began well below baseline. Conversely, a decline is not sufficient evidence for production inhibition: stimulated loss predicts the same initial direction at baseline.
The canonical forms also assume that the drug effect is multiplicative on the selected rate. A more complicated model might change both production and loss, add a direct response term, or allow parameters to vary over time. Such a model is not one of the four single-mechanism canonical forms, and it should not be forced into one category merely because its response rises or falls. Model fit and biological interpretation require evidence beyond the direction of the response.
Check for practical warning signs: a proposed inhibition factor that makes a rate negative; an equation in which the purportedly unchanged process also depends on drug; a drug-free model that does not reproduce the stated baseline; or predicted baseline direction inconsistent with the sign of the modified rate. A fractional inhibition reaching 1 corresponds to complete inhibition of that rate, while values above 1 would imply a negative rate in these forms and violate the assumptions. A zero effect predicts no drug-induced departure from baseline in this model.
Worked Example: Separating Equation Classification From Observed Direction
In a hypothetical experiment, the response is below its drug-free baseline when a drug is administered. The candidate model is $$ \frac{dR}{dt}=k_{\mathrm{in}}[1-I(C)]-k_{\mathrm{out}}R, \qquad 0<I(C)<1. $$ Does the candidate equation necessarily predict that the response will fall immediately after administration?
The equation describes inhibition of production: the drug reduces the input, while the loss term remains unchanged. But the response is not at baseline, so substitute the actual \(R\), rather than automatically using the baseline sign: $$ \frac{dR}{dt}=k_{\mathrm{in}}[1-I(C)]-k_{\mathrm{out}}R. $$ Since \(R<R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\), the drug-free net rate \(k_{\mathrm{in}}-k_{\mathrm{out}}R\) is positive. Whether the drug reverses that positive rate depends on its magnitude: $$ \frac{dR}{dt} =\bigl(k_{\mathrm{in}}-k_{\mathrm{out}}R\bigr)-k_{\mathrm{in}}I(C). $$ The response falls only if \(k_{\mathrm{in}}I(C)>k_{\mathrm{in}}-k_{\mathrm{out}}R\). If the inhibition is smaller than this difference, the response still rises, although more slowly than it would without drug. Thus, the equation's mechanism is unambiguously production inhibition, but the observed direction away from baseline is not determined by the mechanism label alone.
How the canonical models differ from alternatives
A direct-response model relates drug concentration to the response itself, rather than making the drug alter a production or loss rate that then governs the response. Such a model may be suitable when the response follows the drug effect without a turnover process represented explicitly. The indirect models are useful when production and loss provide a plausible explanation for delayed response dynamics. An effect-compartment model, in contrast, represents delay between a measured concentration and an effect-site concentration; that is a different source of delay from response turnover.
Among the four indirect forms, the equation is the primary classifier, not the direction of a single observed change. A response increase is compatible with either inhibition of loss or stimulation of production; a decrease is compatible with either inhibition of production or stimulation of loss. Distinguishing these alternatives generally requires the time course and other relevant evidence, not just the sign at one observation. The canonical classification is also a simplification: it assumes one drug-modified process and does not establish that the chosen process is the true biological target.
Check Your Understanding
Answer each question using the rate equation and the drug-free baseline, not the response direction alone.
- For \(\frac{dR}{dt}=k_{\mathrm{in}}[1+S(C)]-k_{\mathrm{out}}R\), identify the mechanism and predict the initial direction at baseline when \(S(C)>0\).
- For \(\frac{dR}{dt}=k_{\mathrm{in}}-k_{\mathrm{out}}[1-I(C)]R\), identify the mechanism and calculate the initial rate at baseline in terms of \(k_{\mathrm{in}}\) and \(I(C)\).
- Why can an observed rise from baseline not distinguish inhibition of loss from stimulation of production?
- A production-inhibition model is applied when \(R<R_0\). Must the response fall immediately? State the condition that determines its direction.
- What is the drug-free baseline \(R_0\) for \(k_{\mathrm{in}}=9\) response units per hour and \(k_{\mathrm{out}}=0.15\) per hour?
Solutions. (1) The drug stimulates production because the factor \(1+S(C)\) multiplies \(k_{\mathrm{in}}\); at baseline the rate is \(+k_{\mathrm{in}}S(C)\), so the response initially rises. (2) This is inhibition of loss. Substitution of \(R_0=k_{\mathrm{in}}/k_{\mathrm{out}}\) gives \(+k_{\mathrm{in}}I(C)\), so the response initially rises when the inhibition is active. (3) Both mechanisms produce a positive rate at baseline, although one reduces loss and the other increases production. (4) No. Its rate is \((k_{\mathrm{in}}-k_{\mathrm{out}}R)-k_{\mathrm{in}}I(C)\); it falls only when \(k_{\mathrm{in}}I(C)>k_{\mathrm{in}}-k_{\mathrm{out}}R\). (5) \(R_0=9/0.15=60.0\) response units.