Tutorials › Quantitative Systems Pharmacology (QSP) › Ligand Depletion and the Finite-Target Binding Balance

Receptor Binding and Target Engagement · Tutorial 194 of 799

Ligand Depletion and the Finite-Target Binding Balance

Learn to solve the one-site equilibrium binding balance using total ligand and receptor concentrations, and to judge when ligand depletion makes the free-ligand approximation unreliable.

Intermediate 11 min read

What You'll Learn

  • Distinguish total ligand from free ligand and bound ligand in a closed binding system.
  • Derive the quadratic mass-balance equation for one-to-one receptor binding at equilibrium.
  • Select the physically meaningful solution and verify it against conservation and concentration bounds.
  • Calculate receptor occupancy when ligand depletion is appreciable.
  • Assess when treating total ligand as free ligand is a reasonable approximation.
  • Recognize the assumptions and diagnostics that limit use of the equilibrium balance.

Why Total Ligand Is Not Always Free Ligand

An equilibrium occupancy formula uses the free ligand concentration, because free ligand is what participates directly in receptor binding. In an experiment or model, however, the quantity specified may instead be the total amount of ligand added. If receptors bind a substantial fraction of that amount, total ligand and free ligand differ. Using the total concentration as though it were free can then overstate receptor occupancy.

The key is to solve binding and conservation together. For a single reversible binding site, equilibrium relates free ligand, free receptor, and ligand–receptor complex. Two mass balances relate those species to the total ligand and receptor concentrations. Combining the equilibrium relation with both balances gives the finite-target binding balance.

Definition: Let \(L\) be free ligand concentration, \(R\) free receptor concentration, and \(C\) concentration of ligand–receptor complex. Let \(L_T\) and \(R_T\) be total ligand and total receptor concentrations. In a closed system with one-to-one binding, \(L_T=L+C\) and \(R_T=R+C\). The dissociation constant is \(K_d=LR/C\) at equilibrium. The fraction of receptors occupied is \(C/R_T\), when \(R_T>0\).

All concentrations and \(K_d\) must use the same concentration units. This balance assumes a homogeneous system, one reversible ligand-binding site per receptor, 1:1 complex formation, no other processes that consume or produce ligand or receptor during equilibration, and enough time for binding to reach equilibrium. It also assumes that the stated totals refer to the same system volume. If these conditions do not apply, the equations may not describe the system.

Deriving the Finite-Target Balance

The mass balances imply \(R=R_T-C\) and \(L=L_T-C\). Substitute these expressions into the equilibrium definition:

$$ K_d=\frac{(L_T-C)(R_T-C)}{C}. $$

Multiplying by \(C\) and collecting terms gives a quadratic in the complex concentration:

$$ C^2-(L_T+R_T+K_d)C+L_TR_T=0. $$

The quadratic formula gives two mathematical roots. The physically admissible root is the smaller one: it is the only root that remains between zero and \(\min(L_T,R_T)\). Thus the bound concentration is

$$ C=\frac{L_T+R_T+K_d-\sqrt{(L_T+R_T+K_d)^2-4L_TR_T}}{2}. $$

The other root is not physically possible because it exceeds at least one of the total concentrations, implying a negative free concentration. Once \(C\) is found, the free concentrations and occupancy follow directly from the balances.

$$ L=L_T-C,\qquad R=R_T-C,\qquad \text{occupancy}=\frac{C}{R_T}. $$

An equivalent expression solves directly for free ligand. Substituting \(C=L_T-L\) and \(R=R_T-L_T+L\) into the equilibrium relation yields

$$ L^2+(R_T+K_d-L_T)L-K_dL_T=0, \qquad L=\frac{L_T-R_T-K_d+\sqrt{(L_T-R_T-K_d)^2+4K_dL_T}}{2}. $$

This root is nonnegative and no greater than \(L_T\). Either form gives the same solution; the complex form can be convenient for occupancy, while the free-ligand form makes the distinction between total and free ligand explicit.

Worked Examples

Worked Example: Ligand Is Strongly Depleted

