Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Repeated Dosing and Target Recovery

Target-Mediated Drug Disposition and PK/PD · Tutorial 257 of 598

Repeated Dosing and Target Recovery

Use dose interval, drug elimination, and target turnover to interpret changing TMDD exposure and target-mediated response under repeated dosing.

Advanced 12 min read

What You'll Learn

  • Write repeated-dose input and state jumps for the full TMDD model
  • Distinguish residual drug accumulation from target depletion and recovery
  • Derive an interval-dependent target recovery recurrence
  • Calculate drug accumulation and target recovery for hypothetical dosing schedules
  • Recognize when simplified recovery calculations do not represent full TMDD dynamics
  • Diagnose misleading interpretations of changing exposure or response

Why Repeated Dosing Requires Two Clocks

After a dose, drug concentration changes as free drug is eliminated, binds target, and may be released from or removed with drug-target complex. Meanwhile, target is synthesized, turns over, and may be depleted through complex formation and internalization. Repeated dosing therefore involves at least two timescales: the time over which drug remains in the body and the time over which the target pool recovers. A concentration that changes from one dose to the next does not, by itself, reveal which process is responsible.

The useful distinction is between drug accumulation, meaning residual drug from earlier doses contributes to concentrations after later doses, and target depletion or recovery, meaning the abundance of available target changes over time. These processes can occur together, and either can affect observed pharmacodynamic response. The equations below use the full TMDD state variables: free drug \(D\), free target \(R\), and drug-target complex \(C\).

Repeated-Dose Full TMDD Model

Definition: Between doses, one standard full TMDD model is
$$ \begin{aligned} \frac{dD}{dt} &= -k_{\mathrm{el}}D-k_{\mathrm{on}}DR+k_{\mathrm{off}}C,\\ \frac{dR}{dt} &= k_{\mathrm{syn}}-k_{\mathrm{deg}}R-k_{\mathrm{on}}DR+k_{\mathrm{off}}C,\\ \frac{dC}{dt} &= k_{\mathrm{on}}DR-(k_{\mathrm{off}}+k_{\mathrm{int}})C. \end{aligned} $$
Here \(k_{\mathrm{el}}\) is the first-order elimination rate constant for free drug; \(k_{\mathrm{on}}\) and \(k_{\mathrm{off}}\) are association and dissociation rate constants; \(k_{\mathrm{int}}\) is the complex internalization rate constant; \(k_{\mathrm{syn}}\) is the target synthesis rate; and \(k_{\mathrm{deg}}\) is the degradation rate constant for free target. Concentrations must be expressed on compatible volume and amount scales, so the binding terms are dimensionally consistent.

The equations assume a well-mixed compartment, constant rate parameters, and the stated processes as the relevant routes of drug and target turnover. They describe concentrations between doses; for an intravenous bolus of amount \(Q\) into a compartment of volume \(V\), the dose produces the jump \(D(t^+)=D(t^-)+Q/V\), while \(R\) and \(C\) are continuous at the instant of dosing. A repeated schedule with interval \(\tau\) applies this jump at each dosing time \(t_j=j\tau\). Other routes require an absorption or input model rather than this instantaneous jump.

Without drug, the baseline target concentration is \(R_0=k_{\mathrm{syn}}/k_{\mathrm{deg}}\). During dosing, free target alone is not the whole target pool: total target is \(R_{\mathrm{tot}}=R+C\). Adding the target and complex equations gives

$$ \frac{dR_{\mathrm{tot}}}{dt} =k_{\mathrm{syn}}-k_{\mathrm{deg}}R-k_{\mathrm{int}}C. $$

Binding and dissociation cancel in this total-pool balance because they redistribute target between free and complexed forms. Synthesis adds target, while degradation of free target and internalization of complex remove it. Consequently, a low free-target measurement may indicate substantial binding rather than depletion of the total target pool. Target depletion should be interpreted using the measured or modeled pool that matches the biological question.

Recovery Between Doses: A Useful Limiting Calculation

When drug and complex have become negligible after a dose, the target equation reduces to \(dR/dt=k_{\mathrm{syn}}-k_{\mathrm{deg}}R\). Solving this linear equation from a post-dose concentration \(R_{\mathrm{post}}\) gives

$$ R(t+\Delta)=R_0+\bigl(R_{\mathrm{post}}-R_0\bigr)e^{-k_{\mathrm{deg}}\Delta}. $$

Thus recovery is governed by target turnover, not simply by the drug elimination half-life. If the target is at a fraction \(r\) of baseline immediately after the dose, its fraction at time \(\Delta\) is \(1-(1-r)e^{-k_{\mathrm{deg}}\Delta}\). This is a limiting approximation: while appreciable drug or complex remains, binding, dissociation, and complex internalization continue to affect the target trajectory.

Worked Example: Target Recovery After One Dose

