Time-Varying Loss and the Survival-Factor Solution

Why Time-Varying Loss Needs Its Own Weighting

In an indirect response model, the measured response changes through production and loss. When drug effect changes the loss rate over time, response present at an earlier time does not simply decay at one fixed rate: its subsequent loss depends on the entire later history of the drug effect. The same history also determines how much of the response produced along the way remains at the observation time.

The integrating-factor solution makes this history explicit. Its central term is a survival factor: the fraction of response present at one time that remains at a later time under the intervening loss rates. Here, “survival” refers to persistence of response in a deterministic turnover model, not a probability that an individual or event survives.

Definition: Let $R(t)$ be a response, $k_{\mathrm{in}}(t)$ its production rate, and $k_{\mathrm{out}}(t)$ its first-order loss rate. The time-varying-loss model is $$ \frac{dR(t)}{dt}=k_{\mathrm{in}}(t)-k_{\mathrm{out}}(t)R(t). $$ For times $s\leq t$, define the survival factor $$ S(s,t)=\exp\left(-\int_s^t k_{\mathrm{out}}(u)\,du\right). $$ It is the fraction of response present at time $s$ that remains at time $t$ when no additional response is added after $s$.

The model assumes that loss at each instant is proportional to the response then present, with proportionality $k_{\mathrm{out}}(t)$. The production and loss rates are specified functions of time (or functions obtained from a specified concentration-effect model and concentration trajectory). For the standard solution, the rates must be integrable over the time interval; $k_{\mathrm{out}}(t)$ is nonnegative for a loss process. The response is a single turnover state: the equation does not include a separate delay compartment or a drug effect acting directly on the measured response.

Deriving the Survival-Factor Solution

Write the equation in linear form as $dR/dt+k_{\mathrm{out}}(t)R=k_{\mathrm{in}}(t)$. The integrating factor from an initial time $t_0$ to $t$ is $\exp(\int_{t_0}^t k_{\mathrm{out}}(u)\,du)$. Multiplying through by it turns the left side into a product derivative:

\frac{d}{dt}\left[ R(t)\exp\left(\int_{t_0}^t k_{\mathrm{out}}(u)\,du\right) \right] = k_{\mathrm{in}}(t)\exp\left(\int_{t_0}^t k_{\mathrm{out}}(u)\,du\right).

Integrate from $t_0$ to $t$, then divide by the integrating factor at $t$. The result is

R(t)=R(t_0)S(t_0,t)+\int_{t_0}^{t} k_{\mathrm{in}}(s)S(s,t)\,ds.

The first term is the initial response multiplied by the fraction that persists to $t$. In the integral, $k_{\mathrm{in}}(s)\,ds$ is response produced during a short interval around $s$; $S(s,t)$ is the fraction of that contribution remaining at $t$. Later production is exposed to less subsequent loss than earlier production. This is the key distinction between a time-varying-loss solution and a calculation that applies one decay rate to all response accumulated over the interval.

Connecting Loss to Concentration

For stimulated loss, a common specification is a baseline loss rate multiplied by a concentration-effect multiplier:

k_{\mathrm{out}}(t) = k_{\mathrm{out},0} \left[ 1+\frac{E_{\max}C(t)}{EC_{50}+C(t)} \right].

Here $k_{\mathrm{out},0}$ is the drug-free loss rate, $C(t)$ is concentration, $E_{\max}$ is the maximum fractional increase in loss, and $EC_{50}$ is the concentration producing half of that maximum increase under this model. The multiplier is one at zero concentration and rises toward $1+E_{\max}$ at high concentration. A different concentration-effect function can be used, provided it specifies a suitable nonnegative loss rate. The survival factor then depends on the integral of the resulting rate, not merely on the concentration at the observation time.

Worked Example: A Concentration Step Followed by Washout

Consider a hypothetical response with constant production $k_{\mathrm{in}}=10$ response units/hour, baseline loss $k_{\mathrm{out},0}=0.20\ \mathrm{hour}^{-1}$, and initial response $R(0)=50$ response units. Let $E_{\max}=2$, $EC_{50}=2\ \mathrm{mg/L}$, and suppose $C(t)=2\ \mathrm{mg/L}$ from hour 0 to hour 1, then $C(t)=0$ from hour 1 to hour 2. Calculate $R(2)$.

