Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Sampling Designs That Separate TMDD Processes

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

Sampling Designs That Separate TMDD Processes

Learn to choose sampling times, dose levels, and target measurements that provide distinct information about binding, internalization, and target turnover in TMDD models.

Advanced 11 min read

What You'll Learn

  • Relate binding, internalization, and target-turnover timescales to informative sampling windows.
  • Explain why early concentrations, troughs, and target measurements contribute different information.
  • Select dose levels that probe below, near, and above the drug–target affinity range.
  • Use simple calculations to propose sampling schedules and assess their assumptions.
  • Recognize design choices that leave TMDD processes difficult to distinguish.

Why Sampling Design Matters in TMDD

A TMDD model can describe several processes that all change measured drug concentrations: drug binds target, drug–target complex is internalized, and target is synthesized and removed. A concentration decline by itself may not reveal which process produced it. Sampling design should therefore give these mechanisms opportunities to make distinguishable predictions, rather than concentrating measurements in a time window where their effects look similar.

The earlier material on identifiability of binding, internalization, and target turnover motivates this design problem: even a correctly specified model may not support separate parameter estimates if the observations contain little information to distinguish the parameters. Here the practical question is where in time, at what exposure levels, and for which analytes to collect data. The schedule should be based on plausible parameter ranges and the intended assay, and should be refined for the route of administration and study constraints.

Mechanisms and Their Informative Time Windows

For a simple shared-compartment model, let \(C\) be free drug concentration, \(R\) free target concentration, and \(B\) drug–target complex concentration. Let drug enter at amount-per-time rate \(u(t)\) into a compartment of volume \(V\). Assume linear elimination of free drug with rate constant \(k_{\mathrm{el}}\), reversible binding with association rate constant \(k_{\mathrm{on}}\) and dissociation rate constant \(k_{\mathrm{off}}\), target synthesis at rate \(k_{\mathrm{syn}}\), free-target degradation at rate \(k_{\mathrm{deg}}\), and irreversible complex internalization at rate \(k_{\mathrm{int}}\). For clarity, assume no target recycling and no complex-mediated return of drug.

Model and assumptions: In this illustrative model, all concentrations refer to the same compartment, and the rate constants are constant. The drug input \(u(t)\) has units of amount per time, so \(u(t)/V\) has units of concentration per time. The model is
$$ \begin{aligned} \frac{dC}{dt} &= \frac{u(t)}{V}-k_{\mathrm{el}}C-k_{\mathrm{on}}CR+k_{\mathrm{off}}B,\\ \frac{dR}{dt} &= k_{\mathrm{syn}}-k_{\mathrm{deg}}R-k_{\mathrm{on}}CR+k_{\mathrm{off}}B,\\ \frac{dB}{dt} &= k_{\mathrm{on}}CR-(k_{\mathrm{off}}+k_{\mathrm{int}})B. \end{aligned} $$
Here \(C\), \(R\), and \(B\) have concentration units; \(k_{\mathrm{on}}\) has units of concentration\(^{-1}\) time\(^{-1}\); and the other rate constants have units of time\(^{-1}\). The synthesis rate has units of concentration per time. This simplified model is a design aid, not a claim that every study has these exact compartments or pathways.

When free drug is approximately constant over a short interval, the binding complex has an approximate relaxation rate \(k_{\mathrm{on}}C+k_{\mathrm{off}}+k_{\mathrm{int}}\), conditional on treating \(C\) as fixed. If internalization is negligible during binding equilibration, the corresponding binding rate is approximately \(k_{\mathrm{on}}C+k_{\mathrm{off}}\). These are local timescale guides: in a full model, \(C\), \(R\), and \(B\) change together. Binding information is generally concentrated early, when drug and target first encounter one another and concentrations can change rapidly.

Internalization acts on complex over a longer interval when \(B\) persists. Its importance depends not only on \(k_{\mathrm{int}}\), but also on competition with dissociation: in this simple model, the probability that a complex is internalized rather than dissociating during one complex episode is \(k_{\mathrm{int}}/(k_{\mathrm{off}}+k_{\mathrm{int}})\). Thus, a late fall in drug concentration alone may be hard to attribute to internalization without early information on binding and measurements that constrain target or complex behavior.

