Deriving Fractional Target Occupancy at Equilibrium

Why Fractional Occupancy Is Useful

A concentration of drug does not, by itself, say what fraction of a target is bound. That fraction depends on the concentration of ligand available to bind and on the binding equilibrium. Fractional occupancy provides a compact link between those quantities: it describes the proportion of binding-competent target molecules that are ligand-bound at equilibrium.

This relationship is useful in mechanism-based PK/PD models because it gives a specific, interpretable connection between free drug concentration and target engagement. It is not a complete model of the resulting pharmacological response: occupancy alone does not determine the size, direction, or time course of a downstream effect.

Definition: Let $R$ denote free target, $L$ free ligand, and $RL$ the ligand-bound target complex. The fractional target occupancy, $\theta$, is the fraction of the binding-competent target pool that is bound: $$ \theta=\frac{[RL]}{[R]+[RL]}. $$ Here brackets denote concentrations. Thus $0\leq\theta\leq1$, and $\theta$ has no units.

The equation in this tutorial applies to a reversible, one-to-one interaction: one ligand molecule binds one target site, and the target site is either free or occupied by that ligand. The calculation uses the free ligand concentration, not automatically the administered dose, a nominal concentration, or the total ligand concentration in a sample.

Mass-Action Equilibrium and the Occupancy Equation

For the binding reaction $R+L\rightleftharpoons RL$, the equilibrium dissociation constant is defined by the ratio of free target and free ligand concentrations to bound complex concentration at equilibrium. In the concentration-based form used here, it has units of concentration.

Definition: For one-to-one binding at equilibrium, the dissociation constant is $$ K_D=\frac{[R][L]}{[RL]}. $$ Equivalently, the bound-to-free-target concentration ratio is $[RL]/[R]=[L]/K_D$. This tutorial treats concentrations as suitable approximations to the corresponding activities, as is standard for sufficiently dilute solutions.

The target pool in the occupancy definition consists of free target plus ligand-bound target. Write its concentration as $[R]+[RL]$. Starting from the mass-action equilibrium relation, divide numerator and denominator in the definition of $\theta$ by $[RL]$. This expresses occupancy in terms of the ratio of free target to bound target.

\begin{aligned} K_D&=\frac{[R][L]}{[RL]}\\ \frac{[R]}{[RL]}&=\frac{K_D}{[L]}\\ \theta &=\frac{[RL]}{[R]+[RL]}\\ &=\frac{1}{[R]/[RL]+1}\\ &=\frac{1}{K_D/[L]+1}\\ &=\frac{[L]}{K_D+[L]}. \end{aligned}

Equilibrium one-to-one occupancy equation: If the stated binding and equilibrium assumptions hold, fractional target occupancy is $$ \theta=\frac{[L]}{K_D+[L]}, $$ where $[L]$ is the free ligand concentration at equilibrium and $K_D$ is the equilibrium dissociation constant for the same ligand-target interaction, expressed in the same concentration units.

The derivation uses the equilibrium ratio and the definition of the target pool; it does not require the total target concentration to be known. If total binding-competent target is known, the bound concentration can then be calculated as $[RL]=\theta([R]+[RL])$. The equation also does not require that ligand depletion be negligible when the free ligand concentration is measured or otherwise known. What it does require is that the concentration supplied as $[L]$ actually represents free ligand at equilibrium.

Reading the Equation and Its Limits

The dimensionless ratio $[L]/K_D$ determines occupancy. If $[L]=K_D$, the denominator is $2K_D$, so $\theta=0.5$. When free ligand is much smaller than $K_D$, occupancy is low; when free ligand is much larger than $K_D$, occupancy approaches, but does not exceed, one. These are limiting behaviors of the equation, not guarantees that the binding model is appropriate in a particular system.

[L]=K_D\ \Longrightarrow\ \theta=\frac{1}{2}, \qquad [L]\ll K_D\ \Longrightarrow\ \theta\approx\frac{[L]}{K_D}, \qquad [L]\gg K_D\ \Longrightarrow\ \theta\approx1.

The same equation can be rearranged to determine the free ligand concentration associated with a desired fractional occupancy. For $0\leq\theta<1$, multiply through by $K_D+[L]$, then collect the terms containing $[L]$:

\theta(K_D+[L])=[L] \quad\Longrightarrow\quad \theta K_D=[L](1-\theta) \quad\Longrightarrow\quad [L]=\frac{\theta}{1-\theta}K_D.

This form makes an important feature explicit: a specified occupancy below one requires a finite free concentration under the model, whereas exact full occupancy is only approached as free ligand increases without bound. It is an equilibrium concentration calculation, not a dosing rule. Translating a required free concentration into a dose would require additional pharmacokinetic and biological information.

Worked Example: Calculating Occupancy From Free Concentration

Suppose a hypothetical target-binding system has $K_D=12$ nM and a measured free ligand concentration of $8.0$ nM at equilibrium. Calculate fractional occupancy. Assume the measured ligand is free, the binding is one-to-one, and the target pool has only free and ligand-bound forms.

The units match, so substitute the values directly:

\theta=\frac{[L]}{K_D+[L]} =\frac{8.0\ {\rm nM}}{12\ {\rm nM}+8.0\ {\rm nM}} =\frac{8.0}{20.0} =0.400.

The predicted fractional occupancy is $0.400$, or $40.0\%$. As a consistency check, the mass-action ratio implies $[RL]/[R]=[L]/K_D=8.0/12=0.667$. A bound-to-free ratio of $0.667$ corresponds to a bound fraction of $0.667/(1+0.667)=0.400$, matching the result.

