Detecting Tolerance During Repeated Opioid Dosing

Why Repeated Dosing Can Change the Exposure–Effect Relationship

After a single dose, a concentration–analgesia relationship can summarize how response changes with exposure. Under repeated dosing, however, the same concentration may be associated with less analgesia later than earlier. This pattern is consistent with tolerance: a reduction in drug effect during continued or repeated exposure. It is not, by itself, proof of tolerance. Concentration measurement, delayed equilibration, changing baseline symptoms, and other time-related processes can also produce an apparent decline in response.

This case uses a deliberately compact, phenomenological model. It is intended to describe an exposure-linked adaptation pattern, not to assert a particular biological mechanism or guide opioid treatment. Concentrations and analgesia values below are hypothetical. The model’s main diagnostic implication is testable: if the tolerance component is useful, it should reduce systematic prediction errors over time without creating a new pattern in the residuals.

Definition: Let $C(t)$ be opioid concentration at time $t$, $E(t)$ be analgesia on a numerical scale, and $T(t)$ be a dimensionless tolerance state. Define the concentration stimulus $S(C)=C/(EC_{50}+C)$, where $EC_{50}>0$ is the concentration giving half of the maximum stimulus. The model is $$ \frac{dT}{dt}=k_{\mathrm{tol}}\{S(C(t))-T(t)\}, \qquad E(t)=E_0+\frac{E_{\max}S(C(t))}{1+T(t)}. $$ Here $k_{\mathrm{tol}}>0$ has units of inverse time, $E_0$ is baseline analgesia, and $E_{\max}$ is a maximum effect increment in the units of $E$. The tolerance state moves toward the current concentration stimulus, and a larger $T$ attenuates the concentration-driven effect through the denominator $1+T$.

The model assumes that concentration is measured or adequately predicted, that the baseline $E_0$ is stable over the interval being analyzed, and that the chosen analgesia scale can reasonably be represented by this continuous response equation. It also assumes a single adaptation rate for development and recovery, and an inhibitory effect of tolerance that follows the specified denominator. These are modeling assumptions to check, not established facts about an individual patient.

Because $0\leq S(C)<1$ for nonnegative concentration, a nonnegative initial $T$ remains nonnegative under this differential equation. At zero concentration, $S(C)=0$, so $T$ decays toward zero and the drug-driven part of the predicted effect vanishes. The model therefore represents reversible adaptation, but it does not permit a separate onset and recovery rate or multiple tolerance processes.

Response to a Constant Concentration

If concentration is held at a constant value $C^\ast$, then $S(C(t))=S^\ast=S(C^\ast)$ and the tolerance equilibrium is $T_{\mathrm{eq}}=S^\ast$. The differential equation becomes $dT/dt=k_{\mathrm{tol}}(T_{\mathrm{eq}}-T)$. Its solution describes how quickly the state approaches the concentration-specific equilibrium.

Derivation: If $T(0)=T_{\mathrm{eq}}$, substitution into the differential equation gives $dT/dt=0$, so $T(t)=T_{\mathrm{eq}}$ for all $t$. If $T(0)\ne T_{\mathrm{eq}}$, separate variables and integrate: $$ \int_{T(0)}^{T(t)}\frac{dT}{T_{\mathrm{eq}}-T} = k_{\mathrm{tol}}\int_0^t du. $$ The left side is $-\ln|T_{\mathrm{eq}}-T(t)|+\ln|T_{\mathrm{eq}}-T(0)|$. Therefore, $$ \ln\left|\frac{T_{\mathrm{eq}}-T(0)}{T_{\mathrm{eq}}-T(t)}\right| =k_{\mathrm{tol}}t, \qquad T(t)=T_{\mathrm{eq}}+\{T(0)-T_{\mathrm{eq}}\}e^{-k_{\mathrm{tol}}t}. $$ This formula also agrees with the constant solution when $T(0)=T_{\mathrm{eq}}$, although the separated integral itself is used only when the initial state differs from equilibrium.

The exponential term shows that the distance from equilibrium is multiplied by $e^{-k_{\mathrm{tol}}t}$. The time to remove half of that distance is $\ln(2)/k_{\mathrm{tol}}$. If concentration changes, the equilibrium changes as well; for a piecewise-constant concentration over an interval of length $\Delta t$, the corresponding update is

T_{\mathrm{next}}=S(C)+\{T_{\mathrm{current}}-S(C)\}e^{-k_{\mathrm{tol}}\Delta t}.

Applying this recurrence across intervals makes explicit how the prior exposure history influences the current prediction. It is an approximation when concentration varies within an interval; with a continuous concentration profile, the differential equation should instead be integrated along that profile.

Worked Example: Tolerance During a Concentration Plateau

Suppose, hypothetically, concentration is held at $4.00$ ng/mL, $EC_{50}=2.00$ ng/mL, $k_{\mathrm{tol}}=0.200$ h$^{-1}$, and $T(0)=0$. Let baseline analgesia be $E_0=5.00$ percentage points and $E_{\max}=60.0$ percentage points. Calculate the predicted analgesia initially and after 6.00 hours.

