TY - JOUR

T1 - Optimal Control Theory and Calculus of Variations in Mathematical Models of Chemotherapy of Malignant Tumors

AU - Subbotina, Nina

AU - Novoselova, Natalia

AU - Krupennikov, Evgenii

N1 - The contribution of Nina Subbotina and Evgenii Krupennikov was performed as part of research conducted in the Ural Mathematical Center with the financial support of the Ministry of Science and Higher Education of the Russian Federation (Agreement number 075-02-2023-913).

PY - 2023

Y1 - 2023

N2 - This paper is devoted to the analysis of mathematical models of chemotherapy for malignant tumors growing according to the Gompertz law or the generalized logistic law. The influence of the therapeutic agent on the tumor dynamics is determined by a therapy function depending on the time-varying concentration of the drug in the patient’s body. The case of a non-monotonic therapy function with two maxima is studied. It reflects the use of two different therapeutic agents. The state variables of the dynamics are the tumor volume and the amount of the therapeutic agent able to suppress malignant cells (concentration of the drug in the body). The treatment protocol (the rate of administration of the therapeutic agent) is the control in the dynamics. The optimal control problem for this models is considered. It is the problem of the construction of treatment protocols that provide the minimal tumor volume at the end of the treatment. The solution of this problem was obtained by the authors in previous works via the optimal control theory. The form of the considered therapy functions provides a specific structure for the optimal controls. The managerial insights of this structure are discussed. In this paper, the structure of the viability set is described for the model according to the generalized logistic law. It is the set of the initial states of the model for which one can find a treatment protocol that guarantees that the tumor volume remains within the prescribed limits throughout the treatment. The description of the viability set’s structure is based on the optimal control theory and the theory of Hamilton–Jacobi equations. An inverse problem of therapy is also considered, namely the problem of reconstruction of the treatment protocol and identification of the unknown parameter of the intensity of the tumor growth. Reconstruction is carried out by processing information about the observations of the tumor volume dynamics and the measurements of the drug concentration in the body. A solution to this problem is obtained through the use of a method based on the calculus of variations. The results of the numerical simulations are presented herein.

AB - This paper is devoted to the analysis of mathematical models of chemotherapy for malignant tumors growing according to the Gompertz law or the generalized logistic law. The influence of the therapeutic agent on the tumor dynamics is determined by a therapy function depending on the time-varying concentration of the drug in the patient’s body. The case of a non-monotonic therapy function with two maxima is studied. It reflects the use of two different therapeutic agents. The state variables of the dynamics are the tumor volume and the amount of the therapeutic agent able to suppress malignant cells (concentration of the drug in the body). The treatment protocol (the rate of administration of the therapeutic agent) is the control in the dynamics. The optimal control problem for this models is considered. It is the problem of the construction of treatment protocols that provide the minimal tumor volume at the end of the treatment. The solution of this problem was obtained by the authors in previous works via the optimal control theory. The form of the considered therapy functions provides a specific structure for the optimal controls. The managerial insights of this structure are discussed. In this paper, the structure of the viability set is described for the model according to the generalized logistic law. It is the set of the initial states of the model for which one can find a treatment protocol that guarantees that the tumor volume remains within the prescribed limits throughout the treatment. The description of the viability set’s structure is based on the optimal control theory and the theory of Hamilton–Jacobi equations. An inverse problem of therapy is also considered, namely the problem of reconstruction of the treatment protocol and identification of the unknown parameter of the intensity of the tumor growth. Reconstruction is carried out by processing information about the observations of the tumor volume dynamics and the measurements of the drug concentration in the body. A solution to this problem is obtained through the use of a method based on the calculus of variations. The results of the numerical simulations are presented herein.

UR - http://www.scopus.com/inward/record.url?partnerID=8YFLogxK&scp=85175044604

UR - https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcAuth=tsmetrics&SrcApp=tsm_test&DestApp=WOS_CPL&DestLinkType=FullRecord&KeyUT=001095078400001

U2 - 10.3390/math11204301

DO - 10.3390/math11204301

M3 - Article

VL - 11

JO - Mathematics

JF - Mathematics

SN - 2227-7390

IS - 20

M1 - 4301

ER -