RobertLampel at 11:30, 4 September 2026

2026-09-04T11:30:42Z

RobertLampel: /* References */

2026-09-03T11:50:00Z

References

← Older revision		Revision as of 11:50, 3 September 2026
Line 84:		Line 84:


			[[Category:MIOCP]]
	[[Category:ODE model]]		[[Category:ODE model]]
	[[Category:Biomedical model]]		[[Category:Biomedical model]]
	[[Category:Singular arc]]		[[Category:Singular arc]]

RobertLampel: RobertLampel moved page Tumor Anti Anginogenesis to Tumor Anti Angiogenesis: Misspelled title

2026-09-03T11:47:54Z

RobertLampel moved page Tumor Anti Anginogenesis to Tumor Anti Angiogenesis: Misspelled title

← Older revision	Revision as of 11:47, 3 September 2026
(No difference)

RobertLampel: Created page with "{{Dimensions |nd = 1 |nx = 3 |nw = 1 }} The '''tumor anti-angiogenesis problem''' is a classical benchmark in biomedical optimal control. This description is taken from [1]. It is based on the mathematical model of Hahnfeldt et al. [2] for tumor growth under angiogenic signalling, as studied from an optimal control perspective by Ledzewicz and Schättler [3]. Angiogenesis is the process by which tumors rec..."

2026-09-03T11:45:23Z

Created page with "{{Dimensions |nd = 1 |nx = 3 |nw = 1 }} The '''tumor anti-angiogenesis problem''' is a classical benchmark in biomedical optimal control. This description is taken from [1]. It is based on the mathematical model of Hahnfeldt et al. [2] for tumor growth under angiogenic signalling, as studied from an optimal control perspective by Ledzewicz and Schättler [3]. Angiogenesis is the process by which tumors rec..."

New page

{{Dimensions
|nd = 1
|nx = 3
|nw = 1
}}

