Next Article in Journal
Certain Solutions of Abel’s Integral Equations on Distribution Spaces via Distributional Gα-Transform
Previous Article in Journal
L(2, 1)-Labeling Halin Graphs with Maximum Degree Eight
Previous Article in Special Issue
The Stochastic Approach for SIR Epidemic Models: Do They Help to Increase Information from Raw Data?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model-Based Regional Control with Anomalous Diffusion of Multi-Drug Combined Cancer Therapy for Volume Predictions

by
Clara Mihaela Ionescu
1,2,† and
Maria Ghita
1,2,3,*,†
1
Research Group on Dynamical Systems and Control, Faculty of Engineering and Architecture, Ghent University, Tech Lane Science Park 125, 9052 Gent, Belgium
2
Flanders Make Core Lab Engineering in Machineries, Intelligence, Robotics and Electromechanics, Tech Lane Science Park 131, 9052 Gent, Belgium
3
Faculty of Medicine and Health Sciences, Antwerp University, 2610 Wilrijk, Belgium
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2023, 15(1), 51; https://doi.org/10.3390/sym15010051
Submission received: 14 November 2022 / Revised: 13 December 2022 / Accepted: 19 December 2022 / Published: 25 December 2022

Abstract

:
Symmetry breaking in the anatomical lung is triggered by tumorigenesis and disrupted by delivering single or multiple drugs to stop the progression of the tumor and treat cancer. In this study, a prior model of combined drug therapy is augmented to introduce tissue heterogeneity when the drug is applied in multi-drug therapy of lung cancer. Patient-related drug resistance and synergy are investigated as a function of diffusion intensity as drug molecules reach the tumor site. The results indicate that diffusion of drug molecules plays an important role next to other factors such as patient sensitivity to the drug and drug synergy effects. We conclude that the minimal model provides meaningful predictions on tumor growth at the intermediate mesoscale level. With such models at hand, it is now possible to employ model-based control algorithms to optimize the dose profiles in terms of time and amount. In this paper, we present a theoretical framework for control employing networked game theory optimality. Specific situations are discussed in terms of finding optimality at Nash equilibrium in relation to patient response and drug synergy effects.

1. Introduction

Every healthy individual cell follows a well-controlled symmetry pattern in the biological processes of conception, proliferation, differentiation, and death (the process of programmed cell death in multi-cellular organisms referred to as apoptosis). However, once developed, the cell clones may suffer a symmetry breaking due to mutational changes that lead to cancer, either genetic or other alterations.
An understanding of the growth and progression of primary and metastatic tumors requires the elucidation of mechanisms that break the symmetry in the lung and provoke tumorigenesis. Hence, the alteration of the respiratory system is characterized by the collective properties exhibited by tumorous cells and their interaction with the microenvironment and delivered treatment. In general, the tumor microenvironment can be described as a discrete stochastic accumulation of tumorous cells by continuous nonlinear displayed stromal cells. By contrast, healthy lung tissue can be described as having a fractal distribution, without disruption in the symmetry between the cells [1]. A recent study modeling different biological conditions of cancer tissue showed that the proliferation of cancer cells alters the tumor dynamics and heterogeneity [2].
The processes that govern tumor growth have been investigated and translated into mathematical models that allow analyzing of the interactions on the cancer site, and predict the evolution of the tissue and treatment outcome across multiple scales [3,4]. Monitoring pre- and post-treatment of tumorous tissue requires recursive alterations to the initially planned therapy profiles, depending on the patient’s therapeutic response [5]. In particular, for lung cancer-diagnosed patients, it has been recently suggested that a cocktail of multi-drug therapy has better outcomes in terms of clinical effect and patient recovery period [6,7]. Other recent investigations in mouse models indicated that such multi-drug therapy with adaptive protocols achieves better outcomes in terms of cell turnover and therapies which used as little drug as possible worked best [8]. However, models are scarce and the main advances are at the molecular scale of tumor dynamic profiles and dose-effect kinetics [9,10]. As part of the breathing dynamic system, the tumor volume and consistency affect respiratory mechanics [11], requiring adaptive tumor radiotherapy localization [12]. Using tissue properties in lungs [13], anomalous diffusion of the drug into the tissue and multi-drug therapy can be combined to enhance the benefits of augmented models to aid the decision making-mechanisms in therapy profiles [14,15].
Recent advances in modeling lung tumor growth dynamics have proposed simplistic yet descriptive compartmental models for characterizing the pharmacokinetics (PK) and pharmacodynamics (PD) of lung cancer therapy in non-small cell lung cancer (NSCLC). The most prevalent therapy is antiangiogenesis and a model has been proposed in [16,17]. In clinical practice, this is enhanced with periods of stereotactic body radiotherapy (SBRT) [18] and in severe cases, additional immunotherapy [19] is applied. We proposed a PKPD model encompassing all three therapies [20,21] and validated it on a set of patients with lung cancer [22]. The PKPD model was first developed and described in detail in [20]. In [22], we calibrated the model with real data measured in a clinical trial including patients with non-small cell lung cancer. The parameters capture the dynamics of the tumor and the results show the changes in patients’ responses to treatment according to tumor shrinkage. The model fitted the behavior of the lung tumors assessed by measuring the tumor volume of the patients before and after treatment from computed tomography images.
Anomalous diffusion is an important factor in describing the efficacy of a drug on the tumor volume dynamics, thus we propose to augment the model with an additional parameter related to the dose-effect intensity. Having this augmented model available at hand, we also evaluate control algorithms for optimizing multiple scenarios of therapy profiles. Distributed agent optimal control strategies often involve the existence of an equilibrium point where all agents converge towards a common objective or a self-centered objective. As such, one may project a theory that multi-agent systems seeking optimality can well approximate the multi-drug therapy protocols seeking optimal tumor volume response in patients.
In this paper, we introduce a coalition-based game optimality control theory to examine the effects of multi-drug optimal profiles on tumor volume dynamics in lung cancer. We endeavor to build a reliable mechanistic modeling framework for describing drug distribution in solid tumors due to fractional diffusion in a locally targeted approach and simulate their effects on tumor growth. We used a general pharmacokinetic-pharmacodynamic (PKPD) modeling in a tumor environment, focusing on the distribution of local targeted anti-cancer therapies and the dose-effect relationship. To evaluate the proposed rationale, we analyzed different examples simulated using both data from the literature and clinical data.

