Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Separating Target Occupancy From Functional Response

Target Binding, Occupancy and Engagement · Tutorial 235 of 598

Separating Target Occupancy From Functional Response

Learn to compare occupancy and functional response on matched scales and times, and to distinguish response amplification from timing and target-state explanations for a mismatch.

Advanced 10 min read

What You'll Learn

  • Define fractional target occupancy and state the assumptions behind the simple equilibrium binding relation.
  • Distinguish proportional occupancy-response models from nonlinear transduction models.
  • Derive and apply a first-order model for delayed functional response after occupancy changes.
  • Diagnose apparent occupancy-response mismatches caused by different target states.
  • Match assay timing, target context, and response scale before interpreting discordance.
  • Recognize when an occupancy-response comparison cannot identify the underlying mechanism.

Why Occupancy and Response Can Disagree

A target-engagement assay asks how much target is bound by a drug. A functional assay asks what biological effect is observed. These measurements are related, but they are not interchangeable: substantial occupancy need not produce a large response, and a large response need not require occupancy of every target molecule. Interpreting a difference between the two requires attention to how binding is translated into function, when each measurement is made, and which target states each assay samples.

The first requirement is to compare like with like. Occupancy and response should, as far as possible, refer to the same biological system, target population, exposure condition, and time point. If one measurement is made in a cell preparation shortly after dosing and the other in a different preparation or much later, their difference may reflect context or timing rather than a failure of the binding model.

Definition: Fractional occupancy, denoted by \(\theta\), is the fraction of a specified target population bound by the ligand at a stated time. Functional response, denoted by \(R\), is a measured downstream outcome, such as enzyme activity or a signaling readout. Occupancy describes binding; response describes function. A model linking them must specify the target population, assay conditions, and time.

Equilibrium Occupancy and Functional Transduction

For a single homogeneous target population at equilibrium with one reversible ligand, the free ligand concentration \(C\) and dissociation constant \(K_D\) give the simple occupancy relation below. This relation assumes that \(C\) is the free concentration at the target site, binding is at equilibrium, the target population is homogeneous, and there are no competing ligands or cooperative effects. When those assumptions do not hold, the simple relation may not describe the measured occupancy.

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

Even when this equation correctly describes binding, it does not determine the functional response. A direct proportional model assumes that response above baseline is proportional to occupancy. A nonlinear transduction model instead allows downstream processes to amplify or attenuate the effect of bound target. For example, define normalized response \(e=(R-R_0)/E_{\max}\), where \(R_0\) is baseline response and \(E_{\max}\) is the response at complete occupancy under the model. The following phenomenological relation is one possible saturating link:

$$ e=\frac{(1+\kappa)\theta}{\kappa+\theta}, \qquad \kappa>0 $$

This relation gives \(e=0\) at zero occupancy and \(e=1\) at complete occupancy. When \(\kappa\) is small, response can be large even at partial occupancy. That pattern is often described as response amplification or receptor reserve. It is a property of the occupancy-to-response relation, not evidence that the occupancy measurement itself is wrong. The equation is a model, not a universal law; a different target and endpoint may require a different link.

Worked Example: High Response at Partial Occupancy

Hypothetical data. At a matched equilibrium time, a calibrated assay estimates \(\theta=0.40\). Baseline response is 10 units, and the response scale has \(E_{\max}=100\) units above baseline. Suppose the transduction parameter is \(\kappa=0.10\). Calculate the predicted response and compare it with a proportional model.

Under the stated saturating model, the normalized response is

$$ e=\frac{(1+\kappa)\theta}{\kappa+\theta} =\frac{(1+0.10)(0.40)}{0.10+0.40} =\frac{0.44}{0.50}=0.88 $$

The predicted response above baseline is \(0.88(100)=88\) units, so the total response is \(10+88=98\) units. Under a proportional model, the response above baseline would instead be \(0.40(100)=40\) units, for a total of 50 units. Thus, the same 40% occupancy corresponds to either prediction depending on the assumed functional link. These predictions are not competing occupancy estimates; they are different models of transduction.

The calculation assumes the measured occupancy and functional response refer to the same target population and condition, and that the specified transduction relation applies. A single point cannot establish that relation. Measurements across several occupancy levels are needed to assess whether the observed pattern is consistent with it.

Timing: Response May Lag Behind Occupancy

