Resolving Hysteresis in a Bronchodilator Case

Why a Bronchodilator Response Can Lag Concentration

After a bronchodilator dose, plasma concentration may rise and then fall, while bronchodilation develops more slowly and persists after concentration has begun to decline. A plot of response against concentration can then trace a loop rather than a single exposure-response curve. The order in which the points occur matters: the same concentration can be associated with different responses on the rising and falling portions of the time course.

The case below uses hypothetical concentration measurements and baseline-corrected bronchodilation responses. The question is not merely whether concentration and response are associated. It is whether the timing pattern is consistent with a delayed response, and whether a simple effect-compartment model can describe that delay.

Hypothetical Case Data and the Loop

Suppose concentration is measured in hypothetical concentration units and bronchodilation is recorded as a baseline-corrected fraction of the maximum possible response. The six measurements are:

Concentration peaks at 2 h, whereas the largest recorded response is at 4 h. The time ordering therefore suggests a lag. To see the loop direction more directly, linearly interpolate between the 4 h and 6 h measurements: at concentration 4.0, the interpolated response on the declining limb is about \((0.431+0.380)/2=0.4055\). At the measured concentration of 4.0 on the rising limb at 1 h, response is 0.226. The larger response at the same approximate concentration on the declining limb is consistent with counter-clockwise hysteresis. This interpolation is illustrative; it does not add an observed data point or account for measurement uncertainty.

Effect-Compartment Model

A standard way to represent a lag between measured concentration and response is to posit an effect-site concentration, \(C_e(t)\), that approaches the measured concentration, \(C(t)\), at a finite rate. The effect-site concentration is a model quantity; it is not necessarily a directly measured concentration in a physical compartment.

A larger \(k_{e0}\) makes \(C_e(t)\) track \(C(t)\) more quickly. As \(k_{e0}\) becomes very large, the model approaches an immediate concentration-response relationship, \(C_e(t)\approx C(t)\). A smaller positive rate produces a greater lag. This model assumes a constant equilibration rate, a stable response relationship, and no additional time-varying influence such as tolerance or a changing baseline. The case calculation also treats the measured concentrations as known and the response residuals as suitable for least-squares fitting.

Calculating Effect-Site Concentrations

Concentrations are observed only at discrete times. For a transparent worked fit, assume concentration is constant within each interval and equal to its value at the interval endpoint. If concentration is \(C_i\) for an interval of length \(\Delta t_i\), solving the differential equation gives:

This is an exact solution under the stated piecewise-constant input assumption. It is not an exact solution for a continuously changing concentration profile. In practice, denser PK samples or a specified interpolation of the concentration-time curve may be preferable.

Fitting the Equilibration Rate

For this illustration, set \(E_{\max}=0.80\) and \(EC_{50}=4.0\) concentration units in advance, and estimate only \(k_{e0}\). This isolates the rate calculation; it does not establish that the fixed response parameters are known in a real analysis. For each candidate rate, recursively calculate \(C_e\), predict response, and calculate the residual sum of squares:

The rate is selected to minimize \(RSS\) over positive values. The comparison uses the same observations and fixed response parameters for every candidate. Least squares is most directly interpreted when response errors are independent with a common variance; repeated measurements may instead have correlated errors, which this simple calculation does not model.

What the Estimated Rate Does and Does Not Mean

The estimated \(k_{e0}\) describes the speed at which modeled effect-site concentration approaches plasma concentration under this model. Its reciprocal is a characteristic time scale, while the time for a first-order effect compartment to close half the gap between \(C_e\) and a fixed \(C\) is \(\ln(2)/k_{e0}\). At \(k_{e0}\approx0.50\ \mathrm{h}^{-1}\), that half-time is approximately \(0.693/0.50=1.39\) h. This is a model-based equilibration half-time, not the half-life of the drug in plasma or a direct measurement of receptor binding.

A finite estimated rate and a better fit than the direct-effect limit support a delayed-effect description of these data. They do not identify the biological location of the delay. The same pattern could arise from an unmeasured active metabolite, measurement timing, a changing baseline, or another time-dependent process. The rate may also depend on assumptions about concentration between sample times and on the fixed values of \(E_{\max}\) and \(EC_{50}\). A real analysis should examine parameter uncertainty, residual patterns, sensitivity to PK interpolation, and whether allowing response parameters to be estimated changes the rate.

Diagnostics, Alternatives, and Common Pitfalls

Useful checks begin with the time-course plot and the residuals, not just the minimized objective. Plot response against time and concentration against time, then connect concentration-response points in chronological order. For a fitted model, inspect observed-minus-predicted residuals over time and concentration. A persistent signed residual pattern, a missed peak, or large deviations at the late samples can indicate a poor input approximation, a misspecified response curve, or a process beyond simple equilibration. Sparse sampling can also make the loop direction or rate difficult to establish.

A direct concentration-response model is the main simpler alternative: it assumes response follows concentration without a lag. It is parsimonious, but cannot represent different responses at the same concentration on rising and falling limbs. A more flexible time-varying response curve can describe a loop, but may not provide a compact, interpretable equilibration rate and can overfit sparse data. The effect-compartment model offers a mechanistic-looking summary with few parameters, but its effect-site state is latent and its rate is conditional on the model assumptions.

• Pitfall: treating a loop as proof of an effect compartment. The wrong reasoning is that counter-clockwise hysteresis uniquely establishes delayed equilibration. The correction is that the loop supports a time-dependent response; alternative biological and measurement explanations must still be considered.
• Pitfall: reading the rate as a plasma elimination rate. \(k_{e0}\) governs modeled equilibration between measured concentration and effect-site concentration, not plasma elimination.
• Pitfall: reporting a precise rate from a small case data set. An optimizer returns a value even with sparse, rounded observations. Report the conditional estimate, assess uncertainty and sensitivity, and avoid interpreting numerical precision as evidential certainty.