From a Reduced TMDD Model to Calculations
A target-mediated drug disposition calculation is useful when a pharmacokinetic model needs to connect antibody concentration to target binding and antibody removal. The key quantities are related but are not interchangeable: free antibody drives binding, complex concentration determines occupancy, and complex internalization contributes an elimination flux. Calculating each quantity explicitly helps reveal which process dominates at a given concentration.
The examples use a reduced quasi-steady-state (QSS) model at a specified free antibody concentration \(D\). They are hypothetical calculations, not observations from a clinical study. The model describes one target pool and one antibody molecule per target molecule in the complex.
This relation follows from treating complex formation as locally quasi-steady: the rate of complex formation, \(k_{\mathrm{on}}DR\), balances its rate of loss, \((k_{\mathrm{off}}+k_{\mathrm{int}})C\). Consequently, \(C=DR/K_{\mathrm{SS}}\). Substituting \(R=R_{\mathrm{tot}}-C\) and solving for \(C\) gives the expression above. Because internalization is included in the loss of complex, \(K_{\mathrm{SS}}\) need not equal the binding equilibrium dissociation constant \(k_{\mathrm{off}}/k_{\mathrm{on}}\).
The calculations below hold \(R_{\mathrm{tot}}\) fixed and assume the QSS relation is appropriate at the concentration being considered. This is a local, reduced description: it does not calculate how target abundance changes over time or resolve an entire dose profile. The limiting-case results from the earlier material on low- and high-concentration TMDD provide useful checks on the arithmetic.
Clearance Convention and Model Assumptions
For the examples, antibody elimination has two components. Baseline elimination is proportional to free antibody concentration, while internalization removes antibody bound in complex. Let \(CL_0\) be baseline clearance in volume per time and \(V\) the volume used to convert the complex internalization rate into an amount-per-time flux. The total elimination flux is
Here \(J\) has units of amount per time when concentrations are expressed as amount per volume. We define the approximate clearance in these examples by dividing the flux by the free antibody concentration:
This is a free-concentration-normalized clearance, with units of volume per time. Its denominator must be kept in mind when interpreting or comparing a calculated value: a clearance normalized to total antibody concentration \(D+C\) would instead be \(J/(D+C)\), and would generally differ. The clearance here is an approximate local quantity, not a claim that concentration is constant or that the system is at whole-body pharmacokinetic steady state.
The calculation assumes that free antibody is the concentration \(D\) supplied to the binding relation; target abundance is fixed over the calculation; the QSS approximation is valid; and the specified baseline and internalization routes account for the elimination being represented. It also assumes the concentration and rate units are compatible. In particular, \(Vk_{\mathrm{int}}C\) and \(CL_0D\) must both be elimination fluxes in the same units.
Worked Calculations Across Three Regimes
Use the following hypothetical parameter set in all three examples: \(R_{\mathrm{tot}}=20\) nM, \(K_{\mathrm{SS}}=10\) nM, \(k_{\mathrm{int}}=0.20\) h\(^{-1}\), \(V=3.0\) L, and \(CL_0=0.30\) L/h. The three free antibody concentrations are \(1\), \(10\), and \(100\) nM. They represent concentrations below, near, and well above \(K_{\mathrm{SS}}\), respectively.
Worked Example: Low Antibody Concentration
At \(D=1.00\) nM, calculate free target, complex, occupancy, elimination flux, and free-concentration-normalized clearance.
The QSS complex concentration is
The remaining free target is \(R=R_{\mathrm{tot}}-C=20.0-1.818=18.18\) nM. Fractional occupancy is \(C/R_{\mathrm{tot}}=1.818/20.0=0.09091\), or \(9.091\%\). This is also \(D/(K_{\mathrm{SS}}+D)=1/11\).
The baseline elimination flux is \(CL_0D=(0.30\ \mathrm{L/h})(1.00\ \mathrm{nM})=0.300\ \mathrm{nM\,L/h}\). The target-mediated flux is
Thus \(J=0.300+1.091=1.391\ \mathrm{nM\,L/h}\), and \(CL_{\mathrm{app,free}}=J/D=1.391\ \mathrm{L/h}\). Although occupancy is only about \(9.1\%\), the target-mediated contribution to the calculated flux exceeds the baseline contribution for these chosen parameters.
Worked Example: Intermediate Antibody Concentration
At \(D=10.0\) nM, the free antibody concentration equals \(K_{\mathrm{SS}}\). The complex concentration is
Therefore \(R=20.0-10.0=10.0\) nM and occupancy is \(10.0/20.0=0.500\), or \(50.0\%\). This equality of free and bound target follows specifically from \(D=K_{\mathrm{SS}}\) in this model; it is not a general property of arbitrary binding systems.
The baseline flux is \((0.30)(10.0)=3.00\ \mathrm{nM\,L/h}\), and the target-mediated flux is \((3.0)(0.20)(10.0)=6.00\ \mathrm{nM\,L/h}\). Hence \(J=9.00\ \mathrm{nM\,L/h}\), and \(CL_{\mathrm{app,free}}=9.00/10.0=0.900\) L/h. If instead the denominator were total antibody concentration, \(D+C=20.0\) nM, the resulting ratio would be \(9.00/20.0=0.450\) L/h. The different values reflect different clearance conventions, not different elimination fluxes.
Worked Example: High Antibody Concentration
At \(D=100\) nM, which is ten times \(K_{\mathrm{SS}}\), the complex concentration is
Free target is \(R=20.0-18.18=1.818\) nM, and occupancy is \(18.18/20.0=0.9091\), or \(90.91\%\). Occupancy is high but not complete: a finite concentration above \(K_{\mathrm{SS}}\) does not make the model’s unbound target exactly zero.
The baseline flux is \((0.30)(100)=30.0\ \mathrm{nM\,L/h}\). The target-mediated flux is \((3.0)(0.20)(18.18)=10.91\ \mathrm{nM\,L/h}\). Thus \(J=40.91\ \mathrm{nM\,L/h}\), and \(CL_{\mathrm{app,free}}=40.91/100=0.4091\) L/h. Clearance is approaching the baseline value of \(0.30\) L/h because the target-mediated flux grows much less than proportionally to free antibody concentration once the target pool is nearly saturated.
Why the Regimes Behave Differently
The occupancy relation makes the transitions precise. When \(D\) is much smaller than \(K_{\mathrm{SS}}\), \(K_{\mathrm{SS}}+D\) is approximately \(K_{\mathrm{SS}}\), so \(C\) is approximately \(R_{\mathrm{tot}}D/K_{\mathrm{SS}}\) and occupancy is approximately \(D/K_{\mathrm{SS}}\). Complex formation and target-mediated elimination flux are then approximately proportional to \(D\). When \(D=K_{\mathrm{SS}}\), occupancy is one half. When \(D\) is much larger than \(K_{\mathrm{SS}}\), \(C\) approaches \(R_{\mathrm{tot}}\) and occupancy approaches one.
The corresponding clearance limits follow directly from the clearance expression. At low concentration, the approximate target-mediated clearance contribution is \(Vk_{\mathrm{int}}R_{\mathrm{tot}}/K_{\mathrm{SS}}\). At high concentration, that contribution tends to zero, so \(CL_{\mathrm{app,free}}\) tends to \(CL_0\). The target-mediated flux, by contrast, approaches \(Vk_{\mathrm{int}}R_{\mathrm{tot}}\): a saturating flux can coexist with a declining clearance contribution because clearance divides flux by \(D\).
For the hypothetical parameters, \(Vk_{\mathrm{int}}R_{\mathrm{tot}}=(3.0)(0.20)(20)=12.0\ \mathrm{nM\,L/h}\). The low-concentration limit of the target-mediated clearance contribution is \(12.0/10.0=1.20\) L/h, making the total low-concentration limit \(1.50\) L/h. The high-concentration limit is \(0.300\) L/h. The worked values, \(1.391\), \(0.900\), and \(0.4091\) L/h, move between those limits in the expected direction.
Checks, Alternatives, and Common Pitfalls
A few model-based checks can catch many implementation errors. For positive \(D\), \(R_{\mathrm{tot}}\), and \(K_{\mathrm{SS}}\), the calculated complex must satisfy \(0\leq C\leq R_{\mathrm{tot}}\), and free target must equal \(R_{\mathrm{tot}}-C\). Occupancy must be between zero and one and increase with \(D\). The target-mediated flux should approach, but not exceed, \(Vk_{\mathrm{int}}R_{\mathrm{tot}}\). With the stated clearance convention, the target-mediated clearance contribution decreases as \(D\) increases, and total clearance approaches \(CL_0\). A calculated complex greater than total target, a negative free target, or clearance increasing indefinitely with concentration signals an algebra, unit, or denominator error.
The reduced QSS calculation is appropriate for a rapid calculation at specified concentration when binding and complex turnover can be represented by a quasi-steady relation and target abundance can be treated as fixed. A full TMDD ordinary differential equation model is preferable when time-varying free target, complex dynamics, target synthesis and degradation, or binding transients are important. The equilibrium binding relation uses \(k_{\mathrm{off}}/k_{\mathrm{on}}\) rather than \(K_{\mathrm{SS}}\); substituting it without justification omits internalization from the QSS complex-loss balance. The Michaelis–Menten reduction discussed in the earlier material on deriving Michaelis–Menten clearance summarizes saturable elimination, but the explicit complex calculation here retains occupancy and free-target information.
Common mistake: “At high antibody concentration the target is saturated, so target-mediated clearance must be high.” Saturation means the complex concentration and its internalization flux approach finite maxima. Under the free-concentration-normalized convention, dividing that bounded flux by an increasingly large \(D\) makes its clearance contribution decrease. Saturated target-mediated flux is not the same as increasing target-mediated clearance.
The numerical result is only as reliable as its assumptions and inputs. A changing target pool, a slow approach to QSS, multiple target compartments, or additional elimination mechanisms can invalidate the simple calculation. Also, an assay reporting total rather than free antibody changes the concentration used for normalization and may change the reported clearance. State which concentration is supplied to the binding relation and which is used in the clearance denominator before comparing results.
Check Your Understanding
Use the hypothetical parameters given above. Answers are included below the questions so that each calculation can be checked.
- At \(D=5.00\) nM, what are \(C\), free target \(R\), and fractional occupancy?
- At \(D=10.0\) nM, calculate the total elimination flux and \(CL_{\mathrm{app,free}}\).
- What is the low-concentration limit of total free-concentration-normalized clearance for these parameters?
- At \(D=100\) nM, is the target-mediated flux at its limiting maximum \(Vk_{\mathrm{int}}R_{\mathrm{tot}}\)? Explain using the calculated value.
Solutions. (1) \(C=(20)(5)/(10+5)=6.667\) nM; \(R=20-6.667=13.33\) nM; occupancy \(=6.667/20=0.3333\), or \(33.33\%\). (2) At \(D=10\) nM, \(C=10\) nM, so \(J=(0.30)(10)+(3.0)(0.20)(10)=3.00+6.00=9.00\ \mathrm{nM\,L/h}\); \(CL_{\mathrm{app,free}}=9.00/10=0.900\) L/h. (3) The limit is \(CL_0+Vk_{\mathrm{int}}R_{\mathrm{tot}}/K_{\mathrm{SS}}=0.30+12.0/10.0=1.50\) L/h. (4) No. The limiting target-mediated flux is \(12.0\ \mathrm{nM\,L/h}\); at \(D=100\) nM it is \(10.91\ \mathrm{nM\,L/h}\), which approaches but remains below the limit because \(C=18.18\) nM is below \(R_{\mathrm{tot}}=20\) nM.