In the absence of drug and complex, the target approaches its baseline \(R_0=k_{\mathrm{syn}}/k_{\mathrm{deg}}\), with recovery governed by \(k_{\mathrm{deg}}\). Consequently, target measurements after drug and complex have become negligible can reveal turnover on its own characteristic timescale. Drug troughs also matter: concentrations near the end of a dosing interval show whether target-mediated loss remains important at low exposure and whether the target has had time to recover.

Deriving a Practical Sampling Principle

A useful design principle follows from the shape of a simple target-recovery trajectory. Suppose drug and complex are negligible after time zero, and target concentration at that time is \(R_i\). The model reduces to \(dR/dt=k_{\mathrm{syn}}-k_{\mathrm{deg}}R\). Since \(R_0=k_{\mathrm{syn}}/k_{\mathrm{deg}}\), its solution is

$$ R(t)=R_0+(R_i-R_0)e^{-k_{\mathrm{deg}}t}. $$

If \(R_0\) and \(R_i\) are known for this calculation, the sensitivity of the predicted concentration to \(k_{\mathrm{deg}}\) is

$$ \frac{\partial R(t)}{\partial k_{\mathrm{deg}}} =-(R_i-R_0)t e^{-k_{\mathrm{deg}}t}. $$

The absolute sensitivity is proportional to \(t e^{-k_{\mathrm{deg}}t}\). Differentiating this factor gives \(e^{-k_{\mathrm{deg}}t}(1-k_{\mathrm{deg}}t)\), which is zero at \(t=1/k_{\mathrm{deg}}\), positive before that time, and negative after it. Thus, under these simplifying assumptions, a target sample near one turnover time can be more informative about \(k_{\mathrm{deg}}\) than a sample immediately after drug removal. This does not prescribe a universal time point: uncertainty in the rate, measurement error, and lingering drug or complex all affect the useful schedule.

The same logic supports spreading measurements across processes. Very early drug concentrations can reveal rapid concentration changes; intermediate concentrations can inform complex formation and loss; troughs describe end-of-interval behavior; and later target measurements can capture recovery. A schedule should include multiple windows rather than trying to infer all rates from one dense cluster.

Worked Example: Placing Samples Around Plausible Timescales

Consider a hypothetical IV study with plausible values \(k_{\mathrm{on}}=0.0200\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\), \(k_{\mathrm{off}}=0.100\ \mathrm{h}^{-1}\), \(k_{\mathrm{int}}=0.0500\ \mathrm{h}^{-1}\), and \(k_{\mathrm{deg}}=0.0100\ \mathrm{h}^{-1}\). Suppose free drug is expected to be about \(10.0\ \mathrm{nM}\) shortly after dosing. Propose sampling windows, without treating the estimates as known truths.

At \(C=10.0\ \mathrm{nM}\), the approximate binding relaxation rate ignoring internalization is \(k_{\mathrm{on}}C+k_{\mathrm{off}}=(0.0200)(10.0)+0.100=0.300\ \mathrm{h}^{-1}\). Its reciprocal is \(1/0.300=3.33\ \mathrm{h}\). Including internalization in the fixed-\(C\) complex rate gives \(0.300+0.0500=0.350\ \mathrm{h}^{-1}\), or a reciprocal of \(2.86\ \mathrm{h}\). These calculations suggest several early samples, for example at 0.25, 1, and 3 hours, rather than a first sample only at 24 hours.

The internalization half-time is \(\ln(2)/k_{\mathrm{int}}=0.693/0.0500=13.9\ \mathrm{h}\). Samples later in the first day and near the next dose can therefore help characterize persistence and trough behavior. The target-turnover half-time is \(0.693/0.0100=69.3\ \mathrm{h}\), suggesting target measurements over multiple days, including a later point around \(1/k_{\mathrm{deg}}=100\ \mathrm{h}\) if feasible.

A reasonable initial schedule might include baseline; 0.25, 1, and 3 hours; 12 and 24 hours; the predose trough; and target measurements at baseline and approximately 48, 96, and 168 hours. These times span the assumed scales but do not guarantee identifiability. If absorption is not instantaneous, very early sampling must also resolve absorption; a design for subcutaneous dosing cannot simply reuse an IV schedule.

