Tutorials › Quantitative Systems Pharmacology (QSP) › Antigen Heterogeneity in Dual-Target Engagement

Bispecific Antibodies and Antibody-Drug Conjugates · Tutorial 471 of 799

Antigen Heterogeneity in Dual-Target Engagement

Represent antigen heterogeneity explicitly to distinguish the fraction of tumor cells reached by a dual-target agent from the fraction predicted to be killed.

Advanced 11 min read

What You'll Learn

  • Define tumor fractions for cells expressing neither target, either target alone, or both.
  • Calculate single-target coverage, dual-target union coverage, and dual-positive overlap.
  • Compute total predicted killing by weighting subpopulation-specific killing probabilities.
  • Identify how overlap between antigen-positive populations changes the benefit of dual targeting.
  • Calculate how treatment changes the composition of the surviving tumor.
  • Recognize when binary phenotypes or independent functional routes are inadequate.

Why Tumor Heterogeneity Changes the Coverage Calculation

A dual-target agent may recognize two tumor antigens, but that does not mean every tumor cell expresses both, or either. A tumor can contain cells expressing antigen A only, antigen B only, both antigens, or neither. These groups matter because the fraction of cells that can be reached is not necessarily the fraction that will be killed. Coverage describes which cells meet a stated antigen criterion; predicted killing also depends on how effectively treatment acts on each group.

The useful starting point is therefore a mixture representation. It makes overlap explicit, avoids treating two antigen prevalences as if they belonged to separate populations, and allows functional response to be calculated for each subpopulation before being combined.

Definition: Let \(f_{ij}\) denote the fraction of tumor cells with antigen A status \(i\) and antigen B status \(j\), where 1 means positive and 0 means negative according to prespecified positivity thresholds. The four fractions are \(f_{00}, f_{10}, f_{01}, f_{11}\), and they satisfy \(f_{ij}\geq 0\) and \(f_{00}+f_{10}+f_{01}+f_{11}=1\). The marginal fractions positive for A and B are \(f_{10}+f_{11}\) and \(f_{01}+f_{11}\), respectively.

Coverage Is a Set-Overlap Quantity

If a single-target treatment can reach cells positive for its target, A-target coverage is \(C_A=f_{10}+f_{11}\), while B-target coverage is \(C_B=f_{01}+f_{11}\). If a dual-target treatment can reach a cell when either target is present, its union coverage is

$$ C_{\mathrm{either}}=f_{10}+f_{01}+f_{11}=1-f_{00}=C_A+C_B-f_{11}. $$

The subtraction of \(f_{11}\) prevents dual-positive cells from being counted twice. If the specified mechanism instead requires both target antigens on the same cell, the relevant coverage is \(C_{\mathrm{both}}=f_{11}\), not \(C_{\mathrm{either}}\). Thus, “dual-target coverage” is incomplete unless the criterion—either target or both targets—is stated.

These formulas assume that the categories represent the tumor population of interest at the time being modeled, that each cell is assigned to one category, and that the positivity thresholds are defined consistently. They describe antigen-based reach, not proof of binding, synapse formation, or cell death. Functional response requires a separate calculation.

Worked Example: Calculating Coverage From Four Subpopulations

Consider a hypothetical tumor mixture with 12% cells negative for both antigens, 28% positive for A only, 22% positive for B only, and 38% positive for both. The fractions sum to \(0.12+0.28+0.22+0.38=1.00\), so they form a complete mixture.

A-target coverage is \(C_A=0.28+0.38=0.66\), or 66%. B-target coverage is \(C_B=0.22+0.38=0.60\), or 60%.

For a treatment that can reach cells expressing either target, \(C_{\mathrm{either}}=0.28+0.22+0.38=0.88\), or 88%. Equivalently, \(C_A+C_B-f_{11}=0.66+0.60-0.38=0.88\). The 38% dual-positive fraction is counted once. If both antigens were required, coverage would instead be \(C_{\mathrm{both}}=0.38\), or 38%.

Why Marginal Antigen Prevalences Are Not Enough

Knowing that a certain proportion of cells are A-positive and a certain proportion are B-positive does not reveal how much those groups overlap. Because overlap determines the size of the union, marginal prevalences alone cannot determine dual-target coverage. This is especially important when data are reported separately for each marker rather than jointly at the single-cell level.

Worked Example: Same Marginals, Different Dual-Target Coverage

Suppose two hypothetical tumor samples each have 60% A-positive cells and 50% B-positive cells. In Sample 1, let \(f_{11}=0.40\). Then \(f_{10}=0.60-0.40=0.20\), \(f_{01}=0.50-0.40=0.10\), and \(f_{00}=1-(0.40+0.20+0.10)=0.30\). Its either-target coverage is \(1-0.30=0.70\), or 70%.

