Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › How Target Abundance Changes Nonlinear Drug Exposure

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

How Target Abundance Changes Nonlinear Drug Exposure

Learn to predict how changing baseline target abundance shifts steady-state free-drug exposure and the input range for nonlinear disposition.

Advanced 12 min read

What You'll Learn

  • Derive the steady-state TMDD exposure equation with target synthesis and turnover.
  • Relate baseline target abundance to the maximum target-mediated elimination rate.
  • Calculate free-drug concentrations under different target abundances.
  • Identify the infusion input at which target-mediated elimination is half-saturated.
  • Distinguish target-abundance effects from changes in target turnover or binding kinetics.
  • Recognize when the model predicts no finite steady-state concentration.

Why Baseline Target Abundance Matters

A drug can be cleared both by nonspecific elimination and by binding to a target followed by complex internalization. When the target is scarce relative to drug exposure, this target-mediated pathway may be close to saturation. When target abundance is higher, more drug can be captured and eliminated through the target pathway before that saturation is reached. The result can be a substantial change in free-drug exposure, even when the drug’s nonspecific clearance and binding kinetics are unchanged.

The comparison in this tutorial is mechanistic: it asks what the model predicts when baseline target abundance is changed while other specified parameters are held fixed. It is not, by itself, a model of between-subject variability. We use a constant infusion to make the relationship between input and steady-state exposure explicit; a bolus dose and its peak and trough concentrations require consideration of time-varying dynamics.

Steady-State TMDD Model and Its Assumptions

Let \(C\) be free-drug concentration, \(R\) free-target concentration, and \(X\) drug-target complex concentration. The input \(u\) is a constant infusion rate expressed per volume, with units of concentration per time. The model includes linear elimination of free drug, target synthesis and degradation, reversible binding, and internalization of the complex:

$$ \frac{dC}{dt}=u-k_{\mathrm{el}}C-k_{\mathrm{on}}CR+k_{\mathrm{off}}X,\qquad \frac{dR}{dt}=k_{\mathrm{syn}}-k_{\mathrm{deg}}R-k_{\mathrm{on}}CR+k_{\mathrm{off}}X,\qquad \frac{dX}{dt}=k_{\mathrm{on}}CR-(k_{\mathrm{off}}+k_{\mathrm{int}})X. $$

Here \(k_{\mathrm{el}}\), \(k_{\mathrm{deg}}\), and \(k_{\mathrm{int}}\) are rate constants; \(k_{\mathrm{on}}\) and \(k_{\mathrm{off}}\) are association and dissociation rate constants; and \(k_{\mathrm{syn}}\) is the target synthesis rate. The baseline target abundance \(T_0\) is the steady free-target concentration in the absence of drug, so \(T_0=k_{\mathrm{syn}}/k_{\mathrm{deg}}\).

Assumptions: The derivation below assumes a single well-mixed compartment, constant infusion, positive \(k_{\mathrm{deg}}\), \(k_{\mathrm{int}}\), and \(k_{\mathrm{on}}\), and a steady state with no target recycling. Free drug has linear nonspecific elimination. The equations treat internalization as removal of both complex and its associated drug. The conclusions concern this model and its assumptions, not every possible target-disposition mechanism.

At steady state, the complex equation gives \(X=k_{\mathrm{on}}CR/(k_{\mathrm{off}}+k_{\mathrm{int}})\). Define \(K_{\mathrm{ss}}=(k_{\mathrm{off}}+k_{\mathrm{int}})/k_{\mathrm{on}}\), so \(X=CR/K_{\mathrm{ss}}\). This relationship follows from the steady-state complex balance; it does not require that binding be at rapid equilibrium.

Deriving the Exposure–Input Relationship

At steady state, the free-target balance can be combined with the complex balance. Binding and dissociation cancel in the sum of the target and complex equations, leaving \(k_{\mathrm{syn}}=k_{\mathrm{deg}}R+k_{\mathrm{int}}X\). Thus the target pool available for complex internalization is replenished by synthesis and depleted by degradation of free target.

$$ R=T_0-\frac{k_{\mathrm{int}}}{k_{\mathrm{deg}}}X,\qquad X=\frac{CR}{K_{\mathrm{ss}}}. $$