2. Materials and Methods

2.1. A Minimalistic PKPD Model of Lung Tumor Growth

A tumor volume growth model has been proposed in [20], calibrated in [22], and validated with real data from patients in [11]. It consists of first-order compartmental PK models with linear dynamics, followed by dose-effect response nonlinear gains. Although this is a well-known Wiener-type model, identification is difficult due to the lack of modalities for persistent excitation of drug dosing profiles [23,24]. The proliferating tumor volume x 1 in mm 3 and the necrotic tumor volume x 2 in mm 3 , give the total tumor volume. The SBRT inhibitor level x 5 in mg/(mL·day) is expressed based on the radiation dose rate administrated u r (mg · day/mL), according to the dosing profiles adequate for NSCLC. Dosing profiles for antiangiogenic u a and immunotherapy u i give the concentrations x 3 and x 4 , respectively, leading to the overall model given by:
x ˙ 1 = ( a n ) x 1 E · x 1 x ˙ 2 = x 2 ( 0 ) + n x 1 + E · x 1 x ˙ 3 = c a x 3 + u a x e ˙ 3 = c a x e 3 + E t a · x 3 x ˙ 4 = c i x 4 + u i x e ˙ 4 = c i x e 4 + E t i · x 4 x ˙ 5 = c r x 5 + u r x e ˙ 5 = c r x e 5 + E t r · x 5
where a denotes the tumor growth rate, n the necrosis rate, and c r is the clearance rate on the Michaelis–Menten kinetics x 1 x 3 E D 50 r + x 3 (mm 3 /day). As the dead cells are naturally eliminated from the body, the necrotic volume does not depend on the past value, hence x 2 ( 0 ) = 0 . The parameters x e i and E t i are the effects of drug concentrations and the synergic effect between tumor cells and chosen therapy. The dose-effect PD model is in fact a nonlinear gain (for a single drug) or a nonlinear surface (for multiple drugs) and can be considered as averaged effects of all therapies (in absence of real data at the molecular level to confirm this hypothesis) for the tumor E t a l l and drug E d a l l interactions:
E t a l l = E t a + E t i + E t r 3 E d a l l = E a i + E a r + E i r 3 E = E t a l l + E d a l l 2
with E t x denoting the interaction between tumor cells and each drug, while E x y denotes the interaction among drugs. When surface models are used to characterize synergic effects among drugs, the effect drug concentrations x e are normalized to their potency, i.e., to their corresponding half-effect concentration C 50 . The combined effects of two drugs U A and U B are considered as a new drug, and expressed as a Hill curve dose–response relationship 3D surface:
E f f e c t = I γ 1 + I γ
with I denoting the interaction term:
I = U n A + U n B + σ U n A · U n B
with U n A = U A C 50 A and U n B = U B C 50 B the normalized drug effect concentrations and C 50 the concentrations at half effect 50%. The term γ denotes the nonlinearity of the surface, which represents how a patient responds to the drug (effectiveness or resistance of therapy). The term σ denotes the degree of synergy present between the drugs. Values for the model coefficients are those reported in [20] (open access).

2.2. Regional Anomalous Diffusion

