Why Clearance Depends on Concentration
In a conventional linear clearance model, a fixed fraction of drug is eliminated per unit time, so the elimination rate rises in proportion to concentration. Target-mediated drug disposition (TMDD) adds a different route: drug binds to a finite pool of target, and some drug-bound target complexes are internalized. When target is not saturated, this route can make a substantial contribution to elimination. As drug concentration rises, more target becomes occupied, but the amount of complex—and therefore the target-mediated elimination rate—cannot increase without limit if the available target pool is finite.
The resulting pattern is often described as saturable clearance: concentration-normalized clearance can be high at low concentrations and lower at higher concentrations. This description requires care. Binding by itself does not eliminate drug, and a saturated target route does not make all elimination linear. Whether total disposition looks approximately linear at high concentrations depends especially on whether a parallel, nonsaturable elimination pathway is present.
From Binding to Elimination
The full TMDD mass-balance equations distinguish free drug, free target, and complex. For the drug and complex states, a useful form is
Here \(u(t)\) is drug input expressed as concentration per unit time. The association term transfers drug from the free state into the complex; dissociation transfers it back. Neither transfer alone is elimination. Internalization, represented by \(k_{\mathrm{int}}C\), removes complex-bound drug from the modeled drug system. Adding the two equations makes that distinction explicit:
The binding and dissociation terms cancel because they only move drug between states. Thus, under these model assumptions, total drug loss occurs through free-drug elimination and complex internalization. This addition is a short mass-balance justification for using \(k_{\mathrm{int}}C\), rather than the binding rate \(k_{\mathrm{on}}DR\), as the target-mediated elimination rate.
If the total target pool at a given time is \(R_T=R+C\), then \(C\leq R_T\). Consequently,
This is an instantaneous upper bound for a specified target pool, not necessarily a fixed capacity over the entire course of treatment. Target synthesis, degradation, recycling, or changes in target abundance can change \(R_T\). The bound also does not imply that the target-mediated rate is always at its maximum at high drug concentration; the complex must form, and target kinetics may be too slow for occupancy to track a rapid concentration change.
What Changes Across the Concentration Range?
When drug concentration is low relative to the target-binding range, a substantial fraction of target may remain free. Association can then increase as \(D\) rises, because the association rate \(k_{\mathrm{on}}DR\) depends on both free drug and available free target. If complex internalization is appreciable, this can create an additional drug-loss route. In this low-concentration region, the target route may contribute strongly to the total elimination rate per unit of free drug.
As concentration increases, a finite target pool can become increasingly occupied. Complex concentration then has less room to increase, so target-mediated elimination grows less rapidly than free-drug concentration and approaches a capacity constrained by the available target pool. The ratio \(k_{\mathrm{int}}C/D\) can therefore decline even while the absolute target-mediated elimination rate is still rising.
At sufficiently high concentrations, an unsaturated linear pathway may dominate total elimination. In that case, total elimination can become approximately proportional to \(D\), even though the target-mediated component remains saturated. This is an approximate description of the combined pathways, not evidence that the target-mediated route itself has become linear.
Worked Example: Separating Binding Flux From Elimination
Consider a hypothetical instant at which \(D=2.0\) nM, \(R=8.0\) nM, and \(C=2.0\) nM. Let \(k_{\mathrm{on}}=0.10\) nM\(^{-1}\)h\(^{-1}\), \(k_{\mathrm{off}}=0.30\) h\(^{-1}\), and \(k_{\mathrm{int}}=0.20\) h\(^{-1}\). Ignore linear elimination for this example. Calculate the association flux and the target-mediated elimination rate.
The association flux is \(k_{\mathrm{on}}DR=(0.10)(2.0)(8.0)=1.60\) nM/h. The target-mediated elimination rate is instead \(k_{\mathrm{int}}C=(0.20)(2.0)=0.400\) nM/h.
The two rates differ because association creates complex, while only internalization removes complex-bound drug from this model. Complex is also dissociating at rate \(k_{\mathrm{off}}C=(0.30)(2.0)=0.600\) nM/h. At this instant, the complex formation rate minus dissociation and internalization is \(1.60-0.600-0.400=0.600\) nM/h, matching the positive value of \(dC/dt\). It would be incorrect to report the association flux of 1.60 nM/h as drug elimination.
Worked Example: Apparent Clearance at Low and High Concentrations
Suppose two hypothetical model snapshots have the following states. In both, the instantaneous total target pool \(R_T=R+C\) is 10 nM. Use \(k_{\mathrm{el}}=0.050\) h\(^{-1}\) and \(k_{\mathrm{int}}=0.40\) h\(^{-1}\). The low-concentration snapshot has \(D=1.0\) nM, \(R=9.0\) nM, and \(C=1.0\) nM. The high-concentration snapshot has \(D=100\) nM, \(R=1.0\) nM, and \(C=9.0\) nM. Find the elimination rates and concentration-based apparent clearances. These are illustrative states, not a claim that they are equilibrium values.
At the low-concentration snapshot, the linear elimination rate is \(k_{\mathrm{el}}D=(0.050)(1.0)=0.050\) nM/h. Target-mediated elimination is \(k_{\mathrm{int}}C=(0.40)(1.0)=0.400\) nM/h. Total elimination is \(0.050+0.400=0.450\) nM/h, so apparent clearance is \(0.450/1.0=0.450\) h\(^{-1}\).
At the high-concentration snapshot, linear elimination is \((0.050)(100)=5.00\) nM/h, and target-mediated elimination is \((0.40)(9.0)=3.60\) nM/h. Total elimination is \(8.60\) nM/h, so apparent clearance is \(8.60/100=0.0860\) h\(^{-1}\). The target-mediated rate is approaching its instantaneous bound of \((0.40)(10)=4.00\) nM/h, while the linear rate continues to rise with concentration. The combined apparent clearance is lower in the high-concentration snapshot; total elimination itself is not lower.
This calculation also shows why “clearance falls” must not be confused with “elimination rate falls.” Clearance is a rate normalized by concentration. Here total elimination rises from 0.450 to 8.60 nM/h while the normalized value falls from 0.450 to 0.0860 h\(^{-1}\).
Why Internalization and Dissociation Matter
A large amount of binding does not, by itself, establish a large elimination contribution. A complex can dissociate and return drug to the free pool, or it can be internalized. To illustrate the competing exits from a complex, suppose that once a complex has formed, its only possible exits are dissociation at rate \(k_{\mathrm{off}}\) and internalization at rate \(k_{\mathrm{int}}\). Under constant rates, the probability that internalization is the first of those two events is
This follows from the exponential waiting-time model: the total exit rate is \(k_{\mathrm{off}}+k_{\mathrm{int}}\), and the internalization event has rate \(k_{\mathrm{int}}\). The expression describes a single complex’s competing exit event; it is not a substitute for the full model, which also includes repeated binding, target turnover, and changing concentrations.
Worked Example: The Same Binding, Different Elimination Potential
Compare two hypothetical complexes with the same association conditions and \(C=5.0\) nM. In case A, \(k_{\mathrm{off}}=0.90\) h\(^{-1}\) and \(k_{\mathrm{int}}=0.10\) h\(^{-1}\). In case B, \(k_{\mathrm{off}}=0.10\) h\(^{-1}\) and \(k_{\mathrm{int}}=0.90\) h\(^{-1}\). Calculate the instantaneous target-mediated elimination rates and the internalization fractions for a single complex.
For case A, the elimination rate is \((0.10)(5.0)=0.500\) nM/h, and the internalization fraction is \(0.10/(0.90+0.10)=0.100\). For case B, the elimination rate is \((0.90)(5.0)=4.50\) nM/h, and the internalization fraction is \(0.90/(0.10+0.90)=0.900\).
The same complex concentration can therefore correspond to very different target-mediated drug loss when internalization and dissociation rates differ. Occupancy or complex abundance alone does not determine the elimination rate; \(k_{\mathrm{int}}\) is also required.
Assumptions, Diagnostics, and Alternatives
The saturation interpretation rests on a mechanistic structure, not merely on observing a curved concentration–time profile. It assumes that the modeled drug and target pools are correctly specified; association and dissociation move drug between free and bound states; internalization removes complex-bound drug as represented; and the relevant rate constants and target abundance adequately describe the time period being interpreted. The examples above assume a well-mixed concentration compartment and treat their state values as known. In data analysis, those states and parameters are estimated, often with uncertainty.
Useful diagnostic questions include whether concentration-normalized elimination is systematically greater at low concentrations, whether the pattern is consistent across dose levels and time, and whether a proposed target pool is large enough to account for the inferred target-mediated loss. Examine model predictions and residuals over the full concentration range, not only near the highest concentrations. A model can fit the high-concentration observations while missing a steep low-concentration decline, or vice versa. If measurements distinguish free drug, total drug, or complex, their definitions must also match the modeled state: a denominator based on total drug is not interchangeable with one based on free \(D\).
A linear-clearance model is a useful alternative when the observed range does not show a reproducible concentration dependence or when the data cannot support target-specific parameters. It is simpler, but cannot represent a finite target-mediated route. An empirical saturable-clearance model can summarize a curved concentration–elimination relationship with fewer mechanistic states, but it does not by itself establish that target binding and internalization caused the curvature. The full TMDD model represents those states and processes explicitly, at the cost of greater data and identifiability demands.
A common mistake is to say, “The target is saturated at high concentration, so total drug elimination must stop increasing.” The correction is that saturation constrains the target-mediated component, \(k_{\mathrm{int}}C\), not every elimination pathway. In the worked comparison, target-mediated elimination approached 4.00 nM/h, while the parallel linear component rose to 5.00 nM/h. Conversely, if there were no parallel linear pathway and the target-mediated route were fully saturated, total elimination could approach a plateau rather than remain proportional to concentration. High-concentration linearity is therefore conditional, not guaranteed by target saturation alone.
Other limitations matter in practice. Target abundance may change through synthesis, degradation, or internalization, so a capacity inferred at one time may not hold later. Slow binding or internalization can create delays between concentration and elimination, making a static concentration–clearance comparison misleading. Sampling, assay choice, and model misspecification can also make a true mechanistic pattern hard to identify. A decline in estimated apparent clearance is evidence of concentration dependence; on its own, it is not proof of TMDD.
Binding and dissociation redistribute drug; internalization and other elimination processes remove drug from the modeled system.
For a specified target pool, \(C\leq R_T\), so target-mediated elimination cannot exceed \(k_{\mathrm{int}}R_T\) at that instant.
Compare \(k_{\mathrm{int}}C\) with \(k_{\mathrm{el}}D\), and state which concentration measure is used when reporting apparent clearance.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- If \(k_{\mathrm{on}}DR\) is large, does that alone establish a large drug-elimination rate? Why or why not?
- For \(R_T=12\) nM and \(k_{\mathrm{int}}=0.25\) h\(^{-1}\), what is the largest possible instantaneous target-mediated elimination rate under the specified-pool bound?
- If total elimination rises from 2.0 to 8.0 nM/h while free concentration rises from 5.0 to 40 nM, what happens to concentration-based apparent clearance?
- Why can target saturation coexist with approximately linear total elimination at high concentrations?
- Does a falling apparent clearance prove that TMDD is the cause? State one alternative explanation or limitation.
Solutions. (1) No. \(k_{\mathrm{on}}DR\) is complex formation, not drug loss; target-mediated elimination is \(k_{\mathrm{int}}C\). (2) The bound is \(k_{\mathrm{int}}R_T=(0.25)(12)=3.00\) nM/h. (3) Clearance changes from \(2.0/5.0=0.400\) h\(^{-1}\) to \(8.0/40=0.200\) h\(^{-1}\), so it halves even though the elimination rate rises. (4) Saturation constrains the target-mediated component, while an unsaturated linear pathway can continue to increase in proportion to free concentration and dominate the total. (5) No. Other concentration-dependent processes, changing target abundance, measurement differences, or model misspecification can produce or obscure such a pattern.