| Peer-Reviewed

A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed

Received: 30 August 2015     Accepted: 22 September 2015     Published: 13 October 2015
Views:       Downloads:
Abstract

The main aim of this paper is to investigate the potential performance improvements gained from a serial versus parallel implementation of the numerical solution to a system of coupled nonlinear PDEs describing tumour-induced angiogenesis. After applying a suitable finite difference scheme, the resulting hybrid continuous discreet model is solved based on a set of cellular rules defining endothelial cell movement towards a tumour. In addition, the model explicitly incorporates the processes of branching, anastomosis and cell proliferation. Parallel implementations are based on the CUDA programming model with a detailed look at efficient thread deployment and memory management. Results show substantial speedups for the CUDA C language against that of conventional high performance languages, such as C++. Such increased performance highlights the potential for simulating more complex mathematical models of tumour dynamics, such as vascularisation networks, tumour invasion and metastasis, leading to the potential for more rapid experimental results for a range of complex cancer models.

Published in European Journal of Biophysics (Volume 3, Issue 5)
DOI 10.11648/j.ejb.20150305.11
Page(s) 27-37
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2015. Published by Science Publishing Group

Keywords

Cancer Modelling, Tumour Angiogenic Factors (TAF), Tumour-Induced Angiogenesis, Anastomoses, Parallel Programming, Compute Unified Device Architecture (CUDA), Graphical Processing Unit (GPU), High Performance Computing (HPC)