A hybrid model for characterizing regional anomalous diffusion at the surface tension between healthy tissue and tumor tissue is proposed. Having patient-specific parameters, this model has also been validated partially to illustrate its clinical relevance. We propose here an aggregation of local diffusion phenomena to be included as part of the dose-effect relationship, i.e., (3). The effect surface depends not only on the γ and σ parameters characterizing drug–patient interaction but also on local effects in molecular binding between drugs and targeted treatment of the cancer tissue.
Let us assume that the effect of (2) is no longer abiding by the classical Fickean diffusion, but anomalous diffusion. When working with the surface model, it is important to understand the relationship between power-law and exponential functions, and this has been in detail discussed in [25]. Combining the theory of molecular residential time and surface interaction with simplifying algebra from [25], leads to introducing an extra parameter in the relation for x 1 and x 2 , namely:
x ˙ 1 = ( a n ) x 1 E ( 1 α ) · x 1 x ˙ 2 = x 2 ( 0 ) + n x 1 + E ( 1 α ) · x 1
whereas α is the coefficient of diffusion, the particular case of α = 1 giving a linear diffusion pattern along the surface (Fickean diffusion). Indeed, we verify that for α = 1 , the surface remains linear along its domain, as depicted in Figure 1.

2.3. Predictive Control Strategy for Multi-Drug Therapy Optimization

Predictive control has been used in SBRT for compensating the breathing pattern effects in 3D volume changes in lung tumor tissue, to minimize towards zero the radiation of healthy tissue around the tumor [26]. The breathing pattern was identified online and the robot arm used for therapy was guided with feedforward compensation using a specially designed disturbance model in Model Predictive Control (MPC) [12]. When multiple objectives are envisaged within changing context of execution, a prioritized optimizing scheme can be used as that proposed in [27], reducing both computational and numerical complexity involved when the Pareto front is used instead.
In this paper, a centralized MPC strategy with the state-space formulation is proposed. The algorithm is derived starting from the velocity-form methodology from [28], subsequently extended to a multivariable system.
Let us consider the PKPD model from (1) described by the process model in the state-space formulation:
x p ( k + 1 ) = A p x p ( k ) + B p u ( k ) y ( k ) = C p x p ( k )
where k is the discrete-time instant, u R n u , y R n y and x p R n x p are the input, output, and state variables, respectively. Using the methodology described in [28], the difference operation is applied on both sides of (6) resulting:
Δ x p ( k + 1 ) = A p Δ x p ( k ) + B p Δ u ( k ) , Δ y ( k + 1 ) = C p A p Δ x p ( k ) + C p B p Δ u ( k )
Note that (7) introduces the increments of the variables x p , u and y, with Δ y ( k + 1 ) = y ( k + 1 ) y ( k ) . A new state variable x ( k ) = Δ x p ( k ) T y ( k ) T is introduced, resulting the augmented model:
Δ x p ( k + 1 ) y ( k + 1 ) x ( k + 1 ) = A p O n y × n x p T C p A p 1 A Δ x p ( k ) y ( k ) x ( k ) + B p C p B p B Δ u ( k ) y ( k ) = O n y × n x p I n y C Δ x p ( k ) y ( k )
which will be used to design the predictive controller. Note that the new input of the state-space model in velocity-form is Δ u ( k ) .
The model (8) can be written in a compressed form as:
x ( k + 1 ) = A x ( k ) + B Δ u ( k ) y ( k ) = C x ( k )
The centralized MPC cost function is defined as:
J ( x ( k ) , Δ U ( k ) ) = R s Y T R s Y + Δ U ( k ) T R Δ U ( k )
depending on the output future predictor
Y = y ( k + 1 | k ) y ( k + N p | k ) T
and the future input sequence
Δ U ( k ) = Δ u ( k | k ) Δ u ( k + N c 1 | k ) T
with N p the prediction horizon and N c the control horizon ( N c N p ). In this paper, for simplicity, we consider N c = N p . The predicted reference trajectory R s R N p assumed constant and equal with the setpoint at time instant k and the input weight matrix has the form R = α I N c , α 0 .
Using the ’velocity-form’ model (9), the prediction of state and input variables can be computed. The centralized Δ U * ( k ) solution is obtained:
Δ U * = ( B ˜ T B ˜ + R ) 1 B ˜ T [ R s A ˜ x ( k ) ] .
Following the receding horizon principle, only the first n u elements from the optimal solution vector are sent to the process [29].

2.4. Multi-Agent Nash Optimality and Coalition Control

