Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Infusion Dose-Exposure Relationships Under TMDD

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

Infusion Dose-Exposure Relationships Under TMDD

Learn to calculate free and total steady-state drug concentrations across infusion rates and interpret the effects of target saturation on exposure and apparent clearance.

Advanced 12 min read

What You'll Learn

  • Derive the infusion mass balance for a reduced TMDD model with a maintained target pool.
  • Solve for steady-state free drug concentration and calculate bound and total drug.
  • Calculate infusion rate from a desired target occupancy or concentration.
  • Distinguish free-concentration-based effective clearance from total-concentration-based apparent clearance.
  • Explain why increasing infusion rates can produce more-than-proportional increases in free drug concentration.
  • Recognize assumptions and diagnostics that limit interpretation of steady-state predictions.

Why Infusion Rate and Exposure Are Not Proportional Under TMDD

With linear elimination, doubling a constant infusion rate eventually doubles the steady-state concentration. Target-mediated drug disposition (TMDD) can break that relationship: when drug binds a finite target pool, the amount of drug removed through target binding and internalization depends on free drug concentration and target occupancy. As the target becomes saturated, this removal pathway contributes less additional clearance per unit increase in free drug.

This tutorial derives and uses the steady-state infusion relationship for a reduced TMDD model. The central distinction is between free drug, which drives binding and nonspecific elimination in this model, and total drug, which includes drug bound in the drug–target complex. The equations describe a maintained target pool and a steady state; they are not a dosing-interval model.

Definition: Let \(C\) be free drug concentration, \(R\) free target concentration, \(RC\) drug–target complex concentration, and \(R_{\mathrm{tot}}=R+RC\) the maintained total target concentration. Let \(u\) be the constant infusion input per unit volume, \(a\) the first-order nonspecific elimination coefficient for free drug, and \(b\) the first-order rate coefficient for irreversible drug loss through complex internalization. The total drug concentration is \(C_{\mathrm{tot}}=C+RC\). Under the quasi-steady binding relation, \(RC=R_{\mathrm{tot}}C/(K+C)\), where \(K\) is the effective binding constant for that reduction.

Concentrations are expressed in nM and \(u\) in nM/h in the examples. Both \(a\) and \(b\) have units of h\(^{-1}\), while \(K\) has units of nM. In a quasi-steady-state reduction, \(K\) commonly represents \((k_{\mathrm{off}}+k_{\mathrm{int}})/k_{\mathrm{on}}\), rather than the equilibrium dissociation constant \(k_{\mathrm{off}}/k_{\mathrm{on}}\). The appropriate value depends on the binding and internalization assumptions used for the reduction.

Steady-State Relationship and Its Assumptions

The infusion rate must balance irreversible drug loss at steady state. Binding itself transfers drug between the free and complex states; it is not a net loss from their sum. In this reduced model, free drug is lost through nonspecific elimination at rate \(aC\), and complexed drug is lost through internalization at rate \(bRC\). Thus the steady-state mass balance is:

$$ u=aC+bRC $$

Substituting the quasi-steady binding relation gives a single equation for free concentration:

$$ u=aC+bR_{\mathrm{tot}}\frac{C}{K+C} $$

The first term increases linearly with \(C\). The second term increases with \(C\), but approaches the finite maximum \(bR_{\mathrm{tot}}\) as the target becomes saturated. Multiplying through by \(K+C\) gives a quadratic equation:

$$ aC^2+(aK+bR_{\mathrm{tot}}-u)C-uK=0 $$

For \(a>0\), \(K>0\), \(R_{\mathrm{tot}}\geq 0\), and \(u>0\), the constant term is negative and the leading coefficient is positive. The two roots therefore have opposite signs, so exactly one root is positive. That root is the physically meaningful free concentration:

$$ C= \frac{u-aK-bR_{\mathrm{tot}}+ \sqrt{(aK+bR_{\mathrm{tot}}-u)^2+4auK}}{2a} $$

Once \(C\) is found, calculate \(RC=R_{\mathrm{tot}}C/(K+C)\), then \(C_{\mathrm{tot}}=C+RC\). For \(u=0\), the steady-state solution is \(C=RC=C_{\mathrm{tot}}=0\), under these assumptions.

