1
|
Nikitich A, Helmlinger G, Peskov K, Bocharov G. Mathematical modeling of endogenous and exogenously administered T cell recirculation in mouse and its application to pharmacokinetic studies of cell therapies. Front Immunol 2024; 15:1357706. [PMID: 38846946 PMCID: PMC11155669 DOI: 10.3389/fimmu.2024.1357706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 03/19/2024] [Indexed: 06/09/2024] Open
Abstract
Introduction In vivo T cell migration has been of interest to scientists for the past 60 years. T cell kinetics are important in the understanding of the immune response to infectious agents. More recently, adoptive T cell therapies have proven to be a most promising approach to treating a wide range of diseases, including autoimmune and cancer diseases, whereby the characterization of cellular kinetics represents an important step towards the prediction of therapeutic efficacy. Methods Here, we developed a physiologically-based pharmacokinetic (PBPK) model that describes endogenous T cell homeostasis and the kinetics of exogenously administered T cells in mouse. Parameter calibration was performed using a nonlinear fixed-effects modeling approach based on published data on T cell kinetics and steady-state levels in different tissues of mice. The Partial Rank Correlation Coefficient (PRCC) method was used to perform a global sensitivity assessment. To estimate the impact of kinetic parameters on exogenously administered T cell dynamics, a local sensitivity analysis was conducted. Results We simulated the model to analyze cellular kinetics following various T cell doses and frequencies of CCR7+ T cells in the population of infused lymphocytes. The model predicted the effects of T cell numbers and of population composition of infused T cells on the resultant concentration of T cells in various organs. For example, a higher percentage of CCR7+ T cells among exogenously administered T lymphocytes led to an augmented accumulation of T cells in the spleen. The model predicted a linear dependence of T cell dynamics on the dose of adoptively transferred T cells. Discussion The mathematical model of T cell migration presented here can be integrated into a multi-scale model of the immune system and be used in a preclinical setting for predicting the distribution of genetically modified T lymphocytes in various organs, following adoptive T cell therapies.
Collapse
Affiliation(s)
- Antonina Nikitich
- Research Center of Model-Informed Drug Development, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS), Moscow, Russia
| | | | - Kirill Peskov
- Research Center of Model-Informed Drug Development, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS), Moscow, Russia
- Modeling and Simulation Decisions FZ - LLC, Dubai, United Arab Emirates
- University of Science and Technology (STU) “Sirius”, Sochi, Russia
| | - Gennady Bocharov
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS), Moscow, Russia
- Institute for Computer Science and Mathematical Modelling, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
- Moscow Center of Fundamental and Applied Mathematics at INM RAS, Moscow, Russia
| |
Collapse
|
2
|
Kirouac DC, Zmurchok C, Morris D. Making drugs from T cells: The quantitative pharmacology of engineered T cell therapeutics. NPJ Syst Biol Appl 2024; 10:31. [PMID: 38499572 PMCID: PMC10948391 DOI: 10.1038/s41540-024-00355-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 02/26/2024] [Indexed: 03/20/2024] Open
Abstract
Engineered T cells have emerged as highly effective treatments for hematological cancers. Hundreds of clinical programs are underway in efforts to expand the efficacy, safety, and applications of this immuno-therapeutic modality. A primary challenge in developing these "living drugs" is the complexity of their pharmacology, as the drug product proliferates, differentiates, traffics between tissues, and evolves through interactions with patient immune systems. Using publicly available clinical data from Chimeric Antigen Receptor (CAR) T cells, we demonstrate how mathematical models can be used to quantify the relationships between product characteristics, patient physiology, pharmacokinetics and clinical outcomes. As scientists work to develop next-generation cell therapy products, mathematical models will be integral for contextualizing data and facilitating the translation of product designs to clinical strategy.
Collapse
Affiliation(s)
- Daniel C Kirouac
- Notch Therapeutics, Vancouver, BC, Canada.
- The University of British Columbia, School of Biomedical Engineering, Vancouver, BC, Canada.
- Metrum Research Group, Tariffville, CT, USA.
| | | | | |
Collapse
|
3
|
Liao YM, Hsu SH, Chiou SS. Harnessing the Transcriptional Signatures of CAR-T-Cells and Leukemia/Lymphoma Using Single-Cell Sequencing Technologies. Int J Mol Sci 2024; 25:2416. [PMID: 38397092 PMCID: PMC10889174 DOI: 10.3390/ijms25042416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 02/02/2024] [Accepted: 02/10/2024] [Indexed: 02/25/2024] Open
Abstract
Chimeric antigen receptor (CAR)-T-cell therapy has greatly improved outcomes for patients with relapsed or refractory hematological malignancies. However, challenges such as treatment resistance, relapse, and severe toxicity still hinder its widespread clinical application. Traditional transcriptome analysis has provided limited insights into the complex transcriptional landscape of both leukemia cells and engineered CAR-T-cells, as well as their interactions within the tumor microenvironment. However, with the advent of single-cell sequencing techniques, a paradigm shift has occurred, providing robust tools to unravel the complexities of these factors. These techniques enable an unbiased analysis of cellular heterogeneity and molecular patterns. These insights are invaluable for precise receptor design, guiding gene-based T-cell modification, and optimizing manufacturing conditions. Consequently, this review utilizes modern single-cell sequencing techniques to clarify the transcriptional intricacies of leukemia cells and CAR-Ts. The aim of this manuscript is to discuss the potential mechanisms that contribute to the clinical failures of CAR-T immunotherapy. We examine the biological characteristics of CAR-Ts, the mechanisms that govern clinical responses, and the intricacies of adverse events. By exploring these aspects, we hope to gain a deeper understanding of CAR-T therapy, which will ultimately lead to improved clinical outcomes and broader therapeutic applications.
Collapse
Affiliation(s)
- Yu-Mei Liao
- Division of Hematology-Oncology, Department of Pediatrics, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan;
| | - Shih-Hsien Hsu
- Graduate Institute of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
- Department of Medical Research, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan
- Research Center for Environmental Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
- Center of Applied Genomics, Kaohsiung Medical University, Kaohsiung 807, Taiwan
| | - Shyh-Shin Chiou
- Division of Hematology-Oncology, Department of Pediatrics, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 807, Taiwan;
- Research Center for Environmental Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
- Center of Applied Genomics, Kaohsiung Medical University, Kaohsiung 807, Taiwan
- Graduate Institute of Clinical Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
| |
Collapse
|
4
|
Joslyn LR, Huang W, Miles D, Hosseini I, Ramanujan S. "Digital twins elucidate critical role of T scm in clinical persistence of TCR-engineered cell therapy". NPJ Syst Biol Appl 2024; 10:11. [PMID: 38278838 PMCID: PMC10817974 DOI: 10.1038/s41540-024-00335-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 01/11/2024] [Indexed: 01/28/2024] Open
Abstract
Despite recent progress in adoptive T cell therapy for cancer, understanding and predicting the kinetics of infused T cells remains a challenge. Multiple factors can impact the distribution, expansion, and decay or persistence of infused T cells in patients. We have developed a novel quantitative systems pharmacology (QSP) model of TCR-transgenic T cell therapy in patients with solid tumors to describe the kinetics of endogenous T cells and multiple memory subsets of engineered T cells after infusion. These T cells undergo lymphodepletion, proliferation, trafficking, differentiation, and apoptosis in blood, lymph nodes, tumor site, and other peripheral tissues. Using the model, we generated patient-matched digital twins that recapitulate the circulating T cell kinetics reported from a clinical trial of TCR-engineered T cells targeting E7 in patients with metastatic HPV-associated epithelial cancers. Analyses of key parameters influencing cell kinetics and differences among digital twins identify stem cell-like memory T cells (Tscm) cells as an important determinant of both expansion and persistence and suggest that Tscm-related differences contribute significantly to the observed variability in cellular kinetics among patients. We simulated in silico clinical trials using digital twins and predict that Tscm enrichment in the infused product improves persistence of the engineered T cells and could enable administration of a lower dose. Finally, we verified the broader relevance of the QSP model, the digital twins, and findings on the importance of Tscm enrichment by predicting kinetics for two patients with pancreatic cancer treated with KRAS G12D targeting T cell therapy. This work offers insight into the key role of Tscm biology on T cell kinetics and provides a quantitative framework to evaluate cellular kinetics for future efforts in the development and clinical application of TCR-engineered T cell therapies.
Collapse
Affiliation(s)
| | - Weize Huang
- Genentech Inc., South San Francisco, CA, USA
| | - Dale Miles
- Genentech Inc., South San Francisco, CA, USA
| | | | | |
Collapse
|
5
|
Arabameri A, Arab S. Understanding the Interplay of CAR-NK Cells and Triple-Negative Breast Cancer: Insights from Computational Modeling. Bull Math Biol 2024; 86:20. [PMID: 38240892 DOI: 10.1007/s11538-023-01247-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2023] [Accepted: 12/14/2023] [Indexed: 01/23/2024]
Abstract
Chimeric antigen receptor (CAR)-engineered natural killer (NK) cells have recently emerged as a promising and safe alternative to CAR-T cells for targeting solid tumors. In the case of triple-negative breast cancer (TNBC), traditional cancer treatments and common immunotherapies have shown limited effectiveness. However, CAR-NK cells have been successfully employed to target epidermal growth factor receptor (EGFR) on TNBC cells, thereby enhancing the efficacy of immunotherapy. The effectiveness of CAR-NK-based immunotherapy is influenced by various factors, including the vaccination dose, vaccination pattern, and tumor immunosuppressive factors in the microenvironment. To gain insights into the dynamics and effects of CAR-NK-based immunotherapy, we propose a computational model based on experimental data and immunological theories. This model integrates an individual-based model that describes the interplay between the tumor and the immune system, along with an ordinary differential equation model that captures the variation of inflammatory cytokines. Computational results obtained from the proposed model shed light on the conditions necessary for initiating an effective anti-tumor response. Furthermore, global sensitivity analysis highlights the issue of low persistence of CAR-NK cells in vivo, which poses a significant challenge for the successful clinical application of these cells. Leveraging the model, we identify the optimal vaccination time, vaccination dose, and time interval between injections for maximizing therapeutic outcomes.
Collapse
Affiliation(s)
- Abazar Arabameri
- Department of Electrical Engineering, University of Zanjan, Zanjan, Iran.
| | - Samaneh Arab
- Department of Tissue Engineering and Applied Cell Sciences, School of Medicine, Semnan University of Medical Sciences, Semnan, Iran
| |
Collapse
|
6
|
Kirouac DC, Zmurchok C, Deyati A, Sicherman J, Bond C, Zandstra PW. Deconvolution of clinical variance in CAR-T cell pharmacology and response. Nat Biotechnol 2023; 41:1606-1617. [PMID: 36849828 PMCID: PMC10635825 DOI: 10.1038/s41587-023-01687-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Accepted: 01/20/2023] [Indexed: 03/01/2023]
Abstract
Chimeric antigen receptor T cell (CAR-T) expansion and persistence vary widely among patients and predict both efficacy and toxicity. However, the mechanisms underlying clinical outcomes and patient variability are poorly defined. In this study, we developed a mathematical description of T cell responses wherein transitions among memory, effector and exhausted T cell states are coordinately regulated by tumor antigen engagement. The model is trained using clinical data from CAR-T products in different hematological malignancies and identifies cell-intrinsic differences in the turnover rate of memory cells and cytotoxic potency of effectors as the primary determinants of clinical response. Using a machine learning workflow, we demonstrate that product-intrinsic differences can accurately predict patient outcomes based on pre-infusion transcriptomes, and additional pharmacological variance arises from cellular interactions with patient tumors. We found that transcriptional signatures outperform T cell immunophenotyping as predictive of clinical response for two CD19-targeted CAR-T products in three indications, enabling a new phase of predictive CAR-T product development.
Collapse
Affiliation(s)
| | | | | | | | - Chris Bond
- Notch Therapeutics, Vancouver, BC, Canada
| | - Peter W Zandstra
- Notch Therapeutics, Vancouver, BC, Canada
- School of Biomedical Engineering and Michael Smith Laboratories, University of British Columbia, Vancouver, BC, Canada
| |
Collapse
|
7
|
Rajakaruna H, Desai M, Das J. PASCAR: a multiscale framework to explore the design space of constitutive and inducible CAR T cells. Life Sci Alliance 2023; 6:e202302171. [PMID: 37507138 PMCID: PMC10387492 DOI: 10.26508/lsa.202302171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 07/08/2023] [Accepted: 07/12/2023] [Indexed: 07/30/2023] Open
Abstract
CAR T cells are engineered to bind and destroy tumor cells by targeting overexpressed surface antigens. However, healthy cells expressing lower abundances of these antigens can also be lysed by CAR T cells. Various CAR T cell designs increase tumor cell elimination, whereas reducing damage to healthy cells. However, these efforts are costly and labor-intensive, constraining systematic exploration of potential hypotheses. We develop a protein abundance structured population dynamic model for CAR T cells (PASCAR), a framework that combines multiscale population dynamic models and multi-objective optimization approaches with data from cytometry and cytotoxicity assays to systematically explore the design space of constitutive and tunable CAR T cells. PASCAR can quantitatively describe in vitro and in vivo results for constitutive and inducible CAR T cells and can successfully predict experiments outside the training data. Our exploration of the CAR design space reveals that optimal CAR affinities in the intermediate range of dissociation constants effectively reduce healthy cell lysis, whereas maintaining high tumor cell-killing rates. Furthermore, our modeling offers guidance for optimizing CAR expressions in synthetic notch CAR T cells. PASCAR can be extended to other CAR immune cells.
Collapse
Affiliation(s)
- Harshana Rajakaruna
- The Steve and Cindy Rasmussen Institute for Genomics, The Abigail Wexner Research Institute, Nationwide Children's Hospital, Columbus, OH, USA
| | - Milie Desai
- Department of Biology, Indian Institute of Science Education and Research, Pune, India
| | - Jayajit Das
- The Steve and Cindy Rasmussen Institute for Genomics, The Abigail Wexner Research Institute, Nationwide Children's Hospital, Columbus, OH, USA
- Department of Pediatrics and Pelotonia Institute for Immuno-Oncology, College of Medicine, Columbus, OH, USA
- Biophysics Program, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
8
|
Salem AM, Mugundu GM, Singh AP. Development of a multiscale mechanistic modeling framework integrating differential cellular kinetics of CAR T-cell subsets and immunophenotypes in cancer patients. CPT Pharmacometrics Syst Pharmacol 2023; 12:1285-1304. [PMID: 37448297 PMCID: PMC10508581 DOI: 10.1002/psp4.13009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 06/26/2023] [Accepted: 06/28/2023] [Indexed: 07/15/2023] Open
Abstract
Chimeric antigen receptor (CAR) T-cell subsets and immunophenotypic composition of the pre-infusion product, as well as their longitudinal changes following infusion, are expected to affect CAR T-cell expansion, persistence, and clinical outcomes. Herein, we sequentially evolved our previously described cellular kinetic-pharmacodynamic (CK-PD) model to incorporate CAR T-cell product-associated attributes by utilizing published preclinical and clinical datasets from two affinity variants (FMC63 and CAT19 scFv) anti-CD19 CAR T-cells. In step 1, a unified cell-level PD model was used to simultaneously characterize the in vitro killing datasets of two CAR T-cells against CD19+ cell lines at varying effector:target ratios. In step 2, an augmented CK-PD model for anti-CD19 CAR T-cells was developed, by integrating CK dataset(s) from two bioanalytical measurements (quantitative polymerase chain reaction and flow cytometry) in patients with cancer. The model described the differential in vivo expansion properties of CAR T-cell subsets. The estimated expansion rate constant was ~1.12-fold higher for CAR+CD8+ cells in comparison to CAR+CD4+ T-cells. In step 3, the model was extended to characterize the disposition of four immunophenotypic populations of CAR T-cells, including stem-cell memory, central memory, effector memory, and effector cells. The model adequately characterized the longitudinal changes in immunophenotypes post anti-CD19 CAR T-cell infusion in pediatric patients with acute lymphocytic leukemia. Polyclonality in the pre-infusion product was identified as a categorical covariate influencing differentiation of immunophenotypes. In the future, this model could be leveraged a priori toward optimizing the composition of CAR T-cell infusion product, and further understand the CK-PD relationship in patients.
Collapse
Affiliation(s)
- Ahmed M. Salem
- Clinical Pharmacology and Modeling, Precision and Translational MedicineOncology Cell Therapy and Therapeutic Area Unit, Takeda PharmaceuticalsCambridgeMassachusettsUSA
- Center for Translational MedicineUniversity of Maryland School of PharmacyBaltimoreMarylandUSA
| | - Ganesh M. Mugundu
- Clinical Pharmacology and Modeling, Precision and Translational MedicineOncology Cell Therapy and Therapeutic Area Unit, Takeda PharmaceuticalsCambridgeMassachusettsUSA
| | - Aman P. Singh
- Clinical Pharmacology and Modeling, Precision and Translational MedicineOncology Cell Therapy and Therapeutic Area Unit, Takeda PharmaceuticalsCambridgeMassachusettsUSA
| |
Collapse
|
9
|
Brummer AB, Xella A, Woodall R, Adhikarla V, Cho H, Gutova M, Brown CE, Rockne RC. Data driven model discovery and interpretation for CAR T-cell killing using sparse identification and latent variables. Front Immunol 2023; 14:1115536. [PMID: 37256133 PMCID: PMC10226275 DOI: 10.3389/fimmu.2023.1115536] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2022] [Accepted: 03/27/2023] [Indexed: 06/01/2023] Open
Abstract
In the development of cell-based cancer therapies, quantitative mathematical models of cellular interactions are instrumental in understanding treatment efficacy. Efforts to validate and interpret mathematical models of cancer cell growth and death hinge first on proposing a precise mathematical model, then analyzing experimental data in the context of the chosen model. In this work, we present the first application of the sparse identification of non-linear dynamics (SINDy) algorithm to a real biological system in order discover cell-cell interaction dynamics in in vitro experimental data, using chimeric antigen receptor (CAR) T-cells and patient-derived glioblastoma cells. By combining the techniques of latent variable analysis and SINDy, we infer key aspects of the interaction dynamics of CAR T-cell populations and cancer. Importantly, we show how the model terms can be interpreted biologically in relation to different CAR T-cell functional responses, single or double CAR T-cell-cancer cell binding models, and density-dependent growth dynamics in either of the CAR T-cell or cancer cell populations. We show how this data-driven model-discovery based approach provides unique insight into CAR T-cell dynamics when compared to an established model-first approach. These results demonstrate the potential for SINDy to improve the implementation and efficacy of CAR T-cell therapy in the clinic through an improved understanding of CAR T-cell dynamics.
Collapse
Affiliation(s)
- Alexander B. Brummer
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
- Department of Physics and Astronomy, College of Charleston, Charleston, SC, United States
| | - Agata Xella
- Department of Hemtaology and Hematopoietic Cell Translation and Immuno-Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| | - Ryan Woodall
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| | - Vikram Adhikarla
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| | - Heyrim Cho
- Department of Mathematics, University of California, Riverside, Riverside, CA, United States
| | - Margarita Gutova
- Department of Stem Cell Biology and Regenerative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| | - Christine E. Brown
- Department of Hemtaology and Hematopoietic Cell Translation and Immuno-Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| | - Russell C. Rockne
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA, United States
| |
Collapse
|
10
|
Paixão EA, Barros LRC, Fassoni AC, Almeida RC. Modeling Patient-Specific CAR-T Cell Dynamics: Multiphasic Kinetics via Phenotypic Differentiation. Cancers (Basel) 2022; 14:cancers14225576. [PMID: 36428671 PMCID: PMC9688514 DOI: 10.3390/cancers14225576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 11/03/2022] [Accepted: 11/07/2022] [Indexed: 11/16/2022] Open
Abstract
Chimeric Antigen Receptor (CAR)-T cell immunotherapy revolutionized cancer treatment and consists of the genetic modification of T lymphocytes with a CAR gene, aiming to increase their ability to recognize and kill antigen-specific tumor cells. The dynamics of CAR-T cell responses in patients present multiphasic kinetics with distribution, expansion, contraction, and persistence phases. The characteristics and duration of each phase depend on the tumor type, the infused product, and patient-specific characteristics. We present a mathematical model that describes the multiphasic CAR-T cell dynamics resulting from the interplay between CAR-T and tumor cells, considering patient and product heterogeneities. The CAR-T cell population is divided into functional (distributed and effector), memory, and exhausted CAR-T cell phenotypes. The model is able to describe the diversity of CAR-T cell dynamical behaviors in different patients and hematological cancers as well as their therapy outcomes. Our results indicate that the joint assessment of the area under the concentration-time curve in the first 28 days and the corresponding fraction of non-exhausted CAR-T cells may be considered a potential marker to classify therapy responses. Overall, the analysis of different CAR-T cell phenotypes can be a key aspect for a better understanding of the whole CAR-T cell dynamics.
Collapse
Affiliation(s)
- Emanuelle A. Paixão
- Graduate Program, Laboratório Nacional de Computação Científica, Petrópolis 25651-075, Brazil
- Correspondence:
| | - Luciana R. C. Barros
- Center for Translational Research in Oncology, Instituto do Câncer do Estado de São Paulo, Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, São Paulo 01246-000, Brazil
| | - Artur C. Fassoni
- Institute for Mathematics and Computer Science, Universidade Federal de Itajubá, Itajubá 37500-903, Brazil
| | - Regina C. Almeida
- Computational Modeling Department, Laboratório Nacional de Computação Científica, Petrópolis 25651-075, Brazil
| |
Collapse
|
11
|
Cellular kinetics: A clinical and computational review of CAR-T cell pharmacology. Adv Drug Deliv Rev 2022; 188:114421. [PMID: 35809868 DOI: 10.1016/j.addr.2022.114421] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 06/28/2022] [Accepted: 07/04/2022] [Indexed: 12/20/2022]
Abstract
To the extent that pharmacokinetics influence the effectiveness of nonliving therapeutics, so too do cellular kinetics influence the efficacy of Chimeric Antigen Receptor (CAR) -T cell therapy. Like conventional therapeutics, CAR-T cell therapies undergo a distribution phase upon administration. Unlike other therapeutics, however, this distribution phase is followed by subsequent phases of expansion, contraction, and persistence. The magnitude and duration of these phases unequivocally influence clinical outcomes. Furthermore, the "pharmacodynamics" of CAR-T cells is truly dynamic, as cells can rapidly become exhausted and lose their therapeutic efficacy. Mathematical models are among the translational tools commonly applied to assess, characterize, and predict the complex cellular kinetics and dynamics of CAR-T cells. Here, we provide a focused review of the cellular kinetics of CAR-T cells, the mechanisms underpinning their complexity, and the mathematical modeling approaches used to interrogate them.
Collapse
|
12
|
Differential Response to Cytotoxic Drugs Explains the Dynamics of Leukemic Cell Death: Insights from Experiments and Mathematical Modeling. Symmetry (Basel) 2022. [DOI: 10.3390/sym14061269] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
This study presents a framework whereby cancer chemotherapy could be improved through collaboration between mathematicians and experimentalists. Following on from our recently published model, we use A20 murine leukemic cells transfected with monomeric red fluorescent proteins cells (mCherry) to compare the simulated and experimental cytotoxicity of two Federal Drug Administration (FDA)-approved anticancer drugs, Cytarabine (Cyt) and Ibrutinib (Ibr) in an in vitro model system of Chronic Lymphocytic Leukemia (CLL). Maximum growth inhibition with Cyt (95%) was reached at an 8-fold lower drug concentration (6.25 μM) than for Ibr (97%, 50 μM). For the proposed ordinary differential equations (ODE) model, a multistep strategy was used to estimate the parameters relevant to the analysis of in vitro experiments testing the effects of different drug concentrations. The simulation results demonstrate that our model correctly predicts the effects of drugs on leukemic cells. To assess the closeness of the fit between the simulations and experimental data, RMSEs for both drugs were calculated (both RMSEs < 0.1). The numerical solutions of the model show a symmetrical dynamical evolution for two drugs with different modes of action. Simulations of the combinatorial effect of Cyt and Ibr showed that their synergism enhanced the cytotoxic effect by 40%. We suggest that this model could predict a more personalized drug dose based on the growth rate of an individual’s cancer cells.
Collapse
|
13
|
Mahasa KJ, Ouifki R, Eladdadi A, Pillis LD. A combination therapy of oncolytic viruses and chimeric antigen receptor T cells: a mathematical model proof-of-concept. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:4429-4457. [PMID: 35430822 DOI: 10.3934/mbe.2022205] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Combining chimeric antigen receptor T (CAR-T) cells with oncolytic viruses (OVs) has recently emerged as a promising treatment approach in preclinical studies that aim to alleviate some of the barriers faced by CAR-T cell therapy. In this study, we address by means of mathematical modeling the main question of whether a single dose or multiple sequential doses of CAR-T cells during the OVs therapy can have a synergetic effect on tumor reduction. To that end, we propose an ordinary differential equations-based model with virus-induced synergism to investigate potential effects of different regimes that could result in efficacious combination therapy against tumor cell populations. Model simulations show that, while the treatment with a single dose of CAR-T cells is inadequate to eliminate all tumor cells, combining the same dose with a single dose of OVs can successfully eliminate the tumor in the absence of virus-induced synergism. However, in the presence of virus-induced synergism, the same combination therapy fails to eliminate the tumor. Furthermore, it is shown that if the intensity of virus-induced synergy and/or virus oncolytic potency is high, then the induced CAR-T cell response can inhibit virus oncolysis. Additionally, the simulations show a more robust synergistic effect on tumor cell reduction when OVs and CAR-T cells are administered simultaneously compared to the combination treatment where CAR-T cells are administered first or after OV injection. Our findings suggest that the combination therapy of CAR-T cells and OVs seems unlikely to be effective if the virus-induced synergistic effects are included when genetically engineering oncolytic viral vectors.
Collapse
Affiliation(s)
- Khaphetsi Joseph Mahasa
- Department of Mathematics and Computer Science, National University of Lesotho, Roma 180, Maseru, Lesotho
| | - Rachid Ouifki
- Department of Mathematics and Applied Mathematics, North-West University, Mafikeng campus, Private Bag X2046, Mmabatho 2735, South Africa
| | | | | |
Collapse
|
14
|
Brummer AB, Yang X, Ma E, Gutova M, Brown CE, Rockne RC. Dose-dependent thresholds of dexamethasone destabilize CAR T-cell treatment efficacy. PLoS Comput Biol 2022; 18:e1009504. [PMID: 35081104 PMCID: PMC8820647 DOI: 10.1371/journal.pcbi.1009504] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Revised: 02/07/2022] [Accepted: 01/12/2022] [Indexed: 12/14/2022] Open
Abstract
Chimeric antigen receptor (CAR) T-cell therapy is potentially an effective targeted immunotherapy for glioblastoma, yet there is presently little known about the efficacy of CAR T-cell treatment when combined with the widely used anti-inflammatory and immunosuppressant glucocorticoid, dexamethasone. Here we present a mathematical model-based analysis of three patient-derived glioblastoma cell lines treated in vitro with CAR T-cells and dexamethasone. Advanced in vitro experimental cell killing assay technologies allow for highly resolved temporal dynamics of tumor cells treated with CAR T-cells and dexamethasone, making this a valuable model system for studying the rich dynamics of nonlinear biological processes with translational applications. We model the system as a nonautonomous, two-species predator-prey interaction of tumor cells and CAR T-cells, with explicit time-dependence in the clearance rate of dexamethasone. Using time as a bifurcation parameter, we show that (1) dexamethasone destabilizes coexistence equilibria between CAR T-cells and tumor cells in a dose-dependent manner and (2) as dexamethasone is cleared from the system, a stable coexistence equilibrium returns in the form of a Hopf bifurcation. With the model fit to experimental data, we demonstrate that high concentrations of dexamethasone antagonizes CAR T-cell efficacy by exhausting, or reducing the activity of CAR T-cells, and by promoting tumor cell growth. Finally, we identify a critical threshold in the ratio of CAR T-cell death to CAR T-cell proliferation rates that predicts eventual treatment success or failure that may be used to guide the dose and timing of CAR T-cell therapy in the presence of dexamethasone in patients. Bioengineering and gene-editing technologies have paved the way for advance immunotherapies that can target patient-specific tumor cells. One of these therapies, chimeric antigen receptor (CAR) T-cell therapy has recently shown promise in treating glioblastoma, an aggressive brain cancer often with poor patient prognosis. Dexamethasone is a commonly prescribed anti-inflammatory medication due to the health complications of tumor associated swelling in the brain. However, the immunosuppressant effects of dexamethasone on the immunotherapeutic CAR T-cells are not well understood. To address this issue, we use mathematical modeling to study in vitro dynamics of dexamethasone and CAR T-cells in three patient-derived glioblastoma cell lines. We find that in each cell line studied there is a threshold of tolerable dexamethasone concentration. Below this threshold, CAR T-cells are successful at eliminating the cancer cells, while above this threshold, dexamethasone critically inhibits CAR T-cell efficacy. Our modeling suggests that in the presence of high dexamethasone reduced CAR T-cell efficacy, or increased exhaustion, can occur and result in CAR T-cell treatment failure.
Collapse
Affiliation(s)
- Alexander B. Brummer
- Department of Computational and Quantitative Medicine, Division of Mathematical Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
- * E-mail: (ABB); (CEB); (RCR)
| | - Xin Yang
- Department of Hematology and Hematopoietic Cell Translation and Immuno-Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
| | - Eric Ma
- Department of Hematology and Hematopoietic Cell Translation and Immuno-Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
| | - Margarita Gutova
- Department of Stem Cell Biology and Regenerative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
| | - Christine E. Brown
- Department of Hematology and Hematopoietic Cell Translation and Immuno-Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
- * E-mail: (ABB); (CEB); (RCR)
| | - Russell C. Rockne
- Department of Computational and Quantitative Medicine, Division of Mathematical Oncology, Beckman Research Institute, City of Hope National Medical Center, Duarte, California, United States of America
- * E-mail: (ABB); (CEB); (RCR)
| |
Collapse
|
15
|
Deng B, Pan J, Liu Z, Liu S, Chen Y, Qu X, Zhang Y, Lin Y, Zhang Y, Yu X, Zhang Z, Niu X, Luan R, Ma M, Li X, Liu T, Wu X, Niu H, Chang AH, Tong C. Peripheral leukemia burden at time of apheresis negatively affects the clinical efficacy of CART19 in refractory or relapsed B-ALL. Mol Ther Methods Clin Dev 2021; 23:633-643. [PMID: 34901308 PMCID: PMC8640733 DOI: 10.1016/j.omtm.2021.10.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Revised: 09/07/2021] [Accepted: 10/26/2021] [Indexed: 12/01/2022]
Abstract
Our previous clinical study achieved complete remission (CR) rates of >90% following chimeric antigen receptor T cells targeting CD19 (CART19) treatment of refractory/relapsed B cell acute lymphoblastic leukemia (r/r B-ALL); however, the influence of the leukemia burden in peripheral blood (PB) blasts remains unclear. Here, we retrospectively analyzed 143 patients treated with CART19 (including 36 patients with PB blasts) to evaluate the effect of peripheral leukemia burden at the time of apheresis. One hundred seventeen patients with high disease burdens achieved 91.5% CR or incomplete count recovery CR and 86.3% minimal residual disease-negative CR, and 26 patients with low disease burdens obtained 96.2% MRD− CR. Collectively, 9 of 36 (25%) patients with PB blasts and 2 of 107 (1.87%) patients without PB blasts did not respond to CART19 therapy. The leukemia burden in PB negatively influenced ex vivo cell characteristics, including the transduction efficiency of CD3+ T cells and their fold expansion, and in vivo cell dynamics, including peak CART19 proportion and absolute count, fold expansion, and persistence duration. Further studies showed that these patients had higher programmed death-1 expression in CART19 products. Our data imply that PB blasts negatively affected CART19 production and the clinical efficacy of CART19 therapy in patients with r/r B-ALL.
Collapse
Affiliation(s)
- Biping Deng
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Jing Pan
- Department of Hematology, Beijing Boren Hospital, Beijing 100070, China
| | - Zhaoli Liu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Shuangyou Liu
- Department of Hematology, Beijing Boren Hospital, Beijing 100070, China
| | - Yunlong Chen
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Xiaomin Qu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Yu'e Zhang
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Yuehui Lin
- Department of Hematology, Beijing Boren Hospital, Beijing 100070, China
| | - Yanlei Zhang
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, Tongji University School of Medicine, Shanghai 200433, China
| | - Xinjian Yu
- Medical Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Zhongxin Zhang
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Xuansha Niu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Rong Luan
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Ming Ma
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Xiaomei Li
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Tingting Liu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Xi'ai Wu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Huan Niu
- Cytology Laboratory, Beijing Boren Hospital, Beijing 100070, China
| | - Alex H. Chang
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, Tongji University School of Medicine, Shanghai 200433, China
- Corresponding author: Alex H. Chang, Clinical Translational Research Center, Shanghai Pulmonary Hospital, Tongji University School of Medicine, Shanghai, China.
| | - Chunrong Tong
- Department of Hematology, Beijing Boren Hospital, Beijing 100070, China
- Corresponding author: Chunrong Tong, Department of Hematology, Beijing Boren Hospital, No. 6, South Zhengwangfen, Fengtai District, Beijing 100070, China.
| |
Collapse
|
16
|
Adhikarla V, Awuah D, Brummer AB, Caserta E, Krishnan A, Pichiorri F, Minnix M, Shively JE, Wong JYC, Wang X, Rockne RC. A Mathematical Modeling Approach for Targeted Radionuclide and Chimeric Antigen Receptor T Cell Combination Therapy. Cancers (Basel) 2021; 13:cancers13205171. [PMID: 34680320 PMCID: PMC8533817 DOI: 10.3390/cancers13205171] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Revised: 09/30/2021] [Accepted: 10/07/2021] [Indexed: 12/25/2022] Open
Abstract
Simple Summary Targeted radionuclide therapy (TRT) and immunotherapy, an example being chimeric antigen receptor T cells (CAR-Ts), represent two potent means of eradicating systemic cancers. Although each one as a monotherapy might have a limited effect, the potency can be increased with a combination of the two therapies. The complications involved in the dosing and scheduling of these therapies make the mathematical modeling of these therapies a suitable solution for designing combination treatment approaches. Here, we investigate a mathematical model for TRT and CAR-T cell combination therapies. Through an analysis of the mathematical model, we find that the tumor proliferation rate is the most important factor affecting the scheduling of TRT and CAR-T cell treatments with faster proliferating tumors requiring a shorter interval between the two therapies. Abstract Targeted radionuclide therapy (TRT) has recently seen a surge in popularity with the use of radionuclides conjugated to small molecules and antibodies. Similarly, immunotherapy also has shown promising results, an example being chimeric antigen receptor T cell (CAR-T) therapy in hematologic malignancies. Moreover, TRT and CAR-T therapies possess unique features that require special consideration when determining how to dose as well as the timing and sequence of combination treatments including the distribution of the TRT dose in the body, the decay rate of the radionuclide, and the proliferation and persistence of the CAR-T cells. These characteristics complicate the additive or synergistic effects of combination therapies and warrant a mathematical treatment that includes these dynamics in relation to the proliferation and clearance rates of the target tumor cells. Here, we combine two previously published mathematical models to explore the effects of dose, timing, and sequencing of TRT and CAR-T cell-based therapies in a multiple myeloma setting. We find that, for a fixed TRT and CAR-T cell dose, the tumor proliferation rate is the most important parameter in determining the best timing of TRT and CAR-T therapies.
Collapse
Affiliation(s)
- Vikram Adhikarla
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA;
- Correspondence: (V.A.); (R.C.R.)
| | - Dennis Awuah
- Department of Hematology & Hematopoietic Cell Transplantation, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA; (D.A.); (A.K.); (X.W.)
| | - Alexander B. Brummer
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA;
| | - Enrico Caserta
- Department of Hematologic Malignancies Translational Science, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA; (E.C.); (F.P.)
| | - Amrita Krishnan
- Department of Hematology & Hematopoietic Cell Transplantation, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA; (D.A.); (A.K.); (X.W.)
| | - Flavia Pichiorri
- Department of Hematologic Malignancies Translational Science, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA; (E.C.); (F.P.)
| | - Megan Minnix
- Department of Molecular Imaging and Therapy, City of Hope National Medical Center, Duarte, CA 91010, USA; (M.M.); (J.E.S.)
| | - John E. Shively
- Department of Molecular Imaging and Therapy, City of Hope National Medical Center, Duarte, CA 91010, USA; (M.M.); (J.E.S.)
| | - Jeffrey Y. C. Wong
- Department of Radiation Oncology, City of Hope National Medical Center, Duarte, CA 91010, USA;
| | - Xiuli Wang
- Department of Hematology & Hematopoietic Cell Transplantation, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA; (D.A.); (A.K.); (X.W.)
| | - Russell C. Rockne
- Division of Mathematical Oncology, Department of Computational and Quantitative Medicine, Beckman Research Institute, City of Hope National Medical Center, Duarte, CA 91010, USA;
- Correspondence: (V.A.); (R.C.R.)
| |
Collapse
|
17
|
Ottesen JT, Andersen M. Potential of Immunotherapies in Treating Hematological Cancer-Infection Comorbidities-A Mathematical Modelling Approach. Cancers (Basel) 2021; 13:3789. [PMID: 34359690 PMCID: PMC8345105 DOI: 10.3390/cancers13153789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 07/08/2021] [Accepted: 07/23/2021] [Indexed: 11/17/2022] Open
Abstract
BACKGROUND The immune system attacks threats like an emerging cancer or infections like COVID-19 but it also plays a role in dealing with autoimmune disease, e.g., inflammatory bowel diseases, and aging. Malignant cells may tend to be eradicated, to appraoch a dormant state or escape the immune system resulting in uncontrolled growth leading to cancer progression. If the immune system is busy fighting a cancer, a severe infection on top of it may compromise the immunoediting and the comorbidity may be too taxing for the immune system to control. METHOD A novel mechanism based computational model coupling a cancer-infection development to the adaptive immune system is presented and analyzed. The model maps the outcome to the underlying physiological mechanisms and agree with numerous evidence based medical observations. RESULTS AND CONCLUSIONS Progression of a cancer and the effect of treatments depend on the cancer size, the level of infection, and on the efficiency of the adaptive immune system. The model exhibits bi-stability, i.e., virtual patient trajectories gravitate towards one of two stable steady states: a dormant state or a full-blown cancer-infection disease state. An infectious threshold curve exists and if infection exceed this separatrix for sufficiently long time the cancer escapes. Thus, early treatment is vital for remission and severe infections may instigate cancer progression. CAR T-cell Immunotherapy may sufficiently control cancer progression back into a dormant state but the therapy significantly gains efficiency in combination with antibiotics or immunomodulation.
Collapse
Affiliation(s)
- Johnny T. Ottesen
- Center for Mathematical Modeling-Human Health and Disease (COMMAND), Roskilde University, 4000 Roskilde, Denmark;
- IMFUFA, Department of Science and Environment, Roskilde University, 4000 Roskilde, Denmark
| | - Morten Andersen
- Center for Mathematical Modeling-Human Health and Disease (COMMAND), Roskilde University, 4000 Roskilde, Denmark;
- IMFUFA, Department of Science and Environment, Roskilde University, 4000 Roskilde, Denmark
| |
Collapse
|
18
|
Martínez-Rubio Á, Chulián S, Blázquez Goñi C, Ramírez Orellana M, Pérez Martínez A, Navarro-Zapata A, Ferreras C, Pérez-García VM, Rosa M. A Mathematical Description of the Bone Marrow Dynamics during CAR T-Cell Therapy in B-Cell Childhood Acute Lymphoblastic Leukemia. Int J Mol Sci 2021; 22:6371. [PMID: 34198713 PMCID: PMC8232108 DOI: 10.3390/ijms22126371] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/09/2021] [Accepted: 06/10/2021] [Indexed: 01/02/2023] Open
Abstract
Chimeric Antigen Receptor (CAR) T-cell therapy has demonstrated high rates of response in recurrent B-cell Acute Lymphoblastic Leukemia in children and young adults. Despite this success, a fraction of patients' experience relapse after treatment. Relapse is often preceded by recovery of healthy B cells, which suggests loss or dysfunction of CAR T-cells in bone marrow. This site is harder to access, and thus is not monitored as frequently as peripheral blood. Understanding the interplay between B cells, leukemic cells, and CAR T-cells in bone marrow is paramount in ascertaining the causes of lack of response. In this paper, we put forward a mathematical model representing the interaction between constantly renewing B cells, CAR T-cells, and leukemic cells in the bone marrow. Our model accounts for the maturation dynamics of B cells and incorporates effector and memory CAR T-cells. The model provides a plausible description of the dynamics of the various cellular compartments in bone marrow after CAR T infusion. After exploration of the parameter space, we found that the dynamics of CAR T product and disease were independent of the dose injected, initial B-cell load, and leukemia burden. We also show theoretically the importance of CAR T product attributes in determining therapy outcome, and have studied a variety of possible response scenarios, including second dosage schemes. We conclude by setting out ideas for the refinement of the model.
Collapse
Affiliation(s)
- Álvaro Martínez-Rubio
- Department of Mathematics, Universidad de Cádiz, Puerto Real, 11510 Cádiz, Spain; (S.C.); (M.R.)
- Biomedical Research and Innovation Institute of Cádiz (INiBICA), Hospital Universitario Puerta del Mar, 11009 Cádiz, Spain;
| | - Salvador Chulián
- Department of Mathematics, Universidad de Cádiz, Puerto Real, 11510 Cádiz, Spain; (S.C.); (M.R.)
- Biomedical Research and Innovation Institute of Cádiz (INiBICA), Hospital Universitario Puerta del Mar, 11009 Cádiz, Spain;
| | - Cristina Blázquez Goñi
- Biomedical Research and Innovation Institute of Cádiz (INiBICA), Hospital Universitario Puerta del Mar, 11009 Cádiz, Spain;
- Department of Pediatric Hematology and Oncology, Hospital de Jerez, 11407 Cádiz, Spain
| | - Manuel Ramírez Orellana
- Department of Paediatric Haematology and Oncology, Instituto Investigación Sanitaria La Princesa, Hospital Infantil Universitario Niño Jesús, 28006 Madrid, Spain;
| | - Antonio Pérez Martínez
- Translational Research in Pediatric Oncology, Hematopoietic Transplantation and Cell Therapy, IdiPAZ, Hospital Universitario La Paz, 28046 Madrid, Spain; (A.P.M.); (A.N.-Z.); (C.F.)
- Pediatric Hemato-Oncology Department, Hospital Universitario La Paz, 28046 Madrid, Spain
| | - Alfonso Navarro-Zapata
- Translational Research in Pediatric Oncology, Hematopoietic Transplantation and Cell Therapy, IdiPAZ, Hospital Universitario La Paz, 28046 Madrid, Spain; (A.P.M.); (A.N.-Z.); (C.F.)
| | - Cristina Ferreras
- Translational Research in Pediatric Oncology, Hematopoietic Transplantation and Cell Therapy, IdiPAZ, Hospital Universitario La Paz, 28046 Madrid, Spain; (A.P.M.); (A.N.-Z.); (C.F.)
| | - Victor M. Pérez-García
- Mathematical Oncology Laboratory (MOLAB), Instituto de Matemática Aplicada a la Ciencia y la Ingeniería, Universidad de Castilla-La Mancha, 13005 Ciudad Real, Spain;
- Departamento de Matemáticas, Escuela Técnica Superior de Ingenieros Industriales, Universidad de Castilla-La Mancha, 13005 Ciudad Real, Spain
| | - María Rosa
- Department of Mathematics, Universidad de Cádiz, Puerto Real, 11510 Cádiz, Spain; (S.C.); (M.R.)
- Biomedical Research and Innovation Institute of Cádiz (INiBICA), Hospital Universitario Puerta del Mar, 11009 Cádiz, Spain;
| |
Collapse
|
19
|
Mueller-Schoell A, Puebla-Osorio N, Michelet R, Green MR, Künkele A, Huisinga W, Strati P, Chasen B, Neelapu SS, Yee C, Kloft C. Early Survival Prediction Framework in CD19-Specific CAR-T Cell Immunotherapy Using a Quantitative Systems Pharmacology Model. Cancers (Basel) 2021; 13:2782. [PMID: 34205020 PMCID: PMC8199881 DOI: 10.3390/cancers13112782] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2021] [Revised: 05/22/2021] [Accepted: 05/28/2021] [Indexed: 12/20/2022] Open
Abstract
Chimeric antigen receptor (CAR)-T cell therapy has revolutionized treatment of relapsed/refractory non-Hodgkin lymphoma (NHL). However, since 36-60% of patients relapse, early response prediction is crucial. We present a novel population quantitative systems pharmacology model, integrating literature knowledge on physiology, immunology, and adoptive cell therapy together with 133 CAR-T cell phenotype, 1943 cytokine, and 48 metabolic tumor measurements. The model well described post-infusion concentrations of four CAR-T cell phenotypes and CD19+ metabolic tumor volume over 3 months after CAR-T cell infusion. Leveraging the model, we identified a low expansion subpopulation with significantly lower CAR-T cell expansion capacities amongst 19 NHL patients. Together with two patient-/therapy-related factors (autologous stem cell transplantation, CD4+/CD8+ T cells), the low expansion subpopulation explained 2/3 of the interindividual variability in the CAR-T cell expansion capacities. Moreover, the low expansion subpopulation had poor prognosis as only 1/4 of the low expansion subpopulation compared to 2/3 of the reference population were still alive after 24 months. We translated the expansion capacities into a clinical composite score (CCS) of 'Maximum naïve CAR-T cell concentrations/Baseline tumor burden' ratio and propose a CCSTN-value > 0.00136 (cells·µL-1·mL-1 as predictor for survival. Once validated in a larger cohort, the model will foster refining survival prediction and solutions to enhance NHL CAR-T cell therapy response.
Collapse
Affiliation(s)
- Anna Mueller-Schoell
- Department of Clinical Pharmacy and Biochemistry, Institute of Pharmacy, Freie Universitaet Berlin, 12169 Berlin, Germany; (A.M.-S.); (R.M.)
- Graduate Research Training Program PharMetrX, 12169 Berlin, Germany
| | - Nahum Puebla-Osorio
- Department of Lymphoma and Myeloma, Division of Cancer Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA; (N.P.-O.); (M.R.G.); (P.S.)
| | - Robin Michelet
- Department of Clinical Pharmacy and Biochemistry, Institute of Pharmacy, Freie Universitaet Berlin, 12169 Berlin, Germany; (A.M.-S.); (R.M.)
| | - Michael R. Green
- Department of Lymphoma and Myeloma, Division of Cancer Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA; (N.P.-O.); (M.R.G.); (P.S.)
| | - Annette Künkele
- Department of Pediatric Oncology and Hematology, Charité–Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt–Universität zu Berlin, Augustenburger Platz 1, 1335 Berlin, Germany;
- German Cancer Consortium (DKTK), Partner Site Berlin, CCC (Campus Mitte), 10178 Berlin, Germany
| | - Wilhelm Huisinga
- Institute of Mathematics, University of Potsdam, 14476 Potsdam, Germany;
| | - Paolo Strati
- Department of Lymphoma and Myeloma, Division of Cancer Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA; (N.P.-O.); (M.R.G.); (P.S.)
| | - Beth Chasen
- Department of Nuclear Medicine, Division of Diagnostic Imaging, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA;
| | - Sattva S. Neelapu
- Department of Lymphoma and Myeloma, Division of Cancer Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA; (N.P.-O.); (M.R.G.); (P.S.)
| | - Cassian Yee
- Department of Melanoma Medical Oncology, UT MD Anderson Cancer Center, Houston, TX 77030, USA
- Department of Immunology, UT MD Anderson Cancer Center, Houston, TX 70030, USA
| | - Charlotte Kloft
- Department of Clinical Pharmacy and Biochemistry, Institute of Pharmacy, Freie Universitaet Berlin, 12169 Berlin, Germany; (A.M.-S.); (R.M.)
| |
Collapse
|
20
|
Nukala U, Rodriguez Messan M, Yogurtcu ON, Wang X, Yang H. A Systematic Review of the Efforts and Hindrances of Modeling and Simulation of CAR T-cell Therapy. AAPS JOURNAL 2021; 23:52. [PMID: 33835308 DOI: 10.1208/s12248-021-00579-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Accepted: 03/06/2021] [Indexed: 01/08/2023]
Abstract
Chimeric antigen receptor (CAR) T-cell therapy is an immunotherapy that has recently become highly instrumental in the fight against life-threatening diseases. A variety of modeling and computational simulation efforts have addressed different aspects of CAR T-cell therapy, including T-cell activation, T- and malignant cell population dynamics, therapeutic cost-effectiveness strategies, and patient survival. In this article, we present a systematic review of those efforts, including mathematical, statistical, and stochastic models employing a wide range of algorithms, from differential equations to machine learning. To the best of our knowledge, this is the first review of all such models studying CAR T-cell therapy. In this review, we provide a detailed summary of the strengths, limitations, methodology, data used, and data gap in currently published models. This information may help in designing and building better models for enhanced prediction and assessment of the benefit-risk balance associated with novel CAR T-cell therapies, as well as with the data need for building such models.
Collapse
Affiliation(s)
- Ujwani Nukala
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Marisabel Rodriguez Messan
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Osman N Yogurtcu
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Xiaofei Wang
- Office of Tissues and Advanced Therapies, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Hong Yang
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA.
| |
Collapse
|
21
|
Kimmel GJ, Locke FL, Altrock PM. The roles of T cell competition and stochastic extinction events in chimeric antigen receptor T cell therapy. Proc Biol Sci 2021; 288:20210229. [PMID: 33757357 PMCID: PMC8059581 DOI: 10.1098/rspb.2021.0229] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Chimeric antigen receptor (CAR) T cell therapy is a remarkably effective immunotherapy that relies on in vivo expansion of engineered CAR T cells, after lymphodepletion (LD) by chemotherapy. The quantitative laws underlying this expansion and subsequent tumour eradication remain unknown. We develop a mathematical model of T cell-tumour cell interactions and demonstrate that expansion can be explained by immune reconstitution dynamics after LD and competition among T cells. CAR T cells rapidly grow and engage tumour cells but experience an emerging growth rate disadvantage compared to normal T cells. Since tumour eradication is deterministically unstable in our model, we define cure as a stochastic event, which, even when likely, can occur at variable times. However, we show that variability in timing is largely determined by patient variability. While cure events impacted by these fluctuations occur early and are narrowly distributed, progression events occur late and are more widely distributed in time. We parameterized our model using population-level CAR T cell and tumour data over time and compare our predictions with progression-free survival rates. We find that therapy could be improved by optimizing the tumour-killing rate and the CAR T cells' ability to adapt, as quantified by their carrying capacity. Our tumour extinction model can be leveraged to examine why therapy works in some patients but not others, and to better understand the interplay of deterministic and stochastic effects on outcomes. For example, our model implies that LD before a second CAR T injection is necessary.
Collapse
Affiliation(s)
- Gregory J Kimmel
- Department of Integrated Mathematical Oncology, Research Institute, Tampa, FL, USA
| | - Frederick L Locke
- Department of Blood and Marrow Transplant and Cellular Immunotherapy, Research Institute, Tampa, FL, USA
| | - Philipp M Altrock
- Department of Integrated Mathematical Oncology, Research Institute, Tampa, FL, USA.,Department of Blood and Marrow Transplant and Cellular Immunotherapy, Research Institute, Tampa, FL, USA.,Department of Malignant Hematology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| |
Collapse
|
22
|
Arya V, Venkatakrishnan K. Role of Physiologically Based Pharmacokinetic Modeling and Simulation in Enabling Model-Informed Development of Drugs and Biotherapeutics. J Clin Pharmacol 2020; 60 Suppl 1:S7-S11. [PMID: 33205427 DOI: 10.1002/jcph.1770] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2020] [Accepted: 10/05/2020] [Indexed: 12/14/2022]
Affiliation(s)
- Vikram Arya
- Division of Infectious Disease Pharmacology (DIDP), Office of Clinical Pharmacology (OCP), Office of Translational Sciences (OTS), Center for Drug Evaluation and Research, U.S. Food and Drug Administration, Silver Spring, Maryland, USA
| | - Karthik Venkatakrishnan
- EMD Serono Research and Development Institute, Inc. (a business of Merck KGaA, Darmstadt, Germany), Billerica, Massachusetts, USA
| |
Collapse
|