In terms of the dose-effect relationship versus drug resistance, the problem of optimum seeking is both exploitative and explorative, as in reinforcement learning theory. Such learning schemes are highly relevant in network systems with heterogeneous entities, where these entities represent different drug therapy profiles and the agents in the network are the respective drug selection combinations. Finding the optimum implies finding an equilibrium point where the cell turnover outperforms the cell growth in tumor volume. As such, this corresponds broadly to a combination of game theoretic models with learning-based approaches, which we will employ in our analysis. In a decentralized exchange of information, it allows heterogeneous agents to strategically interact with each other, e.g., the choice of drugs affects the degree of synergy effects, and learn to adjust their behaviors, i.e., an adaptive therapy protocol strategy. Information here refers to the structure used to model the knowledge the players in the game acquire and the history of their decision-effect when they make the decisions for their next move. We introduce here some concepts used hereafter:
  • Players are the participants in a game, in competition against each other. In our context these are the different multi-drug selections and protocols competing for the best patient outcome;
  • Actions of a player, denoting here the drug profiles and timeline administered to the patient;
  • Information in game theory refers to acquiring knowledge about the game, skills, and forecasting of move effects in finding optimality; in our context, this refers to the knowledge of how the patient responds to the drug profile both past and forecasted in optimum seeking algorithms;
  • Strategy refers to the association between a player’s move and the information available at that moment; this is fairly similar in our context denoting the controller’s optimal solution-seeking protocol and can be cooperative or non-cooperative, static or adaptive, etc;
  • Utility (or reward) is part of the optimization cost variable and for our case, this is the minimal amount of drug which maximizes the patient outcome, i.e., minimizes a relative ratio between volume growth and cell death rate.
It is necessary to explicitly represent the dynamic nature of the game theory parallelism to the multi-drug decision system, as it evolves over a period of time, i.e., the active treatment period in the patient. The current state of the tumor volume specifies the current situation of the dynamic game (dose-effect relationships), including the set of players (choice of drug cocktail), actions available to them (drug profiles expressed in amount and time interval dosage), and their utilities at this time (relative tumor volume reduction). As an example, a subclass of Markov games with multiagent sequential decision-making under uncertainties has been discussed in [30]. The decision-making process is based on a reinforcement learning (RL) principle, where the future choices of the actions are shaped by feedback of a reward function, in our case this being the therapeutic effect on the patient. The gradient play is most relevant here as it indicates a convergence of the RL scheme towards a Nash equilibrium in dynamic environments. The asymptotic behavior of such systems has been broadly discussed in [30].
When analyzing systems with limited resources, noncooperative games have been proven to be good candidates for reaching Nash equilibrium (NE). To minimize the risk for drug resistance and side effects of drug therapy in patients and improve the quality of life in cancer-treated patients, it is desirable to minimize the amount of drug or intensity of radiation profiles. If we have multiple players denoting multiple therapy profile strategies, reaching NE implies solving the problem of finding the best strategy for one player, given all other players move at optimality solution. That is, it needs to determine the actions that players should take to achieve the best outcome in response to other players’ actions. In terms of finding the best treatment protocols, finding NE implies having an adaptive protocol strategy. If the drug selection remains constant, it represents that the number of players is constant. If the drug selection also varies, then it represents that the number of players in the game changes as well. The coalition consists of multiple agents collectively acting as a virtual player to minimize a coalition cost function, defined as the sum of all agents’ local cost functions, as represented in Figure 2.
To analyze the number of feasible solutions and their convergence for biological applications one may employ the theory of Lyapunov analysis [31], particularly useful in designing continuous-time distributed NE-seeking algorithms. In this context, Figure 2 depicts a non-cooperative game with gradient play for an average consensus collecting the sum of all local cost functions. A convergence analysis for NE seeking in N-coalition games has been discussed in [32].
From a clinical perspective, the selection of drugs to be used in multi-drug combination therapy in cancer patients is crucial, because of their interconnected effects described by drug synergy and drug resistance. This translates to the choice of agents within a coalition set. Aggregative games are a special subclass of non-cooperative games where the decision process of each agent depends on the aggregate effect of all agents in the coalition. In this case, convergence is based on the monotonicity of convex functions, i.e., a gradient descent in tumor volume effect. Solutions for center-based NE seeking in such population games are presented in [33].
In this paper, we investigate the effect of a distributed NE-seeking problem with various coalition profiles, whereas the agents remain the same. In particular, we employ the model from (1) for a set of three drugs, with coalitions defined as cases in the next section.

3. Results

The following is a summary of the obtained results for various sets of the domain of feasible solutions within the coalition game optimality described in the prior section. The following settings were applied: sampling time of 1 day, optimality calculated per day with a recurrent prediction horizon of 7 days, fixed antiangiogenesis therapy as 0.171 mg/mL per day single dose weekly. Single-patient model settings were considered. A supporting immunotherapy profile has been administered as 0.2 mg/mL per day, a single dose weekly. The corresponding concentrations are depicted in Figure 3.
Three coalitions have been evaluated, for which the corresponding optimal profiles are given in Figure 4.
The results of the three coalitions are given in Figure 5, Figure 6 and Figure 7, respectively. The results suggest the model-based control methodology presented here has clinical relevance in analyzing the effect on short-term tumor growth dynamics.