Dose Levels and What to Measure

Time sampling alone is not enough if every dose produces nearly the same degree of target occupancy. A useful dose-ranging design includes exposures below, near, and above the affinity scale. Under a rapid-equilibrium approximation with total target \(T=R+B\), negligible internalization during equilibration, and fixed free drug \(C\), \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\), and

$$ B=\frac{TC}{K_D+C}, \qquad \frac{B}{T}=\frac{C}{K_D+C}. $$

This relation is a design guide, not a substitute for fitting a full dynamic model. It shows why observations only at \(C\gg K_D\) may be uninformative about the exact affinity: occupancy is close to one across a broad range of \(K_D\). Measurements spanning the transition around \(K_D\), alongside lower exposures, can better reveal the concentration dependence. Because TMDD makes dose-to-concentration relationships nonlinear, dose levels should be selected using plausible predicted exposures, not by assuming that doubling dose doubles concentration.

Worked Example: Choosing Exposure Levels Around Affinity

Suppose a hypothetical candidate has \(K_D=5.00\ \mathrm{nM}\), and a dose-ranging design is expected to produce free-drug concentrations of 1.00, 5.00, and 25.0 nM near the relevant observation window. Use the rapid-equilibrium relation to compare expected bound fractions.

At \(1.00\ \mathrm{nM}\), \(B/T=1.00/(5.00+1.00)=0.167\), or 16.7%. At \(5.00\ \mathrm{nM}\), \(B/T=5.00/(5.00+5.00)=0.500\), or 50.0%. At \(25.0\ \mathrm{nM}\), \(B/T=25.0/(5.00+25.0)=0.833\), or 83.3%.

The three exposures therefore probe low, transitional, and high occupancy. If all subjects instead had concentrations near 25 nM, the observed occupancy would be high throughout and provide less contrast for locating \(K_D\). In an actual TMDD study, these are hypothetical exposure targets: dose selection must account for nonlinear disposition, variability, safety, and feasible assay ranges.

The analyte matters as much as the time. Free drug, total drug, free target, and total target do not provide interchangeable observations. Early free-drug measurements can show the changing concentration available for binding. Total drug may include drug in complex, depending on the assay. A target assay that measures only free target can reflect occupancy as well as changes in target abundance; a total-target measure can help distinguish target binding from target depletion if the assay and model definitions align. Collecting paired drug and target observations at selected times can be more informative than collecting many drug concentrations alone.

Worked Example: Scheduling Target-Recovery Measurements

In a hypothetical study, suppose baseline target is \(R_0=10.0\ \mathrm{nM}\), target is \(R_i=4.00\ \mathrm{nM}\) at the point when drug and complex become negligible, and \(k_{\mathrm{deg}}=0.0100\ \mathrm{h}^{-1}\). Calculate the expected target concentration 24 and 72 hours later under the recovery model.

At 24 hours, \(R(24)=10.0+(4.00-10.0)e^{-(0.0100)(24)}=10.0-6.00e^{-0.240}\). Since \(e^{-0.240}=0.7866\), \(R(24)=10.0-4.720=5.280\ \mathrm{nM}\).

At 72 hours, \(R(72)=10.0-6.00e^{-(0.0100)(72)}=10.0-6.00e^{-0.720}\). Since \(e^{-0.720}=0.4868\), \(R(72)=10.0-2.921=7.079\ \mathrm{nM}\). The later observation shows more recovery and helps characterize the time course. This calculation is valid only if drug and complex are negligible throughout recovery; if either remains appreciable, the full target–drug dynamics must be considered.

Building and Comparing Candidate Designs

A practical design can be assembled in stages. First, specify plausible ranges for binding, internalization, turnover, elimination, and absorption parameters. Second, translate those ranges into early, intermediate, trough, and recovery windows. Third, select dose levels whose predicted exposures cover a useful range of target occupancy. Finally, match each measurement to the process it can inform and check feasibility, assay sensitivity, and burden on participants.