If the hypothetical total binding-competent target concentration were $0.200$ nM, then the bound concentration would be $0.400\times0.200=0.0800$ nM, and free target would be $0.200-0.0800=0.120$ nM. These values also satisfy the defining relation: $K_D=(0.120\ {\rm nM})(8.0\ {\rm nM})/(0.0800\ {\rm nM})=12.0$ nM.

Worked Example: Finding the Free Concentration for a Target Occupancy

In a hypothetical system with $K_D=4.0$ nM, what free ligand concentration corresponds to $75.0\%$ occupancy under the one-to-one equilibrium model? Assume the equilibrium dissociation constant and the target definition are appropriate for the system.

Write the desired fraction as $\theta=0.750$ and use the rearranged equation:

[L]=\frac{\theta}{1-\theta}K_D =\frac{0.750}{1-0.750}(4.0\ {\rm nM}) =\frac{0.750}{0.250}(4.0\ {\rm nM}) =3(4.0\ {\rm nM}) =12\ {\rm nM}.

Check by substituting $12$ nM into the occupancy equation:

\theta=\frac{12\ {\rm nM}}{4.0\ {\rm nM}+12\ {\rm nM}} =\frac{12}{16} =0.750.

Thus $12$ nM free ligand, not necessarily $12$ nM total ligand or a particular administered dose, corresponds to $75.0\%$ occupancy under the stated assumptions.

Assumptions, Diagnostics, and Alternatives

Before interpreting a calculated value as target engagement, check that the model matches the biological and measurement context. The one-to-one equation assumes a single effective class of independent binding sites, one ligand per site, reversible binding, and equilibrium. It also assumes that the target pool used to define occupancy is correctly represented by free plus ligand-bound target. If some target is unavailable or is in another state, the calculated fraction may not describe occupancy of all target present.

Equilibrium is a substantive condition: the free ligand and bound complex must have had enough time to reach the equilibrium relation being used. A sample taken during changing exposure may not satisfy it. This tutorial's equation does not describe how quickly occupancy changes, nor does it predict occupancy during a transient. Similarly, if the system has multiple ligands, cooperative binding, irreversible binding, or distinct target populations with different binding behavior, the simple one-to-one relation may not capture the observed pattern.

Useful diagnostics include checking that all concentrations use the same units, verifying that $[L]$ is free concentration at the relevant site and time, and comparing predicted occupancy with an independent occupancy measurement when available. A predicted value outside $0$ to $1$, or one that changes substantially when units are converted consistently, signals a calculation or implementation error. Persistent disagreement with measured occupancy can instead indicate a wrong free concentration, lack of equilibrium, assay limitations, or an inadequate binding model; it is not resolved by treating the formula as universally applicable.

A direct occupancy assay and an equilibrium calculation answer related but different questions. A validated assay estimates occupancy from observed binding or a surrogate measurement, subject to its own measurement assumptions. The equation predicts occupancy from free ligand concentration and $K_D$, conditional on the binding model. Agreement can support the model's adequacy in that setting, while disagreement is a reason to investigate inputs and assumptions. Neither approach alone establishes that occupancy causes a particular clinical response.

Worked Example: Why Total Concentration Is Not a Substitute for Free Concentration

Consider a hypothetical equilibrium sample for which an independent measurement reports $10.0$ nM free ligand, while the total ligand concentration is $100$ nM. Let $K_D=10.0$ nM. Assume the equilibrium one-to-one model applies and use the measured free concentration.

The appropriate prediction is:

\theta=\frac{10.0\ {\rm nM}}{10.0\ {\rm nM}+10.0\ {\rm nM}} =\frac{10.0}{20.0} =0.500.

The predicted occupancy is $50.0\%$. Substituting total concentration instead would give $100/(10+100)=0.909$, or $90.9\%$. That is not a second valid prediction from this equation: it uses the wrong input. The formula requires free ligand concentration. The difference illustrates why the concentration definition must be checked before calculating occupancy; total concentration should not be used as a shortcut for free concentration.

Common Pitfall

Wrong reasoning: “The ligand concentration equals $K_D$, so the target is fully occupied.” This confuses the dissociation constant with a concentration that saturates every binding site. At $[L]=K_D$, the equation gives $\theta=K_D/(K_D+K_D)=0.500$, not $1$.

Correction: $K_D$ is the free concentration that gives half occupancy in this one-to-one equilibrium model. Occupancy increases as free ligand increases, but exact occupancy of $1$ is a limiting value, not a value reached at a finite concentration. Also check whether the stated concentration is free: an administered or total concentration cannot simply be inserted in place of $[L]$.

Key takeaway: For reversible one-to-one binding at equilibrium, use $\theta=[L]/(K_D+[L])$, with free ligand concentration and $K_D$ in matching units. The result is a model-based fraction of the defined target pool; its interpretation depends on the binding assumptions and on having the correct free concentration.

Check Your Understanding

Answer each question using the one-to-one equilibrium model where its assumptions are stated. Short solutions follow the questions.

For $K_D=6.0$ nM and free ligand concentration $6.0$ nM, what is the predicted fractional occupancy?
For $K_D=2.0$ nM, calculate the free ligand concentration required for $80.0\%$ occupancy.
A calculation gives $\theta=1.20$. Is that a possible fractional occupancy under this model? Give one likely source of the problem.
Why is it not valid to substitute total ligand concentration for free ligand concentration without justification?

Solutions. (1) $\theta=6.0/(6.0+6.0)=0.500$, or $50.0\%$. (2) $[L]=[0.800/(1-0.800)](2.0\ {\rm nM})=4(2.0\ {\rm nM})=8.0$ nM. (3) No. The equation gives a value between zero and one for nonnegative free concentration and positive $K_D$; check the formula, units, and inputs. (4) The equation is defined in terms of free ligand at equilibrium. Total concentration may differ, so substituting it can produce a prediction that does not represent occupancy at the supplied free concentration.