Tutorials › Quantitative Systems Pharmacology (QSP) › Receptor Reserve and Occupancy-Response Curves

Receptor Binding and Target Engagement · Tutorial 204 of 799

Receptor Reserve and Occupancy-Response Curves

Learn how operational efficacy links receptor occupancy to response, allowing strong effects below full occupancy and apparent potency to change without altered affinity.

Intermediate 11 min read

What You'll Learn

  • Define receptor reserve using an operational occupancy-response model.
  • Relate free agonist concentration, equilibrium occupancy, and functional response.
  • Derive the agonist’s attainable response maximum and apparent EC50.
  • Compare systems with different operational efficacy but the same binding affinity.
  • Predict how receptor loss changes apparent potency and maximal response.
  • Recognize model assumptions, diagnostic warning signs, and interpretation pitfalls.

Why Occupancy and Response Can Diverge

A receptor-occupancy curve describes the fraction of receptors bound by a ligand. A functional response curve describes what the cell does. These curves need not rise together: a pathway may amplify the signal from occupied receptors, so that a substantial response occurs while many receptors remain unoccupied. This capacity for a response close to the attainable maximum at submaximal occupancy is called receptor reserve.

The earlier material on mapping occupancy to a proximal functional response introduced the distinction between binding and downstream response. Here, an operational model makes that distinction quantitative. It shows how receptor reserve affects apparent potency and why a functional concentration-response curve alone does not establish binding affinity.

Definition: Receptor reserve is the capacity of a receptor-effector system to produce a response near its attainable maximum without occupying all available receptors. In the operational model below, reserve is associated with a large transduction parameter \(\tau\): response rises steeply relative to occupancy when \(\tau\) is large. Reserve is a property of the ligand-receptor-effector system, not a claim that a particular fixed number of receptors is always “spare.”

The Operational Occupancy-Response Model

Assume a reversible, one-site interaction at equilibrium, with free agonist concentration \(C\), dissociation constant \(K_D\), and fractional occupancy \(\theta\). The equilibrium occupancy is

$$ \theta=\frac{C}{K_D+C}. $$

For the response relationship, let \(E_{\mathrm{sys,max}}\) denote the system’s maximum response scale, \(\tau\) the dimensionless operational transduction parameter, and \(n\) the transducer slope parameter. Use a baseline-subtracted response \(E\). The operational model is

$$ \frac{E}{E_{\mathrm{sys,max}}} = \frac{(\tau\theta)^n}{1+(\tau\theta)^n}. $$

The model assumes that occupancy drives a monotone transduction process and that the response is measured under conditions where occupancy and response are adequately represented by their steady-state values. It also assumes a homogeneous receptor population for this single-site representation, controlled free agonist concentration, and no material change in receptor availability during measurement. The model is an approximation: feedback, desensitization, multiple receptor populations, or time-varying exposure may require a richer description.

The parameter \(\tau\) summarizes operational efficacy in the system. It can reflect agonist efficacy together with receptor abundance and coupling efficiency; it should not be interpreted as an affinity constant. The parameter \(n\) controls the steepness of transduction from occupancy to response. It is not, by itself, evidence of cooperative ligand binding.

Key result: At full occupancy, \(\theta=1\), the model predicts \(E/E_{\mathrm{sys,max}}=\tau^n/(1+\tau^n)\). For finite \(\tau\), this is below one. If \(\tau\) is large, however, a response close to this agonist’s attainable maximum can occur at an occupancy well below one.

Worked Example: A Near-Maximal Response at Partial Occupancy

Suppose a hypothetical agonist has \(K_D=10\) nM, \(\tau=9\), and \(n=1\), in a system with \(E_{\mathrm{sys,max}}=100\) response units. Find the occupancy and response at \(C=9\) nM, then compare that response with the agonist’s attainable maximum.

Occupancy. Using the equilibrium occupancy equation,

$$ \theta=\frac{9}{10+9}=\frac{9}{19}=0.4737. $$