Suppose a hypothetical target has a baseline concentration of \(100\ \mathrm{nM}\) and a free-target turnover half-life of 7 days. A sufficiently brief, strong exposure leaves the free target at \(20\ \mathrm{nM}\); assume drug and complex are negligible by day 7. Estimate free target on day 7 after the dose using the recovery approximation.

The degradation rate is \(k_{\mathrm{deg}}=\ln(2)/7=0.09902\ \mathrm{day}^{-1}\). With \(\Delta=7\) days, \(e^{-k_{\mathrm{deg}}\Delta}=e^{-\ln(2)}=0.5000\). Therefore,

$$ R(7)=100+(20-100)(0.5000)=60.0\ \mathrm{nM}. $$

The target has recovered from \(20\ \mathrm{nM}\) to \(60.0\ \mathrm{nM}\), not to baseline. The calculation assumes the post-dose value is free target, target synthesis and degradation remain constant, and drug-mediated processes are negligible during the recovery interval. If complex persists, this exponential calculation alone is not an adequate description.

Repeated Target Depletion and Interval Effects

A simple recurrence illustrates how incomplete recovery can accumulate across a schedule. Suppose, as an explicit approximation, each dose instantaneously removes the same fraction \(f\) of the free target present immediately before dosing. Between doses, assume the drug and complex become negligible and target recovers according to the preceding exponential solution. Let \(R_n\) denote free target immediately before dose \(n\), and let \(a=e^{-k_{\mathrm{deg}}\tau}\). The post-dose value is \((1-f)R_n\), so the next pre-dose value is

$$ R_{n+1}=R_0+\bigl((1-f)R_n-R_0\bigr)a. $$

This recurrence is not a replacement for the full TMDD model. It is a transparent way to isolate the effect of target turnover and interval under the stated pulse-depletion approximation. If it reaches a stable repeating value \(R_*\), set \(R_{n+1}=R_n=R_*\) and solve:

$$ R_*=\frac{R_0(1-a)}{1-a(1-f)}. $$

For \(0\leq a<1\) and \(0\leq f\leq1\), the denominator is positive. The recurrence contracts deviations from its stable value by the factor \(a(1-f)\), whose magnitude is less than one, so it approaches that value under these assumptions. A shorter interval makes \(a\) larger: less recovery occurs before the next dose. This does not establish that a shorter interval always produces a particular exposure or response in a full TMDD system, because dose size and ongoing drug-target dynamics also matter.

Worked Example: Comparing Target Recovery at Two Intervals

Consider a hypothetical target with \(R_0=100\ \mathrm{nM}\), a 7-day turnover half-life, and an assumed per-dose depletion fraction \(f=0.65\). Compare the stable pre-dose target concentration for intervals of 7 and 21 days using the recurrence. The assumptions are that each dose has the same fractional pulse effect and drug and complex are negligible between doses.

For a 7-day interval, \(a=e^{-(\ln2/7)7}=0.5000\). Then

$$ R_*=\frac{100(1-0.5000)}{1-(0.5000)(1-0.65)} =\frac{50.00}{0.8250} =60.61\ \mathrm{nM}. $$

Immediately after a dose, the corresponding target concentration is \((1-0.65)(60.61)=21.21\ \mathrm{nM}\). For a 21-day interval, \(a=e^{-3\ln2}=0.1250\), so

$$ R_*=\frac{100(1-0.1250)}{1-(0.1250)(0.35)} =\frac{87.50}{0.95625} =91.50\ \mathrm{nM}. $$

The predicted pre-dose target is higher with the longer interval because more turnover-mediated recovery occurs between doses. These numbers describe the simplified recurrence, not a fitted TMDD model or a clinical prediction. A full model may yield different depletion fractions at different doses because free drug, target, and complex concentrations change together.

Drug Accumulation Is a Separate Calculation

A common first approximation to drug accumulation is the linear one-compartment model with instantaneous doses, first-order elimination, constant volume, and no target-mediated disposition. Under those assumptions, the fraction of a concentration increment remaining one interval later is \(q=e^{-k_{\mathrm{el}}\tau}\). At periodic steady state, the concentration immediately after a dose is the sum of the current increment and all surviving earlier increments:

$$ C_{\mathrm{post,ss}}=\Delta C(1+q+q^2+\cdots)=\frac{\Delta C}{1-e^{-k_{\mathrm{el}}\tau}}, $$

where \(\Delta C=Q/V\). The accumulation ratio relative to the concentration increment from one dose is therefore \(1/(1-e^{-k_{\mathrm{el}}\tau})\). This result is a benchmark, not a general repeated-dose TMDD formula. In TMDD, target binding and complex disposition can make drug clearance dependent on concentration and on the evolving target pool.

Worked Example: Estimating Linear Drug Accumulation