Assumptions: The infusion is constant and a steady state has been reached; the effective target pool \(R_{\mathrm{tot}}\) is maintained; binding is represented by the stated quasi-steady relation; nonspecific elimination acts on free drug at coefficient \(a\); and internalization is a net drug-loss pathway at coefficient \(b\). The model omits, among other possibilities, elimination of complex by other routes and changes in target abundance. If those processes matter, this reduced steady-state relationship may not describe the system adequately.

The binding relation follows directly from the reduced complex balance. If formation is \(k_{\mathrm{on}}CR\) and complex disappearance is \((k_{\mathrm{off}}+k_{\mathrm{int}})RC\), the quasi-steady condition sets these rates equal:

$$ k_{\mathrm{on}}CR=(k_{\mathrm{off}}+k_{\mathrm{int}})RC \quad\Longrightarrow\quad RC=\frac{CR}{K} $$

Since \(R_{\mathrm{tot}}=R+RC\), substitute \(R=R_{\mathrm{tot}}-RC\) into \(RC=CR/K\). Rearranging gives \(RC(K+C)=R_{\mathrm{tot}}C\), hence \(RC=R_{\mathrm{tot}}C/(K+C)\). This justification depends on the maintained-pool reduction; a full TMDD model with changing target abundance requires its own steady-state target balance.

Worked Example: Free and Total Drug at a Low Infusion Rate

Consider a hypothetical reduced TMDD model with \(a=0.10\ \mathrm{h}^{-1}\), \(b=0.20\ \mathrm{h}^{-1}\), \(R_{\mathrm{tot}}=5.0\ \mathrm{nM}\), \(K=0.20\ \mathrm{nM}\), and \(u=0.050\ \mathrm{nM/h}\). Calculate steady-state free and total drug concentrations.

Substitution into the quadratic gives:

$$ 0.10C^2+(0.10(0.20)+0.20(5.0)-0.050)C-(0.050)(0.20)=0 $$

Thus \(0.10C^2+0.970C-0.0100=0\). Using the positive root, \(C=[-0.970+\sqrt{(0.970)^2+4(0.10)(0.0100)}]/[2(0.10)] =0.0102984\ \mathrm{nM}\) (rounded).

Using that unrounded root in the binding relation gives \(RC=5.0(0.0102984)/(0.20+0.0102984)=0.244851\ \mathrm{nM}\) (rounded). Therefore \(C_{\mathrm{tot}}=0.0102984+0.244851=0.255149\ \mathrm{nM}\), with the sum calculated from unrounded values. The target occupancy is \(RC/R_{\mathrm{tot}}=0.244851/5.0=0.0489702\), or about \(4.90\%\).

Although free concentration is small, much of the drug present is in complex: the complex concentration is about 23.8 times the free concentration. That distinction matters when comparing predicted concentrations with an assay that measures total rather than free drug.

How Saturation Changes Dose Proportionality

The free concentration response to a change in infusion rate can be examined by differentiating the steady-state input function \(u(C)=aC+bR_{\mathrm{tot}}C/(K+C)\):

$$ \frac{du}{dC}=a+\frac{bR_{\mathrm{tot}}K}{(K+C)^2} $$

The local proportionality of free concentration to infusion rate is summarized by the elasticity \(d\ln C/d\ln u\). Differentiating gives:

$$ \frac{d\ln C}{d\ln u} =\frac{u}{C(du/dC)} =\frac{a+bR_{\mathrm{tot}}/(K+C)} {a+bR_{\mathrm{tot}}K/(K+C)^2} $$

For positive \(C\) and a nonzero target-mediated pathway, \(K/(K+C)^2<1/(K+C)\). Therefore the numerator exceeds the denominator and the elasticity is greater than 1: locally, a percentage increase in infusion rate produces a larger percentage increase in free concentration. At low concentrations, binding-mediated removal can be substantial relative to free drug. At high concentrations, the complex approaches its ceiling \(R_{\mathrm{tot}}\), so further increases in input are increasingly balanced by nonspecific elimination, while the free concentration must rise.

This result concerns free concentration. Total concentration also includes a bounded complex pool. Consequently, total concentration need not change by the same factor as free concentration, especially while the bound pool is filling. Dose proportionality should be assessed using the concentration measure relevant to the question or assay.

Worked Example: Comparing Three Infusion Rates

