Tutorials › Quantitative Systems Pharmacology (QSP) › Receptor Subtypes and Mixed-Affinity Occupancy

Receptor Binding and Target Engagement · Tutorial 197 of 799

Receptor Subtypes and Mixed-Affinity Occupancy

Learn to calculate total equilibrium occupancy for receptor subtypes and understand why a single apparent affinity can depend on the concentration range used.

Intermediate 11 min read

What You'll Learn

  • Calculate fractional occupancy separately for each receptor subtype.
  • Combine subtype occupancies using receptor-abundance weights.
  • Distinguish total occupied receptors from fractional total occupancy.
  • Derive the low- and high-concentration limits of a mixed-affinity population.
  • Explain why one apparent dissociation constant may misrepresent a mixture.
  • Check assumptions and diagnose when a pooled occupancy estimate is misleading.

Why Receptor Subtypes Change the Occupancy Calculation

A ligand can encounter more than one receptor subtype in a tissue. If those subtypes have different affinities, they will not generally have the same fractional occupancy at a given free-ligand concentration. The total number of occupied receptors therefore depends on both the occupancy of each subtype and how abundant that subtype is.

The introductory material on fractional occupancy establishes the one-site equilibrium equation. Here, apply that equation separately to each subtype, then combine the results using subtype abundance. This distinction matters: an abundant, lower-affinity subtype can contribute more occupied receptors than a scarce, higher-affinity subtype, even when its own fractional occupancy is lower.

Definition: For receptor subtype \(i\), let \(R_i\) be its total receptor abundance, \(K_{d,i}\) its equilibrium dissociation constant, and \(C\) the shared free-ligand concentration. Under one-site equilibrium binding, the subtype's fractional occupancy is \(\theta_i=C/(K_{d,i}+C)\). The fractional occupancy of the full receptor population is the abundance-weighted mean \(\theta_{\mathrm{tot}}=\sum_i w_i\theta_i\), where \(w_i=R_i/\sum_j R_j\). The total amount of ligand-bound receptor is \(B_{\mathrm{tot}}=\sum_i R_i\theta_i\).

A dissociation constant has concentration units; a smaller \(K_{d,i}\) indicates higher affinity under this equilibrium model. Each weight \(w_i\) is nonnegative and the weights sum to one. Consequently, \(\theta_{\mathrm{tot}}\) is between zero and one. It is not the unweighted mean of the subtype occupancies unless all subtypes have equal abundance.

The calculation assumes that each subtype has a single class of independent binding sites, that the ligand is at equilibrium, and that the free concentration \(C\) is known and common to the subtypes. It also assumes that the specified \(R_i\) values represent the receptor population being summarized. If binding substantially depletes ligand, total ligand concentration cannot simply be substituted for \(C\); the finite-target binding balance must first be used to determine free concentration. The formulas here describe equilibrium occupancy, not a time-dependent approach to equilibrium.

Calculating Occupancy Across Subtypes

For each subtype, calculate \(\theta_i\) from its own \(K_{d,i}\), then multiply by \(R_i\) to obtain its bound receptor amount. Add those amounts to get \(B_{\mathrm{tot}}\). To report a fraction for the combined population, divide by the total receptor abundance. This order of operations keeps the abundance weighting explicit and helps avoid confusing receptor amount with receptor fraction.

Worked Example: A Scarce High-Affinity Subtype

Suppose a hypothetical cell population has 20 pmol of subtype A and 80 pmol of subtype B. Their dissociation constants are \(K_{d,A}=1\) nM and \(K_{d,B}=100\) nM. At a free-ligand concentration of \(C=1\) nM, calculate each occupancy and the total fractional occupancy.

For subtype A, \(\theta_A=C/(K_{d,A}+C)=1/(1+1)=0.5\). Its bound amount is \(20\text{ pmol}\times0.5=10\) pmol.

For subtype B, \(\theta_B=1/(100+1)=1/101=0.009901\). Its bound amount is \(80\text{ pmol}\times1/101=0.7921\) pmol.

Thus \(B_{\mathrm{tot}}=10+0.7921=10.7921\) pmol, and total receptor abundance is \(20+80=100\) pmol. Therefore, \(\theta_{\mathrm{tot}}=10.7921/100=0.1079\), or about 10.79%. Although half of subtype A is occupied, the whole population is only about 10.79% occupied because subtype B is much more abundant and has low occupancy at this concentration.