The concentration stimulus and equilibrium tolerance are $S^\ast=4.00/(2.00+4.00)=0.6667$ and $T_{\mathrm{eq}}=0.6667$. Initially, with $T(0)=0$, $E(0)=5.00+60.0(0.6667)/(1+0)=45.0$ percentage points.

At 6.00 hours, $T(6)=0.6667+(0-0.6667)e^{-(0.200)(6.00)}$. Since $e^{-1.20}=0.3012$, $T(6)=0.6667-0.2008=0.4659$ (rounding the state to four significant figures). Thus, $E(6)=5.00+60.0(0.6667)/(1+0.4660)=5.00+40.00/1.466=32.29$ percentage points.

At unchanged concentration, the model predicts a reduction from 45.0 to 32.29 percentage points as the tolerance state develops. At equilibrium, the predicted effect would be $5.00+40.00/(1+0.6667)=29.00$ percentage points. These are model-based calculations under the stated constant-concentration assumption, not clinical predictions.

Worked Repeated-Dose Pattern and Time-Dependent Errors

For a repeated-dose record, a useful first comparison is between the model with $T(t)$ and a concentration-only model with $T=0$. If the latter overpredicts later analgesia at comparable concentrations, while the tolerance model removes that pattern, the added state may capture adaptation. The comparison is informative only if concentrations, timing, baseline, and measurement quality have also been examined.

Worked Example: Updating Tolerance Across Repeated-Dose Intervals

Consider a hypothetical patient-level dataset summarized by interval concentrations. Assume the concentration is approximated as constant within each 4.00-hour interval, alternating between $4.00$ and $2.00$ ng/mL. Use $EC_{50}=2.00$ ng/mL, $k_{\mathrm{tol}}=0.200$ h$^{-1}$, $T(0)=0$, $E_0=5.00$, and $E_{\max}=60.0$ percentage points. The hypothetical analgesia observations at the ends of the intervals are 35, 26, 29, and 25 percentage points, respectively. Calculate predictions and residuals, defining residual as observed minus predicted.

For every 4.00-hour interval, $e^{-k_{\mathrm{tol}}\Delta t}=e^{-0.800}=0.4493$. The stimuli at concentrations 4.00 and 2.00 ng/mL are $0.6667$ and $0.5000$, respectively. Start with $T(0)=0$, then update: $T_4=0.6667+(0-0.6667)(0.4493)=0.3671$; $T_8=0.5000+(0.3671-0.5000)(0.4493)=0.4403$; $T_{12}=0.6667+(0.4403-0.6667)(0.4493)=0.5650$; and $T_{16}=0.5000+(0.5650-0.5000)(0.4493)=0.5292$.

At 4 hours, the interval concentration is 4.00 ng/mL, so the prediction is $5.00+40.00/(1+0.3671)=34.26$; its residual is $35-34.26=0.74$ percentage points. At 8 hours, the prediction is $5.00+30.00/(1+0.4403)=25.83$; its residual is $26-25.83=0.17$ percentage points. At 12 hours, it is $5.00+40.00/(1+0.5650)=30.56$, giving $29-30.56=-1.56$ percentage points. At 16 hours, it is $5.00+30.00/(1+0.5292)=24.62$, giving $25-24.62=0.38$ percentage points.

For comparison, the concentration-only model predicts 45.0 percentage points whenever concentration is 4.00 ng/mL and 35.0 percentage points whenever it is 2.00 ng/mL. Its residuals are therefore $-10.0$, $-9.00$, $-16.0$, and $-10.0$ percentage points. In this small hypothetical series, the concentration-only model consistently overpredicts analgesia, whereas the tolerance model’s residuals are smaller and do not show the same one-direction pattern. Four observations cannot establish that tolerance is the true explanation; the example illustrates a diagnostic comparison, not a formal model-selection result.

Interpreting the Pattern and Checking Alternatives

A time-dependent decline in observed effect at comparable concentration is compatible with tolerance, but it has competing explanations. Concentration may have been sampled too sparsely to represent peak or trough exposure; incorrect timing can make later observations appear to have lower response at a given recorded concentration. Delayed equilibration between plasma concentration and the effect site can also produce a changing concentration–effect relationship. A changing pain stimulus, rescue medication, disease progression, or drift in the measurement scale can change analgesia without a tolerance process.

The key diagnostic is not whether adding $T(t)$ improves a single summary measure. Inspect predictions and residuals against time, concentration, dose interval, and prior exposure. A successful tolerance component should reduce a systematic time pattern while leaving residuals without obvious structure. Residuals that remain increasingly negative over time may indicate that adaptation is stronger or otherwise different from the chosen model. Residuals that change sign around dose times may point to inadequately represented concentration profiles or delayed effect. A large improvement driven by one observation should be treated cautiously.

Check the exposure record.
Confirm units, dose and sample times, and whether the concentration input represents the interval adequately. If predicted concentrations are used, uncertainty in them can affect the apparent PD pattern.