Suppose a hypothetical closed system has \(L_T=5\) nM, \(R_T=50\) nM, and \(K_d=5\) nM. Find the free ligand, complex, free receptor, and receptor occupancy at equilibrium.

Using the free-ligand root,

$$ L=\frac{5-50-5+\sqrt{(5-50-5)^2+4(5)(5)}}{2} =\frac{-50+\sqrt{2600}}{2} \approx 0.4951\ \mathrm{nM}. $$

The radicand is \(2600\), since \((-50)^2+4(5)(5)=2500+100\). The complex and free receptor concentrations are

$$ C=L_T-L=5-0.4951=4.5049\ \mathrm{nM}, \qquad R=R_T-C=50-4.5049=45.4951\ \mathrm{nM}. $$

Therefore,

$$ \text{occupancy}=\frac{4.5049}{50}=0.09010\approx 9.010\%. $$

The balances check: \(0.4951+4.5049=5.0000\) nM ligand and \(45.4951+4.5049=50.0000\) nM receptor. In contrast, treating \(L_T=5\) nM as free ligand would give occupancy \(5/(5+5)=0.5\) and predict \(25\) nM bound receptor. That bound amount exceeds the \(5\) nM total ligand, so the approximation is plainly impossible here.

Worked Example: A Ligand Excess Makes the Approximation Closer

Consider hypothetical totals \(L_T=100\) nM and \(R_T=10\) nM with \(K_d=10\) nM. Calculate the exact equilibrium and compare it with the approximation \(L\approx L_T\).

The exact free ligand is

$$ L=\frac{100-10-10+\sqrt{(100-10-10)^2+4(10)(100)}}{2} =\frac{80+\sqrt{10400}}{2} \approx 90.9902\ \mathrm{nM}. $$

Hence \(C=100-90.9902=9.0098\) nM and \(R=10-9.0098=0.9902\) nM. The exact occupancy is

$$ \frac{C}{R_T}=\frac{9.0098}{10}=0.9010, $$

or about \(90.10\%\). With the free-ligand approximation, the one-site equilibrium formula gives occupancy \(100/(10+100)=0.9091\), or \(90.91\%\), and predicted complex concentration \(10(0.9091)=9.091\) nM. Here the approximation is fairly close: it overstates complex by about \(0.0812\) nM, or about \(0.90\%\) relative to the exact complex. The exact solution still matters if that difference is important to the application.

Worked Example: Equal Total Ligand and Receptor

For hypothetical \(L_T=R_T=20\) nM and \(K_d=5\) nM, calculate the equilibrium occupancy. This case illustrates that neither total concentration can automatically be treated as free ligand.

The free-ligand solution is

$$ L=\frac{20-20-5+\sqrt{(20-20-5)^2+4(5)(20)}}{2} =\frac{-5+\sqrt{425}}{2} \approx 7.8078\ \mathrm{nM}. $$

Thus \(C=20-7.8078=12.1922\) nM and \(R=20-12.1922=7.8078\) nM. Occupancy is

$$ \frac{C}{R_T}=\frac{12.1922}{20}=0.6096, $$

or approximately \(60.96\%\). As a check, \(LR/C\approx(7.8078)(7.8078)/12.1922=5.000\) nM, matching \(K_d\), and each total equals its corresponding free concentration plus \(C\).

When the Free-Ligand Approximation Is Useful

If depletion is negligible, then \(C\) is small relative to \(L_T\), so \(L=L_T-C\approx L_T\). Substituting \(L_T\) for free ligand in the equilibrium occupancy relation gives

$$ \text{occupancy}\approx\frac{L_T}{K_d+L_T}, \qquad C\approx R_T\frac{L_T}{K_d+L_T}. $$

This approximation is justified by a small bound fraction of the total ligand, not merely by \(L_T\) being larger than \(K_d\). A ligand concentration can be well above \(K_d\) and still be depleted if the amount of receptor is large relative to available ligand. Conversely, when receptor is small compared with ligand, depletion may be modest. There is no universal cutoff for acceptable error: compare approximate and exact results, or judge the error against the precision required for the model or decision.

