Rapid-Binding Equilibrium Approximation

Why a Rapid-Binding Approximation Is Useful

In the full target-mediated drug disposition (TMDD) model, free drug, free target, and drug–target complex change together. When binding and dissociation occur much faster than drug disposition and target turnover, the short binding transient may be too brief to matter for the slower concentration profile. The rapid-binding equilibrium approximation uses that separation of timescales to replace the binding dynamics with an algebraic relation.

The key parameter is the equilibrium dissociation constant \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}\). It describes how free drug and free target are apportioned between unbound and complexed forms under the approximation. This is a statement about the binding subsystem; it does not mean that the entire TMDD system is at equilibrium.

Definition: Let \(C\) denote free drug concentration, \(R\) free target concentration, and \(RC\) drug–target complex concentration. The rapid-binding equilibrium approximation assumes association and dissociation are fast relative to target turnover and drug disposition, so that \(k_{\mathrm{on}}CR\approx k_{\mathrm{off}}RC\). With \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}\), this gives \(RC\approx CR/K_d\).

Deriving the Binding Relation

Consider the mass-action TMDD equations between doses, with input rate \(u(t)\) included for generality:

\begin{aligned} \frac{dC}{dt} &= u(t)-k_{\mathrm{el}}C-k_{\mathrm{on}}CR+k_{\mathrm{off}}RC,\\ \frac{dR}{dt} &= k_{\mathrm{syn}}-k_{\mathrm{deg}}R-k_{\mathrm{on}}CR+k_{\mathrm{off}}RC,\\ \frac{dRC}{dt} &= k_{\mathrm{on}}CR-(k_{\mathrm{off}}+k_{\mathrm{int}})RC. \end{aligned}

Here \(k_{\mathrm{el}}\) is the first-order elimination rate constant for free drug; \(k_{\mathrm{syn}}\) and \(k_{\mathrm{deg}}\) describe target synthesis and degradation; and \(k_{\mathrm{int}}\) is the internalization rate constant for complex. The association rate constant \(k_{\mathrm{on}}\) has units of concentration\(^{-1}\) time\(^{-1}\), while \(k_{\mathrm{off}}\), \(k_{\mathrm{int}}\), \(k_{\mathrm{el}}\), and \(k_{\mathrm{deg}}\) have units of time\(^{-1}\). These equations describe a dynamic system: in general, \(dRC/dt\) is not zero.

For the rapid-binding approximation, focus on the association and dissociation contributions to change in \(RC\). Their characteristic rates are of order \(k_{\mathrm{on}}(C+R)\) and \(k_{\mathrm{off}}\), whereas target turnover, complex internalization, and drug disposition act on the slower timescale. In particular, the approximation requires complex internalization not to appreciably disrupt the fast binding balance. To leading order on the slow timescale, association and dissociation balance:

k_{\mathrm{on}}CR\approx k_{\mathrm{off}}RC.

Dividing by \(k_{\mathrm{on}}\) and defining \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}\) gives \(CR\approx K_dRC\), or \(RC\approx CR/K_d\). The approximation is therefore about the fast binding subsystem, not a claim that synthesis, elimination, or internalization stop.

It is often convenient to define total target \(R_{\mathrm{tot}}=R+RC\), counting target both free and bound. Substituting \(R=R_{\mathrm{tot}}-RC\) into \(RC\approx CR/K_d\) gives:

RC\approx \frac{R_{\mathrm{tot}}C}{K_d+C}, \qquad R\approx \frac{R_{\mathrm{tot}}K_d}{K_d+C}.

The bound fraction of target is consequently \(RC/R_{\mathrm{tot}}\approx C/(K_d+C)\). When \(C=K_d\), the approximation predicts half of the total target is bound. The equations remain time-varying because \(C\) and \(R_{\mathrm{tot}}\) may change over time.

Worked Calculations

Worked Example: Calculating \(K_d\) and Target Occupancy

Suppose a hypothetical system has \(k_{\mathrm{off}}=0.60\ \mathrm{h}^{-1}\) and \(k_{\mathrm{on}}=0.20\ \mathrm{L}\,\mathrm{nmol}^{-1}\,\mathrm{h}^{-1}\). At one time point, free drug is \(C=10.0\ \mathrm{nmol/L}\), and total target is \(R_{\mathrm{tot}}=5.00\ \mathrm{nmol/L}\). Under rapid binding, calculate \(K_d\), \(RC\), and \(R\).

