N.L. Grigorenko, E.N. Khailov, E.V. Grigorieva, A.D. Klimenkova. Optimal strategies of CAR T-Cell therapy in the treatment of leukemia within the Lotka–Volterra predator–prey model ... P. 43-58

A controlled mathematical model of leukemia treatment is considered. The model is based on the three-dimensional Lotka–Volterra predator–prey model, which describes a recently developed leukemia treatment technology called Chimeric Antigen Receptor (CAR) T-cell therapy, and is given on a fixed time interval by a system of four differential equations. The equations describe the interaction between populations of healthy and cancer cells, CAR T-cells, and cytokines. The CAR T-cells act as predators, while healthy and cancer cells act as prey. The CAR T-cell therapy leads to serious side-effects associated with the rapid growth of cytokines, and therefore their dynamics is also included in the model. The model also contains two bounded control functions reflecting the intensity of the therapy (the first control) and the intensity of administration of drugs that suppress the activity of the immune system (the second control). We study the problem of minimizing the objective function related to the number of cancer and healthy cells, as well as cytokines, both at the final moment of a given time interval and during this entire interval. The Pontryagin maximum principle is applied for the analysis of the problem; it is used to establish the bang-bang nature of an optimal first control and to estimate the number of its switchings. It is shown that an optimal second control is a constant function on the entire time interval. The BOCOP-2.2.1 environment is used for the numerical analysis of the problem. The results of numerical calculations are presented, demonstrating various types of optimal protocols for CAR-T therapy.

Keywords: leukemia, nonlinear control system, optimal control, Pontryagin maximum principle, bang-bang control, switching function, generalized Rolle theorem

Received March 25, 2021

Revised May 17, 2021

Accepted June 21, 2021

Nikolai Leont’evich Grigorenko, Dr. Phys.-Math. Sci., Prof., Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, 119992, Russia, e-mail: grigor@cs.msu.ru

Evgenii Nikolaevich Khailov, Cand. Sci. (Phys.-Math.), Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, 119992, Russia, e-mail: khailov@cs.msu.su

Ellina Valer’evna Grigorieva, Cand. Sci. (Phys.-Math.), Prof., Department of Mathematics and Computer Sciences, Texas Woman’s University, TX 76204, USA, e-mail: egrigorieva@mail.twu.edu

Anna Dmitrievna Klimenkova, undergraduate student, Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, 119992, Russia, e-mail: klimenkovaad@mail.ru

REFERENCES

1.   Heymach J. [et all.] Clinical cancer advances 2018: annual report on progress against cancer from the American Society of Clinical Oncology. J. Clin. Oncol., 2018, vol. 36, no. 10, pp. 1020–1044. doi: 10.1200/JCO.2017.77.0446 

2.   June C.H., Sadelain M. Chimeric antigen receptor therapy. N. Engl. J. Med., 2018, vol. 379, no. 1, pp. 64–73. doi: 10.1056/NEJMra1706169 

3.   Hopkins B., Tucker M., Pan Y., Fang N., Huang Z. A model-based investigation of cytokine storm for T-cell therapy. IFAC-PapersOnLine, 2018, vol. 51, no. 19, pp. 76–79. doi: 10.1016/j.ifacol.2018.09.039 

4.   Barros L.R.C., Rodrigues B.J., Almeida R.C. CAR-T cell goes on a mathematical model. J. Cell Immunol., 2020, vol. 2, no. 1, pp. 31–37. doi: 10.33696/immunology.2.016 

5.   Konstorum A., Vella A.T., Adler A.J., Laubenbacher R.C. Addressing current challenges in cancer immunotherapy with mathematical and computational modelling. J. Roy. Soc. Interface, 2017, vol. 14, no. 131, art. ID: 20170150, 10 p. doi: 10.1098/rsif.2017.0150 

6.   Valentinuzzi D., Jeraj R. Computational modelling of modern cancer immunotherapy. Phys. Med. Biol., 2020, vol. 65, no. 24, art. ID: 24TR01, 22 p. doi: 10.1088/1361-6560/abc3fc 

7.   Leon-Triana O., Sabir S., Calvo G.F., Belmonte-Beitia J., Chulian S., Martinez-Rubio A., Rosa M., Perez-Martinez A., Ramirez-Orellana M., Perez-Garcia V.M. CAR T-cell in B-cell acute lymphoblastic leukaemia: insights from mathematical models. Commun. Nonlinear Sci., 2021, vol. 94, art. ID: 105570, doi: 10.1016/j.cnsns.2020.105570 

8.   Perez-Garcia V.M., Leon-Triana O., Rosa M., Perez-Martinez A. CAR T cells for T-cell leukemias: insights from mathematical models. arXiv:2004.14291 [q-bio.TO] 26 Apr 2020, 20 p.

9.   Sahoo P. [et al.] Mathematical deconvolution of CAR T-cell proliferation and exhaustion from real-time killing assay data. J. Roy. Soc. Interface, 2019, vol. 17, no. 162, art. ID: 20190734, 10 p. doi: 10.1098/rsif.2019.0734 

