Discrete and Continuous Dynamical Systems - B
May 2019 , Volume 24 , Issue 5
New horizons in optimal control and biomathematics: A tribute to Urszula Ledzewicz, Helmut Maurer and Heinz Schaettler
Select all articles
In this paper, a non-trivial generalization of a mathematical model put forward in [
This study investigates how optimal control theory may be used to delay the onset of chemotherapy resistance in tumours. An optimal control problem with simple tumour dynamics and an objective functional explicitly penalising drug resistant tumour phenotype is formulated. It is shown that for biologically relevant parameters the system has a single globally attracting positive steady state. The existence of singular arc is then investigated analytically under a very general form of the resistance penalty in the objective functional. It is shown that the singular controls are of order one and that they satisfy Legendre-Clebsch condition in a subset of the domain. A gradient method for solving the proposed optimal control problem is then used to find the control minimising the objective. The optimal control is found to consist of three intervals: full dose, singular and full dose. The singular part of the control is essential in delaying the onset of drug resistance.
Morphogenesis, the shaping of an organism, is a complex biological process accomplished through an well organized interplay between growth, differentiation and cell movement.It is still today one of the major outstanding problems in the biological sciences. Pattern formation has been well-addressed in the literature with the development of many mathematical models including the famous reaction-diffusion ones. We here take a different approach, introducing a controlled cellular automaton in order to model the signal molecules, known as growth factors, that convey information from one cell to another during an organism's development and help maintain the viability of the adult. This control represents extracellular structures that have been associated with the regulation of stem cell proliferation and are called fractones. In this paper we introduce two co-evolving automata, one describing the perturbed diffusion of growth factors and one accounting for the rules of basic cellular functions (proliferation, differentiation, migration and apoptosis). Fractones are introduced as an external input to control the shaping of multi-cellular organisms; we analyze their influence on the emerging shape. We illustrate our theory with 2 and 3 dimensional simulations. This work presents the foundation upon which to develop cellular automata as a tool to simulate the morphodynamics in embryonic development.
Nonlinear dynamics of a reaction-diffusion equation with delay is studied with numerical simulations in 1D and 2D cases. Homogeneous in space solutions can manifest time oscillations with period doubling bifurcations and transition to chaos. Transition between two regions with homogeneous oscillations is provided by quasi-waves, propagating solutions without regular structure and often with complex aperiodic oscillations. Dynamics of space dependent solutions is described by a combination of various waves, e.g., bistable, monostable, periodic and quasi-waves.
Control systems involving unknown parameters appear a natural framework for applications in which the model design has to take into account various uncertainties. In these circumstances the performance criterion can be given in terms of an average cost, providing a paradigm which differs from the more traditional minimax or robust optimization criteria. In this paper, we provide necessary optimality conditions for a nonrestrictive class of optimal control problems in which unknown parameters intervene in the dynamics, the cost function and the right end-point constraint. An important feature of our results is that we allow the unknown parameters belonging to a mere complete separable metric space (not necessarily compact).
In the Hes1 gene expression system the protein (present as dimers) bounds to the promoter of its own DNA blocking transcription of its mRNA. This negative feedback leads to an oscillatory behavior, which is observed experimentally. Classical mathematical model of this system consists of two ordinary differential equations with discrete time delay in the term reflecting transcription. However, transcription takes place in the nucleus while translation occurs in the cytoplasm. This means that the delay present in the system is larger than transcription time. Moreover, in reality it is not discrete but distributed around some mean value. In this paper we present the model of the Hes1 gene expression system and discuss similarities and differences between the model with discrete and distributed delays. It turns out that in the case of distributed delays the region of stability of the steady state is larger than in the case of discrete delay. We also derive conditions that guarantee stability of the steady state for particular delay distributions.
We study mathematical properties of a model describing growth of primary brain tumours called low-grade gliomas (LGGs) and their response to chemotherapy. The motivation for considering this particular type of cancer is its large impact on society. LGGs affect mainly young adults and eventually result in death, despite the tumour growth rate being slow. The model studied consists of two non-autonomous ordinary differential equations and is a generalised version of the model proposed by Bogdańska et al. (Math. Biosci. 2017). We discuss the stability of stationary states, prove global stability of tumour-free steady state and, in some cases, justify the existence of periodic solutions. Assuming that chemotherapy effectiveness remains constant in time, we provide analytical estimates and calculate minimal doses of the drug that should eliminate the tumour for particular patients with LGGs.
In this paper we propose a model of tuberculosis (TB) transmission in a heterogeneous population consisting of two different subpopulations, like homeless and non-homeless people. We use the criss-cross model to describe the illness dynamics. This criss-cross model is based on the simple SIS model with constant inflow into both subpopulations and bilinear transmission function. We find conditions for the existence and local stability of stationary states (disease-free and endemic) and fit the model to epidemic data from Warmian-Masurian Province of Poland. Basic reproduction number
We give a new proof of the maximum principle for optimal control problems with running state constraints. The proof uses the so-called method of
In the present paper we propose and study a simple model of collagen remodeling occurring in latter stage of tendon healing process. The model is an integro-differential equation describing the possibility of an alignment of collagen fibers in a finite time. We show that the solutions may either exist globally in time or blow-up in a finite time depending on initial data. The latter behavior can be related to the healing of injury without the scar formation in a finite time: a full alignment of collagen fibers. We believe that the present model is an essential ingredient of the full description of collagen remodeling.
We discuss and compare numerical methods to solve singular optimal control problems by the direct method. Our discussion is illustrated by an Autonomous Underwater Vehicle (AUV) problem with state constraints. For this problem, we test four different approaches to solve numerically our problem via the direct method. After discretizing the optimal control problem we solve the resulting optimization problem with (ⅰ) A Mathematical Programming Language (
Most clinical trials with combination therapy fail. One of the reasons is that not enough forethought is given to the interaction between the different agents, as well as the potential negative side-effects that may arise in the combined therapy. In the present paper we consider a generic cancer model with combination therapy consisting of chemotherapy agent
We consider a nonlinear system of differential equations describing a process of the psoriasis treatment. Its phase variables are the concentrations of T-lymphocytes, keratinocytes and dendritic cells. A scalar bounded control is introduced into this system to reflect the medication dosage. For such a control system, on a given time interval the minimization problem of the Bolza type functional is stated. Its terminal term is the concentration of keratinocytes at the terminal time, and its integral term is the product of the non-negative weighting coefficient with the total cost of the psoriasis treatment. This cost is linear in the control and proportional to the concentration of keratinocytes. For the analysis of such a problem, the Pontryagin maximum principle is used. As a result of this analysis, it is shown that if the weighting coefficient is zero, then the corresponding optimal control can contain a singular arc. We establish that it is a chattering control, and therefore does not make much sense as a type of a medical treatment. If the weighting coefficient is positive, then the corresponding optimal control is bang-bang, and it can be presented as a type of psoriasis treatment. In addition, when this coefficient tends to zero, such optimal controls can be considered as chattering approximations. Therefore, the convergence of these optimal controls, the corresponding optimal solutions of the original system, and the minimum values of the functional are studied. The obtained theoretical results are illustrated by numerical calculations and the corresponding conclusions are made.
In the paper, we derive a maximum principle for a Bolza problem described by an integro-differential equation of Volterra type. We use the Dubovitskii-Milyutin approach.
We give answer to an open question by proving a sufficient optimality condition for state-linear optimal control problems with time delays in state and control variables. In the proof of our main result, we transform a delayed state-linear optimal control problem to an equivalent non-delayed problem. This allows us to use a well-known theorem that ensures a sufficient optimality condition for non-delayed state-linear optimal control problems. An example is given in order to illustrate the obtained result.
We analyze the structure of optimal protocols for a mathematical model of tumor anti-angiogenic treatment. The control represents the concentration of the agent and we consider the problem to administer an a priori given total amount of agents in order to achieve a minimum tumor volume/maximum tumor reduction. In earlier work, this problem was studied with a log-kill type pharmacodynamic model for drug effects which does not account for saturation of the drug concentration. Here we study the effect of incorporating a Michaelis-Menten (MM) or
This article addresses the problem of controlling a constrained, continuous–time, nonlinear system through Model Predictive Control (MPC). In particular, we focus on methods to efficiently and accurately solve the underlying optimal control problem (OCP). In the numerical solution of a nonlinear OCP, some form of discretization must be used at some stage. There are, however, benefits in postponing the discretization process and maintain a continuous-time model until a later stage. This is because that way we can exploit additional freedom to select the number and the location of the discretization node points.We propose an adaptive time–mesh refinement (AMR) algorithm that iteratively finds an adequate time–mesh satisfying a pre–defined bound on the local error estimate of the obtained trajectories. The algorithm provides a time–dependent stopping criterion, enabling us to impose higher accuracy in the initial parts of the receding horizon, which are more relevant to MPC. Additionally, we analyze the conditions to guarantee closed–loop stability of the MPC framework using the AMR algorithm. The numerical results show that the proposed AMR strategy can obtain solutions as fast as methods using a coarse equidistant–spaced mesh and, on the other hand, as accurate as methods using a fine equidistant–spaced mesh. Therefore, the OCP can be solved, and the MPC law obtained, faster and/or more accurately than with discrete-time MPC schemes using equidistant–spaced meshes.
We consider a generational and continuous-time two-phase model of the cell cycle. The first model is given by a stochastic operator, and the second by a piecewise deterministic Markov process. In the second case we also introduce a stochastic semigroup which describes the evolution of densities of the process. We study long-time behaviour of these models. In particular we prove theorems on asymptotic stability and sweeping. We also show the relations between both models.
We analyze mathematical models for cancer chemotherapy under tumor heterogeneity as optimal control problems. Tumor heterogeneity is incorporated into the model through a potentially large number of states which represent sub-populations of tumor cells with different chemotherapeutic sensitivities. In this paper, a Michaelis-Menten type functional form is used to model the pharmacodynamic effects of the drug concentrations. In the objective, a weighted average of the tumor volume and the total amounts of drugs administered (taken as an indirect measurement for the side effects) is minimized. Mathematically, incorporating a Michaelis-Menten form creates a nonlinear structure in the controls with partial convexity properties in the Hamiltonian function for the optimal control problem. As a result, optimal controls are continuous and this fact can be utilized to set up an efficient numerical computation of extremals. Second order Jacobi type necessary and sufficient conditions for optimality are formulated that allow to verify the strong local optimality of numerically computed extremal controlled trajectories. Examples which illustrate the cases of both strong locally optimal and non-optimal extremals are given.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]