- 1Applied BioMath, Concord, MA, United States
- 2Differentia Biotech, San Francisco, CA, United States
Chimeric antigen receptor T (CAR T) cell therapy has shown remarkable success in treating various leukemias and lymphomas. Cellular kinetic (CK) and pharmacodynamic (PD) behavior of CAR T cell therapy is distinct from other therapies due to its living nature. CAR T CK is typically characterized by an exponential expansion driven by target binding, fast initial decline (contraction), and slow long-term decline (persistence). Due to the dependence of CK on target binding, CK and PD of CAR T therapies are inherently and bidirectionally linked. In this work, we develop a semi-mechanistic model of CAR T CK/PD, incorporating molecular-scale binding, T cell dynamics with multiple phenotypes, and tumor growth and killing. We calibrate this model to published CK and PD data for a CD19-targeting CAR T cell therapy. Using sensitivity analysis, we explore variability in response due to patient- and drug-specific properties. We further explore the impact of tumor characteristics on CAR T-cell expansion and efficacy through individual- and population-level parameter scans.
1 Introduction
Chimeric antigen receptor (CAR) T-cells are T-cells engineered to produce CARs which recognize and bind to a tumor antigen. In CAR T-cell therapy, a patient’s T-cells are extracted and isolated, re-engineered to express a specific CAR, expanded ex vivo, and then infused back into the patient. Six such therapies have been approved for treating a variety of blood cancers (Chen et al., 2023). These therapies have been shown to produce long-lasting response and superior response rates to alternative treatments (Melenhorst et al., 2022; Sermer et al., 2020). As a result of the individualized nature of CAR T manufacturing, the contents of the dosed product will vary from patient to patient. Further, CAR T cellular kinetic behavior is distinct from other therapies due to its “living” nature; it is typically characterized by an exponential expansion, fast initial decline (contraction), and slow long-term decline (persistence). Additionally, cellular kinetics (CK) is not as well-studied as pharmacokinetics for more traditional drugs Chaudhury et al. (2020). Interactions between CAR T-cells and tumor cells are complex since tumor expansion has a significant impact on CAR T-cell expansion. Furthermore, much is still unknown about the workings of CAR T-cells in the body and there is not a standard monitoring process. Modeling can shed light on CAR-T cell CK/PD and inform future studies by mechanistically linking CAR T-cell doses to tumor cell growth and determining optimal drug properties to achieve efficacy and safety. Furthermore, patient characteristics can be incorporated into the model to provide individualized dose predictions and guide patient and indication selection.
Modeling and simulation has been used to understand CAR T-cell dynamics and efficacy (see, for example, reviews by Chaudhury et al. (2020); Nukala et al. (2021)) and the impact of preconditioning (Owens and Bozic, 2021). Until recently, the three distinct phases of CAR T cellular kinetics and the impact of different CD4+ and CD8+ T cell phenotypes had not been mechanistically described. Previous modeling work had captured CAR T cellular kinetics either empirically (Stein et al., 2019), mechanistically but without multiple phenotypes, (Singh et al., 2020), or mechanistically with effector/memory phenotypes but without separating CD4+ and CD8+ T cells (Hardiansyah and Ng, 2019). Recent work by Salem et al. (2023) has incorporated all of these features, developing a mechanistic model incorporating binding-driven CAR T-cell expansion and activity for multiple CD4+ and CD8+ T-cell phenotypes to match clinical data from multiple trials. Further analysis of such models will be useful to understand system behavior, inform engineering of CAR T-cells, and understand variability in patient populations. In particular, sensitivity analysis provides understanding of the key mechanisms driving expansion and efficacy.
Here, we present a semi-mechanistic cellular kinetic-pharmacodynamic (CK-PD) model for CAR T-cell therapy of B-cell non-Hodgkin lymphoma (NHL). Our model includes CD8+ and CD4+ naive, effector, and memory T-cell phenotypes, binding of CARs to their target antigen CD19, binding-driven activation and expansion of T-cells, T-cell death and conversion to memory cells, and binding-driven killing of B cells by CD8+ effector cells. We demonstrate the ability of the model to capture published human CAR T-cell CK and PD data, and perform sensitivity analysis to understand key model features and predict the impact of variability in patient, tumor, and drug characteristics.
2 Methods
2.1 Data
The model was informed by and benchmarked to published human CAR T-cellular kinetics, B cell percentage, and clinical response data from a phase I clinical trial with IM19 CAR T-cells for 13 relapsed or refractory NHL patients (Ying et al., 2021). The CK data and the B cell aplasia data were both digitized using WebPlotDigitizer (Rohatgi, 2022). Two days prior to CAR T-cell infusion, patients were pre-treated with fludarabine and cyclophosphamide for 3 days to deplete endogenous lymphocytes. IM19 CAR T-cells were dosed by weight at
2.2 Model structure
The model consists of a single compartment representing the blood. CAR T-cell and B cell populations and their corresponding receptor burdens are modeled explicitly. CAR T-cells, a fraction of which are CD8+ and the remainder CD4+, are dosed directly into the blood. All CAR T-cells are assumed to be naive at the time of dosing. CARs on both CD8+ and CD4+ CAR T-cells can bind to CD19. CD8+ and CD4+ naive CAR T-cells are activated at a rate proportional to the fraction of CAR that is bound to CD19. Activated CAR T-cells then proliferate and become effector cells. CD8+ effector CAR T-cells can then kill B cells at a rate proportional to the fraction of CARs on CD8+ effector T cells that are bound to CD19. We assume that CD4+ effector CAR T-cells do not kill B cells as we focus only on direct effects (Alizadeh et al., 2023). Effector CAR T-cells either die or become memory cells. Memory CAR T-cells have a longer lifespan than effector cells, but do not participate in B cell killing. The model is intended to describe the initial response to CAR T therapy and therefore does not include any mechanisms for re-activation of memory cells. A diagram of the model reactions is shown in Figure 1. A more detailed description of the model equations is given below, where all states are in units of nmol.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g001.jpg)
Figure 1. Diagram showing interactions represented in the model. CAR-T cells are dosed as part CD8+, part CD4+. Drug product cells are activated by binding to CD19 on malignant B cells. Activated cells replicate and become effector cells. CD8+ effector cells kill B cells. Effector CAR T-cells either die or become memory cells.
2.3 Cell state equations
Infused CAR T cells are dosed directly into the blood (
Tumor cells are able to divide at a rate
Endogenous lymphocytes are produced at a zeroth order rate
2.4 Receptor equations
In addition to the cell-scale dynamics described above, molecular-scale dynamics are explicitly accounted for in the model. Receptor equations for CAR and CD19 are written such that the total receptor densities (CAR per T cell and CD19 per tumor cell) remain constant, as determined by a receptor per cell (RPC) parameter. Receptors are synthesized at a rate
CARs on infused CAR T cells undergo synthesis and internalization, binding/unbinding with CD19, conversion to an activated state, loss from cell death, and release from CAR:CD19 complex when a tumor cell dies. The equations for CD8+ infused CAR and CAR:CD19 complex are shown below; equations for CD4+ CAR are similar.
Since activated CAR T cells in the model are an intermediate state for the purposes of expansion and do not interact with the tumor, we do not include binding of activated CARs to CD19 in the model. CARs on activated CAR T cells follow the cellular kinetics (i.e., activation of infused CAR, division, and differentiation) so that CAR per T cell remains constant.
CARs on effector and memory CAR T cells are synthesized and internalized, bind/unbind with CD19, undergo effector-to-memory conversion, are lost through T cell death, and are released from CAR:CD19 complex when a tumor cell dies. These equations are given below for CD8+ effector and memory CAR; equations for CD4+ CARs are similar.
2.5 Model parameterization
Values for the majority of the model parameters were inferred from literature as described below. The rest of the parameters were fit to individual patient data from Ying et al. (2021), described below.
2.3.1 Patient and tumor properties
The total blood volume was estimated to be 5L, based on the average human adult (Sharma and Sharma, 2018). The concentration of endogenous lymphocytes was assumed to be
2.3.2 CAR T-cell properties
The CAR internalization half-life was estimated to be 6 h based on in vitro measurements for other CD19-targeting CAR T-cells Li et al. (2020). The mean activation time (that is, the time between binding to antigen and the start of cell proliferation) was estimated to be 18 h for CD8+ CAR T-cells (Henrickson et al., 2008; Cui and Kaech, 2010) and 36 h for CD4+ CAR T-cells (Kaech et al., 2002). Average lifespans for memory CAR T-cells were estimated to be 180 days for CD8+ and 240 days for CD4+ (Borghans et al., 2018). The CD4:CD8 ratio of the CAR T-cells for each patient were taken from Ying et al. (2021), and all infused CAR T-cells were assumed to be viable. We assumed expression levels of 12,700 CARs per T-cell for both CD8+ and CD4+ cells based on a published average estimate for a HMW-MAA-specific CAR on CD8+ T cells (Anikeeva et al., 2021). We assumed that CARs bind to CD19 with an affinity of 1 nM based on reported affinities for high affinity CAR T variants (Jayaraman et al., 2020), with a binding on-rate of 0.001/nM/s.
Remaining model parameters, namely, the number of divisions per T-cell, time per T-cell division, drug product and effector cell lifespan, memory cell fraction, and initial tumor burden, were fit to data as described in the following subsection.
2.3.3 Calibration and benchmarking
Considerable variability in CAR T-cell expansion and efficacy is present in the data. To describe individual variability in CK, the following parameters were fit to individual CK trajectories: initial tumor burden, the number of divisions for activated T cells, and the fraction of effector cells that become memory cells. The time per T cell division and drug product and effector cell life spans were fit globally to all patient data. Optimization was performed using a Python-based trust region optimization method. Additionally, the percentage of B cells out of total cells in the model was calibrated to B cell aplasia data by tuning the number of endogenous lymphocytes in the model within a small, biologically reasonable range such that the mean and range of model outputs captured the general trend observed in the data. The rate of tumor cell division was also tuned to match the observed rebound in B cell aplasia data.
2.4 Model simulation and analysis
The model was implemented and simulations were performed with Applied BioMath’s proprietary QSP modeling platform. Analysis and plotting were performed with Python version 3.11.8.
Global sensitivity analysis (GSA) was evaluated using two methods: Sobol indices estimated via the Fourier Amplitude Sensitivity Test (FAST), implemented using SALib (Herman and Usher, 2017; Iwanaga et al., 2022), and partial rank correlation coefficients (PRCC), implemented using Pingouin (Vallat, 2018).
In the GSA, model parameters for which we had individual data or fitted values (body weight, fraction of CD8+ CAR T cells, initial tumor burden, number of CAR T cell divisions, and fraction of memory cells) were varied across the full range of individual values. Where possible, published ranges for individual parameters were used. CD19 expression was varied from 1,500 to 16,825 receptors per cell based on a published range for mantle cell lymphoma (D’Arena et al., 2000). Tumor doubling time was varied from 24 h to 30 days, based on the range reported in Roesch et al. (2014). Binding affinity was varied from 0.32 to 14.3 nM based on the range of values for CD19 CARs reported in Jayaraman et al. (2020). Remaining model parameters were varied 2-fold up and down nominal values. All parameters were sampled from a log-uniform distribution within their respective ranges, with a sample size of 5,000. Simulations were initialized with a
To explore temporal dynamics and explore the impact of tumor characteristics, we performed one-at-a-time scans of tumor division time and CD19 expression per cell. The model was simulated for specific patients as well as for the full patient population using different values of these parameters, while keeping other model parameters fixed.
3 Results
3.1 CK fitting and PD benchmarking
CAR T-cell trajectories vary widely from patient to patient. Our model was developed and calibrated to capture the typical phases of CAR T-cell CK as well as the variability between patients through fitting a combination of patient-specific and global parameter values. Results of optimization of CAR T-cell concentration to clinical data are shown in Figure 2. Figure 2A shows the average trajectory and full range across all 13 patients and Figure 2C shows each fitted patient simulation and data. The model adequately describes the overall behavior of the data despite the significant variability between patients as well as within each patient data set. The full table of final parameter values can be found in the Supplementary Material.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g002.jpg)
Figure 2. Model calibration and benchmarking results. (A) Patient population simulations for CAR T-cell CK. Black line indicates average model fit and shaded region represent the full range of individual trajectories. Points represent data, with colors representing different patients. (B) Patient population simulations for B cell aplasia. (C) Individual patient CAR T-cell CK data and simulations. Each panel represents an individual patient, the ID of which is labeled at the top of each panel.
Endogenous lymphocyte concentration and tumor doubling time were hand-tuned to a small degree to match measurements of B cells as a percentage of total lymphocytes, a measure of the efficacy of the CAR T-cells. Due to challenges with digitization, B cell aplasia data from only 6 of the 13 patients were distinguishable and are shown in Figure 2. The average and range of model simulations for all patient parameterizations capture the general trend of the data well and spans the variability between patients.
3.2 Sensitivity analysis
To explore the impact of model parameters on Cmax and efficacy, we performed Global Sensitivity Analysis (GSA). Results for Cmax are shown in Figure 3A. Sobol indices (first order and total order) and PRCC values are shown for all model parameters that had a p-value less than 0.05 and ranked in the top ten parameters for at least one measure of global sensitivity. The number of CAR T-cell divisions upon activation contributes to more than 80% of the variability in Cmax, which is a far greater contribution than any of the other parameters. The next most influential parameters are tumor growth rate, initial tumor burden, and mAb-CD19 binding affinity, which drive expansion through CAR-antigen interactions. CAR T cell life spans and CD19 expression are also influential. The ordering of parameters is roughly consistent between first order Sobol index, total order Sobol index, and PRCC. However, total order Sobol indices are generally at least two-fold larger than first order Sobol indices, indicating that there are interactions between parameters.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g003.jpg)
Figure 3. Global sensitivity results (Sobol indices and PRCC) for Cmax and tumor burden at day 20 post-treatment. Only parameters with a p-value less than 0.05 and that rank in the top ten for at least one measure of sensitivity are shown.
GSA results for tumor burden at day 20 are shown in Figure 3B. The most influential parameters are the tumor division time, number of T cell divisions, binding affinity, initial tumor burden, and CD19 expression is also influential. Tumor- and binding-related parameters are comparably influential on efficacy as the number of T cell divisions. This is in contrast to the results for Cmax, where the number of T cell divisions was by far the most influential parameter. This indicates that while expansion and efficacy are often correlated, patient properties such as tumor growth rate, initial tumor burden, and CD19 expression are more important for driving efficacy than they are for driving expansion. This is because CAR-CD19 interactions are required for both expansion and tumor cell killing.
3.3 Effects of tumor properties on CAR T expansion and efficacy
To investigate potential mechanisms related to patient to patient variability in response, we evaluated the effects of B cell division time and CD19 expression on B cells. These two parameters, which were informed by literature and not varied in Figure 2C were shown to be influential parameters by the GSA. We first focus on two patient parameterizations (F0104 and F0110) which had distinct CK and tumor growth profiles. F0104 has a typical CK profile consisting of expansion, contraction, and persistence, paired with a clear reduction in tumor growth, while F0110 had continued tumor growth and less defined expansion and contraction phases. Figure 4 shows simulations of CAR T-cell concentration and tumor dynamics for these two patients, scanning over both parameters. Parameters are varied 10-fold up and down from nominal values to explore a wide range of system behaviors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g004.jpg)
Figure 4. Results of scanning key model parameters representing patient characteristics for two individual patient parameterization. Parameters were scanned up to
Scanning over B cell division time, shown in Figure 4A, revealed qualitatively different behavior between the two patients. For patient F0104, the model predicts that a B cell division time corresponds to a more gradual contraction, resulting in a greater concentration of CAR T-cells over time. The division time does not significantly impact the Cmax. However, for patient F0110, the slope of the contraction phase is relatively consistent across division times but the Cmax increases with faster division times. For both patients, the greater expansion of CAR T-cells is not sufficient to reverse tumor cell growth. A faster B cell division time results in more tumor growth regardless of CAR T-cell concentration for the parameter range scanned. Within the first 10–15 days, there is an acute reduction in tumor cells in response to initial CAR T-cell expansion for the fastest tumor cell division time, 1.6 days. However, this effect is transient and the faster B cell division time results in faster rebound of the tumor. For the slower tumor cell division times, CAR T expansion does reduce the tumor size; this combined with the generally slower tumor growth results in slow tumor growth in the longer term.
Figure 4B shows the results of varying CD19 expression on B cells. Higher CD19 expression leads to additional binding to CAR T-cells and subsequent activation, increasing CAR T-cell expansion. This looks different for each of the patient parameterizations; simulated CK for patient F0110 shows greater sensitivity to CD19 expression compared to that of patient F0104. For F0104, higher CD19 expression leads to faster expansion and faster contraction, causing a sharper peak in the CAR T CK. Higher CD19 expression also leads to greater long-term persistence of CAR T-cells. Patient F0110 exhibits greater expansion and persistence with varying levels of CD19 expression, with no evident contraction phase.
Examining the individually fit parameters for F0104 and F0110 sheds light on the unique behaviors of both the CAR T-cells and tumor cells between patient simulations. Patient F0104 has a smaller fraction of effector CAR T-cells that become memory cells, a larger initial tumor burden, and a slightly higher number of CAR T-cell divisions upon activation compared to patient F0110. This leads to greater expansion (and therefore greater efficacy) of the CAR T-cells for patient F0104, but potentially less persistence. For patient F0110, the lower expansion and greater memory cell formation leads to no clear contraction phase in the CK. The corresponding tumor growth curves show no impact of treatment except a small reduction in tumor growth rate at the highest receptor expression level scanned.
To assess the behavior of CAR T-cells and tumor growth on a population level, the parameter value for either B cell division time or CD19 expression was updated one at a time for each patient. These parameters were varied ranges described in literature: Roesch et al. (2014) report NHL doubling times from 24 h to 30 days, and D’Arena et al. (2000) report a standard deviation for CD19 expression of about 1/3 the mean value. The minimum value for CD19 expression reported was much lower (10-fold lower than the mean value), which we also include in the parameter scan. Figure 5 shows the mean and standard deviation across all patients for CK, tumor cell count, and tumor fold change. Overall, the same patterns described above in the patient-specific scans hold true: faster B cell division times yield more CAR T-cell expansion and greater tumor growth, and higher CD19 expression leads to more CAR T-cells and improved tumor cell killing. Within the physiological ranges tested, B cell division time has an impact on both CK and tumor cell growth by close to an order of magnitude, on average. Notably, the rate of tumor regrowth is similar for all tumor doubling times, indicating that the increased persistence of CAR T-cells does counteract the increased tumor growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g005.jpg)
Figure 5. Results of scanning key model parameters representing patient characteristics across all patient parameterizations. Parameters were scanned across ranges consistent with values reported in the literature. Simulations of CAR T-cell concentration, total B cells, and B cell fold change from initial are shown for scans of (A) B cell time per division and (B) CD19 receptors per cell. Mean and one standard deviation of all patient results are shown.
The range of reported CD19 expression is quite varied, and the model predicts that this parameter could have a significant impact on treatment efficacy. Between the maximum and minimum values scanned, within the range of reported values, there is about an order of magnitude difference in the CAR T cell Cmax. Furthermore, for the lowest CD19 RPC, there is essentially no tumor growth inhibition. The three higher RPC values do show inhibition, with a reduction from baseline of up to 10x.
3.4 Exploratory analysis: memory cell killing
To show how the model can explore questions about both individual and population-level dynamics, we performed simulations to understand the potential impact of memory cell killing. In the nominal simulations, we assume that memory cells do not kill tumor cells. For this analysis, we compare the nominal simulations against those in which memory cell killing has the same killing capacity as activated cells. Figure 6 shows the results for both the population level and individual trajectories.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1380018/fsysb-04-1380018-HTML/image_m/fsysb-04-1380018-g006.jpg)
Figure 6. Results of exploring the impact of memory cell killing. Simulations in which memory cells have the same killing capacity as active cells are compared to the nominal case in which memory cells do not kill tumor cells. Simulations of CAR T-cell concentration, total B cells, and B cell fold change from initial are shown in (A) where mean and one standard deviation of all patient results are compared. B cell fold change from initial is shown in (B); multicolored lines represent individual patient trajectories, black curve represents the mean across all patients, and the gray region represents one standard deviation.
First observing the population-level dynamics in Figure 6A, the impact on memory cell killing is observed only in the tumor, not in CAR T CK. Furthermore, the model predicts that any difference is observed after 50 days. This makes sense due to the delayed appearance of memory cells and the subsequent growth of the memory cell population - thus, memory cell killing is predicted to have a small overall impact on reducing tumor growth during the terminal phase of CAR T expansion. Since the exact memory cell populations may vary from patient to patient, the impact of memory cell killing may also be observed on a patient level, shown in Figure 6B. It is evident that some patients show little impact of memory cell killing, sch as F0125 and F0126. On the other hand, tumor growth in patients F0107, F0111, F0123 nearly plateaus as compared to the nominal parameterization which has linear growth.
Although there is insufficient data to inform the true activity of memory CAR-T cells in the model, this hypothetical analysis shows the ability of the model to differentiate the impact of treatment on individual patients as compared to a population-level impact. In these simulations, a moderate population-level effect was the result of an aggregated variety of patient effects, from no impact to a signficiant impact. Furthermore, the model shows in what populations and at what times the impact of these changes might be observed.
4 Discussion
Our mechanistic modeling approach incorporating molecular-scale and cell-scale dynamics successfully captured CAR T CK-PD and revealed key system behaviors. Mechanistic modeling is necessary to capture the interplay of target engagement, T cell expansion, and tumor cell killing. CAR T-cell therapy is distinct from other therapeutics in that CK and PD are inter-dependent. This dependency is demonstrated in our model by the sensitivity of Cmax to tumor and binding parameters.
Global sensitivity analysis revealed that both drug-specific and patient-specific properties can potentially explain variability in response to CAR T therapy. The most influential drug-specific properties are the number of divisions per activated CAR T-cell and the binding of the CAR for CD19. The number of divisions for activated cells is the most influential factor for peak CAR T-cell expansion, and was also highly influential for tumor killing. This number of divisions could potentially be increased through further engineering or refining of manufacturing processes for the CAR T product, for example, through selection for naive cells (Arcangeli et al., 2022). Importantly, while we classify this as a drug-specific property, this could be variable across patients since the CAR T-cells are manufactured from the patients’ own cells. Thus, individual variability in the number of T cell divisions can also contribute to observed variability in CK and efficacy. Binding affinity also impacted both CAR T-cell expansion and efficacy, which could be improved in engineering of the CAR.
While our modeling suggests that CAR T expansion is driven primarily by number of divisions, global sensitivity analysis shows that tumor properties such as CD19 expression and growth rate are comparatively more influential in driving efficacy. Tumor growth rate was also highly influential on CK. This demonstrates two things: (1) while expansion often correlates with efficacy, expansion itself is not necessarily sufficient for tumor shrinkage, and (2) variability in patient characteristics will lead to significant variability in both exposure and response. Modeling provides insight into this variability and can be used to inform patient, target, and indication selection.
In individual- and population-level model simulations, we observed that although a faster tumor growth rate corresponds to increased CAR T-cell expansion and distinct CK profiles, this increased expansion is often not enough to control the faster-growing tumor. This implies that drug characteristics may need to be modified in order to target more aggressive tumors. Notably, while there is little to no predicted tumor shrinkage for this CAR T with faster growing tumors for most patients, treatment is still effective in slowing tumor growth both short- and long-term, providing a benefit to patients. Increased target expression drives both increased expansion and stronger tumor killing. Patients with low target expression may be poor candidates for this type of treatment due to poor expansion and little anti-tumor activity, leading to lack of response. This also suggests that target expression should be a key consideration in both target and indication selection, while balancing toxicity concerns. Furthermore, the power of individualized parameterizations of the model was demonstrated in the memory cell killing exploratory analysis. Although the population-level simulation showed a small overall reduction in tumor growth, some patient trajectories showed signficant reduction while others showed nearly no impact. Although this was a hypothetical exploration due insufficient data, these simulations demonstrate that modeling can have a large impact on understanding individual patient dynamics.
In the future, this model and analysis could help drive decisions in CAR T-cell design, manufacturing, patient selection, patient-specific dose selection, and efficacious dose selection for novel CAR Ts. This model could be further refined by adding other T cell phenotypes, cytokines, immune cells types, and additional reactions such as re-activation of memory cells. One key limitation of the current work is lack of direct measurements of tumor burden over time to inform efficacy. Rather, we relied of B cell aplasia data and assumptions about the native immune population to estimate tumor reduction. Additional efficacy data would help to better constrain the model and may allow for individualized efficacy modeling. Additional patient-specific data such as CD19 expression could enable individualized predictions of efficacy/response through a digital twin approach. Another limitation of this model is that it does not account for effects of CD4+ T cells on tumor cell killing. This model could also be extended to study other targets and indications, including solid tumors for which there are currently no approved CAR T-cell therapies.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7518381/.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
SM: Conceptualization, Formal Analysis, Visualization, Writing–original draft. SG: Data curation, Investigation, Writing–original draft. KS: Conceptualization, Supervision, Writing–review and editing. MR: Formal Analysis, Methodology, Supervision, Visualization, Writing–original draft.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We thank Raibatak Das for developing the first iteration of this model.
Conflict of interest
Authors SM, SG, MR are employed by Applied BioMath. Author KS is employed by Differentia Bio.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fsysb.2024.1380018/full#supplementary-material
References
Alizadeh, D., Wang, D., and Brown, C. E. (2023). Uncovering the role of cd4+ car t cells in cancer immunotherapy. Cancer Res. 83, 2813–2815. doi:10.1158/0008-5472.CAN-23-1948
Anikeeva, N., Panteleev, S., Mazzanti, N., Terai, M., Sato, T., and Sykulev, Y. (2021). Efficient killing of tumor cells by CAR-T cells requires greater number of engaged CARs than TCRs. J. Biol. Chem. 297, 101033. doi:10.1016/j.jbc.2021.101033
Arcangeli, S., Bove, C., Mezzanotte, C., Camisa, B., Falcone, L., Manfredi, F., et al. (2022). Car t cell manufacturing from naive/stem memory t lymphocytes enhances antitumor responses while curtailing cytokine release syndrome. J. Clin. investigation 132, e150807. doi:10.1172/JCI150807
Borghans, J. A., Tesselaar, K., and de Boer, R. J. (2018). Current best estimates for the average lifespans of mouse and human leukocytes: reviewing two decades of deuterium-labeling experiments. Immunol. Rev. 285, 233–248. doi:10.1111/imr.12693
Chaudhury, A., Zhu, X., Chu, L., Goliaei, A., June, C. H., Kearns, J. D., et al. (2020). Chimeric antigen receptor T cell therapies: a review of cellular kinetic-pharmacodynamic modeling approaches. J. Clin. Pharmacol. 60, S147-S159–S159. doi:10.1002/jcph.1691
Chen, Y.-J., Abila, B., and Mostafa Kamel, Y. (2023). CAR-T: what is next? Cancers (Basel) 15, 663. doi:10.3390/cancers15030663
Cui, W., and Kaech, S. M. (2010). Generation of effector CD8+ T cells and their conversion to memory T cells. Immunol. Rev. 236, 151–166. doi:10.1111/j.1600-065X.2010.00926.x
D’Arena, G., Musto, P., Cascavilla, N., Dell’Olio, M., Di Renzo, N., and Carotenuto, M. (2000). Quantitative flow cytometry for the differential diagnosis of leukemic B-cell chronic lymphoproliferative disorders. Am. J. Hematol. 64, 275–281. doi:10.1002/1096-8652(200008)64:4<275::aid-ajh7>3.0.co;2-y
Das, D. K., Gupta, S. K., Datta, B. N., and Sharma, S. C. (1991). Working formulation of the non-Hodgkin’s lymphomas: a study of cell size and mitotic indices in cytologic subtypes. Diagn. Cytopathol. 7, 499–503. doi:10.1002/dc.2840070511
Du, X., Beers, R., FitzGerald, D. J., and Pastan, I. (2008). Differential cellular internalization of anti-CD19 and -CD22 immunotoxins results in different cytotoxic activity. Cancer Res. 68, 6300–6305. doi:10.1158/0008-5472.CAN-08-0461
Hakim, F. T., Memon, S. A., Cepeda, R., Jones, E. C., Chow, C. K., Kasten-Sportes, C., et al. (2005). Age-dependent incidence, time course, and consequences of thymic renewal in adults. J. Clin. investigation 115, 930–939. doi:10.1172/JCI22492
Hardiansyah, D., and Ng, C. M. (2019). Quantitative systems pharmacology model of chimeric antigen receptor T-cell therapy. Clin. Transl. Sci. 12, 343–349. doi:10.1111/cts.12636
Henrickson, S. E., Mempel, T. R., Mazo, I. B., Liu, B., Artyomov, M. N., Zheng, H., et al. (2008). T cell sensing of antigen dose governs interactive behavior with dendritic cells and sets a threshold for T cell activation. Nat. Immunol. 9, 282–291. doi:10.1038/ni1559
Herman, J., and Usher, W. (2017). SALib: an open-source python library for sensitivity analysis. J. Open Source Softw. 2, 97. doi:10.21105/joss.00097
Iwanaga, T., Usher, W., and Herman, J. (2022). Toward SALib 2.0: advancing the accessibility and interpretability of global sensitivity analyses. Socio-Environmental Syst. Model. 4, 18155. doi:10.18174/sesmo.18155
Jayaraman, J., Mellody, M. P., Hou, A. J., Desai, R. P., Fung, A. W., Pham, A. H. T., et al. (2020). CAR-T design: elements and their synergistic function. EBioMedicine 58, 102931. doi:10.1016/j.ebiom.2020.102931
Kaech, S. M., Wherry, E. J., and Ahmed, R. (2002). Effector and memory T-cell differentiation: implications for vaccine development. Nat. Rev. Immunol. 2, 251–262. doi:10.1038/nri778
Li, W., Qiu, S., Chen, J., Jiang, S., Chen, W., Jiang, J., et al. (2020). Chimeric antigen receptor designed to prevent ubiquitination and downregulation showed durable antitumor efficacy. Immunity 53, 456–470. doi:10.1016/j.immuni.2020.07.011
Malik-Chaudhry, H. K., Prabhakar, K., Ugamraj, H. S., Boudreau, A. A., Buelow, B., Dang, K., et al. (2021). TNB-486 induces potent tumor cell cytotoxicity coupled with low cytokine release in preclinical models of B-NHL. MAbs. 13 (1), e1890411-1–e1890411-15. doi:10.1080/19420862.2021.1890411
Melenhorst, J. J., Chen, G. M., Wang, M., Porter, D. L., Chen, C., Collins, M. A., et al. (2022). Decade-long leukaemia remissions with persistence of CD4+ CAR T cells. Nature 602, 503–509. doi:10.1038/s41586-021-04390-6
Nukala, U., Rodriguez Messan, M., Yogurtcu, O. N., Wang, X., and Yang, H. (2021). A systematic review of the efforts and hindrances of modeling and simulation of CAR T-cell therapy. AAPS J. 23, 52–20. doi:10.1208/s12248-021-00579-9
Owens, K., and Bozic, I. (2021). Modeling CAR T-cell therapy with patient preconditioning. Bull. Math. Biol. 83, 42–36. doi:10.1007/s11538-021-00869-5
Press, O. W., Eary, J. F., Appelbaum, F. R., Martin, P. J., Badger, C. C., Nelp, W. B., et al. (1993). Radiolabeled-antibody therapy of B-cell lymphoma with autologous bone marrow support. N. Engl. J. Med. 329, 1219–1224. doi:10.1056/NEJM199310213291702
Roesch, K., Hasenclever, D., and Scholz, M. (2014). Modelling lymphoma therapy and outcome. Bull. Math. Biol. 76, 401–430. doi:10.1007/s11538-013-9925-3
Salem, A. M., Mugundu, G. M., and Singh, A. P. (2023). Development of a multiscale mechanistic modeling framework integrating differential cellular kinetics of car t-cell subsets and immunophenotypes in cancer patients. CPT Pharmacometrics and Syst. Pharmacol. 12, 1285–1304. doi:10.1002/psp4.13009
Sermer, D., Batlevi, C., Palomba, M. L., Shah, G., Lin, R. J., Perales, M.-A., et al. (2020). Outcomes in patients with DLBCL treated with commercial CAR T cells compared with alternate therapies. Blood Adv. 4, 4669–4678. doi:10.1182/bloodadvances.2020002118
Sieber, T., Schoeler, D., Ringel, F., Pascu, M., and Schriever, F. (2003). Selective internalization of monoclonal antibodies by B-cell chronic lymphocytic leukaemia cells. Br. J. Haematol. 121, 458–461. doi:10.1046/j.1365-2141.2003.04305.x
Singh, A. P., Zheng, X., Lin-Schmidt, X., Chen, W., Carpenter, T. J., Zong, A., et al. (2020). Development of a quantitative relationship between CAR-affinity, antigen abundance, tumor cell depletion and CAR-T cell expansion using a multiscale systems PK-PD model. MAbs. 12 (1), e1688616-1–e1688616-21. doi:10.1080/19420862.2019.1688616
Spiegel, J. Y., Patel, S., Muffly, L., Hossain, N. M., Oak, J., Baird, J. H., et al. (2021). Car t cells with dual targeting of cd19 and cd22 in adult patients with recurrent or refractory b cell malignancies: a phase 1 trial. Nat. Med. 27, 1419–1431. doi:10.1038/s41591-021-01436-0
Stein, A. M., Grupp, S. A., Levine, J. E., Laetsch, T. W., Pulsipher, M. A., Boyer, M. W., et al. (2019). Tisagenlecleucel model-based cellular kinetic analysis of chimeric antigen receptor–T cells. CPT pharmacometrics and Syst. Pharmacol. 8, 285–295. doi:10.1002/psp4.12388
Vallat, R. (2018). Pingouin: statistics in python. J. Open Source Softw. 3, 1026. doi:10.21105/joss.01026
Ying, Z., He, T., Wang, X., Zheng, W., Lin, N., Tu, M., et al. (2021). Distribution of chimeric antigen receptor-modified T cells against CD19 in B-cell malignancies. BMC cancer 21, 198–207. doi:10.1186/s12885-021-07934-1
Keywords: CAR T-cell therapy, mechanistic modeling, quantitative systems pharmacology, sensitivity analysis, cellular kinetics
Citation: Minucci S, Gruver S, Subramanian K and Renardy M (2024) A multi-scale semi-mechanistic CK/PD model for CAR T-cell therapy. Front. Syst. Biol. 4:1380018. doi: 10.3389/fsysb.2024.1380018
Received: 31 January 2024; Accepted: 15 August 2024;
Published: 29 August 2024.
Edited by:
Shayn Peirce-Cottler, University of Virginia, United StatesReviewed by:
Elisa Domínguez-Hüttinger, National Autonomous University of Mexico, MexicoRada Amin, University of Nebraska-Lincoln, United States
Copyright © 2024 Minucci, Gruver, Subramanian and Renardy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sarah Minucci, sarah.minucci@appliedbiomath.com