From Occupancy to a Functional Readout
Receptor occupancy describes binding; it does not, by itself, specify the size of a functional response. A measured fraction of receptors bound might correspond to a proportional change in a proximal signal, a response that levels off, or little detectable change until occupancy passes a threshold. A model connecting these quantities makes that extra assumption explicit.
Here, the input is the occupancy supplied by a binding model, and the output is a proximal functional readout, such as a signaling measurement. The earlier material on receptor occupancy and agonist and antagonist efficacy distinguishes binding from efficacy; this tutorial focuses on the mathematical mapping between an occupancy value and a readout. The mappings below are alternatives, not successive steps that every receptor must follow.
A static mapping assumes that the readout can be related to the occupancy at the same time point without explicitly modelling a delay or additional time-varying states. It also assumes that the reported occupancy is a suitable summary for the receptors contributing to the measured readout. For the simple forms below, assume one effective receptor population, a consistent readout scale, and occupancy measured or estimated on the same basis across observations. If receptor populations, cell types, or measurement conditions differ, a single mapping may not describe them all.
Three Candidate Mappings
A linear mapping assumes a constant change in readout per unit change in occupancy:
Here \(Y_0\) is the predicted readout at zero occupancy, in readout units, and \(a\) is the change in readout for a full unit increase in fractional occupancy, also in readout units. A positive \(a\) gives an increasing mapping; a negative \(a\) gives a decreasing one. Linearity is a useful local approximation, but it should not be extrapolated outside the occupancy range for which it is plausible.
A saturable mapping allows the incremental response to diminish as occupancy increases:
In this form, \(A\) is the asymptotic response increment of the mathematical function, in readout units, and \(K_\theta\) is a positive occupancy-scale parameter, dimensionless when occupancy is a fraction. At \(\theta=K_\theta\), the increment is \(A/2\). Because physical occupancy is restricted to \(0\leq\theta\leq1\), the asymptotic increment \(A\) need not be reached: if \(K_\theta\) is large compared with 1, the function may not visibly plateau in the feasible occupancy range.
A thresholded linear mapping allows a baseline region with no occupancy-related increment, followed by a response that increases above a threshold:
The threshold \(\theta_c\) is between 0 and 1, and \(g\) is the increase in readout per unit occupancy above that threshold. This mapping is continuous at \(\theta_c\), but its slope changes there. A flat predicted region means no change in the modelled readout below threshold; it does not prove that receptors below that occupancy are biologically inactive. In particular, a threshold imposed by an assay’s detection limit is not necessarily a biological threshold.
Worked Example: Comparing All Three Mappings
Suppose the input occupancy is \(\theta=0.20\) or \(\theta=0.60\), and the proximal readout is reported in hypothetical signal units. Compare three specified models: a linear mapping with \(Y_0=2\) and \(a=10\); a saturable mapping with \(Y_0=2\), \(A=12\), and \(K_\theta=0.30\); and a thresholded linear mapping with \(Y_0=2\), \(g=15\), and \(\theta_c=0.25\). These values are illustrative, not observations from a real study.
At \(\theta=0.20\), the linear prediction is \(Y=2+10(0.20)=4.00\) signal units. The saturable prediction is \(Y=2+12[0.20/(0.30+0.20)]=2+12(0.40)=6.80\) signal units. Since \(0.20<0.25\), the thresholded model gives \(Y=2+15\max(0,0.20-0.25)=2.00\) signal units.
At \(\theta=0.60\), the linear prediction is \(Y=2+10(0.60)=8.00\) signal units. The saturable prediction is \(Y=2+12[0.60/(0.30+0.60)]=2+12(2/3)=10.0\) signal units. The thresholded prediction is \(Y=2+15(0.60-0.25)=2+5.25=7.25\) signal units.
The models make different predictions at the same occupancy because their shapes and parameters differ. These calculations compare specified mappings; they do not establish which shape best describes a biological system.
Why a Saturating Shape Can Arise
One simple justification for a saturating mapping is a downstream component that is activated by occupied receptors and has a finite total amount. Let \(q\) be the amount of that component in an activated state and \(Q_{\mathrm{tot}}\) its total amount. Suppose activation occurs at a rate proportional to occupancy and the amount still inactive, while deactivation occurs at a rate proportional to the activated amount:
For this illustrative derivation, assume occupancy is held constant long enough for the downstream component to reach equilibrium, the rate constants are positive, and this one state is sufficient to describe the proximal readout. At equilibrium, set \(dq/dt=0\) and solve:
If \(Y=Y_0+Aq/Q_{\mathrm{tot}}\), this gives the saturable form with \(K_\theta=k_{\mathrm{off}}/k_{\mathrm{on}}\). This derivation shows how a finite downstream pool can generate a saturating relationship; it does not show that every proximal response uses this mechanism. The static formula also relies on equilibration. If occupancy changes faster than the downstream state can respond, this equilibrium mapping may misrepresent the time course.
The linear mapping can approximate the saturable form at low occupancy relative to \(K_\theta\). When \(\theta/K_\theta\) is small,
Thus, data confined to a narrow low-occupancy range may look approximately linear even when a saturating mapping is plausible. That similarity makes it difficult to distinguish the forms without informative observations over a wider occupancy range.
Worked Example: Applying the Low-Occupancy Approximation
Use the saturable parameters \(Y_0=2\) signal units, \(A=12\) signal units, and \(K_\theta=0.30\). At occupancy \(\theta=0.03\), the exact prediction is \(Y=2+12[0.03/(0.30+0.03)]\).
The exact fraction is \(0.03/0.33=0.09091\), so \(Y=2+12(0.09091)=3.091\) signal units. The low-occupancy approximation gives \(Y\approx2+(12/0.30)(0.03)=2+1.20=3.20\) signal units. The approximation is 0.1091 signal units higher than the exact prediction. The discrepancy is modest here because \(\theta/K_\theta=0.10\), but the approximation becomes less reliable as occupancy approaches or exceeds \(K_\theta\).
Using Data to Assess a Mapping
A response mapping is most informative when occupancy and the proximal readout are paired under comparable conditions. Before interpreting a fitted or specified mapping, check that the occupancy values cover a useful range and that readout units, baseline conditions, and timing are consistent. The mathematical assumptions are not confirmed merely because a formula can be evaluated.
Some practical diagnostics are:
- Check the range and endpoints. Occupancy should be between 0 and 1 when expressed as a fraction. Compare predicted and observed readouts near zero occupancy and at the highest measured occupancies.
- Look for systematic shape mismatch. Residuals that change direction across occupancy—for example, repeated underprediction at intermediate values and overprediction at high values—can indicate that a straight line is missing curvature. A plateau-like pattern or a delayed rise can likewise be inconsistent with a proposed shape.
- Check the meaning of the threshold. A flat observed region may reflect limited measurement sensitivity, sparse sampling, or a true threshold-like process. The mapping alone cannot distinguish these explanations.
- Check the time assumption. If the same occupancy is associated with different readouts at different times, a static mapping may be inadequate. A downstream delay or changing system state may matter, rather than a different occupancy-to-response shape alone.
- Check parameter support. If data never approach a plateau, an estimated saturation scale or asymptote may be weakly constrained. If measurements do not bracket a proposed threshold, its location is likewise difficult to determine from those data.
The linear form is simple and may be suitable over a restricted range. The saturable form is appropriate when diminishing increments are expected or supported by observations. The thresholded form is useful when a no-increment region is a defensible approximation. More flexible functions are possible, but added flexibility requires enough informative data to distinguish it from simpler alternatives. None of these mappings, without further evidence, identifies the underlying cellular mechanism.
Worked Example: Inferring a Linear Mapping from Two Points
In a hypothetical calibration set, occupancy \(0.10\) is paired with a readout of 5.0 signal units, and occupancy \(0.50\) with 13.0 signal units. Assume the linear form \(Y=Y_0+a\theta\) is suitable over this range, and treat both occupancy values as known accurately for this calculation.
The slope is the change in readout divided by the change in occupancy: \(a=(13.0-5.0)/(0.50-0.10)=8.0/0.40=20.0\) signal units per unit fractional occupancy. Substitute the first point to find the baseline: \(Y_0=5.0-(20.0)(0.10)=3.0\) signal units. The resulting mapping is \(Y=3.0+20.0\theta\).
At \(\theta=0.30\), the prediction is \(Y=3.0+(20.0)(0.30)=9.0\) signal units. The two points determine a straight line, but do not test whether linearity is adequate; additional observations are needed to assess the shape and the model’s assumptions.
Common Mistakes and Limitations
- Equating occupancy with response. A fraction bound is not a signal value. The mapping and its units must be stated before translating one into the other.
- Calling \(A\) the response at full occupancy. In the saturable formula, \(A\) is the limiting increment as \(\theta\) grows without bound. Since occupancy cannot exceed 1, the predicted increment at full occupancy is \(A/(K_\theta+1)\), not generally \(A\).
- Interpreting a threshold as proof of biological inactivity. A thresholded mapping may approximate a detection floor or limited data resolution. The functional relationship alone does not establish a molecular switch.
- Extrapolating beyond the data. A line can predict implausible values outside its supported range, while a fitted saturation or threshold can be poorly identified if observations do not cover the relevant part of the occupancy scale.
- Ignoring time and population differences. A single static mapping may fail when the readout lags occupancy or when the reported occupancy averages populations with different responses.
These mappings are compact ways to state assumptions about signal transmission from occupied receptors to a proximal readout. Their value is not that one must always be chosen, but that each makes a distinct, testable shape assumption. Interpretation should stay within the measured occupancy range and the biological and measurement conditions used to define the readout.
Check Your Understanding
Answer each question using the definitions and assumptions above. Short solutions follow.
- For \(Y=4+8\theta\), what is the predicted readout at occupancy \(0.25\), and what are the units of the slope if \(Y\) is measured in signal units?
- For \(Y=1+10\theta/(0.20+\theta)\), what is the response increment at \(\theta=0.20\)?
- For a thresholded mapping with \(\theta_c=0.40\), \(Y_0=3\), and \(g=5\), what are the predicted readouts at occupancies 0.30 and 0.70?
- Why can observations only at low occupancy make linear and saturable mappings hard to distinguish?
Solutions. 1. \(Y=4+8(0.25)=6\) signal units; the slope is 8 signal units per unit fractional occupancy. 2. The increment is \(10[0.20/(0.20+0.20)]=10(0.5)=5\) readout units, so the total predicted readout is 6. 3. At 0.30, occupancy is below threshold, so \(Y=3\). At 0.70, \(Y=3+5(0.70-0.40)=4.5\). 4. At low occupancy relative to \(K_\theta\), the saturable form is approximately linear, \(Y\approx Y_0+(A/K_\theta)\theta\), so the observed range may not reveal curvature.