10.   Nichelatti M. A mathematical model for the chimeric antigen receptor T cell (CAR-T) therapy as a Lotka–Volterra system. J. Math. Stat. Res., 2020, vol. 2, no. 3, art. ID: 136, 4 p. doi: 10.36266/JMSR/136 

11.   Khailov E.N., Klimenkova A.D., Korobeinikov A. Optimal control for anti-cancer therapy. In: Extended abstracts spring 2018. A. Korobeinikov, M. Caubergh, T. Lazaro, J. Sardanyes (eds). Basel: Birkhauser, 2019, pp. 35–43, Ser. Trends in mathematics; vol. 11. doi: 10.1007/978-3-030-25261-8_6 

12.   Grigorenko N.L., Khailov E.N., Klimenkova A.D., Korobeinikov A. Program and positional control strategies for the Lotka – Volterra competition model. In: Stability, Control and Differential Games. Proceedings of the International Conference “Stability, Control, Differential Games” (SCDG2019). A. Tarasyev, V. Maksimov, T. Filippova (eds). Cham: Springer Nature, 2020, pp. 39–49. doi: 10.1007/978-3-030-42831-0_4 

13.   Grigorenko N.L., Khailov E.N., Grigorieva E.V., Klimenkova A.D. Optimal strategies in the treatment of cancers in the Lotka – Volterra mathematical model of competition. Trudy Inst. Mat. i Mekh. UrO RAN, 2020, vol. 26, no. 1, pp. 71–88 (in Russian). doi: 10.21538/0134-4889-2020-26-1-71-88 

14.   Mostolizadeh R., Afsharnezhad Z., Marciniak-Czochra A. Mathematical model of chimeric anti-gene receptor (CAR) T cell therapy with presence of cytokine. Numer. Algebr. Control Optim., 2018, vol. 8, no. 1, pp. 63–80. doi: 10.3934/naco.2018004 

15.   Khailov E., Grigorieva E., Klimenkova A. Optimal CAR T-cell immunotherapy strategies for a leukemia treatment model. Games, 2020, vol. 11, no. 4, art. ID: 53, 26 p. doi: 10.3390/g11040053 

16.   Pillis L.G., Radunskaya A. A mathematical tumor model with immune resistance and drug therapy: an optimal control approach. J. Theoretical Medicine, 2001, vol. 3, art. ID: 318436, pp. 79–100. doi: 10.1080/10273660108833067 

17.   Pillis L.G., Radunskaya A. The dynamics of an optimally controlled tumor model: a case study. Math. Comput. Model., 2003, vol. 37, no. 11, pp. 1221–1244. doi: 10.1016/S0895-7177(03)00133-X 

18.   Bazykin A.D. Nonlinear dynamics of interacting populations. Singapore: World Scientific, 1998, 216 p. doi: 10.1142/2284 . Original Russian text published in Bazykin A.D. Nelineinaya dinamika vzaimodeistvuyushchikh populyatsii. Moscow; Izhevsk: Institut komp’yuternykh issledovanii, 2003, 368 p.

19.   Lee E.B., Markus L. Foundations of optimal control theory. N Y; London; Sydney: John Wiley & Sons, Inc., 1967, 576 p. Translated to Russian under the title Osnovy teorii optimal’nogo upravleniya. Moscow: Nauka Publ., 1972, 576 p. ISBN: 0471522635 .

20.   Vasil’ev F.P. Metody optimizatsii [Optimization methods]. Moscow: Factorial Press, 2002, 824 p.

21.   Hartman P. Ordinary differential equations. N Y: Wiley, 1964, 612 p. Translated to Russian under the title Obyknovennye differentsial’nye uravneniya. Moscow: Mir Publ., 1970, 720 p.

22.   Dmitruk A.V. A generalized estimate on the number of zeros for solutions of a class of linear differential equations. SIAM J. Control Optim., 1992, vol. 30, no. 5, pp. 1087–1091. doi: 10.1137/0330057 

23.   Kuzenkov O.A., Ryabova E.A. Matematicheskoe modelirovanie protsessov otbora [Mathematical simulation of selection processes]. Nizhnii Novgorod: Izdatel’stvo Nizhnegorodskogo universiteta, 2007, 324 p. ISBN: 978-5-91326-011-6 .

24.   Bonnans F., Martinon P., Giorgi D., Grelard V., Maindrault S., Tissot O., Liu J. BOCOP 2.2.1 — user guide [e-resource]. August 8, 2019. Available on: http://bocop.org 

Cite this article as: N.L. Grigorenko, E.N. Khailov, E.V. Grigorieva, A.D. Klimenkova. Optimal strategies of CAR T-Cell therapy in the treatment of leukemia within the Lotka–Volterra predator–prey model. Trudy Instituta Matematiki i Mekhaniki UrO RAN, 2021, vol. 27, no. 3, pp. 43–58.