Travelling wave solutions of the reaction-diffusion mathematical model of glioblastoma growth: An Abel equation based approach

Tiberiu Harko; Man Kwong Mak

doi:10.3934/mbe.2015.12.41

Article Contents

2015, Volume 12, Issue 1: 41-69. Doi: 10.3934/mbe.2015.12.41

This issue Previous Article A double age-structured model of the co-infection of tuberculosis and HIV Next Article Dynamics of competitive systems with a single common limiting factor

Travelling wave solutions of the reaction-diffusion mathematical model of glioblastoma growth: An Abel equation based approach

Tiberiu Harko^1, and
Man Kwong Mak^2,

1.
Department of Mathematics, University College London, Gower Street, London WC1E 6BT
2.
Department of Computing and Information Management, Hong Kong Institute of Vocational Education, Chai Wan, Hong Kong

Received: April 30, 2014

Revised: August 31, 2014

Published: November 2014

Abstract Related Papers Cited by

Abstract

We consider quasi-stationary (travelling wave type) solutions to a nonlinear reaction-diffusion equation with arbitrary, autonomous coefficients, describing the evolution of glioblastomas, aggressive primary brain tumors that are characterized by extensive infiltration into the brain and are highly resistant to treatment. The second order nonlinear equation describing the glioblastoma growth through travelling waves can be reduced to a first order Abel type equation. By using the integrability conditions for the Abel equation several classes of exact travelling wave solutions of the general reaction-diffusion equation that describes glioblastoma growth are obtained, corresponding to different forms of the product of the diffusion and reaction functions. The solutions are obtained by using the Chiellini lemma and the Lemke transformation, respectively, and the corresponding equations represent generalizations of the classical Fisher--Kolmogorov equation. The biological implications of two classes of solutions are also investigated by using both numerical and semi-analytical methods for realistic values of the biological parameters.

Keywords:

Mathematics Subject Classification: Primary: 35C07, 34A05; Secondary: 92C37.

Citation:

Related Papers

Cited by

References