First, the dissociation constant is \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}=(0.60\ \mathrm{h}^{-1})/(0.20\ \mathrm{L}\,\mathrm{nmol}^{-1}\,\mathrm{h}^{-1})=3.00\ \mathrm{nmol/L}\). The units are concentration, as required.

Using the rapid-binding relation, \(RC\approx R_{\mathrm{tot}}C/(K_d+C)=(5.00)(10.0)/(3.00+10.0)=50.0/13.0=3.846\ \mathrm{nmol/L}\). The free target is \(R\approx R_{\mathrm{tot}}-RC=5.00-3.846=1.154\ \mathrm{nmol/L}\). Thus the predicted bound fraction is \(3.846/5.00=0.7692\), or about 76.9%. These are approximate values conditional on the rapid-binding assumption.

Worked Example: Comparing \(K_d\) With the Quasi-Steady-State Constant

Suppose a hypothetical model has \(k_{\mathrm{off}}=0.60\ \mathrm{h}^{-1}\), \(k_{\mathrm{on}}=0.20\ \mathrm{L}\,\mathrm{nmol}^{-1}\,\mathrm{h}^{-1}\), and \(k_{\mathrm{int}}=0.20\ \mathrm{h}^{-1}\). At a given time, let \(C=3.00\ \mathrm{nmol/L}\) and \(R_{\mathrm{tot}}=12.0\ \mathrm{nmol/L}\). Compare the rapid-binding prediction with the quasi-steady-state (QSS) binding relation.

Rapid binding uses \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}=0.60/0.20=3.00\ \mathrm{nmol/L}\). Therefore, \(RC\approx 12.0(3.00)/(3.00+3.00)=6.00\ \mathrm{nmol/L}\).

In the QSS approximation, setting \(dRC/dt\approx0\) in the full complex equation gives \(k_{\mathrm{on}}CR\approx(k_{\mathrm{off}}+k_{\mathrm{int}})RC\). Its corresponding constant is \(K_{\mathrm{SS}}=(k_{\mathrm{off}}+k_{\mathrm{int}})/k_{\mathrm{on}}=(0.60+0.20)/0.20=4.00\ \mathrm{nmol/L}\). Thus the QSS relation gives \(RC\approx12.0(3.00)/(4.00+3.00)=36.0/7.00=5.143\ \mathrm{nmol/L}\).

The results differ because the QSS balance includes loss of complex through internalization, whereas the rapid-binding balance neglects that loss relative to association and dissociation. Here internalization is one-third of \(k_{\mathrm{off}}\), so the two constants are not close; the numerical difference warns against treating the approximations as interchangeable. It does not, by itself, establish which approximation best describes a particular time course.

Worked Example: Selecting a Physically Constrained Root

Suppose a hypothetical assay inversion yields the equation \(C^2-16C+40=0\), with \(C\) in \(\mathrm{nmol/L}\). An independent measurement establishes that the relevant concentration is greater than \(8.00\ \mathrm{nmol/L}\). Find the admissible value and use it to calculate rapid-binding occupancy when \(K_d=4.00\ \mathrm{nmol/L}\) and \(R_{\mathrm{tot}}=16.0\ \mathrm{nmol/L}\).

By the quadratic formula, \(C=(16\pm\sqrt{16^2-4(1)(40)})/2=(16\pm\sqrt{96})/2\). Since \(\sqrt{96}=9.798\), the two roots are \((16+9.798)/2=12.899\) and \((16-9.798)/2=3.101\ \mathrm{nmol/L}\). The independently specified condition \(C>8.00\ \mathrm{nmol/L}\) selects \(C=12.9\ \mathrm{nmol/L}\) to three significant figures.

The rapid-binding complex concentration is then \(RC\approx16.0(12.9)/(4.00+12.9)=206.4/16.9=12.21\ \mathrm{nmol/L}\). Free target is \(R\approx16.0-12.21=3.79\ \mathrm{nmol/L}\), and the bound fraction is \(12.9/16.9\approx0.7633\), or 76.3%. The quadratic alone has two roots; the independent concentration constraint, not the binding approximation, determines which root is admissible in this example.

Assumptions, Diagnostics, and Alternatives

The approximation is most defensible when binding relaxes much faster than the quantities that move the system away from binding equilibrium. A useful qualitative comparison is between binding rates such as \(k_{\mathrm{off}}\) and \(k_{\mathrm{on}}(C+R)\), and the slower rates \(k_{\mathrm{el}}\), \(k_{\mathrm{deg}}\), and \(k_{\mathrm{int}}\), together with the rates of change induced by input and target synthesis. There is no universal numerical cutoff: the relevant comparison depends on concentrations and on the time interval being modeled.

