Sampling Times That Reveal Turnover Rates

Why Sampling Time Matters

An indirect response does not generally follow a change in drug effect instantly. The response reflects a balance between its production and loss, so its rate of approach to a new level contains information about turnover. But that information is time-dependent: measurements taken only very early or only after the response has nearly settled can reveal little about the rate. Sampling should capture the response while it is changing, as well as its subsequent recovery when feasible.

The guidance here is conditional on a specified indirect-response model. For the basic case, drug changes the input rate while the loss rate remains constant. A drug effect that changes the loss rate has a different on-treatment rate of approach; that rate should not be mistaken for the baseline turnover rate.

Step Changes: The Onset and Recovery Curves

Suppose production is \(k_{\text{in}}\), loss is first-order, and the response is at baseline steady state before drug is given. Then \(R_0=k_{\text{in}}/k_{\text{out}}\). If a constant drug effect changes production at time zero, while \(k_{\text{out}}\) stays constant, the response moves toward a new equilibrium \(R_{\text{ss}}\):

Writing \(A=R_{\text{ss}}-R_0\), the response displacement from baseline during onset is \(\Delta R(t)=A(1-e^{-k_{\text{out}}t})\). If the effect is removed at time \(T\), and the system returns to the same baseline with the same loss rate, its displacement \(D\) then decays as \(D(T+s)=D(T)e^{-k_{\text{out}}s}\). These expressions assume an abrupt change to a constant effect, a constant first-order loss rate, and no baseline drift during the interval.

Onset and recovery provide complementary views. During onset, an early sample detects the initial direction and slope, a sample around \(\tau\) captures substantial curvature, and a later sample helps establish how close the response is to its new level. During recovery, samples at separated times show the exponential decline directly. A sample only at the peak or plateau gives little information about the rate of approach; samples on both sides of a peak can help establish when the response changes direction.

Why the Time Constant Is Informative

A short sensitivity calculation explains why sampling around \(\tau\) is useful. If the equilibrium displacement \(A\) is fixed, the sensitivity of the onset displacement to the loss rate is \(\partial\Delta R/\partial k_{\text{out}}=At e^{-k_{\text{out}}t}\). Its magnitude is small near time zero because little time has elapsed, and it becomes small at late times because the exponential term has decayed. Differentiating with respect to time shows that this sensitivity is largest at \(t=1/k_{\text{out}}\):

The same time point maximizes the absolute sensitivity of an exponential recovery displacement to \(k_{\text{out}}\), since \(\left|\partial(D_0e^{-k_{\text{out}}t})/\partial k_{\text{out}}\right|=|D_0|t e^{-k_{\text{out}}t}\). For nonzero \(D_0\), this magnitude is maximized at \(t=1/k_{\text{out}}\); if \(D_0=0\), it is zero at every time. This is a local sensitivity result, not a guarantee that one sample at \(\tau\) identifies the rate in a fitted model. If \(A\) or \(D_0\) is unknown, it must also be estimated; repeated measurements at distinct times help distinguish the rate from the response amplitude.

The expected fractions also give a practical scale for a schedule. Onset has completed 63.2%, 86.5%, and 95.0% of its total change at \(\tau\), \(2\tau\), and \(3\tau\), respectively. Recovery has 36.8%, 13.5%, and 5.0% of its initial displacement remaining at those times. Very early and very late observations still have roles, but an informative design should ordinarily include observations in the changing region, rather than concentrating all samples at either extreme.

Planning Around a Peak

For a time-varying drug effect, the response is driven by both the changing input and turnover. A peak is the time at which the response stops rising and begins to fall; it is not necessarily the time of peak drug concentration. For a production-driven response, the instantaneous change in response is the drug-modified production minus the loss at the current response. At a smooth peak, those two rates balance. Thus, a peak observation alone does not generally identify \(k_{\text{out}}\): the effect size and the exposure profile also influence the balance.

A practical schedule therefore aims to observe the response before the expected peak, near it, and after it. Samples on the rising side show the lag and onset; measurements near the peak help locate the change in direction; observations on the falling side or after drug effect wanes show recovery. The response time constant supplies a first guide to spacing: if prior information suggests \(\tau\) is several hours, observations separated by only a few minutes are unlikely to span meaningful turnover, while a schedule that ends well before \(\tau\) may miss much of the rate-sensitive curvature.

Assumptions, Diagnostics, and Alternatives

Sampling guidance depends on knowing what process governs each interval. If the drug effect changes loss rather than production, the on-treatment exponential rate is the drug-modified loss rate, not necessarily baseline \(k_{\text{out}}\). If concentration changes continuously, the response is not generally the step solution above. Use the appropriate time-varying model and, when needed, evaluate candidate sampling schedules against its predicted response and parameter sensitivities.

A useful diagnostic is whether the observed onset or recovery displacement is compatible with the assumed exponential pattern. Recovery ratios over equal time intervals should be similar when the same first-order rate applies: \(D(t+\Delta)/D(t)=e^{-k_{\text{out}}\Delta}\). Strongly inconsistent ratios can indicate measurement noise, a drifting baseline, incorrect timing, a changing drug effect, or a misspecified turnover mechanism. A measured peak that falls between widely spaced samples is poorly localized; denser sampling around the expected turning point can help, although noise may make the largest observed value differ from the true peak.

A two-point recovery calculation is simple and transparent, but sensitive to errors in either observation and to baseline error. A joint fit to onset, peak-region, and recovery observations can use more of the profile and estimate unknown amplitudes, but it relies on the model and its assumptions. A rich, nearly uniform schedule may be more robust when the time constant is uncertain, at the cost of more observations. A targeted schedule can be efficient when plausible rates and the effect profile are known; it should still span enough time to distinguish rate from amplitude.

Common mistake: “The response is largest at the peak, so that is the best time to estimate \(k_{\text{out}}\).” The peak height depends on effect magnitude and exposure as well as turnover, and a single peak value does not show how fast the response approached or left that peak. Instead, sample the rising and falling phases, include times informed by the expected time constant, and interpret the peak in the context of the full input-response model.