Use the same hypothetical parameters as in the previous example. Calculate free and total concentrations at \(u=0.50\) and \(2.00\ \mathrm{nM/h}\), then compare the results with \(u=0.050\ \mathrm{nM/h}\).

For \(u=0.50\ \mathrm{nM/h}\), the quadratic is \(0.10C^2+0.520C-0.100=0\). Its positive root is \(C=0.185678\ \mathrm{nM}\). Using the unrounded root, \(RC=5.0C/(0.20+C)=2.40716\ \mathrm{nM}\), and \(C_{\mathrm{tot}}=C+RC=2.59284\ \mathrm{nM}\). The mass balance checks: \(0.10(0.185678)+0.20(2.40716)=0.500000\ \mathrm{nM/h}\), to the shown precision.

For \(u=2.00\ \mathrm{nM/h}\), the quadratic is \(0.10C^2-0.980C-0.400=0\). Its positive root is \(C=10.1924\ \mathrm{nM}\). Using the unrounded root, \(RC=5.0C/(0.20+C)=4.90378\ \mathrm{nM}\), and \(C_{\mathrm{tot}}=15.0962\ \mathrm{nM}\). The input check is \(0.10(10.1924)+0.20(4.90378)=2.00000\ \mathrm{nM/h}\), to the shown precision.

Increasing the infusion tenfold, from \(0.050\) to \(0.50\ \mathrm{nM/h}\), increases free concentration about 18.0-fold, but total concentration about 10.2-fold. Increasing the infusion fourfold, from \(0.50\) to \(2.00\ \mathrm{nM/h}\), increases free concentration about 54.9-fold and total concentration about 5.82-fold. The large increase in free concentration at the upper rate occurs as the complex approaches its maximum of \(5.0\ \mathrm{nM}\); the calculated occupancy is about \(4.90378/5.0=0.980756\), or \(98.1\%\).

These comparisons illustrate the model’s nonlinear prediction, not a universal numerical pattern. The magnitude depends on the parameters, the concentration range, and whether exposure is measured as free or total drug.

Clearance Depends on Which Concentration Is Used

At steady state, two useful ratios can be reported, but they answer different questions. The input-to-free-concentration ratio is \(u/C\); from the mass balance, it equals \(a+bRC/C\). It includes the concentration-dependent contribution of target-mediated loss relative to free drug. The input-to-total-concentration ratio is \(u/C_{\mathrm{tot}}\), a total-concentration-based apparent clearance. It is not generally equal to either \(a\) or \(b\), because drug is distributed between free and bound forms and the elimination pathways act on different states.

For the three hypothetical infusion rates, the ratios are approximately:

Infusion \(u\) (nM/h)Free \(C\) (nM)Total \(C_{\mathrm{tot}}\) (nM)\(u/C\) (h\(^{-1}\))\(u/C_{\mathrm{tot}}\) (h\(^{-1}\))
0.0500.01029840.2551494.8550.1960
0.500.1856782.592842.6930.1928
2.0010.192415.09620.19620.1325

The ratios change with infusion rate, so a single constant “clearance” cannot summarize every concentration measure and range in this example. Always state which concentration is in the denominator and whether the reported quantity is a model coefficient, an effective ratio, or an apparent clearance based on total measured drug.

Worked Example: Infusion Rate for a Desired Occupancy

Suppose the desired target occupancy is \(80\%\) in the same hypothetical model. Under the binding relation, occupancy is \(RC/R_{\mathrm{tot}}=C/(K+C)\). Solve for the free concentration:

$$ 0.80=\frac{C}{0.20+C} \quad\Longrightarrow\quad 0.80(0.20+C)=C \quad\Longrightarrow\quad C=0.80\ \mathrm{nM} $$

The corresponding complex concentration is \(RC=0.80(5.0)=4.0\ \mathrm{nM}\). The required steady-state infusion input follows from the mass balance: \(u=aC+bRC=(0.10)(0.80)+(0.20)(4.0)=0.080+0.800=0.880\ \mathrm{nM/h}\).

Thus the model predicts that \(0.880\ \mathrm{nM/h}\) is required to maintain \(80\%\) occupancy at steady state under the stated assumptions. This calculation does not specify how quickly that occupancy is reached.

Model Choice, Diagnostics, and Common Pitfalls