The same calculation can be written directly in terms of abundance weights:

$$ \theta_{\mathrm{tot}}(C)=\sum_i \frac{R_i}{\sum_j R_j}\frac{C}{K_{d,i}+C}. $$

This is a sum of one-site occupancy curves, not generally another one-site occupancy curve. The distinction becomes important when a mixed population is summarized by fitting or reporting a single apparent dissociation constant.

Why One Apparent Affinity Can Mislead

For a single receptor population, fractional occupancy follows a hyperbola: \(\theta=C/(K_d+C)\). A mixture also rises with concentration and approaches full occupancy under the assumptions above, but its curve is a weighted sum of hyperbolas with different \(K_d\) values. Its shape can therefore differ from that of any single one-site population.

One way to see the problem is to ask what single \(K_d\) would reproduce the mixture at one chosen concentration. Rearranging the one-site equation gives \(K_{d,\mathrm{app}}=C(1-\theta_{\mathrm{tot}})/\theta_{\mathrm{tot}}\). This value is a concentration-specific summary: it need not reproduce the mixture at other concentrations.

Worked Example: Apparent Affinity Depends on Concentration

Use the hypothetical mixture from the previous example, with 20% subtype A (\(K_{d,A}=1\) nM) and 80% subtype B (\(K_{d,B}=100\) nM). Calculate its total occupancy at \(C=10\) nM and compare the single-site apparent \(K_d\) at 1 nM and 10 nM.

At \(C=10\) nM, subtype A occupancy is \(10/(1+10)=10/11=0.9091\). Subtype B occupancy is \(10/(100+10)=10/110=0.09091\). The mixture occupancy is \(0.2(0.9091)+0.8(0.09091)=0.2545\), or about 25.45%.

At \(C=1\) nM, the previous calculation gave \(\theta_{\mathrm{tot}}=0.1079\). The single-site value that matches that point is \(K_{d,\mathrm{app}}=1(1-0.1079)/0.1079=8.27\) nM.

At \(C=10\) nM, the value that matches the mixture is \(K_{d,\mathrm{app}}=10(1-0.2545)/0.2545=29.29\) nM. These two apparent values differ substantially. Neither is the unique affinity of the mixture; each merely matches the mixture at its own concentration.

The dependence can also be understood from the low- and high-concentration limits. For small positive \(C\), each subtype has \(\theta_i\) approximately \(C/K_{d,i}\), so the initial slope of the mixed curve depends on the abundance-weighted reciprocals of the dissociation constants. At high concentration, the unoccupied fraction of subtype \(i\) is \(K_{d,i}/(K_{d,i}+C)\), approximately \(K_{d,i}/C\). These limits imply two different candidate summaries for a single-site fit.

$$ \theta_{\mathrm{tot}}(C)\approx C\sum_i\frac{w_i}{K_{d,i}} \quad (C\ \text{small}),\qquad 1-\theta_{\mathrm{tot}}(C)\approx\frac{\sum_i w_iK_{d,i}}{C} \quad (C\ \text{large}). $$

Matching the low-concentration slope with a single hyperbola would require \(K_{d,\mathrm{low}}=1/\sum_i(w_i/K_{d,i})\), the weighted harmonic mean. Matching the high-concentration unoccupied fraction would require \(K_{d,\mathrm{high}}=\sum_iw_iK_{d,i}\), the weighted arithmetic mean. For positive weights, the weighted arithmetic mean is at least as large as the weighted harmonic mean, with equality only when all positively weighted \(K_{d,i}\) values are equal. Thus, unless all represented subtypes have the same affinity, one single \(K_d\) cannot in general match both ends of the curve.

Worked Example: Comparing Low- and High-Concentration Summaries

Consider a hypothetical population with 25% subtype A (\(K_{d,A}=2\) nM) and 75% subtype B (\(K_{d,B}=20\) nM). Calculate the single-site \(K_d\) suggested by the initial slope and the one suggested by the high-concentration limit.

The low-concentration summary is the weighted harmonic mean: \(K_{d,\mathrm{low}}=1/[0.25/2+0.75/20]=1/(0.125+0.0375)=1/0.1625=6.154\) nM.

The high-concentration summary is the weighted arithmetic mean: \(K_{d,\mathrm{high}}=0.25(2)+0.75(20)=0.5+15=15.50\) nM.

The summaries differ by more than a factor of two. A single-site approximation could still be useful over a specified concentration range, but its fitted value would depend on that range and on the fitting criterion. It should not be interpreted automatically as the dissociation constant of a physical receptor subtype.