Substitute the second relation into the first and solve for \(R\). For positive rates, define \(K_m=k_{\mathrm{deg}}K_{\mathrm{ss}}/k_{\mathrm{int}}\). The resulting expressions show how baseline abundance enters the steady-state target pool and the internalization rate:

$$ R=\frac{T_0K_m}{K_m+C},\qquad k_{\mathrm{int}}X=\frac{k_{\mathrm{deg}}T_0C}{K_m+C} =\frac{k_{\mathrm{syn}}C}{K_m+C}. $$

The drug balance at steady state is \(u=k_{\mathrm{el}}C+k_{\mathrm{int}}X\). Therefore the input–exposure relationship is:

$$ u=k_{\mathrm{el}}C+\frac{k_{\mathrm{syn}}C}{K_m+C}, \qquad V_{\max,\mathrm{TM}}=k_{\mathrm{syn}}=k_{\mathrm{deg}}T_0. $$

The maximum target-mediated elimination rate in this turnover model is \(k_{\mathrm{syn}}\), not \(k_{\mathrm{int}}T_0\). At high free-drug concentration, the available target is increasingly depleted into complex, and the internalization rate approaches the rate at which target is supplied. This distinction matters when target abundance is changed by altering synthesis while degradation remains fixed.

Uniqueness and the Nonlinear Input Range

For \(k_{\mathrm{el}}>0\), the right-hand side of the steady-state balance is zero at \(C=0\), increases continuously with \(C\), and tends to infinity as \(C\) increases. Its derivative is \(k_{\mathrm{el}}+k_{\mathrm{syn}}K_m/(K_m+C)^2\), which is positive. Consequently, for every positive input \(u\), there is exactly one positive steady-state concentration.

Multiplying the balance by \(K_m+C\) gives a quadratic that can be used to calculate that concentration:

$$ k_{\mathrm{el}}C^2+ \left(k_{\mathrm{el}}K_m+k_{\mathrm{syn}}-u\right)C-uK_m=0. $$

The positive root is obtained by the quadratic formula. A useful landmark is \(C=K_m\), where target-mediated elimination is half its maximum. Substitution into the balance gives \(u=k_{\mathrm{el}}K_m+k_{\mathrm{syn}}/2\). Thus, if \(k_{\mathrm{deg}}\) and binding and internalization parameters are held fixed, increasing \(T_0\) increases \(k_{\mathrm{syn}}\) and shifts this input landmark upward. The concentration landmark \(K_m\), however, does not change under those particular conditions.

Interpretation: “The nonlinear range shifts to higher input” means that a larger constant infusion rate is needed to reach the same degree of target-pathway saturation when target abundance is higher and the other stated parameters are fixed. It does not mean that concentration is proportional to dose throughout the nonlinear range.

Worked Example: Comparing Free-Drug Exposure at Two Target Abundances

Consider hypothetical steady-state parameters \(k_{\mathrm{el}}=0.100\ \mathrm{h}^{-1}\), \(k_{\mathrm{deg}}=0.0200\ \mathrm{h}^{-1}\), \(K_{\mathrm{ss}}=5.00\ \mathrm{nM}\), \(k_{\mathrm{int}}=0.100\ \mathrm{h}^{-1}\), and \(u=0.500\ \mathrm{nM/h}\). Then \(K_m=(0.0200)(5.00)/(0.100)=1.00\ \mathrm{nM}\). Compare \(T_0=50.0\ \mathrm{nM}\) and \(T_0=100\ \mathrm{nM}\), with \(k_{\mathrm{deg}}\) fixed.

For \(T_0=50.0\ \mathrm{nM}\), \(k_{\mathrm{syn}}=k_{\mathrm{deg}}T_0=1.00\ \mathrm{nM/h}\). The balance is \(0.100C+C/(1.00+C)=0.500\). The positive quadratic root is:

$$ C=\frac{0.500-0.100-1.00+ \sqrt{(0.100+1.00-0.500)^2+4(0.100)(0.500)(1.00)}}{2(0.100)} =0.7417\ \mathrm{nM}. $$

Checking the balance with the unrounded concentration, nonspecific elimination is \(0.100(0.7417)=0.07417\ \mathrm{nM/h}\), and target-mediated elimination is approximately \(1.00(0.7417)/(1.00+0.7417)=0.4258\ \mathrm{nM/h}\). Their sum is \(0.5000\ \mathrm{nM/h}\), up to rounding.

