Why Binding Kinetics Matter for MABEL
An affinity value summarizes the balance between association and dissociation at equilibrium. It does not, by itself, say how quickly a drug occupies its target after exposure begins, or how quickly occupancy falls after concentration declines. Those distinctions matter when target engagement is assessed over a finite period, when drug concentrations change appreciably, or when occupancy persists after free drug has fallen.
For an occupancy-based MABEL estimate, an equilibrium calculation can be adequate if binding is close to equilibrium at the relevant concentrations and times. If not, the concentration associated with a chosen occupancy threshold depends on the time course of exposure as well as on affinity. The kinetic calculation is therefore not automatically “more conservative”: slow association can reduce early occupancy, while slow dissociation can prolong occupancy during washout.
From Binding Rates to Occupancy Over Time
Consider a target with total concentration \(R_{\mathrm{tot}}\), of which \(B\) is drug-bound. Assume reversible one-to-one binding, a homogeneous accessible target pool, and a free-drug concentration \(C\) that is constant over the interval being analyzed. Assume also that binding is the only process changing the bound-target amount: for example, target synthesis, internalization, and degradation are not represented. The mass-action model is
With \(\theta=B/R_{\mathrm{tot}}\), the corresponding fractional-occupancy equation is
At equilibrium, the derivative is zero. Solving gives the familiar affinity-only expression. For a receptor initially unoccupied, the time-dependent solution follows from solving the linear differential equation.
The equilibrium occupancy depends on \(K_D\), whereas the rate of approach to equilibrium depends on \(k_{\mathrm{on}}C+k_{\mathrm{off}}\). The characteristic time scale is \(1/(k_{\mathrm{on}}C+k_{\mathrm{off}})\); after several such time scales at constant concentration, occupancy is close to its equilibrium value. If concentration varies with time, the constant-concentration solution cannot simply be applied at each instant as though equilibrium were immediate. Instead, occupancy follows the changing concentration with a lag.
This expression makes explicit that occupancy at time \(t\) reflects both current concentration and earlier exposure. In practice, it may be evaluated numerically using a predicted concentration-time profile.
Worked Example: Early Occupancy Below the Affinity Prediction
Suppose \(k_{\mathrm{on}}=1.00\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\), \(k_{\mathrm{off}}=0.200\ \mathrm{h}^{-1}\), and free-drug concentration is held at \(0.0630\ \mathrm{nM}\). Starting from \(\theta(0)=0\), calculate occupancy after 2.00 hours and compare it with the equilibrium prediction.
The equilibrium dissociation constant is \(K_D=0.200/1.00=0.200\ \mathrm{nM}\). At the stated concentration, the equilibrium occupancy is
The rate constant for approach to equilibrium is \(k_{\mathrm{on}}C+k_{\mathrm{off}}=(1.00)(0.0630)+0.200=0.263\ \mathrm{h}^{-1}\). Thus, at 2.00 hours,
Occupancy is about 9.80% at 2.00 hours, even though the eventual equilibrium occupancy is about 23.95%. The equilibrium formula is not wrong; it answers what occupancy would be after sufficient time at the maintained concentration, rather than what occupancy has been reached by this particular time.
Same Affinity, Different Time Courses
Two compounds can have the same \(K_D\) but different values of both rate constants. In that case, their equilibrium occupancy at a given free concentration is the same, but their rates of association and dissociation can differ substantially. This is important when the MABEL assessment concerns early exposure, a short infusion, a fluctuating concentration profile, or occupancy after drug concentrations decline.
Worked Example: Comparing Fast and Slow Binding at the Same Affinity
Compare two hypothetical compounds at a constant concentration of \(1.00\ \mathrm{nM}\), starting from unoccupied target. Compound A has \(k_{\mathrm{on}}=0.100\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\) and \(k_{\mathrm{off}}=1.00\ \mathrm{h}^{-1}\). Compound B has \(k_{\mathrm{on}}=0.00100\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\) and \(k_{\mathrm{off}}=0.0100\ \mathrm{h}^{-1}\).
Both have \(K_D=10.0\ \mathrm{nM}\), so both have equilibrium occupancy \(1.00/(1.00+10.0)=0.09091\), or 9.091%. After 2.00 hours, Compound A has \(k_{\mathrm{on}}C+k_{\mathrm{off}}=1.10\ \mathrm{h}^{-1}\), giving
For Compound B, the corresponding rate is \(0.0110\ \mathrm{h}^{-1}\), so
At 2.00 hours, predicted occupancy is about 8.08% for A and 0.198% for B, despite identical affinity. Compound B approaches the same eventual occupancy much more slowly. For a sufficiently long exposure, this early difference diminishes; for a short exposure, it may be central to the interpretation.
The same rate constants also govern binding after concentration falls. When \(C=0\), the model reduces to \(d\theta/dt=-k_{\mathrm{off}}\theta\), so occupancy decays exponentially with half-life \(\ln(2)/k_{\mathrm{off}}\), provided the other model assumptions remain reasonable. A small \(k_{\mathrm{off}}\) can therefore sustain occupancy during washout, even if free-drug concentration is low.
Worked Example: A Brief Exposure and Occupancy During Washout
A hypothetical compound has \(k_{\mathrm{on}}=0.100\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\) and \(k_{\mathrm{off}}=1.00\ \mathrm{h}^{-1}\). Free concentration is \(10.0\ \mathrm{nM}\) for 0.500 hours and then falls to zero. Starting from \(\theta(0)=0\), calculate occupancy at the end of the exposure and one hour after washout begins.
During exposure, \(K_D=10.0\ \mathrm{nM}\), \(\theta_{\mathrm{eq}}=10.0/(10.0+10.0)=0.500\), and \(k_{\mathrm{on}}C+k_{\mathrm{off}}=2.00\ \mathrm{h}^{-1}\). Therefore,
After concentration falls to zero, occupancy decays with \(k_{\mathrm{off}}=1.00\ \mathrm{h}^{-1}\). One hour later,
Thus occupancy is about 31.6% at the end of the brief exposure and remains about 11.6% one hour after washout. The calculation assumes concentration falls immediately to zero; a gradual decline would require using the full time-varying concentration model.
Using Kinetics in an Occupancy-Based MABEL Estimate
A MABEL estimate requires a pharmacologically justified low-effect level, not just a binding calculation. For illustration, suppose prior mechanistic evidence supports using 10% target occupancy at a defined observation time as the low-effect threshold. The kinetic calculation can then find the exposure that reaches that threshold under a specified concentration-time scenario. The threshold and observation time must be scientifically justified for the product; neither 10% nor any particular time is universal.
Worked Example: A Time-Specific Occupancy Threshold
For Compound B from the comparison example, let \(K_D=10.0\ \mathrm{nM}\), \(k_{\mathrm{on}}=0.00100\ \mathrm{nM}^{-1}\mathrm{h}^{-1}\), and \(k_{\mathrm{off}}=0.0100\ \mathrm{h}^{-1}\). Assume constant free concentration for 2.00 hours and initially unoccupied target. Find the concentration that gives 10.0% occupancy at 2.00 hours. Compare with the equilibrium-only calculation.
At equilibrium, 10.0% occupancy corresponds to
That concentration is not predicted to produce 10.0% occupancy within two hours for this slowly associating compound. The time-specific calculation is
Solving this equation numerically gives approximately \(C=53.2\ \mathrm{nM}\). Substitution checks the result: the equilibrium fraction at that concentration is \(53.2/(53.2+10.0)=0.8418\), and the approach rate is \(0.00100(53.2)+0.0100=0.0632\ \mathrm{h}^{-1}\). Hence,
The equilibrium-only concentration is much lower because it describes eventual occupancy, not occupancy reached by two hours from an initially unoccupied state. If, purely hypothetically, a specified human pharmacokinetic scenario maintained \(0.0200\ \mathrm{nM}\) per \(\mathrm{\mu g/kg}\) of dose over this interval, the time-specific concentration would map to \(53.2/0.0200=2660\ \mathrm{\mu g/kg}\). This is a model-based exposure-to-dose conversion for the stated scenario, not a clinical dose recommendation. The result depends strongly on the assumed concentration profile, kinetic parameters, and occupancy threshold.
Assumptions, Diagnostics, and Alternatives
A kinetic occupancy estimate is only as credible as the binding and exposure inputs. The simple model assumes reversible one-to-one binding, constant kinetic parameters, a known free concentration at the target site, and no important changes in target abundance or accessibility over the modeled interval. It also treats the measured concentration as representative of concentration at the target. These assumptions may be poor for heterogeneous tissues, multivalent biologics, cooperative binding, target turnover, or concentration-dependent internalization.
Useful diagnostics include checking whether occupancy measurements over time are consistent with the predicted rise and washout, whether the sampling schedule can resolve the estimated rate scale, and whether the assumed target-site concentration profile is plausible. A single equilibrium binding measurement may constrain \(K_D\) but cannot uniquely determine \(k_{\mathrm{on}}\) and \(k_{\mathrm{off}}\). Rate estimates therefore need suitable kinetic data, with attention to assay conditions and uncertainty. If observed occupancy rises or falls systematically faster or slower than predicted, the model may be missing processes or using unrepresentative parameters.
| Approach | What it represents | When it can be useful | Main limitation |
|---|---|---|---|
| Affinity-only occupancy | \(\theta=C/(C+K_D)\) at equilibrium | Binding is near equilibrium over the relevant exposure interval | Does not predict lag or persistence |
| Kinetic occupancy | Occupancy changes according to \(k_{\mathrm{on}}\), \(k_{\mathrm{off}}\), and \(C(t)\) | Timing, brief exposure, or washout affects interpretation | Requires rate and concentration-time inputs and model assumptions |
| Direct functional concentration-response | Measured biological response over tested concentrations and assay conditions | A relevant functional assay directly informs the low-effect region | Assay context and time course may not represent human target engagement |
The affinity-only approach is not a competing physical theory; it is the equilibrium limit of the kinetic model. A full kinetic model is also not automatically preferable if its parameters are poorly identified. Where data permit, compare predicted occupancy with direct time-course measurements and with functional evidence. Propagate uncertainty in binding rates and human exposure rather than presenting a single threshold as exact.
Binding kinetics matter to MABEL because the relevant early pharmacology may depend on when target engagement occurs and how long it persists, not only on the eventual fraction bound. A defensible estimate states the occupancy or functional threshold, the time window, the concentration profile, and the kinetic assumptions. It then checks whether available data support those choices before converting an exposure target into a dose.
Check Your Understanding
Use the kinetic model and assumptions in this tutorial. Solutions are included below the questions.
- If \(k_{\mathrm{off}}\) doubles while \(k_{\mathrm{on}}\) remains fixed, what happens to \(K_D\) and, generally, to equilibrium occupancy at fixed concentration?
- A compound has \(k_{\mathrm{off}}=0.200\ \mathrm{h}^{-1}\). What is its occupancy half-life after free concentration falls to zero?
- Two compounds have the same \(K_D\), but one has smaller \(k_{\mathrm{on}}\) and proportionally smaller \(k_{\mathrm{off}}\). Which has greater equilibrium occupancy at the same concentration, and which tends to approach equilibrium more slowly?
- Why can an equilibrium concentration threshold be inappropriate when defining a low-occupancy threshold two hours after exposure begins?
Solutions. (1) \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\), so doubling \(k_{\mathrm{off}}\) doubles \(K_D\); at fixed \(C\), \(C/(C+K_D)\) decreases. (2) The half-life is \(\ln(2)/k_{\mathrm{off}}=0.6931/0.200=3.47\ \mathrm{h}\). (3) Their equilibrium occupancies are equal because their \(K_D\) values are equal. The compound with smaller rate constants generally approaches equilibrium more slowly, as indicated by the smaller \(k_{\mathrm{on}}C+k_{\mathrm{off}}\). (4) The equilibrium calculation describes eventual occupancy under maintained concentration. At two hours, occupancy may still be changing, so the threshold must be evaluated with the time-dependent binding equation and the specified exposure history.