The Time-Dependent Occupancy Equation

Why Occupancy Changes Over Time

An equilibrium occupancy formula describes a limiting state, but a receptor population may begin far from that state. After free ligand is introduced, ligand-receptor complexes form through association and disappear through dissociation. The balance between these two processes determines how quickly occupancy changes and where it tends under constant conditions.

This tutorial derives the time-dependent equation for a single reversible binding site when the free-ligand concentration is held constant. The central result is an exponential approach from the initial occupancy toward the equilibrium occupancy. Its derivation also provides useful checks: the solution must start at the specified initial value, remain between zero and one when the initial occupancy is physically valid, and approach the equilibrium level.

Model and Assumptions

Let \(R\) denote unbound receptor and \(RL\) the ligand-bound receptor complex. Let \(R_{\mathrm{tot}}\) be the total concentration of binding-competent receptor, \(L\) the free-ligand concentration, and \(\theta(t)\) the fraction of receptors occupied at time \(t\). Thus \(\theta(t)=[RL](t)/R_{\mathrm{tot}}\), and the unoccupied fraction is \(1-\theta(t)\).

The model assumes a single homogeneous class of independent binding sites; one ligand molecule binds each site; association and dissociation follow mass-action kinetics; and the total receptor concentration is constant. It also assumes that free-ligand concentration \(L\) is constant over the time interval being modeled. In particular, treating \(L\) as constant requires either that it is controlled externally or that binding does not appreciably deplete it. The model does not include receptor synthesis, loss, internalization, cooperative binding, or competing ligands.

When the occupied concentration is \([RL]\), the association flux is \(k_{\mathrm{on}}L[R]\), and the dissociation flux is \(k_{\mathrm{off}}[RL]\). Since \([R]=R_{\mathrm{tot}}-[RL]\), the net change in occupied receptor is association minus dissociation. Dividing by the constant \(R_{\mathrm{tot}}\) converts this balance into an equation for fractional occupancy:

The product \(k_{\mathrm{on}}L\) has units of time\(^{-1}\), as does \(k_{\mathrm{off}}\), so both terms on the right have units of occupancy per time. Define \(k_{\mathrm{obs}}=k_{\mathrm{on}}L+k_{\mathrm{off}}\), which is positive for positive rates and nonnegative ligand concentration. The initial-value problem is therefore a linear differential equation:

Deriving the Time-Dependent Solution

The constant value that makes the derivative zero is \(\theta_{\infty}=k_{\mathrm{on}}L/k_{\mathrm{obs}}\). This is the limiting occupancy under the stated constant conditions. Subtract this constant from the time-varying occupancy by defining \(u(t)=\theta(t)-\theta_{\infty}\). Substitution into the differential equation gives \(du/dt=-k_{\mathrm{obs}}u\).

An integrating factor solves this equation without requiring division by \(u\), so the derivation also covers the case \(u(0)=0\). Multiply both sides by \(e^{k_{\mathrm{obs}}t}\); the left side is the derivative of a product:

Thus \(e^{k_{\mathrm{obs}}t}u(t)\) is constant and equals \(u(0)=\theta_0-\theta_{\infty}\). This yields the solution for every initial occupancy, including \(\theta_0=\theta_{\infty}\), when \(u(t)\) is identically zero.

The form of \(\theta_{\infty}\) agrees with the one-site equilibrium occupancy equation: writing \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\), it is \(L/(L+K_D)\) when \(k_{\mathrm{on}}>0\). The solution makes the time course explicit as well. The deviation from the limiting value is the initial deviation multiplied by \(e^{-k_{\mathrm{obs}}t}\). If occupancy begins below the limit, it rises toward it; if it begins above the limit, it falls toward it. Under these assumptions it does not overshoot.

How Association and Dissociation Shape the Time Course

The limiting occupancy and the speed of approach are determined by related but distinct combinations of the rates. At a fixed \(L\), the ratio of association to total rate contributions sets \(\theta_{\infty}\), while their sum \(k_{\mathrm{obs}}\) sets how rapidly the initial difference shrinks. Consequently, two rate pairs can produce the same limiting occupancy but different time courses.

Checks, Alternatives, and Limitations

The equation provides several practical checks on calculations and model use:

• Initial condition: setting \(t=0\) in the solution gives \(\theta(0)=\theta_{\infty}+(\theta_0-\theta_{\infty})=\theta_0\).
• Limiting behavior: for \(k_{\mathrm{obs}}>0\), \(e^{-k_{\mathrm{obs}}t}\) tends to zero as time increases, so \(\theta(t)\) tends to \(\theta_{\infty}\).
• Direction of change: the derivative at any time is \(d\theta/dt=k_{\mathrm{obs}}(\theta_{\infty}-\theta)\). Its sign points toward the limiting occupancy.
• Physical bounds: if \(k_{\mathrm{on}}L\geq0\), \(k_{\mathrm{off}}\geq0\), and \(k_{\mathrm{obs}}>0\), then \(0\leq\theta_{\infty}\leq1\). If \(0\leq\theta_0\leq1\), the solution is a weighted average of \(\theta_0\) and \(\theta_{\infty}\), so it remains in that interval.

For this constant-condition model, the analytical solution is preferable to numerical integration: it is exact and makes the effects of initial occupancy and rates transparent. Numerical integration can be useful when the model is expanded to include additional states or processes, but it should reproduce these checks when applied to this simple equation. An equilibrium-only calculation gives \(\theta_{\infty}\) but not the transient occupancy at a specified time. It is therefore not a substitute when the time course itself is the question.

A common mistake is to use \(\theta_{\infty}\) as though it were the occupancy at every time after ligand exposure. That is only correct if the system starts at the limiting occupancy or has had sufficient time to approach it. Another is to interpret \(k_{\mathrm{off}}\) alone as the speed of the full occupancy response: under constant free ligand, the exponential rate is \(k_{\mathrm{on}}L+k_{\mathrm{off}}\), not \(k_{\mathrm{off}}\) by itself.

The formula can fail as a description of a real time course if free ligand is not constant, binding appreciably depletes ligand, receptor abundance changes, or the receptor population has multiple binding behaviors. In such cases the stated assumptions are not met, and a different model may be required. The solution is also a deterministic prediction for the specified parameters and initial condition; it does not itself represent measurement error or uncertainty in those inputs.