In Sample 2, let \(f_{11}=0.20\). Then \(f_{10}=0.40\), \(f_{01}=0.30\), and \(f_{00}=0.10\). Its either-target coverage is \(1-0.10=0.90\), or 90%.

Both samples have identical marginal prevalences, but their dual-target union coverages differ by 20 percentage points because their overlap differs. A model built only from the two marginal prevalences cannot distinguish these cases. Joint antigen measurement, or an explicit assumption about overlap, is needed.

From Coverage to Predicted Killing

Let \(k_{ij}\) be the probability that a cell in subpopulation \(ij\) is killed over a specified treatment interval and under specified exposure conditions. The \(k_{ij}\) values may be estimated from suitable functional data or supplied by a mechanistic model. They need not be the same across antigen-positive groups: antigen density, accessibility, signaling, and other cell properties may differ. The total predicted fraction killed is the mixture-weighted average:

$$ K_{\mathrm{total}}=f_{00}k_{00}+f_{10}k_{10}+f_{01}k_{01}+f_{11}k_{11}. $$

This result assumes the \(f_{ij}\) values describe the starting population for the interval, the \(k_{ij}\) values apply to that interval and exposure, and each starting cell belongs to exactly one subpopulation. It does not require killing outcomes in different cells to be independent; it uses their expected proportions. If antigen status or exposure changes materially during treatment, a time-varying or state-transition model is needed instead of applying one fixed mixture for the whole course.

A useful illustrative assumption is that each antigen-specific functional route kills a positive cell with a specified probability and that the two routes act independently on a dual-positive cell. If their probabilities are \(e_A\) and \(e_B\), the probability that neither route kills is \((1-e_A)(1-e_B)\). Therefore, the dual-positive killing probability under this assumption is \(k_{11}=1-(1-e_A)(1-e_B)\). This is a modeling choice, not a general property of bispecific agents: shared signaling limits, correlated responses, or other interactions can make the independence assumption inappropriate.

Worked Example: Predicting Killing Across a Heterogeneous Mixture

Use the hypothetical mixture above: \(f_{00}=0.12\), \(f_{10}=0.28\), \(f_{01}=0.22\), and \(f_{11}=0.38\). Suppose that under a stated exposure, the A-mediated route kills 50% of A-positive cells and the B-mediated route kills 40% of B-positive cells. Assume cells negative for both are not killed by this mechanism, and assume the two routes act independently on dual-positive cells.

The subpopulation killing probabilities are \(k_{00}=0\), \(k_{10}=0.50\), and \(k_{01}=0.40\). For dual-positive cells, \(k_{11}=1-(1-0.50)(1-0.40)=1-(0.50)(0.60)=0.70\).

Substituting into the mixture equation gives \(K_{\mathrm{total}}=(0.12)(0)+(0.28)(0.50)+(0.22)(0.40)+(0.38)(0.70)\). The contributions are \(0+0.140+0.088+0.266=0.494\). The model therefore predicts that 49.4% of the starting tumor cells are killed over the specified interval.

For comparison, an A-only treatment with the same assumed A-route activity has killing probabilities 0 for \(f_{00}\) and \(f_{01}\), and 0.50 for \(f_{10}\) and \(f_{11}\). Its predicted total killing is \((0.28+0.38)(0.50)=0.330\), or 33.0%. The difference in this calculation is 16.4 percentage points. It follows from the stated mixture and response assumptions; it is not a universal estimate of the benefit of dual targeting.

A Short Justification: Why the Mixture Average Applies

Imagine \(N\) starting cells, with approximately \(Nf_{ij}\) cells in category \(ij\). If a fraction \(k_{ij}\) of that category is killed, the expected number killed from the category is \(Nf_{ij}k_{ij}\). Adding across the mutually exclusive categories gives the expected total killed:

$$ E[\text{cells killed}]=N\sum_{i,j}f_{ij}k_{ij}. $$

Dividing by \(N\) gives \(K_{\mathrm{total}}=\sum_{i,j}f_{ij}k_{ij}\). The calculation is a weighted average because each category contributes in proportion to its starting share. It would be incorrect to average the four \(k_{ij}\) values equally unless all four subpopulations had equal fractions.

Worked Example: Interpreting the Surviving Mixture

Continue the previous hypothetical calculation. The survival probabilities are \(1-k_{ij}\): 1.00 for double-negative cells, 0.50 for A-only cells, 0.60 for B-only cells, and 0.30 for dual-positive cells. The surviving fractions of the original tumor are therefore \(0.12\), \(0.28(0.50)=0.140\), \(0.22(0.60)=0.132\), and \(0.38(0.30)=0.114\).

The total surviving fraction is \(0.12+0.140+0.132+0.114=0.506\), consistent with \(1-K_{\mathrm{total}}=1-0.494=0.506\). Among survivors, the double-negative fraction is \(0.12/0.506=0.237\), or 23.7%; the A-only fraction is \(0.140/0.506=27.7\%\); the B-only fraction is \(0.132/0.506=26.1\%\); and the dual-positive fraction is \(0.114/0.506=22.5\%\).