A hypothetical drug has a terminal half-life of 5 days and is given as an intravenous bolus every 7 days. Estimate the accumulation ratio under the linear, one-compartment assumptions, and explain what this calculation does not establish about target recovery.

The elimination rate is \(k_{\mathrm{el}}=\ln(2)/5=0.1386\ \mathrm{day}^{-1}\). The fraction remaining after 7 days is

$$ q=e^{-0.1386(7)}=0.3789. $$

Thus the steady-state accumulation ratio is

$$ \frac{1}{1-q}=\frac{1}{1-0.3789}=1.610. $$

Under this simplified model, the post-dose steady-state concentration increment is 1.610 times the increment from a single dose. This does not mean target has recovered by a particular amount: the calculation contains no target synthesis, target turnover, binding, or complex internalization. Nor should the ratio be assumed to describe actual TMDD exposure without checking whether linear disposition is a reasonable approximation over the concentration range.

Reading Exposure and Response Together

Interpreting repeated-dose profiles requires tracking the states relevant to the observation. Free drug concentration, total drug concentration, free target, total target, and occupancy are not interchangeable. A response linked to occupancy can change because drug exposure changed, because target abundance changed, or because both changed. The introductory material on linking target occupancy to pharmacodynamic response establishes why response interpretation must follow the specified exposure-response relationship; the repeated-dose question adds time-varying target and drug states.

For example, lower free drug after a later dose does not prove that drug accumulation is absent: higher target availability can increase binding and target-mediated removal, while residual drug can still be present. Conversely, a higher trough concentration does not prove target recovery. The concentration profile and target trajectory must be evaluated jointly, with the observation model stated clearly. When the full model is used, numerical integration across every dose and the intervals between them preserves the changing state carried into the next dose.

Diagnostics, Alternatives, and Limitations

Useful checks start with model predictions at matched times in the dosing cycle. Compare pre-dose and post-dose free drug, free target, complex, and total target across successive cycles. If these states have not become periodic, a steady repeating profile has not yet been reached. Check whether the observed changes are consistent with the specified dose jumps, synthesis and turnover, and the model’s concentration definitions. Residual patterns or mismatches between free and total measurements may point to an inadequate observation model or missing disposition process, rather than proving target depletion.

The limiting exponential recovery calculation is appropriate only after drug-mediated processes are negligible and the free target returns toward a stable baseline under constant synthesis and degradation. The pulse-depletion recurrence adds the stronger assumption that every dose removes the same fraction of the target available before that dose. The linear accumulation ratio assumes no target-mediated clearance and first-order elimination. These are useful comparisons because each isolates a process, but none should be substituted for a full TMDD simulation when binding and complex persist through the interval.

A model can also be numerically plausible yet weakly supported by data. Sparse sampling may fail to distinguish rapid changes in drug concentration from slower target recovery; concentration assays may measure different combinations of free and bound species. Parameter uncertainty in target turnover or complex internalization can substantially alter predicted recovery. Treat numerical simulations as conditional on the model, parameters, dosing history, and measured species—not as direct evidence that one mechanism has been uniquely identified.

Common pitfall: “The trough concentration increased, so the target must have recovered.” This reasoning confuses drug accumulation with target turnover. A trough is a drug-concentration observation; target recovery must be assessed from target states or measurements and their turnover dynamics. In TMDD, changing target availability can itself alter drug disposition.

Check Your Understanding

Answer each question using the stated assumptions. Brief solutions follow the questions.

  1. If \(R_0=80\ \mathrm{nM}\), the target half-life is 4 days, and free target is \(16\ \mathrm{nM}\) after a dose, what is the predicted free target 4 days later if drug and complex are negligible?
  2. In the pulse-depletion recurrence, what does \(a=e^{-k_{\mathrm{deg}}\tau}\) represent, and how does it change when the interval becomes longer?
  3. A linear accumulation calculation predicts a ratio of 1.5. Does this determine the fraction of target recovered before the next dose? Why or why not?
  4. Why can free target decrease while total target does not decrease by the same amount?
  5. Name two assumptions required for the simple exponential recovery equation to apply between doses.

Solutions. (1) The target half-life equals the interval, so the unrecovered difference halves: \(R=80+(16-80)(0.5)=48.0\ \mathrm{nM}\). (2) It is the fraction of a deviation from baseline that remains after one interval under first-order target turnover. A longer interval makes \(a\) smaller, allowing more recovery. (3) No. The ratio concerns drug concentration under linear elimination and contains no target turnover or target concentration. (4) Target can move from free to complexed form through binding; this redistribution lowers \(R\) while increasing \(C\), without immediately removing the same amount of total target. (5) Among the required assumptions are negligible drug and complex during the interval, constant \(k_{\mathrm{syn}}\) and \(k_{\mathrm{deg}}\), and a stable baseline \(R_0=k_{\mathrm{syn}}/k_{\mathrm{deg}}\).

Clinical Biostats