A common comparison error is to treat an occupancy measurement and a functional response as if both changed immediately. Binding can change faster than downstream signaling, enzyme regulation, or a measured biomarker. To describe a simple delay, suppose that occupancy changes quickly to a fixed level and remains there, and that this level implies a constant equilibrium response \(R_{\mathrm{eq}}\). Assume the observed response approaches \(R_{\mathrm{eq}}\) according to a first-order process with time constant \(\tau>0\):

$$ \frac{dR}{dt}=\frac{R_{\mathrm{eq}}-R}{\tau} $$

This model assumes a constant equilibrium response during the interval, a single time constant, and no additional delay or feedback process. It is a useful diagnostic model, not a complete description of every biological response. The solution is

$$ R(t)=R_{\mathrm{eq}}+\bigl(R(0)-R_{\mathrm{eq}}\bigr)e^{-t/\tau} $$

To justify the result without dividing by a quantity that could be zero, let \(D(t)=R(t)-R_{\mathrm{eq}}\). Since \(R_{\mathrm{eq}}\) is constant, the differential equation gives \(D'(t)=-D(t)/\tau\). Multiplying by the integrating factor \(e^{t/\tau}\) yields

$$ \frac{d}{dt}\left(e^{t/\tau}D(t)\right) =e^{t/\tau}\left(D'(t)+\frac{D(t)}{\tau}\right)=0 $$

Therefore \(e^{t/\tau}D(t)\) is constant and equals \(D(0)\), which gives the stated solution. If \(D(0)=0\), then \(R(0)=R_{\mathrm{eq}}\) and the solution is \(R(t)=R_{\mathrm{eq}}\) for all \(t\); the derivation remains valid in this case. If \(D(0)\ne 0\), the same expression shows the deviation from equilibrium shrinks exponentially, with the magnitude reduced to \(e^{-1}\) of its initial value after one time constant.

Worked Example: Occupancy Changes Before the Response

Hypothetical data. After a step increase in drug concentration, occupancy rapidly reaches and remains at a level that implies \(R_{\mathrm{eq}}=80\) response units. The response at the moment of the step is \(R(0)=20\) units, and the time constant is \(\tau=10\) minutes. Predict the response at 10 and 30 minutes.

At 10 minutes,

$$ R(10)=80+(20-80)e^{-10/10} =80-60e^{-1} \approx 80-60(0.3679) =57.93\ \text{units} $$

At 30 minutes,

$$ R(30)=80+(20-80)e^{-30/10} =80-60e^{-3} \approx 80-60(0.04979) =77.01\ \text{units} $$

Occupancy is already at its new level in this example, while response is still below its equilibrium value at both times. The response approaches 80 units as time increases. A low early response therefore does not, by itself, show that occupancy failed to occur or that the equilibrium occupancy-response relationship is incorrect. The prediction depends on the assumed step change, constant \(R_{\mathrm{eq}}\), and single time constant; if occupancy continues to change, this constant-equilibrium solution does not apply directly.

Target State: The Measured Population May Not Be the Functional Population

A binding assay may count target molecules that are not all functionally equivalent. Some may be in a signaling-competent state, while others are inactive, differently localized, or otherwise unable to contribute to the measured endpoint. If ligand binding occurs in both groups, aggregate occupancy can differ substantially from occupancy of the functionally competent group. A mismatch is then informative about the target population, but it cannot be interpreted simply as a difference in drug potency.

The following example separates these quantities. It assumes that the binding assay counts all target molecules equally, that the stated occupancy fractions apply to the indicated groups, and that response is proportional to occupancy of competent target only. These assumptions are deliberately explicit; a real assay may require a more detailed model.

Worked Example: Aggregate Binding Includes Inactive Target

Hypothetical data. A sample contains 100 target units: 60 are signaling-competent and 40 are not. The drug occupies 50% of competent target and 90% of noncompetent target. The binding assay reports total occupied target as a fraction of all 100 units. The functional assay is calibrated so that occupying all 60 competent units gives a 100-unit response, with response proportional to occupied competent target.

The number of occupied competent target units is \(60(0.50)=30\). The number of occupied noncompetent units is \(40(0.90)=36\). Thus aggregate occupancy is

$$ \theta_{\mathrm{aggregate}} =\frac{30+36}{100} =0.66 $$

The binding assay therefore reports 66% occupancy. But only 30 of the 60 competent units are occupied, so the response under the stated calibration is

$$ R=100\left(\frac{30}{60}\right)=50\ \text{response units} $$

The 66% aggregate occupancy and 50-unit response are consistent with the assumed target-state model. They do not imply a measurement error. If the binding assay selectively measured competent target, its occupancy estimate would be 50%, not 66%. Determining which interpretation is appropriate requires knowing what target states each assay detects.

Practical Comparison and Diagnostics

Several models can be considered, but their assumptions differ. The most useful first step is to check assay alignment and inspect the pattern of discrepancy before choosing a more complex explanation.

ApproachKey assumptionPattern it can explainMain limitation
Direct proportional linkResponse above baseline scales linearly with occupancyA roughly proportional response across the observed rangeMay miss amplification or saturation
Nonlinear transduction linkA specified nonlinear function maps occupancy to responseLarge response at partial occupancy or a curved occupancy-response patternThe function is not identified by one matched measurement
Dynamic response modelResponse approaches a defined equilibrium with a specified time courseResponse changing after occupancy has already changedA single time constant may not capture feedback or multiple delays
Target-state stratificationAssays differ in which target states they count or those states differ in functionAggregate binding that does not track the competent target fractionRequires evidence about target states and assay selectivity

A practical diagnostic sequence is to first verify that occupancy and response are expressed on well-defined scales and refer to the same biological setting. Next, compare their sampling times: a systematic response lag after a rapid occupancy change supports investigating dynamics. Then ask whether the binding and functional assays sample the same target states. Finally, assess whether a nonlinear transduction relation is plausible across multiple exposure or occupancy levels. No single discrepancy pattern proves one mechanism, because timing, target state, and nonlinear transduction can coexist.

Warning signs include comparing occupancy at peak exposure with a delayed response measured later, using total target in one assay and competent target in another, or drawing a general occupancy-response curve from only one condition. If response and occupancy diverge in a way that changes with time, additional time-matched measurements can test a lag explanation. If the mismatch persists at stable conditions, state-selective assays or independent evidence about target competence may be informative. If response is consistently large at partial occupancy across matched conditions, amplification is one possible explanation, but it should be evaluated against alternatives rather than assumed.

Common Pitfall and Limitations

Common mistake: “The response is only half of its maximum, so occupancy must be 50%.” This conclusion assumes a direct proportional link and matched timing and target populations. If transduction is nonlinear, the response can be high at partial occupancy; if response is delayed or binding assays count nonfunctional target, the two measurements can differ for other reasons. Estimate occupancy from an appropriate binding model and interpret response through an explicitly stated functional and temporal model.

These comparisons have limits. A good fit to a nonlinear link does not establish the biological mechanism behind that curve. A fitted time constant summarizes the response dynamics under the stated model; it does not prove that the underlying system has only one process. Likewise, aggregate occupancy alone cannot reveal how binding is distributed across functional target states. Mismatches are therefore best treated as prompts to test assumptions and collect measurements that distinguish explanations, not as automatic evidence for or against target engagement.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. At equilibrium, what assumptions are required to use \(\theta=C/(K_D+C)\), and what does \(\theta\) represent?
  2. A matched assay estimates \(\theta=0.25\). Under \(e=(1+\kappa)\theta/(\kappa+\theta)\) with \(\kappa=0.05\), calculate \(e\).
  3. Under the first-order response model, \(R_{\mathrm{eq}}=60\) units, \(R(0)=60\) units, and \(\tau=5\) minutes. What is \(R(t)\)?
  4. A binding assay reports 80% aggregate occupancy, but the functional endpoint is weak. Give two distinct explanations that should be checked before concluding that the binding result is wrong.

Solutions. (1) The simple relation assumes one homogeneous target population, one reversible ligand, free concentration at the target site, and equilibrium without competition or cooperativity. \(\theta\) is the fraction of that specified target population bound. (2) Substitution gives \(e=(1.05)(0.25)/(0.05+0.25)=0.2625/0.30=0.875\), or 87.5% of the modeled maximum response. (3) Since \(R(0)=R_{\mathrm{eq}}\), the response remains at equilibrium: \(R(t)=60\) units for all \(t\). (4) Response may lag occupancy, or the binding assay may count target states that do not contribute to the functional endpoint. Nonlinear or attenuated transduction is another possibility. These explanations require different evidence and should not be treated as interchangeable.

Clinical Biostats