4. Discussion

Carcinogenesis is progressively developing by acquiring specific collective characteristics of tumorous cells in order to persist within the tissue. Since all cell populations have an intrinsic self-renewal capacity, neoplastic cells will also rapidly increase, developing their blood supply for facilitating nutrient demands. The concepts promoting tumor growth and progression are important due to their involvement in multi-drug therapeutic applications, where drugs are developed to interfere with each of these capabilities of tumor growth. The formation of new multi-drug treatment options to target specific morphological and functional abnormal properties is only possible via earlier simulations of predicted optimal responses [34,35].
The augmented model provided here allows investigations into regional anomalous diffusion and dose-effect variability in patients undergoing multi-drug therapy for lung cancer. In fact, the principles governing the equations of the proposed model are generally applicable to any other kind of tumor tissue growth analysis. As such, the variability observed in terms of diffusion is rather small, given the limited scale of analysis and lack of micro-scale modeling of molecular binding patterns [36]. However, we foresee the potential of this model in being complementary to other meso- and micro-scale model analyses for a more comprehensive investigation. Residual times of drug molecules in sub-diffusive environment α 1 are not included since they cannot be validated. Diffusion can be validated through MRI data and observing patterns of enhanced colored bio-markers in the various areas of the tissue.
The projection between multi-drug optimal predictive control and multi-agent coalition (multi players) networked games is indeed novel in this application domain. A stability analysis following searchability and convergence to Nash equilibrium has not been performed. For this, more information is necessary, as to the clinically feasible knowledge domain of search for valid solutions while monotonicity must be guaranteed for convergence. As the research is in its infancy, the limited clinical data available does not allow the validation of such models.
From a modeling point of view, distributed games over networks including Nash equilibrium convergence are defined as stochastic problems, and as such, stochastic models may better fit the task. Stochastic models have been recently largely employed to characterize pandemic evolution [37]. When using multi-agent networked coalitions, this becomes a problem of decision-making of multiple self-interested decision-makers, where uncertainty and risk can be explicitly formulated under the Nash equilibrium problem [38]. However, in this control methodology, one must be aware of the difference between anticipative versus adaptive models for optimal search algorithms. In the case of the best treatment profile of combined drug therapy [39,40], the delineation between the two concepts is not yet defined, creating confusion in terminology. We expect that model-based predictive control plays the anticipative role [41], but the adaptation of the model itself to the new patient set of conditions must be integrated into a recurrent optimal search solution [42,43].
Using the conceptual framework that comprises the integral components of most forms of cancer, we envision significant advances in the development of mathematical modeling approaches in cancer research. The power of mathematics lies in its capacity to quantitatively describe the driving mechanisms investigated above, despite the complexity of cancer. Through mathematical formalism, we can address an abstraction of the underlying biological hypotheses, evaluate assumptions, investigate alternative solutions and make predictions that have to be tested and validated in multiple experiments. Individualized and targeted anti-cancer therapy requires the use of population-derived models despite the manual effort based on the clinician’s level of experience. Introducing quantitative approaches to predict tumor dynamics and patients’ responses based on mathematical models and prior data of the patient will provide model-driven predictions in the treatment planning process of cancer [44].

5. Conclusions

This study set out to prove the use of PKPD computational models and fractional tools to characterize the physiological processes of distribution and diffusion of drugs in the tumor. The variety of biological properties of cancer tissue provokes heterogeneous disparity of diffusion in space and time. We have proposed a mathematical approach to enhance drug absorption in the primary tumor and characterize the tumor dynamics under drug action. The biological events in healthy cancer tissue are influenced by microscale processes, taking into account the spatial distribution and the range of time. Moreover, the relation between multi-drug optimal predictive control and multi-agent coalition networked games brings novelty to this topic, anticipating the patient’s response to treatment and allowing its simulation. The principal theoretical implication of this study is that the proposed mathematical formulations are capable to define the macroscale properties of the anomalous tumor dynamics, capturing spatio-temporal anomalous drug diffusion, and identifying the specific response to therapy.

Author Contributions