Compare time patterns.
Plot observed and predicted analgesia over time, and examine observed-minus-predicted residuals by time and concentration for each candidate model.

Challenge the interpretation.
Ask whether delayed drug effect, baseline change, measurement drift, or other time-varying influences can explain the same residual pattern. Do not label a decline as tolerance solely because it follows repeated dosing.

Check adequacy and scope.
Look for implausible parameter estimates, poor prediction at low exposure, and systematic residual structure. State explicitly the exposure range and dosing conditions over which the model has been assessed.

The alternatives answer different questions. A concentration-only model is simpler and is appropriate when the response is adequately explained by current concentration with no remaining time trend. An effect-compartment model represents a delay between measured concentration and effect; it can explain hysteresis without requiring a loss of sensitivity. A tolerance model instead makes the effect at a given concentration depend on prior exposure through $T(t)$. A time-varying baseline model may be more appropriate when the underlying symptom level changes independently of drug exposure. These explanations can be difficult to distinguish from sparse concentration and response data, so model comparison should be guided by sampling design and the expected time course, not just by a better numerical fit.

Worked Example: Quantifying the Effect of Continued Exposure

Use the plateau assumptions from the first example: $C=4.00$ ng/mL, $EC_{50}=2.00$ ng/mL, $k_{\mathrm{tol}}=0.200$ h$^{-1}$, $T(0)=0$, $E_0=5.00$, and $E_{\max}=60.0$ percentage points. Compare predicted analgesia after 12.0 hours of continued exposure with the concentration-only prediction.

The tolerance state is $T(12)=0.6667+(0-0.6667)e^{-(0.200)(12)}$. Since $e^{-2.40}=0.09072$, $T(12)=0.6667-0.06048=0.6062$. The tolerance-model effect is $E(12)=5.00+40.00/(1+0.6062)=5.00+40.00/1.6062=29.90$ percentage points.

The concentration-only model sets $T=0$, so at the same concentration it predicts $5.00+40.00=45.0$ percentage points. The difference is $29.90-45.0=-15.1$ percentage points: under the stated model, continued exposure is associated with a 15.1-percentage-point lower predicted effect at 12 hours than the no-tolerance prediction. This contrast quantifies the consequence of the assumed adaptation component; it does not establish a causal clinical effect.

Limitations and a Common Pitfall

This simple model uses one state, one rate constant, and one empirical attenuation function. It assumes the tolerance state responds to current concentration stimulus and that the same rate governs its increase and recovery. Real adaptation may depend on dose history in a more complex way, may operate on multiple timescales, or may not be captured by plasma concentration alone. The model also treats $E_0$ as fixed and may be inappropriate for a bounded, ordinal, or otherwise non-continuous response scale without a suitable observation model. Parameter estimates can be weakly identified when concentration changes little, when sampling is sparse, or when tolerance and delayed effect produce similar patterns.

Common pitfall: “Analgesia fell after repeated doses, so tolerance has been demonstrated.” This reasoning mistakes a time sequence for a unique mechanism. The correction is to test whether the decline remains after accounting for concentration and timing, compare plausible alternatives such as delayed equilibration and baseline drift, and inspect time-dependent prediction errors. Even a better-fitting tolerance model supports only the adequacy of that description for the observed data; it does not, by itself, prove a biological mechanism.

Key takeaway: A tolerance state makes predicted effect depend on both current concentration and exposure history. Its value is clearest when it explains a systematic time pattern in response at comparable concentrations, while diagnostics and alternative explanations remain part of the interpretation.

Check Your Understanding

Answer each question using the model and assumptions stated above. Solutions follow the questions.

For a constant concentration, what is the equilibrium tolerance state $T_{\mathrm{eq}}$?
If the initial tolerance state equals $T_{\mathrm{eq}}$, what is $T(t)$, and why must this case be handled separately in the separation-of-variables derivation?
At $C=2.00$ ng/mL and $EC_{50}=2.00$ ng/mL, calculate $S(C)$. If $T=0.500$, $E_0=5.00$, and $E_{\max}=60.0$ percentage points, calculate predicted analgesia.
Name one residual pattern that may suggest the tolerance model is inadequate, and one alternative explanation for declining analgesia over time.

Solutions: (1) $T_{\mathrm{eq}}=S(C)=C/(EC_{50}+C)$, the stimulus at that constant concentration. (2) $T(t)=T_{\mathrm{eq}}$, because the differential equation has zero derivative there; the separated integrand $1/(T_{\mathrm{eq}}-T)$ is singular at that initial value, so the integral derivation applies only when $T(0)\ne T_{\mathrm{eq}}$. (3) $S=2.00/(2.00+2.00)=0.500$. The prediction is $5.00+60.0(0.500)/(1+0.500)=5.00+30.0/1.500=25.0$ percentage points. (4) Residuals that become increasingly negative over time may indicate inadequately modeled adaptation; delayed concentration-to-effect equilibration, changing baseline symptoms, or exposure measurement error could also explain a declining response.