For \(T_0=100\ \mathrm{nM}\), \(k_{\mathrm{syn}}=2.00\ \mathrm{nM/h}\). The positive root of \(0.100C+2.00C/(1.00+C)=0.500\) is \(C=0.3066\ \mathrm{nM}\). The higher target abundance therefore predicts lower free-drug concentration at the same infusion input. In this model, the higher synthesis rate supports a larger target-mediated elimination flux at that input.

Worked Example: Locating the Half-Saturation Input

Use the same hypothetical kinetic parameters as above. At \(C=K_m=1.00\ \mathrm{nM}\), the target-mediated elimination rate is exactly half of \(k_{\mathrm{syn}}\). For \(T_0=50.0\ \mathrm{nM}\), the input at this concentration is:

$$ u=k_{\mathrm{el}}K_m+\frac{k_{\mathrm{syn}}}{2} =(0.100)(1.00)+\frac{1.00}{2} =0.600\ \mathrm{nM/h}. $$

For \(T_0=100\ \mathrm{nM}\), \(k_{\mathrm{syn}}=2.00\ \mathrm{nM/h}\), so the corresponding input is \(0.100+2.00/2=1.10\ \mathrm{nM/h}\). Both cases have the same concentration landmark, \(1.00\ \mathrm{nM}\), because \(K_m\) was held fixed. But the higher-abundance case requires a larger input to reach it. This calculation identifies one clear saturation landmark, not a universal boundary separating linear from nonlinear disposition.

Worked Example: Exposure Above the Target-Mediated Capacity

Retain the parameters above, set \(u=3.00\ \mathrm{nM/h}\), and compare \(T_0=50.0\ \mathrm{nM}\) with \(T_0=100\ \mathrm{nM}\). For the lower abundance, \(k_{\mathrm{syn}}=1.00\ \mathrm{nM/h}\), and the positive root is:

$$ C=\frac{3.00-0.100-1.00+ \sqrt{(0.100+1.00-3.00)^2+4(0.100)(3.00)(1.00)}}{0.200} =20.47\ \mathrm{nM}. $$

At this concentration, nonspecific elimination is approximately \(2.047\ \mathrm{nM/h}\), and target-mediated elimination is approximately \(20.47/(1.00+20.47)=0.9534\ \mathrm{nM/h}\). The rates sum to \(3.000\ \mathrm{nM/h}\), allowing for rounding.

For \(T_0=100\ \mathrm{nM}\), \(k_{\mathrm{syn}}=2.00\ \mathrm{nM/h}\). Solving \(0.100C+2.00C/(1.00+C)=3.00\) gives \(C=11.59\ \mathrm{nM}\). The target pathway is close to its maximum in both cases, but the increased target abundance still lowers free-drug exposure at this input. With \(k_{\mathrm{el}}>0\), input above the target-mediated maximum does not prevent a finite steady state: the excess input is balanced by increasing nonspecific elimination.

Diagnostics, Alternatives, and Limitations

A useful check is to calculate the two elimination components separately and verify that they sum to the imposed input. Also verify that the computed free-target and complex concentrations are nonnegative and satisfy the steady-state target balance. A mismatch can arise from rounding, inconsistent units, or mixing the maximum rate from a fixed-pool approximation with the turnover model used here.

The equation shows two distinct ways that a changed \(T_0\) could accompany a changed exposure curve. If \(T_0\) is increased by increasing \(k_{\mathrm{syn}}\) while \(k_{\mathrm{deg}}\) and the binding and internalization rates are fixed, \(k_{\mathrm{syn}}\) and the target-mediated capacity increase, while \(K_m\) stays fixed. If \(T_0\) changes because \(k_{\mathrm{deg}}\) changes while synthesis is held fixed, then \(k_{\mathrm{syn}}=k_{\mathrm{deg}}T_0\) stays fixed but \(K_m=k_{\mathrm{deg}}K_{\mathrm{ss}}/k_{\mathrm{int}}\) changes. It is therefore not enough to know that baseline abundance differs; the mechanism of that difference matters.

