When Equilibrium Occupancy Is a Valid Approximation

Why Compare Timescales?

The equilibrium occupancy equation gives the fraction of receptors bound when free-ligand concentration is held constant long enough for binding to relax. In a pharmacology model, however, free ligand may change before that relaxation is complete. Using equilibrium occupancy at every time point is then an approximation: it treats the bound fraction as if it responded immediately to the current ligand concentration.

The central check is whether binding relaxes much faster than free-ligand concentration changes. This is a timescale comparison, not a claim that equilibrium is reached exactly. Whether the approximation is useful also depends on the occupancy error that can be tolerated for the intended model or decision.

This definition assumes one reversible binding site, constant rate constants, a stable receptor pool, and free ligand that is prescribed rather than depleted by binding. Occupancy is the fraction of available receptors that are bound, denoted \(\theta\), with \(0\leq\theta\leq1\). The dissociation constant is \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\), and the equilibrium occupancy at fixed \(L\) is \(\theta_{\mathrm{eq}}=L/(L+K_D)\). The relaxation time is not the same as the dissociation residence time \(1/k_{\mathrm{off}}\): association contributes to relaxation whenever ligand is present.

What the Relaxation Time Tells Us

For a constant free-ligand concentration, the time-dependent one-site equation is \(d\theta/dt=k_{\mathrm{on}}L(1-\theta)-k_{\mathrm{off}}\theta\). Regrouping its terms gives a direct description of the approach to equilibrium:

The solution from an initial occupancy \(\theta(0)\) is

This solution shows why \(\tau_{\mathrm{relax}}\) is a characteristic, rather than a time to exact equilibrium. After one relaxation time, the remaining difference from equilibrium is \(e^{-1}\), or about 36.8%, of its initial value. After three relaxation times it is about 5.0%; after five it is about 0.67%. The percentages apply to the initial occupancy gap, not to occupancy itself.

Comparing Binding and Ligand-Change Timescales

A useful local measure of how quickly a positive free-ligand concentration changes is its relative rate of change, \(|dL/dt|/L\). Its reciprocal, \(\tau_L=L/|dL/dt|\), is the time required for a change comparable in size to the current concentration if the instantaneous rate persisted. Compare this with the binding relaxation rate using the dimensionless ratio

When \(\varepsilon\) is much smaller than one, binding is fast relative to the local fractional change in ligand. This supports an equilibrium approximation, subject to an error tolerance. If the ratio is not small, a kinetic occupancy model is generally needed to represent the delay in binding. The ratio is a screening measure, not a universal cutoff: a model requiring accuracy within 0.01 occupancy units needs a stricter check than one for which a 0.1-unit discrepancy is immaterial.

The relative-change timescale is undefined at \(L=0\) and can become uninformative near zero. In those situations, inspect the predicted change in equilibrium occupancy directly, or compare the absolute ligand change over a binding relaxation time with a relevant concentration scale. One useful local check is how much the equilibrium occupancy itself would change during one relaxation time.

Justifying the Slow-Change Approximation

The timescale argument can be made more precise by allowing \(L\) to vary while keeping the one-site assumptions. At each time, define the instantaneous equilibrium value \(\theta_{\mathrm{eq}}(t)=L(t)/(L(t)+K_D)\), and the tracking error \(e(t)=\theta(t)-\theta_{\mathrm{eq}}(t)\). Substituting \(\theta=e+\theta_{\mathrm{eq}}\) into the kinetic equation gives

The error is pulled toward zero at rate \(\lambda(t)\), while changes in the instantaneous equilibrium value continually create tracking error. The exact solution from time \(0\) is

The first term is the decaying error present at the start. The second is the accumulated effect of equilibrium occupancy moving while binding catches up. This expression gives the justification for the comparison: when \(\theta_{\mathrm{eq}}\) changes slowly relative to relaxation, each recent contribution is small and older contributions decay.

For a local, slowly varying situation, the error is approximately \(e\approx-(d\theta_{\mathrm{eq}}/dt)/\lambda\). This approximation is most useful when the rates are changing gradually over several relaxation times; it is not an exact formula for arbitrary input. Since \(d\theta_{\mathrm{eq}}/dt=K_D(dL/dt)/(L+K_D)^2\), this error estimate also incorporates how sensitively occupancy responds at the current concentration. For example, changing \(L\) rapidly at a concentration where equilibrium occupancy is already near one may produce less occupancy error than the same absolute concentration change near the steep part of the occupancy curve.

Choosing Between Equilibrium and Kinetic Occupancy

The equilibrium approximation replaces dynamic occupancy with the algebraic value \(L/(L+K_D)\). It is computationally simple and appropriate when binding relaxation is fast enough for the required accuracy. The kinetic equation retains occupancy as a state and represents its finite response time; it is preferable when ligand changes rapidly, when initial occupancy is far from equilibrium, or when a stated tolerance is not met. A useful workflow is:

The comparison assumes the free concentration is known independently of receptor binding. If ligand depletion makes free concentration depend on bound ligand, or if receptor availability changes, this simple check does not establish that the equilibrium one-site approximation is adequate. Likewise, multiple binding states, time-varying rate constants, transport delays, and downstream response delays require additional model structure. The approximation concerns receptor occupancy only; it does not guarantee that a downstream pharmacological response is at equilibrium.

Diagnostics and Common Pitfalls

• Using residence time as relaxation time. The residence time \(1/k_{\mathrm{off}}\) describes dissociation after binding. At fixed positive ligand concentration, the relaxation rate is \(k_{\mathrm{on}}L+k_{\mathrm{off}}\), so association also affects adjustment toward equilibrium.
• Calling a small timescale ratio proof of accuracy. The ratio is a screening diagnostic. An initial occupancy gap can persist for several relaxation times, and acceptable error depends on the model's purpose. Check transient error and state a tolerance.
• Confusing movement with remaining error. Under a constant step, the occupancy change from its starting value is the initial equilibrium gap multiplied by \(1-e^{-\lambda t}\). The remaining distance to equilibrium is the gap multiplied by \(e^{-\lambda t}\). These are different quantities.
• Using relative ligand change near zero. The ratio \(|dL/dt|/L\) becomes undefined at zero and can be very large near zero. Inspect the expected occupancy change or use a meaningful absolute concentration scale instead.
• Assuming a good occupancy approximation validates the whole model. This test does not address ligand depletion, receptor turnover, changing binding properties, or delays between occupancy and response.

In practice, compare the two approaches at representative concentrations and transition rates when the occupancy tolerance is important. If the equilibrium and kinetic calculations differ materially over the times or concentrations of interest, retain kinetic occupancy rather than interpreting equilibrium occupancy as an instantaneous response.