Thus about 47.37% of receptors are occupied. With \(n=1\), the response is

$$ E=100\frac{9(0.4737)}{1+9(0.4737)} =100\frac{4.263}{5.263} =81.00\ \text{response units}. $$

At saturating agonist concentration, occupancy approaches one and the model’s response approaches \(100(9)/(1+9)=90.00\) response units. Therefore, the response at 47.37% occupancy is \(81/90=0.9000\), or 90.00% of this agonist’s attainable maximum. This illustrates reserve: a near-maximal response does not require full occupancy.

Deriving the Functional Maximum and Apparent Potency

The occupancy equation approaches one as \(C\) becomes very large. Substituting \(\theta=1\) in the response model gives the agonist’s attainable response maximum:

$$ E_{\max,\mathrm{agonist}} = E_{\mathrm{sys,max}}\frac{\tau^n}{1+\tau^n}. $$

This maximum equals the system scale only in the limiting case of very large \(\tau\). A finite-\(\tau\) agonist may therefore have a lower response plateau, even when its receptors are saturated. A large \(\tau\) allows the response to approach the system scale at less than full occupancy.

Define functional \(EC_{50}\) as the agonist concentration that produces half of its own attainable maximum, not half of \(E_{\mathrm{sys,max}}\). To derive it, let \(x=(\tau\theta)^n\). At the functional \(EC_{50}\), the response fraction of the system scale is half its saturating value:

$$ \frac{x}{1+x} = \frac{1}{2}\frac{\tau^n}{1+\tau^n}. $$

Multiplying both sides by \(2(1+x)(1+\tau^n)\) and collecting terms gives \(x(2+\tau^n)=\tau^n\). Thus \(\tau\theta=[\tau^n/(2+\tau^n)]^{1/n}\), and \(\theta=1/(2+\tau^n)^{1/n}\). Since \(C=K_D\theta/(1-\theta)\), substitution yields

$$ EC_{50} = \frac{K_D}{(2+\tau^n)^{1/n}-1}. $$

For \(n=1\), this reduces to \(EC_{50}=K_D/(1+\tau)\). Accordingly, increasing \(\tau\) lowers the functional \(EC_{50}\) even if \(K_D\) does not change. This is an increase in apparent functional potency, not an increase in binding affinity.

Worked Example: Same Affinity, Different Apparent Potency

Consider two hypothetical systems with the same agonist and \(K_D=10\) nM, both with \(n=1\) and \(E_{\mathrm{sys,max}}=100\) response units. System A has \(\tau=1\), and System B has \(\tau=9\). Calculate each attainable maximum and functional \(EC_{50}\).

For System A,

$$ E_{\max,\mathrm{agonist}}=100\frac{1}{1+1}=50.00\ \text{response units}, \qquad EC_{50}=\frac{10}{1+1}=5.000\ \text{nM}. $$

For System B,

$$ E_{\max,\mathrm{agonist}}=100\frac{9}{1+9}=90.00\ \text{response units}, \qquad EC_{50}=\frac{10}{1+9}=1.000\ \text{nM}. $$

The functional \(EC_{50}\) is fivefold lower in System B, although the binding affinity parameter is \(10\) nM in both systems. The difference follows from the change in \(\tau\), which represents a change in operational efficacy, receptor abundance, coupling, or some combination—not a change in \(K_D\).

How Receptor Availability Changes the Curve

If receptor number falls while agonist affinity and downstream conditions remain unchanged, operational transduction capacity can fall as well. In this simple model, represent that change as a reduction in \(\tau\). The consequences are a higher functional \(EC_{50}\) and a lower attainable response maximum. The direction and size of those changes depend on the model and conditions; they are not a universal numerical rule for every biological system.

Worked Example: Reduced Receptor Availability

Suppose a hypothetical system has \(K_D=10\) nM, \(n=1\), and \(E_{\mathrm{sys,max}}=100\) response units. Before receptor loss, \(\tau=9\). After receptor availability is reduced, suppose \(\tau=1.8\). Calculate the new maximum and functional \(EC_{50}\), and compare them with the original values.

