A System of Coupled Nonlinear Partial Differential Equations Describing Avascular Tumour Growth Are Solved Numerically Using Parallel Programming to Assess Computational Speedup
Computational Biology and Bioinformatics
Volume 3, Issue 5, October 2015, Pages: 65-73
Received: Aug. 3, 2015;
Accepted: Aug. 13, 2015;
Published: Aug. 26, 2015
Views 3756 Downloads 120
Paul M. Darbyshire, Department of Computational Biophysics, Algenet Cancer Research, Nottingham, UK
The challenging issues of cancer prevention and cure lie in the need for a more detailed knowledge of the internal processes and mechanisms of tumour growth. We present a mathematical model of avascular tumour growth formulated in a system of coupled nonlinear PDEs. The interaction between the surrounding tissue and cell motility of the developing tumour are also included to more realistic replicate an in-vivo environment. The mathematical model is solved using finite difference methods and implemented in the C programming language. The CUDA programming framework is then introduced to allow a parallelisation of the sequential C implementation. Results show a dramatic Speedup of around 26x that of conventional implementations in C. Such increased computational efficiency clearly highlights the possibility of improvements in the numerical simulation of more complex mathematical models of 2D and 3D tumour growth, such as angiogenesis and vascularisation. Parallelisation of such models can greatly facilitate researchers, clinicians and oncologists by performing time-saving in-silico experiments that have the potential to highlight new cancer treatments and therapies without the need for the use of valuable resources associated with excessive pre-clinical trials.
Paul M. Darbyshire,
A System of Coupled Nonlinear Partial Differential Equations Describing Avascular Tumour Growth Are Solved Numerically Using Parallel Programming to Assess Computational Speedup, Computational Biology and Bioinformatics.
Vol. 3, No. 5,
2015, pp. 65-73.
Araujo, R. and McElwain, D. A history of the study of solid tumor growth: the contribution of mathematical modelling. Bulletin of Mathematical Biology, 66. 2004.
Sutherland, R.M. and Durand, R. E. Growth and cellular characteristics of multicell spheroids. Recent Results in Cancer Research 95, 24-49. 1984.
Sutherland, R. M., Sordat, B., Bamat, J., Gabbert, H., Bourrat, B., Mueller-Klieser,W. Oxygenation and differentiation in multicellular spheroids of human colon carcinoma. Cancer Research, Vol. 46, 5320–5329. 1986.
Sutherland, R. M. Cell and environment interaction in tumour microregions: the multicell spheroid model. Science, Vol. 240, 177-184. 1988.
Freyer, J. P. and Schor, P. L. Regrowth of cells from multicell tumour spheroids. Cell and Tissue Kinetics, 20, 249. 1987.
Orme, M. and Chaplain M.A.J. A mathematical model of vascular tumor growth and invasion. Mathematical and Computational Modelling, 23. 1996.
Sherratt, J. A. Wave front propagation in a competition equation with a new motility term modelling contact inhibition between cell populations. Proceedings of the Royal Society of London, A456, 2365–2386. 2000.
Sherratt, J. A and Chaplain M.A.J. A new mathematical model for avascular tumour growth. Journal of Mathematical Biology, Vol. 43, pp291–312. 2001.
Huttenlocher, A., Lakonishok, M., Kinder, M., Wu, S., Truong, T., Knudsen, K. A. and Horwitz, A. F. Integrin and cadherin synergy regulates contact inhibition of migration and motile activity. Journal of Cell Biology 141, 515-526. 1998.
McElwain, D. L. S. and Pettet, G. J. Cell migration in multicell spheroids: swimming against the tide. Bulletin of Mathematical Biology, 55, 655–674. 1993.
Nvidia Corporation. CUDA C programming guide. Version 6.0. 2014.
Amdahl, G. M. Validity of the Single Processor Approach to Achieving Large-Scale Computing Capabilities. AFIPS Conference Proceedings (30): 483–485. 1967.