The initial dual-positive group made up 38% of the tumor but 22.5% of the predicted survivors. These survivor proportions are conditional on the modelled treatment interval and assumptions. They illustrate how differential killing can reshape a mixture; they do not establish how the tumor will evolve over longer periods.

Model Choices, Diagnostics, and Limitations

A four-category model is transparent and often useful when measurements support clear positive and negative groups. It is not the only option. If antigen abundance varies continuously, replacing each group with a single binary status can discard important information. A continuous-density model represents each cell by its A and B antigen levels and averages a response function over their joint distribution. In notation, if \(g(a,b)\) is the joint density of antigen levels and \(k(a,b)\) is the predicted killing probability at those levels, then total killing is \(\int\int k(a,b)g(a,b)\,da\,db\). This approach requires data or justified assumptions about the joint distribution and response function; separate average antigen levels do not generally specify either one.

There is also a difference between an empirical subpopulation model and a mechanistic binding-state model. The empirical approach uses \(k_{ij}\) values directly and can be practical when functional data are available. A mechanistic model may represent binding, bridging, or signaling states and derive \(k_{ij}\) from those processes. The state-based and synapse-signaling concepts provide a basis for such extensions, but they do not remove the need to represent the mixture of tumor phenotypes. More detailed models introduce more parameters and need corresponding measurements for identification.

Useful diagnostics begin with bookkeeping. Check that all four fractions are nonnegative and sum to one, and that the marginals agree with the joint fractions. Compare model-predicted killing with observed responses separately for antigen-defined subgroups where possible, rather than validating only against an overall average. Large systematic errors in one subgroup can be hidden by a plausible total. Sensitivity analysis over uncertain overlap and killing probabilities is also informative: if the predicted benefit of dual targeting changes substantially across plausible values, the conclusion depends on information the current data do not determine.

Other warning signs include assigning all antigen-positive cells the same response despite substantial density variation, treating marginal antigen prevalences as independent without evidence, and applying a fixed mixture after treatment has selected among the subpopulations. Positivity thresholds also need care: assay noise or threshold changes can move cells between categories and alter calculated coverage. When measurements are noisy or the boundaries are biologically gradual, a continuous or probabilistic classification may be more appropriate.

The independence calculation for two functional routes is particularly easy to overinterpret. It assumes that, conditional on the cell being dual-positive and on the specified exposure, the two route-specific killing events combine as independent opportunities. If both routes depend on the same limiting downstream process, then \(1-(1-e_A)(1-e_B)\) may overpredict killing. In that case, use an interaction model or estimate \(k_{11}\) from data rather than imposing the independence formula.

Common pitfall: Adding A-positive and B-positive percentages and calling the result dual-target coverage double-counts cells positive for both antigens. For example, 66% A-positive plus 60% B-positive is not 126% coverage. Subtract the overlap: \(66\%+60\%-38\%=88\%\). Also keep coverage distinct from predicted killing: cells may meet an antigen criterion without being killed.
Key takeaway: Represent the joint antigen mixture first. Calculate coverage from the union or intersection required by the mechanism, then predict total killing by weighting each subpopulation’s response by its starting fraction. Report the assumptions that connect antigen status to function.

Check Your Understanding

Work each question, then compare with the solutions below.

  1. A hypothetical mixture has \(f_{00}=0.20\), \(f_{10}=0.30\), \(f_{01}=0.10\), and \(f_{11}=0.40\). What are A coverage, B coverage, and either-target coverage?
  2. Why do A-positive and B-positive marginal fractions alone fail to determine either-target coverage?
  3. For the mixture in Question 1, suppose \(k_{00}=0\), \(k_{10}=0.40\), \(k_{01}=0.50\), and \(k_{11}=0.60\). Calculate total predicted killing.
  4. Give one reason the independent-route formula for dual-positive killing might not apply.

Solutions. 1. A coverage is \(0.30+0.40=0.70\), or 70%. B coverage is \(0.10+0.40=0.50\), or 50%. Either-target coverage is \(0.30+0.10+0.40=0.80\), or 80%. 2. The marginal fractions do not specify the overlap \(f_{11}\), which must be subtracted when forming the union: \(C_{\mathrm{either}}=C_A+C_B-f_{11}\). 3. \(K_{\mathrm{total}}=(0.20)(0)+(0.30)(0.40)+(0.10)(0.50)+(0.40)(0.60)=0+0.120+0.050+0.240=0.410\), or 41.0%. 4. The routes may share a limiting downstream process or have correlated functional outcomes, so their combined probability need not equal \(1-(1-e_A)(1-e_B)\).

Clinical Biostats