Before the change, the values are \(90.00\) response units and \(10/(1+9)=1.000\) nM. After the change,

$$ E_{\max,\mathrm{agonist}}=100\frac{1.8}{1+1.8} =64.29\ \text{response units}, \qquad EC_{50}=\frac{10}{1+1.8} =3.571\ \text{nM}. $$

The model therefore predicts a lower response plateau and a higher functional \(EC_{50}\), while \(K_D\) remains \(10\) nM. This pattern is consistent with reduced operational reserve under the stated assumptions. It does not, by itself, prove receptor loss: other changes to coupling or pathway state could also alter \(\tau\).

Diagnostics, Alternatives, and Limitations

When using this model to interpret data, check that the measurements support its assumptions. A concentration-response curve should be measured over a range broad enough to identify its plateau; otherwise, the maximum and \(EC_{50}\) may be poorly determined. If occupancy is available, compare measured responses at matched occupancies. Systematic differences from the model’s monotone occupancy-response relation may indicate time dependence, feedback, receptor heterogeneity, or another missing process. Also check that free agonist concentration is being represented appropriately; if ligand depletion is material, nominal concentration need not equal free concentration.

The operational model is not the only way to summarize these data. An empirical Hill concentration-response curve can describe a response plateau, \(EC_{50}\), and slope without explicitly representing occupancy or transduction. It may be useful for summarizing measured responses, but its fitted \(EC_{50}\) does not on its own identify \(K_D\) or receptor reserve. A model that assumes response is directly proportional to occupancy is simpler, but cannot represent a strongly amplified response at low occupancy followed by a saturating functional response. Choose the model that matches the question and the available measurements; added mechanistic parameters are not automatically identifiable from a single response curve.

The operational model also has important limits. \(\tau\) and \(K_D\) can be difficult to estimate separately from response data alone. Changes in receptor availability, coupling, ligand efficacy, or assay conditions may produce similar functional curve changes. If receptor number or occupancy is measured independently, those data can help distinguish explanations. For time-dependent systems, equilibrium occupancy and a steady-state response may not describe the measurements; a dynamic model may be needed instead.

  • Do not equate occupancy with response. The model includes a transduction step between them.
  • Do not call a lower \(EC_{50}\) stronger binding. Functional potency depends on both binding and operational transduction.
  • Do not confuse the two maxima. \(E_{\mathrm{sys,max}}\) is a response scale; \(E_{\max,\mathrm{agonist}}\) is the attainable plateau for the agonist in the modeled system.
  • Do not infer reserve from a curve shift alone. Changes to receptor availability, coupling, or other system properties may be confounded.
Key takeaway: In the operational model, \(\tau\) links occupancy to response. A high \(\tau\) can produce a near-maximal response at submaximal occupancy and a lower functional \(EC_{50}\), without changing binding affinity \(K_D\).

Check Your Understanding

Answer each question, then compare with the solutions.

  1. With \(n=1\), \(K_D=12\) nM, and \(\tau=3\), what is the model-predicted functional \(EC_{50}\)?
  2. For \(n=1\) and \(\tau=4\), what fraction of \(E_{\mathrm{sys,max}}\) is the agonist’s attainable response maximum?
  3. Two systems have the same \(K_D\), but one has a lower functional \(EC_{50}\). Does this alone show that its binding affinity is higher?
  4. Name two biological or experimental features that could make the equilibrium operational model inadequate.

Solutions. (1) \(EC_{50}=K_D/(1+\tau)=12/(1+3)=3.000\) nM. (2) The fraction is \(\tau/(1+\tau)=4/5=0.8000\), or 80.00% of \(E_{\mathrm{sys,max}}\). (3) No. A change in \(\tau\) can lower functional \(EC_{50}\) without changing \(K_D\). (4) Examples include material ligand depletion, time-varying occupancy or response, desensitization, feedback, or a heterogeneous receptor population.

Clinical Biostats