The exact finite-target balance is usually the most direct option when both totals are known and the assumptions hold. The free-ligand approximation is simpler and can be adequate when its small-depletion condition is checked. A numerical solution of \(K_d=(L_T-C)(R_T-C)/C\) can also be used, but it solves the same equilibrium balance and should still be checked against \(0\leq C\leq\min(L_T,R_T)\). If equilibrium is not reached, an equilibrium calculation is not a substitute for a time-dependent binding model.

Diagnostics, Failure Modes, and Limits

After calculating a solution, verify both conservation equations and the concentration bounds. For nonnegative totals and positive \(K_d\), a valid solution must have \(0\leq C\leq\min(L_T,R_T)\), \(0\leq L\leq L_T\), and \(0\leq R\leq R_T\). It must also satisfy \(LR/C=K_d\) when \(C>0\). Failure of any check usually signals an algebra, unit, or root-selection error. Small rounding differences are expected, but substantial mass-balance discrepancies are not.

A particularly useful warning sign is an approximate bound concentration greater than \(L_T\). That violates conservation and directly demonstrates that treating total ligand as free is unsuitable. Less dramatic depletion can still matter: compare the exact and approximate values for the quantities used downstream, rather than judging only whether the approximate occupancy looks plausible.

  • Using total ligand in place of free ligand without checking depletion. The occupancy formula requires free ligand; correct the analysis by solving the finite-target balance when bound ligand is appreciable.
  • Choosing the larger quadratic root. It can imply a negative free concentration; select the root within the physical concentration bounds.
  • Mixing units. For example, combining ligand in nM with \(K_d\) in \(\mu\)M changes the numerical result; convert all concentrations to a common unit first.
  • Applying equilibrium to a transient situation. If the system has not approached equilibrium or totals change during binding, the static balance does not predict the time course.

The model also omits binding to other molecules, multiple receptor classes, spatial gradients, and processes that change ligand or receptor amounts. Such processes can make the simple totals incomplete. The finite-target balance corrects one specific issue—depletion by a single 1:1 receptor interaction—not every possible difference between nominal dose and biologically available ligand.

Key takeaway: When ligand depletion may matter, use the total ligand and receptor mass balances together with the one-site equilibrium relation. Recover free ligand from \(L=L_T-C\), and verify the result through conservation, physical bounds, and the equilibrium definition.

Check Your Understanding

Work each problem using the finite-target balance. Short solutions follow the questions.

  1. A system has \(L_T=12\) nM, \(R_T=4\) nM, and \(K_d=8\) nM. Write the quadratic equation for \(C\), and state the allowable range for its physical root.
  2. For \(L_T=2\) nM, \(R_T=20\) nM, and \(K_d=2\) nM, explain why the free-ligand approximation is suspect before calculating an exact value.
  3. If a calculation returns \(C=15\) nM when \(L_T=10\) nM and \(R_T=30\) nM, identify the diagnostic failure.
  4. State the condition that makes \(L\approx L_T\) reasonable, and explain why \(L_T>K_d\) alone is not sufficient.

Solutions. For question 1, \(C^2-(12+4+8)C+(12)(4)=0\), or \(C^2-24C+48=0\). The physical root must lie between 0 and 4 nM. For question 2, at most 2 nM ligand is available while receptor totals 20 nM; binding can consume a large fraction of ligand, so setting free ligand equal to 2 nM requires checking. The exact root is \(L=[2-20-2+\sqrt{(-20)^2+4(2)(2)}]/2=(-20+\sqrt{416})/2\approx0.1980\) nM, showing strong depletion. For question 3, \(C\) exceeds \(L_T\), violating ligand conservation and the required bound \(C\leq\min(L_T,R_T)\). For question 4, the approximation requires \(C/L_T\) to be small enough for the intended use; \(L_T>K_d\) says nothing by itself about how much ligand the finite receptor pool binds.

Clinical Biostats