[1]	M. S. Ahluwalia, J. de Groot, W. M. Liu and C. L. Gladson, Targeting SRC in glioblastoma tumors and brain metastases: Rationale and preclinical studies, Cancer Letters, 298 (2010), 139-149.doi: 10.1016/j.canlet.2010.08.014.
[2]	E. C. Alvord, Jr., Patterns of growth of gliomas, American Journal of Neuroradiology, 16 (1995), 1013-1017.
[3]	M. Bellini, R. Deza and N. Giovanbattista, Exact travelling annular waves in generalized reaction-diffusion equations, Physics Letters A, 232 (1997), 200-206.doi: 10.1016/S0375-9601(97)00360-5.
[4]	T. Brugarino, Exact solutions of nonlinear differential equations using the Abelian equation of the first type, Il Nuovo Cimento, 119 (2004), 975-982.
[5]	A. Chiellini, Sull'integrazione dell'equazione differenziale $y^{'} + Py^ 2 +Qy^ 3 =0$, Bollettino dell Unione Matematica Italiana, 10 (1931), 301-307.
[6]	T. Demuth and M. E. Berens, Molecular mechanisms of glioma cell migration and invasion, Journal of Neuro-Oncology, 70 (2004), 217-228.doi: 10.1007/s11060-004-2751-6.
[7]	S. E. Eikenberry, T. Sankar, M. C. Preul, E. J. Kostelich, C. J. Thalhauser and Y. Kuang, Virtual glioblastoma: Growth, migration and treatment in a three-dimensional mathematical model, 2009, Cell Prolif., 42 (2009), 511-528.
[8]	R. A. Fisher, The Genetical Theory of Natural Selection, Oxford, Oxford University Press, 1999.
[9]	R. A. Fisher, The wave of advance of advantageous genes, Ann. Eugenics, 7 (1937), 355-369.doi: 10.1111/j.1469-1809.1937.tb02153.x.
[10]	N. Follain, J.-M. Valleton, L. Lebrun, B. Alexandre, P. Schaetzel, M. Metayer and S. Marais, Simulation of kinetic curves in mass transfer phenomena for a concentration-dependent diffusion coefficient in polymer membranes, Journal of Membrane Science, 349 (2010), 195-207.doi: 10.1016/j.memsci.2009.11.044.
[11]	P. Gerlee and S. Nelander, The Impact of Phenotypic Switching on Glioblastoma Growth and Invasion, PLOS Computational Biology, 8 (2012), 1002556, 12pp.doi: 10.1371/journal.pcbi.1002556.
[12]	P. Gerlee and S. Nelander, Travelling wave analysis of a mathematical model of glioblastoma growth, preprint, arXiv:1305.5036.
[13]	S. G. Giatili and G. S. Stamatakos, A detailed numerical treatment of the boundary conditions imposed by the skull on a diffusion-reaction model of glioma tumor growth. Clinical validation aspects, Applied Mathematics and Computation, 218 (2012), 8779-8799.doi: 10.1016/j.amc.2012.02.036.
[14]	B. H. Gilding and R. Kersner, The Characterization of Reaction-Convection-Diffusion Processes by Travelling Waves, Journal of differential equations, 124 (1996), 27-79.doi: 10.1006/jdeq.1996.0002.
[15]	B. H. Gilding and R. Kersner, Travelling Waves in Nonlinear Diffusion Convection Reaction, Birkhäuser, 2004.doi: 10.1007/978-3-0348-7964-4.
[16]	V. V. Golubev, Vorlesungen über Differentialgleichungen im Komplexen, Deutsch. Verlag Wissenschaft., Berlin, 1958.
[17]	T. Harko and M. K. Mak, Relativistic dissipative cosmological models and Abel differential equation, Computers and Mathematics with Applications, 46 (2003), 849-853.doi: 10.1016/S0898-1221(03)90147-7.
[18]	T. Harko, F. S. N. Lobo and M. K. Mak, Exact analytical solutions of the Susceptible-Infected- Recovered (SIR) epidemic model and of the SIR model with equal deaths and births, Applied Mathematics and Computation, 236 (2014), 184-194.doi: 10.1016/j.amc.2014.03.030.
[19]	T. Harko, F. S. N. Lobo and M. K. Mak, A Chiellini type integrability condition for the generalized first kind Abel differential equation, Universal Journal of Applied Mathematics, 1 (2013), 101-104.
[20]	T. Harko, F. S. N. Lobo and M. K. Mak, A class of exact solutions of the Liénard type ordinary non-linear differential equation, Journal of Engineering Mathematics, accepted for publication, arXiv:1302.0836.
[21]	H. L. P. Harpold, E. C. Alvord, Jr. and K. R. Swanson, The Evolution of Mathematical Modeling of Glioma Proliferation and Invasion, J Neuropathol Exp Neurol, 66 (2007), 1-9.doi: 10.1097/nen.0b013e31802d9000.
[22]	H. Hatzikirou, L. Brusch, K. Schaller, M. Simon and A. Deutsch, Prediction of traveling front behavior in a lattice-gas cellular automaton model for tumor invasion, Computers and Mathematics with Applications, 59 (2010), 2326-2339.doi: 10.1016/j.camwa.2009.08.041.
[23]	H. Hatzikirou, D. Basanta, M. Simon, K. Schaller and A. Deutsch, "Go or Grow": the key to the emergence of invasion in tumour progression?, Mathematical Medicine and Biology, 29 (2012), 49-65.doi: 10.1093/imammb/dqq011.
[24]	P. Hinow, P. Gerlee, L. J. McCawley, V. Quaranta, M. Ciobanu, S. Wang, J. M. Graham, B. P. Ayati, J. Claridge, K. R. Swanson, M. Loveless and A. R. A. Anderson, A Spatial Model of Tumor-Host Interaction: Application of Chemotherapy, Math. Biosci. Eng., 6 (2009), 521-546.doi: 10.3934/mbe.2009.6.521.
[25]	S. Jbabdi, E. Mandonnet, H. Duffau, L. Capelle, K. R. Swanson, M. Pélégrini-Issac, R. Guillevin and H. Benali, Simulation of anisotropic growth of low-grade gliomas using diffusion tensor imaging, Magnetic Resonance in Medicine, 54 (2005), 616-624.doi: 10.1002/mrm.20625.
[26]	S. Joannès, L. Mazé and A. R. Bunsell, A concentration-dependent diffusion coefficient model for water sorption in composite, Composite Structures, 108 (2014), 111-118.
[27]	D. R. Johnson and B. P. O'Neill, Glioblastoma survival in the United States before and during the temozolomide era, Journal of Neurooncology, 107 (2012), 359-364.doi: 10.1007/s11060-011-0749-4.
[28]	M. Kinoshita, N. Hashimoto, T. Goto, N. Kagawa, H. Kishima, S. Izumoto, H. Tanaka, N. Fujita and T. Yoshimine, NeuroImage, 43 (2008), 29-35.
[29]	E. Kamke, Differentialgleichungen: Lösungsmethoden und Lösungen, Chelsea, New York, 1959.
[30]	A. Kolmogorov, I. Petrovskii and N. Piscounov, Étude de l'équation de la diffusion avec croissance de la quantité de matiére et son application a un problème biologique, Bulletin de l'Université d'État à Moscou, Série Internationale, Section A Mathématiques et Mécanique, 1 (1937), 1-25.
[31]	H. Lemke, Über eine von R. Liouville untersuchte Differentialgleichung erster Ordnung, Sitzungs. Berl. Math. Ges., 18 (1920), 26-31.
[32]	C. K. N. Li, The glucose distribution in 9L rat brain multicell tumour spheroids and its effect on cell necrosis, Cancer, 50 (1982), 2066-2073.
[33]	D. N. Louis, H. Ohgaki, O. D. Wiestler, W. K. Cavenee, P. C. Burger, A. Jouvet, B. W. Scheithauer and P. Kleihues, The 2007 WHO Classification of Tumours of the Central Nervous System, Acta Neuropathol., 114 (2007), 97-109.doi: 10.1007/s00401-007-0243-4.
[34]	M. K. Mak, H. W. Chan and T. Harko, Solutions generating technique for Abel type non-linear ordinary differential equations, Computers and Mathematics with Applications, 41 (2001), 1395-1401.doi: 10.1016/S0898-1221(01)00104-3.
[35]	M. K. Mak and T. Harko, New method for generating general solution of Abel differential equations, Computers and Mathematics with Applications, 43 (2002), 91-94.doi: 10.1016/S0898-1221(01)00274-7.
[36]	S. C. Mancas and H. C. Rosu, Integrable dissipative nonlinear second order differential equations via factorizations and Abel equations, Physics Letters, A 377 (2013), 1234-1238.doi: 10.1016/j.physleta.2013.04.024.
[37]	S. C. Mancas and H. C. Rosu, Integrable Ermakov-Pinney equations with nonlinear Chiellini 'damping', preprint, arXiv:1301.3567.
[38]	S. C. Mancas, G. Spradlin and H. Khanal, Weierstrass traveling wave solutions for dissipative Benjamin, Bona, and Mahony (BBM) equation, Journal of Mathematical Physics, 54 (2013), 081502, 15pp.doi: 10.1063/1.4817342.
[39]	E. Mandonnet, J.-Y. Delattre, M.-L. Tanguy, K. R. Swanson, A. F. Carpentier, H. Duffau, P. Cornu, R. Van Effenterre, E. C. Alvord. Jr. and L. Capelle, Continuous growth of mean tumor diameter in a subset of grade II gliomas, Annals of Neurology, 53 (2003), 524-528.
[40]	D. C. Markham, M. J. Simpson, P. K. Maini, E. A. Gaffney and R. E. Baker, Comparing methods for modelling spreading cell fronts, Journal of Theoretical Biology, 353 (2014), 95-103.doi: 10.1016/j.jtbi.2014.02.023.
[41]	A. Martínez-González, G. F. Calvo, L. A. Pérez Romasanta and V. M. Pérez-García, Hypoxic Cell Waves Around Necrotic Cores in Glioblastoma: A Biomathematical Model and Its Therapeutic Implications, Bulletin of Mathematical Biology, 74 (2012), 2875-2896.doi: 10.1007/s11538-012-9786-1.
[42]	M. Marusic, Z. Bajzer, J. P. Freyer and S. Vuk-Palovic, Analysis of growth of multicellular tumour spheroids by mathematical models, Cell Prolif, 27 (1994), 73-94.
[43]	A. Matzavinos and M. A. J. Chaplain, Traveling-wave analysis of a model of the immune response to cancer, 2004, C. R. Biologies, 327 (2004), 995-1008.
[44]	E. Mehrara, E. Forssell-Aronsson, H. Ahlman and P. Bernhardt, Quantitative analysis of tumor growth rate and changes in tumor marker level: Specific growth rate versus doubling time, Acta Oncologica, 48 (2009), 591-597.doi: 10.1080/02841860802616736.
[45]	J. D. Murray, On travelling wave solutions in a model for the Belousov-Zhabotinskii reaction, Journal of Theoretical Biology, 56 (1976), 329-353.doi: 10.1016/S0022-5193(76)80078-1.
[46]	J. D. Murray, Mathematical Biology, 3rd ed, New York, NY: Springer- Verlag, 2002.doi: 10.1007/b98869.
[47]	B. Perthame, M. Tang and N. Vauchelet, Traveling wave solution of the Hele-Shaw model of tumor growth with nutrient, Math. Models Methods Appl. Sci., 24 (2014), 2601-2626. arXiv:1401.3649.doi: 10.1142/S0218202514500316.
[48]	A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, Chapman & Hall/CRC, Boca Raton London New York Washington, D.C., 2003.
[49]	R. Rockne, E. C. Alvord Jr., J. K. Rockhill and K. R. Swanson, A mathematical model for brain tumor response to radiation therapy, J. Math. Biol., 58 (2009), 561-578.doi: 10.1007/s00285-008-0219-6.
[50]	A. M. Samsonov, On exact quasistationary solutions to a nonlinear reaction-diffusion equation, Physics Letters A, 245 (1998), 527-536.doi: 10.1016/S0375-9601(98)00458-7.
[51]	D. L. Silbergeld, R. C. Rostomily and E. C. Alvord Jr., The cause of death in patients with glioblastoma is multifactorial: clinical factors and autopsy findings in 117 cases of supratentorial glioblastoma in adults, Journal of Neuro-Oncology, 10 (1991), 179-185.doi: 10.1007/BF00146880.
[52]	H. E. Skipper, F. M. Schabel and H. H. Lloyd, Dose-response and tumor cell repopulation rate in chemotherapeutic trials, Adv. Cancer Chemother., 1 (1979), 205-253.
[53]	A. Solovyova, P. Schuck, L. Costenaro and C. Ebel, Non-Ideality by Sedimentation Velocity of Halophilic Malate Dehydrogenase in Complex Solvents, Biophysical Journal, 81 (2001), 1868-1880.doi: 10.1016/S0006-3495(01)75838-9.
[54]	A. M. Stein, T. Demuth, D. Mobley, M. Berens and L. M. Sander, A Mathematical Model of Glioblastoma Tumor Spheroid Invasion in a Three-Dimensional In Vitro Experiment, Biophysical Journal, 92 (2007), 356-365.doi: 10.1529/biophysj.106.093468.
[55]	K. R. Swanson, C. Bridgea, J. D. Murray and E. C. Alvord Jr., Virtual and real brain tumors: Using mathematical modeling to quantify glioma growth and invasion, Journal of the Neurological Sciences, 216 (2003), 1-10.doi: 10.1016/j.jns.2003.06.001.
[56]	K. R. Swanson and E. C. Alvord Jr., The concept of gliomas as a travelling wave - application of a mathematical model to high - and low-grade gliomas, Can J Neurol Sci, 29 (2002), p395.
[57]	P. Tracqui, G. C. Cruywagen, D. E. Woodward, G. T. Bartooll, J. D. Murray and E. C. Alvord, Jr., A mathematical model of glioma growth: The effect of chemotherapy on spatio-temporal growth, Cell Prolif., 28 (1995), 17-31.doi: 10.1111/j.1365-2184.1995.tb00036.x.
[58]	P. Y. Wen and S. Kesari, Malignant Gliomas in Adults, The New England Journal of Medicine, 359 (2008), 492-507.doi: 10.1056/NEJMra0708126.
[59]	G. B. West, J. H. Brown and B. J. Enquist, A general model for ontogenetic growth, Nature, 413 (2001), 628-631.doi: 10.1038/35098076.
[60]	D. E. Woodward, J. Cook, P. Tracqui, G. C. Cruywagen, J. D. Murray and E. C. Alvord, Jr., A mathematical model of glioma growth: The effect of extent of surgical resection, Cell. Prolif., 29 (1996), 269-288.doi: 10.1111/j.1365-2184.1996.tb01580.x.
[61]	E. D. Yorke, Z. Fuks, L. Norton, V. Whitmore and C. C. Ling, Modeling the Development of Metastases from Primary and Locally Recurrent Tumors: Comparison with a Clinical Data Base for Prostatic Cancer, Cancer Research, 53 (1993), 2987-2993.