Conceptualization, C.M.I. and M.G.; methodology, C.M.I.; software, M.G.; validation, M.G.; formal analysis, C.M.I.; investigation, M.G.; resources, M.G.; writing—original draft preparation, C.M.I. and M.G.; writing—review and editing, C.M.I. and M.G.; visualization, M.G.; supervision, C.M.I.; project administration, C.M.I.; funding acquisition, C.M.I. and M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Special Research Fund of Ghent University, project grant number 01J01619 and doctoral fellowship number 01D15919 (M.G.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
MPCModel Predictive Control
NENash Equilibrium
NSCLCNon-Small Cell Lung Cancer
PKPDPharmacokinetic-Pharmacodynamic
RLReinforcement Learning
SBRTStereotactic Body Radiation Therapy

References

  1. Hassan, Z.; Raza, N.; Abdel-Aty, A.-H.; Zakarya, M.; Rahman, R.U.; Yasmeen, A.; Muse, A.H.; Mahmoud, E.E. New fractal soliton solutions and sensitivity visualization for double-chain DNA model. J. Funct. Spaces 2022, 2022, 2297866. [Google Scholar] [CrossRef]
  2. Amoddeo, A. Indirect contributions to tumor dynamics in the first stage of the avascular phase. Symmetry 2020, 12, 1546. [Google Scholar] [CrossRef]
  3. Axenie, C.; Bauer, R.; Martínez, M.R. The multiple dimensions of networks in cancer: A perspective. Symmetry 2021, 13, 1559. [Google Scholar] [CrossRef]
  4. Go, J. Mathematical analysis for the effects of medicine supplies to a solid Tumor. Symmetry 2021, 13, 1988. [Google Scholar] [CrossRef]
  5. Cacace, F.; Cusimano, V.; Palumbo, P. Optimal impulsive control with application to antiangiogenic tumor therapy. IEEE Trans. Control Syst. Technol. 2020, 28, 106–117. [Google Scholar] [CrossRef]
  6. Huang, W.; Chen, L.; Kang, L.; Jin, M.; Sun, P.; Xin, X.; Gao, Z.; Bae, Y.H. Nanomedicine-based combination anticancer therapy between nucleic acids and small molecular drugs. Adv. Drug Deliv. Rev. 2017, 115, 82–97. [Google Scholar] [CrossRef]
  7. Eisenstein, M. New lung-cancer drugs extend survival times. Nature 2020, 587, S10–S12. [Google Scholar] [CrossRef]
  8. Thomas, D.S.; Cisneros, L.H.; Anderson, A.R.A.; Maley, C.C. In silico investigations of multi-drug adaptive therapy protocols. Cancers 2022, 14, 2699. [Google Scholar] [CrossRef]
  9. Rudin, C.M.; Poitier, J.T.; Byers, L.A.; Dive, C.; Dowlati, A.; George, J.; Heymach, J.V.; Johnson, J.E.; Lehman, J.M.; MacPherson, D.; et al. Molecular subtypes of small cell lung cancer: A synthesis of human and mouse model data. Nat. Rev. Cancer 2019, 19, 289–297. [Google Scholar] [CrossRef]
  10. Karachaliou, N.; Pilloto, S.; Lazzari, C.; Bria, E.; De Marini, F.; Rosell, R. Cellular and molecular biology of small cell lung cancer. an overview. Trans. Lung Cancer Res. 2016, 5, 2–15. [Google Scholar]
  11. Ghita, M.; Billiet, C.; Copot, D.; Verellen, D.; Ionescu, C.M. Parameterisation of respiratory impedance in lung cancer patients from forced oscillation lung function test. IEEE Trans. Biomed. Eng. 2022. [Google Scholar] [CrossRef]
  12. Ionescu, C.M.; Copot, C.; Verellen, D. Motion compensation for robotic lung tumour radiotherapy in remote locations: A personalised medicine approach. Acta Astronaut. 2017, 132, 59–66. [Google Scholar] [CrossRef]
  13. Ionescu, C.M. The Human Respiratory System: An Analysis of the Interplay between Anatomy, Structure, Breathing and Fractal Dynamics; Series in BioEngineering; Springer: London, UK, 2013. [Google Scholar]
  14. Ionescu, C.; Kelly, J.F. Fractional calculus for respiratory mechanics: Power law impedance, viscoelasticity, and tissue heterogeneity. Chaos Solitons Fractals 2017, 102, 433–440. [Google Scholar] [CrossRef]
  15. Ionescu, C.; Lopes, A.; Copot, D.; Machado, J.A.T.; Bates, J.H.T. The role of fractional calculus in modeling biological phenomena: A review. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 141–159. [Google Scholar] [CrossRef]
  16. Drexler, D.A.; Sapi, J.; Kovacs, L. Modeling of tumor growth incorporating the effects of necrosis and the effect of bevacizumab. Complexity 2017, 5985031. [Google Scholar] [CrossRef] [Green Version]
  17. Sapi, J.; Kovacs, L.; Drexler, D.A.; Kocsis, P.; Gajari, D.; Sapi, Z. Tumor volume estimation and quasi-continuous administration for most effective bevacizumab therapy. PLoS ONE 2015, 10, e0142190. [Google Scholar] [CrossRef] [Green Version]
  18. Prezzano, K.M.; Ma, S.J.; Hermann, G.M.; Rivers, C.I.; Gomez-Suescun, J.A.; Singh, A.K. Stereotactic body radiation therapy for non-small cell lung cancer: A review. World J Clin Oncol 2019, 10, 14–27. [Google Scholar] [CrossRef]
  19. Shields, M.D.; Marin-Acevedo, J.A.; Pellini, B. Immunotherapy for advanced non-small cell lung cancer: A decade of progress. Am. Soc. Clin. Oncol. Educ. Book 2021, 41, e105–e127. [Google Scholar] [CrossRef]
  20. Ionescu, C.M.; Ghita, M.; Copot, D.; Derom, E.; Verellen, D. A minimal PKPD interaction model for evaluating synergy effects of combined NSCLC therapies. J. Clin. Med. 2020, 9, 1832. [Google Scholar] [CrossRef]
  21. Ghita, M.; Copot, D.; Billiet, C.; Verellen, D. Ionescu, C.M. Lung cancer dynamics using fractional order impedance modeling on a mimicked lung tumor setup. J. Adv. Res. 2021, 32, 61–71. [Google Scholar] [CrossRef]
  22. Ghita, M.; Billiet, C.; Copot, D.; Verellen, D.; Ionescu, C.M. Model calibration of pharmacokineticpharmacodynamic lung tumour dynamics for anticancer therapies. J. Clin. Med. 2022, 11(4), 1006. [Google Scholar] [CrossRef] [PubMed]
  23. Haryanto, A.; Hong, K.-S. Maximum likelihood identification of Wiener-Hammerstein models. Mech. Syst. Signal Process 2003, 41, 54–70. [Google Scholar] [CrossRef]
  24. Shaikh, A.H.; Barbe, K. Study of Random Forest to Identify Wiener-Hammerstein System. IEEE Trans. Instrum. Meas. 2021, 70, 1–12. [Google Scholar] [CrossRef]
  25. Ionescu, C.M. A computationally efficient Hill curve adaptation strategy during continuous monitoring of dose-effect relation in anaesthesia. Nonlinear Dyn. 2018, 92, 843–852. [Google Scholar] [CrossRef]
  26. Ghita, M.; Copot, D.; De Keyser, R.; Ionescu, C.M. Pharmaco-impedance modelling for lung cancer therapy with predictive control. In Proceedings of the 9th International Conference on Systems and Control (ICSC), Caen, France, 24–26 November 2021; pp. 423–428. [Google Scholar]
  27. Ionescu, C.; Cajo Diaz, R.A.; Zhao, S.; Ghita, M.; Ghita, M.; Copot, D. A low computational cost, prioritized, multi-objective optimization procedure for predictive control towards cyber physical systems. IEEE Access 2020, 8, 128152–128166. [Google Scholar] [CrossRef]
  28. Wang, L. Control System Design and Implementation Using MATLAB; Springer: London, UK, 2009. [Google Scholar]
  29. Rossiter, J.A. A First Course in Predictive Control, 2nd ed.; CRC Press: Cambridge, UK, 2022. [Google Scholar]
  30. Li, T.; Peng, G.; Zhy, Q.; Basar, T. The confluence of networks, games and learning a game-theoretic framework for multiagent decision making over network. IEEE Control Syst. Mag. 2022, 42, 35–67. [Google Scholar] [CrossRef]
  31. Xu, C.; Farman, M.; Hasan, A.; Akgül, A.; Zakarya, M.; Albalawi, W.; Park, C. Lyapunov stability and wave analysis of COVID-19 omicron variant of real data with fractional operator. Alex. Eng. J. 2022, 61, 11787–11802. [Google Scholar] [CrossRef]
  32. Hu, G.; Pang, Y.; Sun, C.; Hong, Y. Distributed Nash equilibrium seeking: Continuous-time control-theoretic approaches. IEEE Control Syst. Mag. 2022, 42, 68–86. [Google Scholar] [CrossRef]
  33. Belgioioso, G.; Yi, P.; Grammatico, S.; Pavel, L. Distributed generalized Nash equilibirum seeking: An operator-theoretic perspective. IEEE Control Syst. Mag. 2022, 42, 87–102. [Google Scholar] [CrossRef]
  34. Wu, J.; Tan, Y.; Chen, Z.; Zhao, M. Data decision and drug therapy based on non-small cell lung cancer in a big data medical system in developing countries. Symmetry 2018, 10, 152. [Google Scholar] [CrossRef] [Green Version]
  35. Drexler, D.A.; Ferenci, T.; Furedi, A.; Szakacs, G.; Kovacs, L. Experimental data-driven tumor modeling for chemotherapy. IFACPapersOnline 2020, 53, 16245–16250. [Google Scholar] [CrossRef]
  36. Cacace, F.; Cusimano, V.; Germani, A.; Palumbo, P.; Papi, M. Optimal continuous-discrete linear filter and moment equations for nonlinear diffusions. IEEE Trans. Automat. Contr. 2020, 65, 3961–3976. [Google Scholar] [CrossRef]
  37. Borri, A.; Palumbo, P.; Papa, F. The stochastic approach for SIR epidemic models: Do they help to increase information from raw data? Symmetry 2022, 14, 2330. [Google Scholar] [CrossRef]
  38. Lei, J.; Shanbhag, U.V. Stochastic Nash equilibrium problems: Models, analysis, and algorithms. IEEE Control Syst. Mag. 2022, 42, 103–124. [Google Scholar] [CrossRef]
  39. Capasso, A.; Lang, J.; Pitts, T.M.; Jordan, K.R.; Lieu, C.H.; Diamond, J.R.; Kopetz, S.; Barbee, J.; Peterson, J.; Freed, B.M.; et al. Characterization of immune responses to anti-PD-1 mono- and combination therapy in hematopoietic humanized mice implanted with tumor xenographs. J. Immunother. Cancer 2019, 7, 37. [Google Scholar] [CrossRef]
  40. Sculier, J.P.; Paesmans, M.; Bureau, G.; Dabois, G.; Libert, P.; Van Cutsem, G.V.O.; Berchier, M.C.; Ries, F.; Michel, J. Multiple-drug weekly chemotherapy versus standard combination regimen in small-cell lung cancer: A phase III randomized study conducted by the European Lung Cancer Working Party. J. Clin Oncol. 1993, 11, 1858–1865. [Google Scholar] [CrossRef]
  41. Kovacs, L.; Szeles, A.; Sapi, J.; Drexler, D.A.; Rudas, I.; Harmati, I.; Sapi, Z. Model-based angiogenic inhibition of tumor growth using modern robust control method. Comput. Methods Programs Biomed. 2014, 114, e98–e110. [Google Scholar] [CrossRef] [Green Version]
  42. Almetwally, E.M.; Jawa, T.M.; Sayed-Ahmed, N.; Park, C.; Zakarya, M.; Dey, S. Analysis of unit-Weibull based on progressive type-II censored with optimal scheme. Alex. Eng. J. 2023, 63, 321–338. [Google Scholar] [CrossRef]
  43. Khalid, A.; Rehan, A.; Nisar, K.S.; Abdel-Aty, A.-H.; Zakarya, M. Splines solutions of higher-order BVPs that ARISE in consistent magnetized force field. Fractals 2022, 30, 2240043. [Google Scholar] [CrossRef]
  44. Hahn, J.-O.; Inan, O.T. Physiological closed-loop control in critical care: Opportunities for innovations. Prog. Biomed. Eng. 2022, 4, 033001. [Google Scholar] [CrossRef]
Figure 1. Surface response E for α = 1 .
Figure 1. Surface response E for α = 1 .
Symmetry 15 00051 g001
Figure 2. Concept of coalition game NE optimality problem, with Players representing various therapy protocols, and different drug combination therapy.
Figure 2. Concept of coalition game NE optimality problem, with Players representing various therapy protocols, and different drug combination therapy.
Symmetry 15 00051 g002
Figure 3. Concentration profiles as a result of infusion dose in u a and u i .
Figure 3. Concentration profiles as a result of infusion dose in u a and u i .
Symmetry 15 00051 g003
Figure 4. Dosimetry profiles of SBRT u r for three cases of coalitions.
Figure 4. Dosimetry profiles of SBRT u r for three cases of coalitions.
Symmetry 15 00051 g004
Figure 5. Coalition 1: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Figure 5. Coalition 1: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Symmetry 15 00051 g005
Figure 6. Coalition 2: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Figure 6. Coalition 2: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Symmetry 15 00051 g006
Figure 7. Coalition 3: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Figure 7. Coalition 3: results for various diffusion coefficients in terms of dose-effect surface and residual active tumor volume ( α = 0.08 (left), α = 0.03 (middle), α = 0.007 (right)).
Symmetry 15 00051 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ionescu, C.M.; Ghita, M. Model-Based Regional Control with Anomalous Diffusion of Multi-Drug Combined Cancer Therapy for Volume Predictions. Symmetry 2023, 15, 51. https://doi.org/10.3390/sym15010051

AMA Style

Ionescu CM, Ghita M. Model-Based Regional Control with Anomalous Diffusion of Multi-Drug Combined Cancer Therapy for Volume Predictions. Symmetry. 2023; 15(1):51. https://doi.org/10.3390/sym15010051

Chicago/Turabian Style

Ionescu, Clara Mihaela, and Maria Ghita. 2023. "Model-Based Regional Control with Anomalous Diffusion of Multi-Drug Combined Cancer Therapy for Volume Predictions" Symmetry 15, no. 1: 51. https://doi.org/10.3390/sym15010051

APA Style

Ionescu, C. M., & Ghita, M. (2023). Model-Based Regional Control with Anomalous Diffusion of Multi-Drug Combined Cancer Therapy for Volume Predictions. Symmetry, 15(1), 51. https://doi.org/10.3390/sym15010051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop