From Engager Dose to Cell Killing
A T-cell engager must connect a T cell to a target cell for a productive interaction, but the concentration that occupies each binding arm is not necessarily the concentration that maximizes cell bridging. At high engager concentrations, the two cell types can become separately occupied by different engager molecules, leaving fewer available opportunities to form a bridge. A useful quantitative model should therefore distinguish arm occupancy, bridge abundance, and killing, rather than treating them as interchangeable measurements.
The model here is a deliberately reduced, hypothetical example. It compares two tumor-cell populations, one with low and one with high target-antigen density, as free engager concentration changes. Its equations are not a complete molecular binding model; their purpose is to make assumptions and predictions explicit.
A Competition-Adjusted Bridging Model
The occupancy equations are simple saturable binding expressions. To represent the loss of productive bridges when high free drug separately occupies the two cell types, the model includes a competition factor \(K_H/(K_H+C)\). This factor is a phenomenological reduction in bridge formation, not an additional measured arm occupancy. The resulting bridge equation is
For the worked setting, \(K_E=2\) ng/mL, \(K_{A,L}=10\) ng/mL, \(K_{A,H}=3\) ng/mL, and \(K_H=10\) ng/mL. The maximum bridge scales are \(B_{\max,L}=2350.53\) and \(B_{\max,H}=2366.112\) bridge units per target cell. These hypothetical scales represent population-specific capacities; they are not inferred from occupancy alone.
At low concentration, both occupancies are small, so bridge abundance initially increases with dose. At high concentration, the competition factor becomes small. In fact, as \(C\) becomes large, the two occupancy fractions approach one while \(K_H/(K_H+C)\) approaches zero. Thus \(B_j(C)\) eventually declines. This model can produce a bell-shaped bridge response without claiming that every engager has the same response curve.
Worked Example: Occupancy and Bridging at 5 ng/mL
Worked Example: Occupancy and Bridging at 5 ng/mL
Calculate each arm occupancy and the predicted bridge abundance at \(C=5\) ng/mL. Then translate bridge abundance into killing over 24 hours using the hypothetical relationship \(k_j=\gamma B_j\), where \(\gamma=1.00\times10^{-4}\) per bridge unit per hour.
Effector-arm occupancy. Both populations share the same effector-arm parameter: \(\theta_E=5/(2+5)=5/7=0.7143\), or 71.43%.
Target-arm occupancy. For low-density cells, \(\theta_{A,L}=5/(10+5)=0.3333\), or 33.33%. For high-density cells, \(\theta_{A,H}=5/(3+5)=0.6250\), or 62.50%.
Low-density bridges. The competition factor is \(10/(10+5)=0.6667\). Therefore, \(B_L=2350.53(0.714286)(0.333333)(0.666667)=373.1\) bridge units per target cell.
High-density bridges. \(B_H=2366.112(0.714286)(0.625000)(0.666667)=704.2\) bridge units per target cell.
Killing over 24 hours. With a constant hazard, the surviving fraction is \(S_j(t)=\exp(-k_jt)\), so the killed fraction is \(1-S_j(t)\). The low-density hazard is \(k_L=(1.00\times10^{-4})(373.1)=0.03731\) per hour. Its predicted killing is \(100[1-\exp\{-(0.03731)(24.00)\}]=59.16\%\). The high-density hazard is \(k_H=(1.00\times10^{-4})(704.2)=0.07042\) per hour, giving \(100[1-\exp\{-(0.07042)(24.00)\}]=81.55\%\).
The predicted difference is \(81.55-59.16=22.39\) percentage points in this hypothetical setting. It is a model-based contrast between the two populations, not evidence that antigen density alone causes the same difference in a biological experiment.
Why the Dose Response Can Turn Down
The high-dose behavior follows directly from the bridge equation. Taking the logarithmic derivative with respect to \(C\) gives
At very low \(C\), the \(2/C\) term dominates, so bridge abundance rises with concentration. At very high \(C\), the four terms together behave approximately as \(2/C-3/C=-1/C\), which is negative. The bridge abundance must therefore turn downward somewhere between those regions. For the low-density parameters, setting this derivative to zero yields \(C^2-10C-40=0\); the positive solution is \(C=(10+\sqrt{260})/2=13.06\) ng/mL. This is the model's low-density bridge-abundance peak, not a universal optimal clinical dose.
Worked Example: Comparing Three Doses
Compare predicted occupancies and bridges at 1, 5, and 20 ng/mL. At 1 ng/mL, the occupancies are \(\theta_E=1/3=0.3333\), \(\theta_{A,L}=1/11=0.09091\), and \(\theta_{A,H}=1/4=0.2500\); the competition factor is \(10/11=0.9091\). Substitution gives \(B_L=2350.53(0.3333)(0.09091)(0.9091)=64.75\) and \(B_H=2366.112(0.3333)(0.2500)(0.9091)=179.3\) bridge units per target cell.
At 5 ng/mL, the previous calculation gives 373.1 and 704.2 bridge units per target cell. At 20 ng/mL, \(\theta_E=20/22=0.9091\), the target occupancies are \(20/30=0.6667\) and \(20/23=0.8696\), and the competition factor is \(10/30=0.3333\). Thus \(B_L=2350.53(0.9091)(0.6667)(0.3333)=474.9\), while \(B_H=2366.112(0.9091)(0.8696)(0.3333)=623.5\).
The high-density population has greater predicted bridge abundance at all three doses, but its bridge abundance at 20 ng/mL is below its value at 5 ng/mL. For example, the high-density 20 ng/mL hazard is \(0.06235\) per hour, giving \(100[1-\exp\{-(0.06235)(24)\}]\approx77.61\%\) killing. This is below the 81.55% prediction at 5 ng/mL. The model illustrates why high individual-arm occupancy does not guarantee maximal bridging or killing.
From Bridge Abundance to Selective Killing
The killing link used above is a constant-hazard model. If each cell has a constant instantaneous killing hazard \(k_j\), then the probability of remaining alive through time \(t\) is \(S_j(t)=\exp(-k_jt)\), and the fraction killed is \(1-S_j(t)\). This follows from the survival equation \(dS_j/dt=-k_jS_j\) with \(S_j(0)=1\): separating variables and integrating gives \(\log S_j(t)=-k_jt\), hence the exponential expression. The units are consistent because \(k_j\) is per hour and \(t\) is in hours.
In this example, \(k_j=\gamma B_j\) assumes killing hazard is directly proportional to bridge abundance. That is a simplifying assumption, not a general law. It is most defensible over a range where bridge abundance is a useful proxy for productive contacts and where effector-cell availability is not limiting. If killing saturates, an alternative could be \(k_j=k_{\max}B_j/(B_{50}+B_j)\), with \(B_{50}\) the bridge abundance giving half of the maximum hazard. Such a model needs additional parameters and data.
Worked Example: Interpreting Selectivity at a Fixed Time
At 5 ng/mL, the model predicts 59.16% killing in the low-density population and 81.55% in the high-density population by 24 hours. The absolute difference is 22.39 percentage points. The relative ratio of predicted killed fractions is \(81.55/59.16\approx1.378\), so the high-density prediction is about 1.38 times the low-density prediction under these assumptions.
These two summaries answer different questions: percentage-point difference measures the absolute separation in killed fractions, while the ratio compares their relative magnitudes. Neither is a clinical effect estimate. They are conditional predictions from the specified parameter values, time interval, and hazard relationship.
Diagnostics, Alternatives, and Limitations
A model is useful only if its assumptions remain credible for the question being asked. Check whether free concentration is approximately maintained, whether target and effector populations are sufficiently well mixed, and whether measured bridge abundance or killing supports the assumed dose dependence. If observed killing continues to rise at high concentrations while this model predicts a decline, possible explanations include an inappropriate competition factor, changing exposure, or biology absent from the reduced model. A mismatch is a reason to revisit assumptions, not to relabel the prediction as an observation.
An occupancy-only product, \(\theta_E\theta_{A,j}\), is a simpler alternative. It is easier to interpret but approaches one at high concentration and cannot represent the high-dose decline in bridges. A more detailed binding and cell-interaction model could track free drug, singly bound species, bridges, receptor turnover, and cell encounters explicitly. Such a model can represent more mechanisms, but it requires more data and identifiable parameters. The present reduced model is a transparent compromise; its competition factor and calibrated bridge scales should not be mistaken for uniquely established molecular quantities.
Other limitations matter for interpretation. The two populations are each treated as uniform, although actual tumors may contain a distribution of antigen densities. The model does not represent spatial access, effector-to-target ratios, serial killing, cytokine effects, or time-varying exposure. The assumed linear conversion from bridges to killing may fail when effector cells are saturated or when signaling has thresholds. Consequently, the predicted selectivity is conditional on the model structure and parameterization.
Common Pitfall
- Wrong reasoning: “At 20 ng/mL, both binding arms are more occupied than at 5 ng/mL, so there must be more bridges and more killing.”
- Correction: Arm occupancy and bridge abundance are distinct. In this model the competition factor falls from \(10/15=0.6667\) at 5 ng/mL to \(10/30=0.3333\) at 20 ng/mL. The resulting high-density bridge prediction falls from 704.2 to 623.5 units per cell, and the predicted 24-hour killing falls from 81.55% to about 77.61%.
Check Your Understanding
Answer each question using the model and hypothetical parameters in this tutorial. Brief solutions follow the questions.
- At 5 ng/mL, what is the effector-arm occupancy, and what are the low- and high-density target-arm occupancies?
- At 20 ng/mL, why does the high-density population have fewer predicted bridges than at 5 ng/mL even though its target-arm occupancy is higher?
- If a population has a constant killing hazard of \(0.0500\) per hour, what fraction is predicted killed after 12 hours?
- Name one assumption that could make the model's predicted density selectivity unreliable in a heterogeneous tumor.
Solutions. (1) The occupancies are \(5/7=0.7143\), \(5/15=0.3333\), and \(5/8=0.6250\), respectively. (2) The competition factor drops from 0.6667 to 0.3333; its decline more than offsets the increase in the two arm occupancies, so predicted bridges fall from 704.2 to 623.5 per target cell. (3) The killed fraction is \(1-\exp[-(0.0500)(12)]=1-\exp(-0.600)=0.4512\), or 45.12%. (4) Treating all cells in each population as having one antigen density can misrepresent a tumor with a broad density distribution; other valid answers include assuming free concentration is fixed or assuming hazard is proportional to bridges.