Why Irreversible Binding Depends on Time
For a rapidly reversible ligand, target engagement can often be described by association and dissociation, with occupancy changing as free-ligand concentration changes. Covalent or effectively irreversible ligands behave differently: once a target molecule is modified, it does not simply dissociate back into the active pool on the time scale of interest. Consequently, the fraction of target modified can continue to increase during sustained exposure, even when concentration is constant. It can also remain elevated after the drug has left the system.
Here the outcome is the fraction of the target that is covalently modified and functionally inactivated. This fraction is sometimes called irreversible occupancy, but it is not necessarily the fraction of all receptors occupied by any ligand. The model separates irreversible modification from pharmacological efficacy: a modified receptor is assumed to be inactivated, while the magnitude of any downstream response is outside this model.
The recovery term is a phenomenological representation of restoration of the active target pool. It assumes that recovered target replaces modified target so that the total target amount stays constant. It does not specify the biological process responsible for recovery. A later treatment of receptor turnover can represent that process more explicitly.
Model and Assumptions
The model assumes a well-mixed target pool with only two states: active and modified/inactivated. Drug modifies active target with a concentration-dependent rate, and modified target returns to the active pool with first-order rate constant \(k_{\mathrm{rec}}\). Modification is treated as irreversible on the modeled time scale: there is no dissociation term for the modified drug-target complex. The model also assumes that \(C(t)\) is the free concentration at the target site, not necessarily the measured plasma concentration.
A common effective inactivation-rate relationship is \(k_{\mathrm{obs}}(C)=k_{\mathrm{inact}}C/(K_I+C)\), where \(k_{\mathrm{inact}}\) is the maximum inactivation rate and \(K_I\) is a concentration scale. This relationship can arise when drug first forms a reversible encounter complex and then reacts, under conditions where formation of that complex is sufficiently rapid for an effective rate to be used. It is a modeling choice, not a universal law. If the mechanism or data do not support it, \(k_{\mathrm{obs}}\) should be specified or estimated using a more suitable model.
The rate of change in \(q(t)\) is the rate at which active target is modified, minus the rate at which modified target recovers:
For constant concentration \(C\), write \(k_{\mathrm{obs}}=k_{\mathrm{obs}}(C)\). Then \(\frac{dq}{dt}=k_{\mathrm{obs}}-(k_{\mathrm{obs}}+k_{\mathrm{rec}})q\). The constant-exposure solution from initial fraction \(q_0\) is
Here \(q_{\mathrm{ss}}\) is the long-term inactivated fraction under that constant concentration, not necessarily complete inactivation. It is below one when recovery is present. When exposure stops, \(C=0\), so \(k_{\mathrm{obs}}(0)=0\) and \(\frac{dq}{dt}=-k_{\mathrm{rec}}q\). Thus the washout-period solution, starting from \(q_w\), is \(q(t)=q_w e^{-k_{\mathrm{rec}}t}\), where \(t\) now measures time since washout.
Why the Constant-Exposure Solution Has This Form
For constant concentration, define \(\lambda=k_{\mathrm{obs}}+k_{\mathrm{rec}}\). The differential equation is \(dq/dt+\lambda q=k_{\mathrm{obs}}\). Multiplying by the integrating factor \(e^{\lambda t}\) gives \(\frac{d}{dt}(q e^{\lambda t})=k_{\mathrm{obs}}e^{\lambda t}\). Integrating from \(0\) to \(t\) yields \(q(t)e^{\lambda t}-q_0=(k_{\mathrm{obs}}/\lambda)(e^{\lambda t}-1)\). Rearranging gives the stated solution, because \(k_{\mathrm{obs}}/\lambda=q_{\mathrm{ss}}\). This also shows that the difference \(q(t)-q_{\mathrm{ss}}\) decays exponentially at rate \(\lambda\). If both rates are zero, there is no change and \(q(t)=q_0\); the formula involving \(q_{\mathrm{ss}}\) is not used because its denominator is zero.
Worked Example: Inactivation During an Eight-Hour Exposure
Suppose a hypothetical ligand has \(k_{\mathrm{inact}}=0.12\ \mathrm{h}^{-1}\) and \(K_I=2.0\ \mathrm{\mu M}\). Free concentration is held at \(C=6.0\ \mathrm{\mu M}\) for 8 hours, recovery has \(k_{\mathrm{rec}}=0.020\ \mathrm{h}^{-1}\), and initially \(q_0=0\). Calculate the inactivated and active fractions at 8 hours.
The concentration-dependent inactivation rate is \(k_{\mathrm{obs}}=0.12(6.0)/(2.0+6.0)=0.0900\ \mathrm{h}^{-1}\). The assumptions of constant free concentration and constant rates apply by construction in this hypothetical example. The long-term inactivated fraction at this concentration is \(q_{\mathrm{ss}}=0.0900/(0.0900+0.020)=0.8182\).
Using the constant-exposure solution, \(q(8)=0.8182+(0-0.8182)e^{-(0.110)(8)} =0.8182(1-e^{-0.880})\). Since \(e^{-0.880}\approx0.4148\), \(q(8)\approx0.8182(0.5852)=0.4788\). The active fraction is \(1-q(8)=0.5212\). Thus about 47.88% of the target is inactivated after 8 hours, and about 52.12% remains active under these model assumptions. The inactivated fraction is below its long-term value because the exposure has not continued indefinitely.
Exposure Duration, Concentration, and Recovery
If recovery is negligible during exposure, set \(k_{\mathrm{rec}}=0\). When \(k_{\mathrm{obs}}>0\), at constant concentration the inactivated fraction approaches one as \(t\to\infty\), according to \(q(t)=1-(1-q_0)e^{-k_{\mathrm{obs}}t}\). If \(k_{\mathrm{obs}}=0\), then \(q(t)=q_0\). Longer exposure increases inactivation, but the increase becomes progressively smaller as fewer active target molecules remain. If recovery is appreciable, the competition between modification and recovery creates a plateau \(q_{\mathrm{ss}}\). Increasing exposure duration then moves \(q(t)\) toward that plateau rather than toward complete inactivation.
For a desired inactivated fraction \(q^*\), the constant-exposure model can also give the required duration. When \(q_0<q^*<q_{\mathrm{ss}}\),
The target level must be below \(q_{\mathrm{ss}}\) to be reached in finite time from below. The formula is not a way to obtain a finite duration for a target level above the plateau.
Worked Example: Calculating Inactivation and Recovery Across a Pulse
In a hypothetical experiment, the effective inactivation rate during a 4-hour exposure is \(0.150\ \mathrm{h}^{-1}\), \(k_{\mathrm{rec}}=0.030\ \mathrm{h}^{-1}\), and \(q_0=0.200\). Drug concentration then falls to zero for 12 hours. Calculate the inactivated fraction at the end of exposure and at the end of washout.
During exposure, \(q_{\mathrm{ss}}=0.150/(0.150+0.030)=0.8333\), and the combined rate is \(0.180\ \mathrm{h}^{-1}\). Therefore, \(q(4)=0.8333+(0.200-0.8333)e^{-(0.180)(4)} =0.8333-0.6333e^{-0.720}\). Using unrounded values, \(q(4)=0.833333\ldots-0.633333\ldots e^{-0.720}\approx0.5251\).
During washout, there is no further inactivation in this model. The modified fraction decays at the recovery rate: \(q_{\mathrm{washout}}=0.5251e^{-(0.030)(12)} =0.5251e^{-0.360}\). Using unrounded values, \(q_{\mathrm{washout}}\approx(0.833333\ldots-0.633333\ldots e^{-0.720})e^{-0.360}\approx0.3663\). The active fraction after washout is therefore approximately \(1-0.3663=0.6337\). The nonzero residual inactivation reflects finite recovery over the 12-hour interval.
Worked Example: Time to a Desired Inactivation Level
Suppose a hypothetical exposure produces \(k_{\mathrm{obs}}=0.080\ \mathrm{h}^{-1}\), with \(k_{\mathrm{rec}}=0.020\ \mathrm{h}^{-1}\) and initial inactivation \(q_0=0.10\). How long does constant exposure take to reach \(q^*=0.60\)?
The plateau is \(q_{\mathrm{ss}}=0.080/(0.080+0.020)=0.80\), so the desired fraction \(0.60\) is attainable from \(q_0=0.10\). Substitution into the duration equation gives \(t=-\ln[(0.80-0.60)/(0.80-0.10)]/0.10 =-\ln(0.20/0.70)/0.10\). Since \(-\ln(0.2857)\approx1.253\), the required time is approximately \(12.53\ \mathrm{h}\). The answer is a duration under constant exposure at the specified effective rate; a changing concentration would require accounting for the changing rate over time.
Diagnostics, Alternatives, and Limitations
Before interpreting a prediction, check that concentration is free concentration at the target site, that the exposure schedule is represented correctly, and that the rate constants use the same time units as the time variable. The model predicts \(0\leq q(t)\leq1\) when \(0\leq q_0\leq1\) and both rates are nonnegative. A fitted trajectory outside that range, or systematic deviations from the predicted approach to a plateau and subsequent recovery, suggests a coding error or a model that does not capture the data.
A key model check is whether the long-exposure plateau makes biological sense. If the data show continuing inactivation despite a substantial proposed recovery rate, the assumed balance may be wrong. Conversely, if apparent engagement declines during washout much faster or slower than \(e^{-k_{\mathrm{rec}}t}\), a single constant recovery rate may not describe the process. Changes in free concentration, drug loss, multiple target states, incomplete functional inactivation, and variation between target pools are possible explanations; the simple model does not distinguish among them.
Compared with reversible one-site binding, this model has no dissociation of the modified complex. A reversible model instead includes association and dissociation, and its occupancy can decline as ligand concentration falls even without target replacement. For irreversible binding, washout alone stops new modification but does not instantly remove existing modified target; decline requires recovery in the model. The irreversible model is appropriate only when dissociation is negligible over the time scale being studied. If modification is reversible on that time scale, a reversible or explicitly reversible-adduct model is preferable.
The concentration-response expression for \(k_{\mathrm{obs}}\) also has limits. It assumes a particular effective relationship between concentration and chemical inactivation. A changing or depleted free-drug concentration, slow encounter-complex formation, multiple binding modes, or time-varying reactivity can invalidate that relationship. In such settings, use a time-varying \(C(t)\) in the differential equation if the rate law remains appropriate, or choose a mechanistic model with additional states. The model predicts target state, not downstream response; equal inactivation fractions need not imply equal functional effects.
Check Your Understanding
Answer each question, then compare with the solutions below.
- If \(C=K_I\), what is \(k_{\mathrm{obs}}\) under the stated concentration-response relationship?
- If \(q=0.60\) at washout and \(k_{\mathrm{rec}}=0.040\ \mathrm{h}^{-1}\), what is \(q\) 10 hours later?
- At constant exposure, what is \(q_{\mathrm{ss}}\) if \(k_{\mathrm{obs}}=0.12\ \mathrm{h}^{-1}\) and \(k_{\mathrm{rec}}=0.030\ \mathrm{h}^{-1}\)?
- Can a target inactivation fraction of \(0.90\) be reached from below if \(q_{\mathrm{ss}}=0.80\)? Explain.
- Why does the model not predict that inactivated target immediately becomes active when free drug concentration reaches zero?
Solutions. (1) \(k_{\mathrm{obs}}=k_{\mathrm{inact}}K_I/(K_I+K_I)=k_{\mathrm{inact}}/2\). (2) \(q=0.60e^{-(0.040)(10)}=0.60e^{-0.400}\approx0.4022\). (3) \(q_{\mathrm{ss}}=0.12/(0.12+0.030)=0.8000\). (4) No. Under constant exposure the inactivated fraction approaches \(0.80\) from below and cannot reach \(0.90\). (5) Zero concentration stops new modification, but the model removes already modified target only through recovery, at rate \(k_{\mathrm{rec}}\); therefore the inactivated fraction declines over time rather than disappearing instantly.