The derivation also assumes that the mass-action representation is appropriate and that concentrations and rate constants are on a consistent scale. If \(k_{\mathrm{int}}\) is not small relative to the dissociation and association processes, complex is removed before binding can remain close to equilibrium. Likewise, a rapid change in dose or a rapid change in free drug can create a transient during which the algebraic relation is not accurate.

Practical checks include comparing predicted time courses from the reduced relation with solutions of the full mass-action equations over the intended dose range and observation times. Large, systematic discrepancies around dosing changes, during rapid concentration changes, or as target concentrations change indicate that the assumed timescale separation may be inadequate. Parameter estimates that require implausibly fast binding to fit data are another reason to inspect the assumption, although estimates alone do not prove that it fails.

The main alternative at this stage is to retain the full differential equations, which preserves finite binding and dissociation kinetics but requires solving the additional dynamic state. The QSS approximation instead balances complex formation with both dissociation and internalization, producing \(K_{\mathrm{SS}}=(k_{\mathrm{off}}+k_{\mathrm{int}})/k_{\mathrm{on}}\). Rapid binding uses \(K_d=k_{\mathrm{off}}/k_{\mathrm{on}}\) and is appropriate only when internalization does not materially disturb the fast binding balance. Choosing between them is a mechanistic assumption, not a choice between two equivalent formulas.

Common pitfall: Setting \(dRC/dt=0\) and automatically calling the result rapid binding. In the full model that operation gives \(k_{\mathrm{on}}CR=(k_{\mathrm{off}}+k_{\mathrm{int}})RC\), the QSS balance. Rapid binding instead requires association and dissociation to dominate the slower processes, including complex internalization, so its relation is \(k_{\mathrm{on}}CR\approx k_{\mathrm{off}}RC\). Check which terms are being balanced before selecting the constant.

The approximation is valuable because it removes a fast binding transient from a model focused on slower disposition and turnover. Its limitation is the same simplification: it cannot represent a meaningful delay in association or dissociation. It should therefore be understood as a reduced description for a specified regime, not a universal replacement for the full TMDD equations.

Check Your Understanding

Answer each question using the assumptions and definitions above. The solutions are included so that the calculations can be checked directly.

If \(k_{\mathrm{off}}=0.45\ \mathrm{h}^{-1}\) and \(k_{\mathrm{on}}=0.15\ \mathrm{L}\,\mathrm{nmol}^{-1}\,\mathrm{h}^{-1}\), what is \(K_d\)?
Solution: \(K_d=0.45/0.15=3.00\ \mathrm{nmol/L}\).
Under rapid binding, what fraction of target is bound when \(C=2K_d\)?
Solution: \(RC/R_{\mathrm{tot}}\approx C/(K_d+C)=2K_d/(K_d+2K_d)=2/3\), or 66.7%.
For \(C=6.00\ \mathrm{nmol/L}\), \(K_d=2.00\ \mathrm{nmol/L}\), and \(R_{\mathrm{tot}}=9.00\ \mathrm{nmol/L}\), calculate \(RC\) and \(R\).
Solution: \(RC\approx9.00(6.00)/(2.00+6.00)=6.75\ \mathrm{nmol/L}\). Thus \(R\approx9.00-6.75=2.25\ \mathrm{nmol/L}\).
With \(k_{\mathrm{off}}=0.50\ \mathrm{h}^{-1}\), \(k_{\mathrm{int}}=0.30\ \mathrm{h}^{-1}\), and \(k_{\mathrm{on}}=0.20\ \mathrm{L}\,\mathrm{nmol}^{-1}\,\mathrm{h}^{-1}\), calculate \(K_d\) and \(K_{\mathrm{SS}}\). Why do they differ?
Solution: \(K_d=0.50/0.20=2.50\ \mathrm{nmol/L}\), while \(K_{\mathrm{SS}}=(0.50+0.30)/0.20=4.00\ \mathrm{nmol/L}\). QSS includes complex internalization in its balance; rapid binding neglects that loss relative to the fast binding processes.
Why may the rapid-binding relation be unreliable immediately after a sharp dose change?
Solution: Free drug may change rapidly, so association and dissociation may need time to restore their balance. If that transient is not short relative to the relevant observation timescale, the algebraic equilibrium relation may not describe the changing complex concentration.

Rapid-Binding Equilibrium Approximation

What You'll Learn