The '''tumor anti-angiogenesis problem''' is a classical benchmark in biomedical optimal control. This description is taken from [[#OCPjl | [1]]].

It is based on the mathematical model of Hahnfeldt et al. [[#Hahnfeldt|[2]]] for tumor growth under angiogenic signalling, as studied from an optimal control perspective by Ledzewicz and Schättler [[#Ledzewicz|[3]]].
Angiogenesis is the process by which tumors recruit a network of blood vessels to sustain their growth; anti-angiogenic therapy targets this process by administering inhibitor agents.
The two biological states are the tumor volume and the carrying capacity of the vascular endothelium (which sets the maximum supportable tumor size).
A third state tracks the cumulative inhibitor dose, enforcing a total-dose constraint.
The control is the instantaneous inhibitor dose rate.
The objective is to minimise the tumor volume at the final time.
The problem is known to exhibit a rich optimal control structure, including bang–bang arcs, singular arcs, and boundary arcs corresponding to the dose constraint.

== Mathematical formulation ==

<math>
\begin{array}{lll}
\displaystyle \min_{u} && p(t_\mathrm{f}) \\[6pt]
\text{subject to} \\
\quad \dot{p}(t) &=& -\xi\, p(t)\,\ln\!\left(\dfrac{p(t)}{q(t)}\right), \\[6pt]
\quad \dot{q}(t) &=& b\, p(t) - \mu\, q(t) - \left(d\, p(t)^{2/3} + e\, u(t)\right) q(t), \\[6pt]
\quad \dot{y}(t) &=& u(t), \\
\quad (p(0),\, q(0),\, y(0)) &=& (p_0,\, q_0,\, 0), \\
\quad y(t_\mathrm{f}) &\leq& A, \\
\quad u(t) &\in& [0,\, a] \quad \forall\, t \in [0, t_\mathrm{f}],
\end{array}
</math>


where <math>p(t)</math> is the tumor volume (mm3), <math>q(t)</math> is the endothelial carrying capacity (mm3), and <math>y(t)</math> is the cumulative inhibitor dose (mg/kg) administered up to time <math>t</math>.
The tumor dynamics follow the Gompertz law <math>\dot{p} = -\xi\,p\ln(p/q)</math>, which models logistic-like growth bounded by the carrying capacity <math>q</math>.
The carrying capacity decreases due to natural endothelial cell loss (<math>\mu\,q</math>), inhibitor-induced endothelial cell loss (<math>e\,u\,q</math>), and stimulated shedding proportional to tumor surface area (<math>d\,p^{2/3}\,q</math>), while increasing due to endothelial recruitment by the tumor (<math>b\,p</math>).
The state <math>y(t)</math> enforces the total-dose constraint <math>y(t_\mathrm{f}) \leq A</math> as an isoperimetric inequality.

== Parameters ==
These fixed values are used within the model:

{| border="1" align="center" cellpadding="5" cellspacing="0"
|- bgcolor=#c7c7c7
! Symbol !! Value !! Description
|-
| align=center | <math>\xi</math> || align=right | 0.084 || Tumor growth rate (day−1)
|-
| align=center | <math>b</math> || align=right | 5.85 || Endothelial birth rate coefficient (day−1)
|-
| align=center | <math>\mu</math> || align=right | 0.02 || Endothelial natural loss rate (day−1)
|-
| align=center | <math>d</math> || align=right | 0.00873 || Inhibition coefficient for endothelial shedding (mm−2 day−1)
|-
| align=center | <math>e</math> || align=right | 0.66 || Inhibitor effectiveness on endothelium (kg mg−1 day−1)
|-
| align=center | <math>a</math> || align=right | 75 || Maximum inhibitor dose rate (mg kg−1 day−1)
|-
| align=center | <math>A</math> || align=right | 15 || Maximum cumulative inhibitor dose (mg kg−1)
|-
| align=center | <math>p_0</math> || align=right | 12000 || Initial tumor volume (mm3)
|-
| align=center | <math>q_0</math> || align=right | 15000 || Initial endothelial carrying capacity (mm3)
|-
| align=center | <math>t_\mathrm{f}</math> || align=right | 1.2 || Horizon of the control problem (days)
|}

== Reference Solutions ==

Here is one local solution to the above control problem.

<gallery caption="Reference solution plots" widths="500px" heights="300px" perrow="1">
Image: Tumor_Anti_Angiogenesis.png | Tumor volume, carrying capacity, cumulative dose, and optimal inhibitor schedule for a reference solution.
</gallery>

== Miscellaneous and Further Reading ==
The equilibrium <math>p^* = q^*</math> corresponds to a tumor at maximum sustainable size.
Without treatment, <math>p(t)</math> converges to this steady state.
Ledzewicz and Schättler [[#Ledzewicz|[3]]] characterise the optimal control using Pontryagin's minimum principle and show that the optimal policy consists of a sequence of full-dose (<math>u = a</math>), no-dose (<math>u = 0</math>), and singular arcs, depending on the initial conditions and remaining dose budget.
This formulation and detailed analysis can be found in [[#OCPjl|[1]]], [[#Hahnfeldt|[2]]], and [[#Ledzewicz|[3]]].

== References ==
[1] Caillau, J.-B., Cots, O., Gergaud, J., & Martinon, P. OptimalControlProblems.jl: a collection of optimal control problems with ODE's in Julia. https://github.com/control-toolbox/OptimalControlProblems.jl 
[2] Hahnfeldt, P., Panigrahy, D., Folkman, J., & Hlatky, L. (1999). Tumor development under angiogenic signaling: A dynamical theory of tumor growth, treatment response, and postvascular dormancy. Cancer Research, 59(19), 4770–4775. 
[3] Ledzewicz, U., & Schättler, H. (2007). Antiangiogenic therapy in cancer treatment as an optimal control problem. SIAM Journal on Control and Optimization, 46(3), 1052–1079. 

[[Category:ODE model]]
[[Category:Biomedical model]]
[[Category:Singular arc]]

Tumor Anti Angiogenesis - Revision history

RobertLampel at 11:30, 4 September 2026

RobertLampel: /* References */

RobertLampel: RobertLampel moved page Tumor Anti Anginogenesis to Tumor Anti Angiogenesis: Misspelled title