Composition, Measurements, and Model Checks

Two tissues exposed to the same free-ligand concentration can have different total occupancy if their subtype proportions differ. Even if every subtype has identical binding properties across the tissues, changing the weights \(w_i\) changes the weighted sum. Conversely, a measured difference in total occupancy does not by itself establish a change in receptor affinity: it could reflect a change in subtype abundance.

The abundance-weighted formula describes the fraction of all receptors that are occupied. An experimental signal may instead weight subtypes differently if detection efficiency differs by subtype. If a signal is proportional to bound receptors with subtype-specific proportionality factors \(q_i\), its total is proportional to \(\sum_i q_iR_i\theta_i\), not necessarily \(\sum_i R_i\theta_i\). The measurement model must therefore be checked before treating a pooled signal as total receptor occupancy.

Useful checks include confirming that all \(R_i\) refer to the same receptor-amount units, that the \(K_{d,i}\) and \(C\) use the same concentration units, and that the weights sum to one. At \(C=0\), the model gives zero occupancy; as \(C\) increases without bound, it approaches one for every subtype with finite positive \(K_{d,i}\). A calculated total occupancy outside the range zero to one signals an arithmetic, weighting, or model-definition error.

Key result: Total equilibrium occupancy is the receptor-abundance-weighted average of subtype-specific occupancies. A one-site apparent \(K_d\) may summarize a mixture over a stated concentration range, but it is not generally a concentration-independent property of the mixed population.

Common Pitfalls and Alternatives

  • Using an unweighted mean. Averaging subtype occupancies equally gives each subtype equal influence, even when their receptor abundances differ. Use \(w_i=R_i/\sum_jR_j\) for the fraction of all receptors occupied.
  • Averaging the dissociation constants first. Substituting an average \(K_d\) into one hyperbola does not generally reproduce the weighted sum of subtype-specific curves. Calculate each \(\theta_i\) first, then combine.
  • Calling a pointwise fit the mixture's true affinity. A single-site \(K_{d,\mathrm{app}}\) inferred at one concentration may not match the curve elsewhere. State the concentration range and model used when reporting an apparent value.
  • Confusing total bound amount with fractional occupancy. \(B_{\mathrm{tot}}\) has receptor-amount units; \(\theta_{\mathrm{tot}}\) is dimensionless and requires division by total receptor abundance.

The most direct approach, when subtype abundances and affinities are available, is to model the subtypes explicitly. A single apparent affinity can be a practical approximation when data are limited, but it compresses information and can hide subtype contributions. Which description is useful depends on the question: predicting total occupied receptor may require the mixture, while a narrowly targeted empirical summary may need only a fit over the observed concentration range. Neither representation resolves uncertainty in the subtype abundances or binding parameters; those inputs should be supported by evidence or treated as uncertain in a broader model.

Check Your Understanding

Answer each question before reading its solution.

  1. A population contains 60 pmol of subtype A and 40 pmol of subtype B. At a stated ligand concentration, their occupancies are 0.20 and 0.80. What is total fractional occupancy?
  2. Why is the unweighted mean of subtype occupancies generally not the fraction of all receptors occupied?
  3. A mixture has positively weighted subtypes with different finite \(K_d\) values. Can one one-site \(K_d\) match both the low- and high-concentration limits? Explain.
  4. For 50% subtype A with \(K_{d,A}=4\) nM and 50% subtype B with \(K_{d,B}=12\) nM, calculate total occupancy at \(C=4\) nM.

Solutions. (1) The weights are 0.60 and 0.40, so \(\theta_{\mathrm{tot}}=0.60(0.20)+0.40(0.80)=0.12+0.32=0.44\), or 44%. (2) The unweighted mean assigns equal contribution to each subtype, whereas the fraction of all receptors occupied weights each subtype by its share of total receptor abundance. (3) No, except when all subtypes with positive abundance have the same \(K_d\): the low-concentration limit requires the weighted harmonic mean, and the high-concentration limit requires the weighted arithmetic mean; these are unequal for differing positive values. (4) Subtype A occupancy is \(4/(4+4)=0.5\); subtype B occupancy is \(4/(12+4)=0.25\). Equal abundance gives \(\theta_{\mathrm{tot}}=0.50(0.50)+0.50(0.25)=0.375\), or 37.5%.

Clinical Biostats