1
Anchor the baseline.
Measure predose drug and target so baseline target abundance and background concentrations are not inferred only from postdose data.
2
Resolve early change.
Place more than one early concentration sample near the anticipated binding or absorption timescale. One early sample cannot establish the shape of a rapid transition.
3
Cover complex loss and trough.
Measure during the period when complex-mediated loss may accumulate and near the dosing-interval end, especially across repeated dosing.
4
Observe target recovery.
Continue target measurements after exposure declines, using times that span plausible turnover rather than stopping at the final PK sample.
5
Check separation, not just density.
Compare model predictions across plausible parameter combinations. A design is weak if different binding, internalization, and turnover values produce nearly identical predictions at all planned observations.

A sparse design may be appropriate when burden is severe, but it risks missing rapid binding and offering too few points to distinguish curvature. A dense early PK design is useful for fast changes but may still provide little turnover information if target is never measured later. A dose-ranging design adds exposure contrast, while a repeated-dose design adds trough and recovery information; neither replaces the need for early observations. The strongest feasible design usually combines complementary windows and analytes rather than maximizing sample count in one window.

Design evaluation can formalize this comparison. For example, local sensitivity analysis asks how much predicted observations change when one parameter changes near a nominal value. If two parameters produce nearly proportional changes across all scheduled observations, their effects may be difficult to separate. A Fisher information calculation can summarize this local information, but its conclusions depend on the assumed model, parameter values, and observation-error model. It cannot rescue a design that omits the time window or analyte needed to reveal a mechanism.

Common Pitfalls and Limitations

  • Sampling only at the peak and trough. Those points can characterize exposure endpoints but may miss the early binding transition and intermediate complex dynamics. Add early and intermediate observations when feasible.
  • Using only high doses. Saturating exposure can obscure the concentration range where occupancy changes with drug. Include exposures near and below the plausible affinity scale when justified.
  • Calling free-target decline target turnover. A free-target assay can fall because target is occupied, because target abundance changes, or both. Interpret it with the assay definition and, where possible, paired total-target or drug measurements.
  • Assuming a late sample is automatically informative. If drug or complex persists, the target-recovery equation does not apply. Confirm that the proposed recovery window corresponds to low exposure under plausible model predictions.
  • Overinterpreting a single schedule. Informative times depend on uncertain parameters, route, assay error, and between-subject variability. Consider a range of plausible values and adapt the schedule when early data update those ranges.

A common mistake is to argue, “Binding is rapid, so the first postdose sample is enough to estimate binding.” Rapid binding can instead occur largely before that sample, leaving only an apparent equilibrium or an early concentration change that is confounded with other processes. The correction is to place multiple samples as early as practical, ensure assay timing and absorption are understood, and recognize that if the transition is faster than feasible sampling, binding-rate constants may not be separately estimable even with a well-executed study.

Key takeaway: Distinguish mechanisms by matching measurements to their time scales and exposure conditions: early concentrations for binding dynamics, intermediate and trough observations for persistent complex effects, and target measurements during recovery for turnover. Dose levels should span the informative occupancy range, and the analyte must match the biological quantity being interpreted.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. If \(k_{\mathrm{deg}}=0.0200\ \mathrm{h}^{-1}\), what are the target-turnover half-time and the approximate time \(1/k_{\mathrm{deg}}\)?
  2. Why might a study with concentrations always far above \(K_D\) provide weak information about affinity?
  3. What additional interpretation is needed when a target assay measures free target rather than total target?
  4. For \(K_D=8.00\ \mathrm{nM}\), calculate the rapid-equilibrium bound fraction at \(C=8.00\ \mathrm{nM}\).
  5. Why can a sample at 24 hours fail to measure target recovery even if target is measured at that time?

Solutions. (1) The half-time is \(\ln(2)/0.0200=34.7\ \mathrm{h}\), and \(1/k_{\mathrm{deg}}=50.0\ \mathrm{h}\). (2) When \(C\gg K_D\), \(B/T=C/(K_D+C)\) is close to one across a range of \(K_D\) values, so occupancy changes little as affinity changes. (3) A decline in free target can reflect occupancy as well as a change in total target abundance; interpret it using the assay definition and, if available, paired total-target or drug data. (4) \(B/T=8.00/(8.00+8.00)=0.500\), or 50.0%. (5) Drug or complex may still be present and influencing target, so the drug-free recovery model may not apply; timing should be judged against plausible exposure and complex persistence.

Clinical Biostats