Oscillatory dynamics in a model of vascular tumour growth--implications for chemotherapy.
BACKGROUND
Investigations of solid tumours suggest that vessel occlusion may occur when increased pressure from the tumour mass is exerted on the vessel walls. Since immature vessels are frequently found in tumours and may be particularly sensitive, such occlusion may impair tumour blood flow and have a negative impact on therapeutic outcome. In order to study the effects that occlusion may have on tumour growth patterns and therapeutic response, in this paper we develop and investigate a continuum model of vascular tumour growth.


RESULTS
By analysing a spatially uniform submodel, we identify regions of parameter space in which the combination of tumour cell proliferation and vessel occlusion give rise to sustained temporal oscillations in the tumour cell population and in the vessel density. Alternatively, if the vessels are assumed to be less prone to collapse, stable steady state solutions are observed. When spatial effects are considered, the pattern of tumour invasion depends on the dynamics of the spatially uniform submodel. If the submodel predicts a stable steady state, then steady travelling waves are observed in the full model, and the system evolves to the same stable steady state behind the invading front. When the submodel yields oscillatory behaviour, the full model produces periodic travelling waves. The stability of the waves (which can be predicted by approximating the system as one of lambda-omega type) dictates whether the waves develop into regular or irregular spatio-temporal oscillations. Simulations of chemotherapy reveal that treatment outcome depends crucially on the underlying tumour growth dynamics. In particular, if the dynamics are oscillatory, then therapeutic efficacy is difficult to assess since the fluctuations in the size of the tumour cell population are enhanced, compared to untreated controls.


CONCLUSIONS
We have developed a mathematical model of vascular tumour growth formulated as a system of partial differential equations (PDEs). Employing a combination of numerical and analytical techniques, we demonstrate how the spatio-temporal dynamics of the untreated tumour may influence its response to chemotherapy.


METHODS
This manuscript was reviewed by Professor Zvia Agur and Professor Marek Kimmel.
Stamper
I J
IJ
Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK.
Owen
M R
MR
Maini
P K
PK
Byrne
H M
HM
D000818:Animals; D006801:Humans; D008962:Models, Theoretical; D009369:Neoplasms / Q000188:drug therapy* / Q000473:pathology*; D009389:Neovascularization, Pathologic
Annu Rev Biomed Eng. 1999;1:241-63
11701489
J Theor Biol. 1997 Nov 7;189(1):63-80
9398504
Science. 2005 Jan 7;307(5706):58-62
15637262
Cancer Cell. 2004 Dec;6(6):529-31
15607955
Mol Pharmacol. 1994 Apr;45(4):649-56
8183243
J Theor Biol. 2003 Jan 21;220(2):201-13
12468292
J Neurol Sci. 2003 Dec 15;216(1):1-10
14607296
Nature. 2004 Feb 19;427(6976):695
14973470
J Theor Biol. 2006 Feb 21;238(4):817-27
16384586
Bull Math Biol. 2005 Jan;67(1):79-99
15691540
Math Med Biol. 2004 Mar;21(1):1-34
15065736
Neoplasia. 2000 Jul-Aug;2(4):325-38
11005567
J Cell Physiol. 1992 May;151(2):386-94
1572910
Cancer Res. 1999 Aug 1;59(15):3776-82
10446995
Int J Radiat Oncol Biol Phys. 2004 Jun 1;59(2):572-8
15145178
Bull Math Biol. 2003 Jul;65(4):609-40
12875336
Cancer Res. 1988 May 15;48(10):2641-58
3282647
Eur J Cancer. 2005 Jan;41(1):159-67
15618001
Bull Math Biol. 2007 Jul;69(5):1747-76
17333419
C R Biol. 2004 Nov;327(11):995-1008
15628222
Angiogenesis. 2002;5(3):203-14
12831061
Cancer Res. 1986 May;46(5):2295-9
3697976