The full dynamic TMDD equations describe transient binding, target recovery, and complex accumulation. The steady-state balance here is an alternative when the question is specifically the long-run infusion input–exposure relationship. A fixed-total-target saturation approximation can instead be useful over a time scale on which target abundance is treated as effectively constant; in that approximation, the maximum internalization rate is often written \(k_{\mathrm{int}}T_{\mathrm{tot}}\). That is not the maximum rate derived for a synthesizing and degrading target at steady state, where the capacity is \(k_{\mathrm{syn}}\). Rapid-equilibrium or Michaelis–Menten approximations may further simplify binding, but they require their own assumptions and should not be substituted without checking the relevant rate relationships.

This steady-state result has limits. It does not predict the timing of the concentration transition after a dose, transient target depletion, or peak and trough concentrations under repeated bolus dosing. Further, if \(k_{\mathrm{el}}=0\), a finite steady state exists only for inputs below the target-mediated capacity when that capacity is positive. Specifically, for \(k_{\mathrm{syn}}>0\), \(u=k_{\mathrm{syn}}C/(K_m+C)\) has a finite solution only for \(0\leq u<k_{\mathrm{syn}}\); at \(u=k_{\mathrm{syn}}\), the rate is approached only as \(C\) tends to infinity. If \(k_{\mathrm{syn}}=0\) as well, the target-mediated capacity is zero: at \(u=0\), every \(C\geq0\) balances when \(k_{\mathrm{el}}=0\), while positive input has no steady-state balance. These edge cases differ from the positive-\(k_{\mathrm{el}}\) examples above.

Common Pitfall

  • Wrong reasoning: “Doubling baseline target abundance must double the target-mediated maximum elimination rate \(k_{\mathrm{int}}T_0\), so the exposure curve follows that capacity.”
  • Correction: In the turnover model used here, steady-state target synthesis is \(k_{\mathrm{syn}}=k_{\mathrm{deg}}T_0\), and the maximum target-mediated elimination rate is \(k_{\mathrm{syn}}\). The expression \(k_{\mathrm{int}}T_0\) belongs to a different fixed-pool approximation and is not generally the steady-state capacity when target is synthesized, degraded, and depleted through complex internalization.
  • Wrong reasoning: “A higher target abundance always shifts the concentration at half-saturation.”
  • Correction: Under the controlled change where \(k_{\mathrm{deg}}\), \(K_{\mathrm{ss}}\), and \(k_{\mathrm{int}}\) stay fixed, \(K_m\) is unchanged. The input needed to reach \(C=K_m\) shifts because \(k_{\mathrm{syn}}\) changes.
Key takeaway: In the steady-state turnover model, target abundance affects free-drug exposure through both the target synthesis-supported elimination capacity and the input needed to reach a given concentration. Interpret the shift only after specifying which biological parameters are held fixed.

Check Your Understanding

Answer each question, then compare with the short solutions below.

  1. If \(k_{\mathrm{deg}}\) is fixed and \(T_0\) doubles by a change in \(k_{\mathrm{syn}}\), what happens to the maximum target-mediated elimination rate?
  2. At what free-drug concentration is target-mediated elimination half of its maximum in the derived steady-state model?
  3. With \(k_{\mathrm{el}}>0\), does input above \(k_{\mathrm{syn}}\) rule out a finite steady-state concentration? Why?
  4. If \(T_0\) changes because \(k_{\mathrm{deg}}\) changes while \(k_{\mathrm{syn}}\) is fixed, which quantities in the exposure relationship may change?
  5. When \(k_{\mathrm{el}}=0\) and \(k_{\mathrm{syn}}>0\), can the model have a finite steady state at \(u=k_{\mathrm{syn}}\)?

Solutions. (1) It doubles, because the maximum is \(k_{\mathrm{syn}}=k_{\mathrm{deg}}T_0\) when \(k_{\mathrm{deg}}\) is fixed. (2) At \(C=K_m\), since \(C/(K_m+C)=1/2\). (3) Yes. Increasing \(C\) raises the linear elimination term without bound, so it can balance any finite positive input. (4) \(K_m=k_{\mathrm{deg}}K_{\mathrm{ss}}/k_{\mathrm{int}}\) may change, while the capacity \(k_{\mathrm{syn}}\) remains fixed. (5) No; the target pathway approaches \(k_{\mathrm{syn}}\) only as \(C\) tends to infinity, so there is no finite solution at equality.

Clinical Biostats