The reduced steady-state equation is useful when the question is the relationship between a constant input and steady-state exposure, and the binding and target-pool assumptions are defensible. It is more mechanistic than a simple Michaelis–Menten elimination approximation because it explicitly represents the bound pool and its occupancy. A Michaelis–Menten form can be a convenient approximation over a specified concentration range, but its parameters may not map uniquely to binding, target turnover, and internalization processes. A full TMDD model is preferable when transient binding, changing target abundance, or other state-specific processes are important. It requires solving the coupled state equations rather than relying on this single steady-state equation.

Several checks help identify when the reduced relationship is being misused. First, confirm that the input is truly constant and that observations represent steady state; a concentration still changing over time cannot be interpreted using this balance. Second, check that the predicted complex does not exceed \(R_{\mathrm{tot}}\), and that calculated \(C\), \(RC\), and \(C_{\mathrm{tot}}\) are nonnegative. Third, verify the mass balance by substituting the calculated concentrations back into \(u=aC+bRC\). Fourth, compare predictions with the concentration definition used for data: an assay measuring total drug should not be compared directly with predicted free drug without an appropriate observation model.

Systematic disagreement across infusion rates can indicate that the maintained-pool or quasi-steady binding assumptions are inadequate, that a modeled elimination pathway is missing, or that the system has not reached steady state. A good fit at one infusion rate alone does not establish that the predicted saturation pattern is reliable at other rates.

  • Mistake: treating binding as drug elimination. Binding moves drug from \(C\) to \(RC\); in the combined free-plus-complex balance, the binding terms cancel. Correction: count irreversible loss pathways, such as the specified complex internalization, in the input balance.
  • Mistake: assuming total drug is the free concentration. When a substantial fraction is bound, \(C_{\mathrm{tot}}=C+RC\) can be much larger than \(C\). Correction: calculate and label both quantities.
  • Mistake: interpreting an apparent clearance ratio as a fixed model parameter. Ratios such as \(u/C\) vary with occupancy in this model. Correction: state the concentration denominator and report the infusion-rate range and units.
  • Mistake: assuming a constant target pool is automatically valid. Target synthesis, degradation, recycling, or abundance changes can alter the relationship. Correction: use the reduced equation only when its target-pool assumption is reasonable.
Key takeaway: In the reduced steady-state TMDD model, \(u=aC+bR_{\mathrm{tot}}C/(K+C)\). The target-mediated term saturates, so increasing infusion can produce a more-than-proportional increase in free concentration. Report free and total exposure separately, and interpret apparent clearance only with its concentration basis specified.

Check Your Understanding

Use the stated reduced model assumptions. Answers follow each question so that the calculations can be checked.

  1. For \(a=0.10\ \mathrm{h}^{-1}\), \(b=0.20\ \mathrm{h}^{-1}\), \(R_{\mathrm{tot}}=5.0\ \mathrm{nM}\), \(K=0.20\ \mathrm{nM}\), and \(u=0.50\ \mathrm{nM/h}\), write the quadratic for \(C\).
    Solution: \(0.10C^2+(0.020+1.00-0.50)C-(0.50)(0.20)=0\), or \(0.10C^2+0.520C-0.100=0\).
  2. If \(C=0.80\ \mathrm{nM}\) and \(K=0.20\ \mathrm{nM}\), what is occupancy under the binding relation?
    Solution: \(C/(K+C)=0.80/(0.20+0.80)=0.80\), or \(80\%\).
  3. Why does the target-mediated term not increase without bound as \(C\) increases?
    Solution: \(RC=R_{\mathrm{tot}}C/(K+C)\) approaches \(R_{\mathrm{tot}}\) as \(C\) becomes large. Therefore \(bRC\) approaches the finite maximum \(bR_{\mathrm{tot}}\).
  4. At steady state, what does \(u/C_{\mathrm{tot}}\) represent, and why is it not generally equal to \(a\)?
    Solution: It is the infusion input divided by total drug concentration, an apparent clearance based on total concentration. The coefficient \(a\) applies only to nonspecific elimination of free drug; total concentration also includes complex, which has a separate loss pathway.
  5. Name one diagnostic check for a calculated steady-state solution and give the equation used.
    Solution: Substitute the calculated concentrations into the mass balance and check that \(u=aC+bRC\), allowing for rounding. Also verify \(0\leq RC\leq R_{\mathrm{tot}}\) and \(C_{\mathrm{tot}}=C+RC\).
Clinical Biostats