The baseline is consistent with these production and loss rates because $10/0.20=50$ response units. During the first hour, the concentration-effect multiplier is $1+2(2)/(2+2)=2$, so $k_{\mathrm{out}}=0.40\ \mathrm{hour}^{-1}$. During the second hour, the multiplier returns to one and $k_{\mathrm{out}}=0.20\ \mathrm{hour}^{-1}$. The rate is constant within each interval, so the linear model has the interval solution $R(t+\Delta)=k_{\mathrm{in}}/k_{\mathrm{out}}+[R(t)-k_{\mathrm{in}}/k_{\mathrm{out}}]e^{-k_{\mathrm{out}}\Delta}$.

At hour 1, substitute the first interval’s values: $$ R(1)=\frac{10}{0.40}+\left(50-\frac{10}{0.40}\right)e^{-0.40(1)} =25+25(0.6703)=41.76\ \text{response units}. $$

During the second hour, the drug-free equilibrium is $10/0.20=50$ response units. Therefore, $$ R(2)=\frac{10}{0.20}+\left(41.76-\frac{10}{0.20}\right)e^{-0.20(1)} =50-8.242(0.8187)=43.25\ \text{response units}. $$

The response has not returned to baseline at hour 2. Although the loss rate returned to its baseline value at hour 1, the response still needed time to recover.

When the Loss Rate Changes Continuously

If $C(t)$ changes continuously, the concentration-effect function generally makes $k_{\mathrm{out}}(t)$ continuous as well. The survival factor still has the same definition, but its integral may require numerical evaluation. The response integral may also require quadrature, or the original differential equation may be solved numerically. These are computational choices for evaluating the same specified model, not reasons to replace the history of loss with the rate at the endpoint.

Worked Example: A Continuously Rising Concentration

Suppose, hypothetically, that $R(0)=50$ response units, $k_{\mathrm{in}}=10$ response units/hour, and $k_{\mathrm{out},0}=0.20\ \mathrm{hour}^{-1}$. From hour 0 to hour 2, let $C(t)=t\ \mathrm{mg/L}$, with time measured in hours, $E_{\max}=1.5$, and $EC_{50}=1\ \mathrm{mg/L}$. Find the approximate response at hour 2.

The resulting loss rate is $0.20[1+1.5t/(1+t)]\ \mathrm{hour}^{-1}$. Its cumulative loss from time 0 to $t$ is $$ K(t)=\int_0^t k_{\mathrm{out}}(u)\,du =0.5t-0.3\ln(1+t). $$ Thus $K(2)=0.6704$, and $S(0,2)=e^{-0.6704}=0.5114$. The initial response contributes $50(0.5114)=25.57$ response units at hour 2.

For production occurring at time $s$, the survival factor to hour 2 is $e^{-[K(2)-K(s)]}$. Using Simpson’s rule with step size $h=0.5$ hour, the values of this integrand at $s=0,0.5,1,1.5,2$ are approximately $0.5114, 0.5815, 0.6849, 0.8226, 1.0000$. Therefore, $$ \int_0^2 e^{-[K(2)-K(s)]}\,ds \approx \frac{0.5}{3} \left[0.5114+1+4(0.5815+0.8226)+2(0.6849)\right] =1.416\ \mathrm{hours}. $$

The surviving production contribution is approximately $10(1.416)=14.16$ response units. Combining it with the initial response contribution gives $R(2)\approx25.57+14.16=39.73$ response units. This value is an approximation because the production integral was evaluated by numerical quadrature; a finer grid can be used to check its numerical accuracy.

Why the Order of Loss Rates Matters

The survival factor for the initial response depends on cumulative loss, $\int k_{\mathrm{out}}(u)\,du$. But the final response also contains contributions produced at every intermediate time. Consequently, two schedules can have the same cumulative loss—and the same survival of the initial response—yet produce different final responses because their production contributions experience different amounts of subsequent loss.

Worked Example: Same Cumulative Loss, Different Final Responses

Use the hypothetical rates $k_{\mathrm{in}}=10$ response units/hour and $R(0)=50$ response units. Compare two two-hour schedules, each with one hour at $k_{\mathrm{out}}=0.40\ \mathrm{hour}^{-1}$ and one hour at $k_{\mathrm{out}}=0.20\ \mathrm{hour}^{-1}$.