References
[1] Arnaoutova, I. and Kleinman, H. K. In vitro angiogenesis: Endothelial cell tube formation on gelled basement membrane extract. Nature Protocols, 5, 628–635. 2010.
[2] Gimbrone MA, Cotran RS, Leapman SB, Folkman J: Tumor growth and neovascularization: An experimental model using the rabbit cornea. Journal of the National Cancer Institute. 52, 413–427. 1974.
[3] Folkman, J. and Haudenschild, C. Angiogenesis in vitro. Nature, 288, 551–556. 1980.
[4] Muthukkaruppan, V. R., Kubai, L., Auerbach, R. Tumorinduced neovascularization in the mouse eye. Journal of the National Cancer Institute, 69, 699–705. 1982.
[5] Jain RK, Schlenger K, Hӧckel M. and Yuan F: Quantitative angiogenesis assays: Progress and problems. Nature Medicine 3: 1203–1208.1997.
[6] Hanahan, D and Folkman, J. Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell, 86, 353–364. 1996.
[7] Albini, A., Tosetti, A. F., Li, W. V., Noonan, D. M. and Li, W. W. Cancer prevention by targeting angiogenesis Nature Reviews Clinical Oncology 9, 498-509. 2012.
[8] Ferrara, N. and Kerbel, R. S. Angiogenesis as a therapeutic target. Nature, 438 967–974. 2005.
[9] Carmeliet, P. Angiogenesis in life, disease and medicine. Nature, 438: 932–936. 2005.
[10] Bouard S. de, Herlin, P. and Christensen, J. G. Antiangio-genic and anti-invasive effects of sunitinib on experimental human glioblastoma. Neuro-Oncology, Vol. 9, No. 4, 412– 423. 2007.
[11] Norden, A. D, Drappatz, J. and Wen P. Y. Novel antiangiogenic therapies for malignant gliomas. The Lancet Neurology, Vol. 7, No. 12, 1152–1160. 2008.
[12] Peirce, S. M. Computational and mathematical modeling of angiogenesis. Microcirculation, 15(8), 739–751. 2008.
[13] M. Scianna, M., Bell. C. and Preziosi L. A review of mathematical models for the formation of vascular networks. Oxford Centre for Collaborative Applied Mathematics. 2012.
[14] Darbyshire, P. M. 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, 65-73. 2015.
[15] Rejniak A. K. and Anderson A.R.A. Hybrid models of tumor growth. Wiley Interdisciplinary Reviews: Systems Biology and Medicine, 3(1), 115–125. 2011.
[16] Paweletz, N. and M. Knierim M. Tumor-related angiogenesis. Critical Reviews in Oncology and Hematology, 9, 197–242. 1989.
[17] Paku, S. and N. Paweletz. First steps of tumor-related angiogenesis. Laboratory Investigation, 65, 334–346. 1991.
[18] Schor S. L., Schor A. M., Brazill G. W. The effects of fibronectin on the migration of human foreskin fibroblasts and Syrian hamster melanoma cells into three-dimensional gels of lattice collagen fibres. Journal of Cell Science, 48, 301–314, 1981.
[19] Bowersox J. C. and Sorgente N. Chemotaxis of aortic endothelial cells in response to fibronectin. Cancer Research 42, 2547–2551. 1982.
[20] Quigley J. P., Lacovara J., and Cramer E. B. The directed migration of B-16 melanoma-cells in response to a haptotactic chemotactic gradient of fibronectin. Journal of Cell Biology 97, A450–451. 1983.
[21] Stokes C. L., Lauffenburger D. A., and Williams S. K. Migration of individual microvessel endothelial cells: stochastic model and parameter measurement. Journal of Cell Science, 99: 419–430. 1991.
[22] Stokes C. L., Rupnick M. A., Williams S. K., and Lauffenburger D. A. Chemotaxis of human microvessel endothelial cells in response to acidic fibroblast growth factor. Laboratory Investigation, 63, 657–668, 1991.
[23] Stokes C. L. and Lauffenburger D. A. Analysis of the roles of microvessel endothelial cell random motility and chemotaxis in angiogenesis. Journal of Theoretical Biology, 152, 377–403. 1991.
[24] Anderson, A. R. A. and Chaplain, M. Continuous and discrete mathematical models of tumour-induced angiogenesis, Bulletin of Mathematical Biology, 60, 857-900. 1998.
[25] Anderson, A., B. D. S. Sleeman, I. M. Young and B. S. Griffiths. Nematode movement along a chemical gradient in a structurally heterogeneous environment: II. Theory. Fundamental and Applied Nematology, 20, 165–172. 1997.
[26] Nvidia Corporation. CUDA C programming guide. Version 6.0. 2014.
[27] Cheng, J., Grossman, M and McKercher, Ty. Professional CUDA C Programming. Wrox. 2014.
[28] Venkatasubramanian, S. and Vuduc, R. W. Tuned and wildly asynchronous stencil kernels for hybrid CPU/GPU systems. In Proceedings of the Association of Computing Machinery International Conference on Supercomputing, New York. 2009.
[29] Amorim, R., Haase, G., Liebmann, M. and Weber dos Santos, R. Comparing CUDA and OpenGL implementations for a Jacobi iteration. In Proceedings of High Performance Computing and Simulation Conference, Berlin. 2009.
[30] Cecilia, J. M., Garcıa, J. M. and Ujaldon, M. CUDA 2D stencil computations for the Jacobi method. Springer, 173-183. 2012.
Cite This Article
  • APA Style

    Paul M. Darbyshire. (2015). A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed. European Journal of Biophysics, 3(5), 27-37. https://doi.org/10.11648/j.ejb.20150305.11

    Copy | Download

    ACS Style

    Paul M. Darbyshire. A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed. Eur. J. Biophys. 2015, 3(5), 27-37. doi: 10.11648/j.ejb.20150305.11

    Copy | Download

    AMA Style

    Paul M. Darbyshire. A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed. Eur J Biophys. 2015;3(5):27-37. doi: 10.11648/j.ejb.20150305.11

    Copy | Download

  • @article{10.11648/j.ejb.20150305.11,
      author = {Paul M. Darbyshire},
      title = {A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed},
      journal = {European Journal of Biophysics},
      volume = {3},
      number = {5},
      pages = {27-37},
      doi = {10.11648/j.ejb.20150305.11},
      url = {https://doi.org/10.11648/j.ejb.20150305.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ejb.20150305.11},
      abstract = {The main aim of this paper is to investigate the potential performance improvements gained from a serial versus parallel implementation of the numerical solution to a system of coupled nonlinear PDEs describing tumour-induced angiogenesis. After applying a suitable finite difference scheme, the resulting hybrid continuous discreet model is solved based on a set of cellular rules defining endothelial cell movement towards a tumour. In addition, the model explicitly incorporates the processes of branching, anastomosis and cell proliferation. Parallel implementations are based on the CUDA programming model with a detailed look at efficient thread deployment and memory management. Results show substantial speedups for the CUDA C language against that of conventional high performance languages, such as C++. Such increased performance highlights the potential for simulating more complex mathematical models of tumour dynamics, such as vascularisation networks, tumour invasion and metastasis, leading to the potential for more rapid experimental results for a range of complex cancer models.},
     year = {2015}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - A Hybrid Continuous-Discrete Model of Tumour-Induced Angiogenesis is Solved Numerically in Parallel and Performance Improvements Analysed
    AU  - Paul M. Darbyshire
    Y1  - 2015/10/13
    PY  - 2015
    N1  - https://doi.org/10.11648/j.ejb.20150305.11
    DO  - 10.11648/j.ejb.20150305.11
    T2  - European Journal of Biophysics
    JF  - European Journal of Biophysics
    JO  - European Journal of Biophysics
    SP  - 27
    EP  - 37
    PB  - Science Publishing Group
    SN  - 2329-1737
    UR  - https://doi.org/10.11648/j.ejb.20150305.11
    AB  - The main aim of this paper is to investigate the potential performance improvements gained from a serial versus parallel implementation of the numerical solution to a system of coupled nonlinear PDEs describing tumour-induced angiogenesis. After applying a suitable finite difference scheme, the resulting hybrid continuous discreet model is solved based on a set of cellular rules defining endothelial cell movement towards a tumour. In addition, the model explicitly incorporates the processes of branching, anastomosis and cell proliferation. Parallel implementations are based on the CUDA programming model with a detailed look at efficient thread deployment and memory management. Results show substantial speedups for the CUDA C language against that of conventional high performance languages, such as C++. Such increased performance highlights the potential for simulating more complex mathematical models of tumour dynamics, such as vascularisation networks, tumour invasion and metastasis, leading to the potential for more rapid experimental results for a range of complex cancer models.
    VL  - 3
    IS  - 5
    ER  - 

    Copy | Download

Author Information
  • Department of Computational Biophysics, Algenet Cancer Research, Nottingham, UK

  • Sections