Why Keep the Two T-Cell Populations Separate?
CAR T cells and endogenous tumor-reactive T cells can both contribute to tumor killing, but they do not recognize tumor in the same way. A CAR recognizes its target through the engineered receptor, whereas an endogenous T cell responds through its native antigen receptor and depends on antigen presentation. Tumor damage caused by CAR T cells may therefore influence endogenous activity indirectly: released tumor material can become available for presentation, and local inflammatory cues can support endogenous T-cell activation. These connections are not a reason to combine both populations into one generic effector compartment.
A useful QSP model preserves the two recognition routes and then represents the links between them. This makes it possible to ask whether a modeled CAR effect is direct killing, indirect recruitment or activation of endogenous cells, or both. The equations below are a deliberately compact representation of those mechanisms, not a claim that all tumors or CAR products behave identically.
A Compact Interaction Model
Assume that CAR T-cell abundance is specified over the period being modeled, or supplied by a separate CAR-cell model. This isolates the interactions of interest without re-deriving CAR activation or exhaustion. Let \(a_C\) and \(a_E\), each between zero and one, represent the fractions of tumor burden effectively recognizable by the CAR and endogenous populations. The saturable killing fluxes are
Here \(k_C\) and \(k_E\) are maximum killing rates per unit of the corresponding T-cell population, and \(H_C\) and \(H_E\) are half-saturation constants expressed in the same tumor-burden units as \(T\). These parameters need not be equal: the populations can differ in recognition, potency, and access to targets. If a population has no recognizable target (\(a_j=0\)), its corresponding killing flux is zero.
A minimal set of state equations is
The tumor equation includes logistic growth, with rate \(r\) and carrying capacity \(T_{\max}\), as well as both killing fluxes. The \(P\) equation is a lumped representation of tumor-derived antigen becoming available for presentation; \(\eta_C\) and \(\eta_E\) allow the two killing routes to contribute differently. It does not imply that CAR T cells themselves present antigen. The \(I\) equation represents a local inflammatory cue, with production coefficients \(\gamma_C\) and \(\gamma_E\) and clearance rate \(\delta_I\). The endogenous T-cell equation includes baseline input \(s_E\), loss at rate \(d_E\), and signal-dependent expansion or activation represented by \(r_E\). \(K_P\) and \(K_I\) are half-saturation constants for presentation and inflammatory signal.
The multiplicative signal terms express an assumption: in this illustrative model, endogenous expansion is strongest when both tumor-antigen presentation and local inflammatory support are present. They are phenomenological factors, not a detailed model of antigen-presenting cells, T-cell priming, trafficking, or cytokine networks. All rate constants and state units must be chosen consistently; for example, if time is measured in days, rates have units per day.
When Can Endogenous T Cells Expand?
The model gives a direct condition for the per-capita change in \(E\). First consider a period with \(s_E=0\). Dividing its equation by \(E>0\) gives
Thus, endogenous cells expand at that instant only if the signal-dependent expansion term exceeds their loss rate. CAR killing alone does not guarantee that condition: killing may produce little effective presentation, little local inflammatory signal, or both. Conversely, if \(s_E>0\), the total \(E\) can increase from input even when its per-capita net growth term is negative; the two contributions should not be confused.
A short derivation clarifies the link between killing and the two mediators. If \(K_C\) and \(K_E\) are approximately constant over a period, each mediator equation has the form \(dX/dt=Q-\delta X\), where \(Q\) is its production flux. At steady state, setting the derivative to zero gives \(X=Q/\delta\). Therefore,
These expressions are useful for examining the mechanism, but they are not a substitute for simulating transient dynamics when mediator clearance is slow or killing changes rapidly. They also show why antigen supply and inflammatory support should be represented separately: their source coefficients and clearance rates can differ.
Worked Example: Separating the Two Killing Contributions
Consider a hypothetical tumor with \(T=2\) tumor-burden units, \(C=1\) CAR-cell unit, and \(E=0.5\) endogenous-cell units. Let \(a_C=1\), \(H_C=2\), and \(k_C=3\) tumor-burden units per CAR-cell unit per day. For the endogenous population, let \(a_E=0.6\), \(H_E=1\), and \(k_E=2\) tumor-burden units per endogenous-cell unit per day.
The CAR killing flux is
tumor-burden units per day. The endogenous killing flux is
tumor-burden units per day. The combined modeled killing is \(1.50+0.545=2.045\) units per day, but the separate terms remain available for calculating their distinct downstream contributions.
For an instantaneous tumor-growth calculation, take \(r=0.4\) per day and \(T_{\max}=10\) tumor-burden units. Logistic growth is \(0.4(2)(1-2/10)=0.640\) units per day, so
tumor-burden units per day at this state. This is a model-predicted instantaneous decline, not a guarantee that the tumor will keep declining as the states change.
Finally, let \(\eta_C=0.8\), \(\eta_E=0.5\), \(\delta_P=0.5\) per day, and \(P=1\). Then \(dP/dt=0.8(1.50)+0.5(0.545)-0.5(1)=0.973\) presentation units per day. The calculation preserves the distinction between the killing sources while allowing both to contribute to presentation.
Worked Example: Checking the Endogenous Expansion Condition
Suppose constant killing fluxes over a sufficiently long hypothetical interval are \(K_C=1\) and \(K_E=0\) tumor-burden units per day. Set \(\eta_C=1\), \(\delta_P=0.5\) per day, \(\gamma_C=0.4\), and \(\delta_I=0.4\) per day. Under the constant-flux steady-state approximation,
Now let \(K_P=0.5\), \(K_I=0.5\), \(r_E=0.75\) per day, \(d_E=0.30\) per day, and \(E=0.4\) endogenous-cell units. The two signal factors are \(2/(0.5+2)=0.8\) and \(1/(0.5+1)=0.667\). The per-capita net rate is
per day. With \(s_E=0\), the predicted instantaneous change is \(dE/dt=0.100(0.4)=0.040\) endogenous-cell units per day. The result depends on both mediator levels; it does not mean that CAR cells have become endogenous tumor-reactive cells.
Why More CAR Killing Does Not Necessarily Mean More Endogenous Activity
The coefficients linking tumor killing to presentation and inflammatory signal determine whether damage couples to endogenous-cell activity. They are not interchangeable. For example, an intervention or tumor context could permit substantial CAR-mediated killing while producing little of the modeled inflammatory cue. The following hypothetical comparison holds killing and antigen-presentation production fixed and changes only inflammatory-signal production.
Worked Example: Same CAR Killing, Different Local Signal
Let \(K_C=2\) tumor-burden units per day, \(K_E=0\), \(\eta_C=1\), \(\delta_P=0.5\) per day, and \(K_P=1\). Then \(P^*=2/0.5=4\), so the presentation factor is \(4/(1+4)=0.8\). Also take \(K_I=0.5\), \(r_E=0.9\) per day, \(d_E=0.35\) per day, \(E=0.3\), and \(s_E=0\).
In scenario A, set \(\gamma_C=0.1\) and \(\delta_I=0.5\) per day. Then \(I^*=0.1(2)/0.5=0.4\), and the inflammatory factor is \(0.4/(0.5+0.4)=0.444\). The endogenous per-capita rate is
per day. Thus \(dE/dt=-0.030(0.3)=-0.0090\) endogenous-cell units per day at that steady mediator state.
In scenario B, retain the same CAR killing and presentation parameters but set \(\gamma_C=0.5\). Now \(I^*=0.5(2)/0.5=2\), and the inflammatory factor is \(2/(0.5+2)=0.8\). The per-capita rate becomes \(0.9(0.8)(0.8)-0.35=0.226\) per day, giving \(dE/dt=0.226(0.3)=0.0678\) endogenous-cell units per day. In this model, the difference arises from the assumed coupling to local signal, not from a change in CAR killing or tumor antigen presentation.
Assumptions, Diagnostics, and Alternatives
The equations are most useful when their simplifying assumptions are visible. They assume that the chosen recognition fractions summarize target availability; each killing flux has a saturable response; antigen presentation and local signal can each be represented by one well-mixed state; and the signal-dependent term is an adequate summary of endogenous-cell activation over the modeled interval. The steady-state calculations further assume that killing fluxes are approximately constant long enough for mediator levels to approach steady state. When these conditions are not credible, simulate the mediator dynamics directly and consider richer states supported by data.
Useful checks include the following:
- Check state and flux units. Each killing flux must have units of tumor burden per time, and each source-minus-clearance equation must have units of its state per time.
- Check the no-recognition limits. Setting \(a_C=0\) must eliminate CAR-mediated killing, and setting \(a_E=0\) must eliminate endogenous-mediated killing.
- Check mediator behavior. With no production, positive \(P\) and \(I\) should decay under positive clearance rates. Under constant nonnegative production, the mediator steady states should be nonnegative.
- Check attribution against data. Tumor shrinkage alone may constrain total killing but may not identify how much came from each population or whether endogenous activity changed. Separate measurements or informative perturbations are needed to distinguish these mechanisms.
A common mistake is to add \(C\) and \(E\) into one effector pool and apply a single killing and activation rule. That may reproduce one tumor trajectory, but it removes the distinction between CAR recognition and endogenous antigen presentation and makes the indirect link from tumor damage to endogenous activity impossible to interpret. The correction is to retain separate states and killing fluxes, then specify which fluxes contribute to each mediator and which signals affect \(E\).
Compared with the compact ODE model here, a more detailed alternative could explicitly track antigen-presenting cells, peptide–MHC complexes, priming in lymphoid tissue, trafficking, or multiple endogenous T-cell states. Such detail is warranted when the data support it and the research question requires it; otherwise, extra states can be weakly identifiable. At the other extreme, a purely empirical model could let \(E\) change with time or tumor burden without an explicit antigen and signal pathway. That may be useful for prediction within a narrow setting, but it is less informative about the mechanism. No formulation removes the need to justify its recognition assumptions, parameter values, and observation model.
This tutorial treats \(I\) only as a local, lumped input to endogenous activity. It does not represent systemic inflammatory feedback, a detailed cytokine network, or a general account of clinical toxicity. Those mechanisms require additional states and evidence rather than being inferred from this local coupling term.
Check Your Understanding
Answer each question using the model definitions above; solutions are included for checking your reasoning.
- If \(a_C=0\) while \(a_E>0\), which killing flux must be zero, and why?
- With constant \(K_C=0.6\), \(K_E=0\), \(\eta_C=0.5\), and \(\delta_P=0.3\) per day, what is \(P^*\)?
- At \(P=1\), \(I=1\), \(K_P=1\), \(K_I=1\), \(r_E=0.8\) per day, and \(d_E=0.25\) per day, is the per-capita endogenous rate positive or negative when \(s_E=0\)?
- Why can combining \(C\) and \(E\) into one effector state obscure the mechanism, even if the combined state predicts tumor decline?
Solutions. 1. \(K_C=0\), because its numerator contains \(a_C T\), which is zero; \(K_E\) need not be zero because its recognition fraction is positive. 2. \(P^*=(0.5(0.6))/0.3=1.0\) presentation unit. 3. Each signal factor is \(1/(1+1)=0.5\), so the per-capita rate is \(0.8(0.5)(0.5)-0.25=0.200-0.250=-0.050\) per day; it is negative. 4. The combined state conceals which receptor-mediated recognition route caused killing and whether tumor damage generated presentation and local signals that changed endogenous activity.