For high loss first, the calculation in the preceding example gives $R(1)=41.76$, followed by $R(2)=43.25$ response units. For low loss first, the response stays at 50 during the first hour because $10-0.20(50)=0$. The high-loss second hour then gives $$ R(2)=\frac{10}{0.40}+\left(50-\frac{10}{0.40}\right)e^{-0.40} =25+25(0.6703)=41.76\ \text{response units}. $$

Both schedules have cumulative loss $0.40(1)+0.20(1)=0.60$, so the initial response contributes $50e^{-0.60}=27.44$ response units in either schedule. The final responses nevertheless differ: production occurring early in the high-loss-first schedule has time to be removed during the later low-loss interval, whereas production in the low-loss-first schedule is exposed to high loss near the end. The survival-factor integral accounts for this timing; cumulative loss alone does not determine the full response.

Assumptions, Diagnostics, and Alternatives

The solution is exact for the stated linear model when its rates and initial response are specified. Its application to data still depends on whether those assumptions are plausible. In particular, the loss term must be proportional to $R(t)$, and the time-varying rate supplied to the calculation must adequately represent the drug effect on loss. A concentration-effect multiplier is not itself evidence that loss is the mechanism; that is a model assumption to assess against the response data and study design.

Useful checks include confirming that concentration, time, and rate units are consistent; verifying that $k_{\mathrm{out}}(t)$ remains nonnegative; checking that the numerical solution approaches the analytic piecewise-constant solution when rates are held constant; and checking that the response changes in the direction implied by the production-loss balance. When $k_{\mathrm{in}}$ is constant and $R$ is above the instantaneous equilibrium $k_{\mathrm{in}}/k_{\mathrm{out}}(t)$, the derivative is negative; when it is below that value, the derivative is positive. With changing rates, this equilibrium is a moving reference, not a claim that the response instantly reaches equilibrium.

Systematic residual patterns, implausible response predictions, or sensitivity to a small change in the assumed concentration trajectory can signal model inadequacy or insufficient data to distinguish mechanisms. The solution does not address measurement error or parameter estimation. Nor does it cover nonlinear loss such as a saturable loss term, multiple response compartments, or an additional delay mechanism; those require a different model formulation. The introductory material on integrating-factor solutions for time-varying drug effects gives the general method, while this application emphasizes the contribution of time-varying loss to survival weighting. The earlier material on convolution solutions for drug-driven production is a useful contrast: the same survival weighting can be used to describe production contributions, but here the changing loss rate determines how each contribution is retained.

Common pitfall: Replacing $k_{\mathrm{out}}(t)$ throughout the interval by its value at the final time, or by its time average, and then applying a constant-rate formula is generally not exact. It ignores when production occurred. Use the integrated survival factor for each contribution, or solve the specified differential equation numerically.

Check Your Understanding

Answer each question, then compare with the solutions below.

If $k_{\mathrm{out}}(t)=0.3\ \mathrm{hour}^{-1}$ for two hours, what fraction of the response present at the start remains after two hours?
With $k_{\mathrm{in}}=8$ response units/hour and $k_{\mathrm{out}}=0.2\ \mathrm{hour}^{-1}$, what is the constant-rate equilibrium response?
For the model $k_{\mathrm{out}}=k_{\mathrm{out},0}[1+E_{\max}C/(EC_{50}+C)]$, what is the loss-rate multiplier when $C=EC_{50}$?
Two loss-rate schedules have the same integral over an interval. Does that guarantee the same final response when production continues throughout the interval? Explain.

Solutions. (1) The survival factor is $e^{-0.3(2)}=e^{-0.6}=0.5488$, so approximately 54.88% remains. (2) At equilibrium, production equals loss: $k_{\mathrm{in}}=k_{\mathrm{out}}R$, so $R=8/0.2=40$ response units. (3) Substituting $C=EC_{50}$ gives $1+E_{\max}/2$, because $C/(EC_{50}+C)=1/2$. (4) No. The schedules give the same survival of the initial response, but production added at different times experiences different subsequent loss; the survival-factor integral can therefore give different final responses.