A multiple time-scale computational model of a tumor and its micro environment.

Experimental evidence suggests that a tumor's environment may be critical to designing successful therapeutic protocols: Modeling interactions between a tumor and its environment could improve our understanding of tumor growth and inform approaches to treatment. This paper describes an efficient, flexible, hybrid cellular automaton-based implementation of numerical solutions to multiple time-scale reaction-diffusion equations, applied to a model of tumor proliferation. The growth and maintenance of cells in our simulation depend on the rate of cellular energy (ATP) metabolized from nearby nutrients such as glucose and oxygen. Nutrient consumption rates are functions of local pH as well as local concentrations of oxygen and other fuels. The diffusion of these nutrients is modeled using a novel variation of random-walk techniques. Furthermore, we detail the effects of three boundary update rules on simulations, describing their effects on computational efficiency and biological realism. Qualitative and quantitative results from simulations provide insight on how tumor growth is affected by various environmental changes such as micro-vessel density or lower pH, both of high interest in current cancer research.

[1]  S. V. Sotirchos,et al.  Variations in tumor cell growth rates and metabolism with oxygen concentration, glucose concentration, and extracellular pH , 1992, Journal of cellular physiology.

[2]  T. Schulz,et al.  Induction of Oxidative Metabolism by Mitochondrial Frataxin Inhibits Cancer Growth , 2006, Journal of Biological Chemistry.

[3]  M. Dewhirst,et al.  Concepts of oxygen transport at the microcirculatory level. , 1998, Seminars in radiation oncology.

[4]  M Gryczyński,et al.  Intratumoral microvessels density and morphometric study of angiogenesis as prognostic factor in laryngeal cancer , 2003 .

[5]  Steven M. Wise,et al.  An adaptive multigrid algorithm for simulating solid tumor growth using mixture models , 2011, Math. Comput. Model..

[6]  Antonio Fasano,et al.  Necrotic core in EMT6/Ro tumour spheroids: Is it caused by an ATP deficit? , 2010, Journal of theoretical biology.

[7]  Kevin J Painter,et al.  From a discrete to a continuous model of biological cell movement. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  Richard P. Hill,et al.  Acute Hypoxia Enhances Spontaneous Lymph Node Metastasis in an Orthotopic Murine Model of Human Cervical Carcinoma , 2004, Cancer Research.

[9]  J S Rasey,et al.  Determination of the radiobiologically hypoxic fraction in multicellular spheroids from data on the uptake of [3H]fluoromisonidazole. , 1995, Radiation research.

[10]  Antonio Fasano,et al.  Regression and regrowth of tumour cords following single-dose anticancer treatment , 2003, Bulletin of mathematical biology.

[11]  D. G. Mallet,et al.  Spatial tumor-immune modeling , 2006 .

[12]  Keigo Endo,et al.  Present role and future prospects of positron emission tomography in clinical oncology , 2006, Cancer science.

[13]  M. Chaplain,et al.  Modelling the role of cell-cell adhesion in the growth and development of carcinomas , 1996 .

[14]  Raissa M. D'Souza,et al.  Dimension-Splitting for Simplifying Diffusion in Lattice-Gas Models , 2001 .

[15]  T S Edgington,et al.  Expression of tissue factor by melanoma cells promotes efficient hematogenous metastasis. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Ami Radunskaya,et al.  A delay differential equation model for tumor growth , 2003, Journal of mathematical biology.

[17]  Tarik Tihan,et al.  The hypoxic response of tumors is dependent on their microenvironment. , 2003, Cancer cell.

[18]  Helen M Byrne,et al.  Mathematical modelling of the use of macrophages as vehicles for drug delivery to hypoxic tumour sites. , 2004, Journal of theoretical biology.

[19]  Antonio Fasano,et al.  ATP Production and Necrosis Formation in a Tumour Spheroid Model , 2007 .

[20]  P. Maini,et al.  A cellular automaton model for tumour growth in inhomogeneous environment. , 2003, Journal of theoretical biology.

[21]  D. Fraker,et al.  An integrated approach to measuring tumor oxygen status using human melanoma xenografts as a model. , 2003, Cancer research.

[22]  Georgios Ch. Sirakoulis,et al.  An Intelligent Cellular Automaton Model for Crowd Evacuation in Fire Spreading Conditions , 2007 .

[23]  L. D. Pillis,et al.  A Validated Mathematical Model of Cell-Mediated Immune Response to Tumor Growth , 2005 .

[24]  Luigi Preziosi,et al.  Cancer Modelling and Simulation , 2003 .

[25]  F Castiglione,et al.  Design and implementation of an immune system simulator , 2001, Comput. Biol. Medicine.

[26]  A. Anderson,et al.  Hybrid models of tumor growth , 2011, Wiley interdisciplinary reviews. Systems biology and medicine.

[27]  P. Maini,et al.  Metabolic changes during carcinogenesis: potential impact on invasiveness. , 2007, Journal of theoretical biology.

[28]  A. Deutsch,et al.  Modeling of self-organized avascular tumor growth with a hybrid cellular automaton. , 2002, In silico biology.

[29]  R. Gillies,et al.  Why do cancers have high aerobic glycolysis? , 2004, Nature Reviews Cancer.

[30]  Ami Radunskaya,et al.  A mathematical tumor model with immune resistance and drug therapy: an optimal control approach , 2001 .

[31]  Alexander Margulis,et al.  Microenvironmental control of premalignant disease: the role of intercellular adhesion in the progression of squamous cell carcinoma. , 2005, Seminars in cancer biology.

[32]  J. Smolle Cellular Automaton Simulation of Tumour Growth – Equivocal Relationships between Simulation Parameters and Morphologic Pattern Features , 1998, Analytical cellular pathology : the journal of the European Society for Analytical Cellular Pathology.

[33]  V. Cristini,et al.  Nonlinear simulation of tumor growth , 2003, Journal of mathematical biology.

[34]  E. T. Gawlinski,et al.  A Cellular Automaton Model of Early Tumor Growth and Invasion: The Effects of Native Tissue Vascularity and Increased Anaerobic Tumor Metabolism , 2001 .

[35]  M. Chaplain,et al.  Continuous and discrete mathematical models of tumor-induced angiogenesis , 1998, Bulletin of mathematical biology.

[36]  K. Berg,et al.  Energy metabolism in human melanoma cells under hypoxic and acidic conditions in vitro. , 1997, British Journal of Cancer.

[37]  K. Rejniak An immersed boundary framework for modelling the growth of individual cells: an application to the early tumour development. , 2007, Journal of theoretical biology.

[38]  M. Guppy,et al.  Cancer metabolism: facts, fantasy, and fiction. , 2004, Biochemical and biophysical research communications.

[39]  Mark A. J. Chaplain,et al.  Mathematical Modelling of Spatio-temporal Phenomena in Tumour Immunology , 2006 .

[40]  A. Giatromanolaki,et al.  Oxygen and glucose consumption in gastrointestinal adenocarcinomas: Correlation with markers of hypoxia, acidity and anaerobic glycolysis , 2006, Cancer science.

[41]  José C. M. Mombach,et al.  The interplay between cell adhesion and environment rigidity in the morphology of tumors , 2003 .

[42]  Norishige Chiba,et al.  Visual simulation of crack pattern based on 3D surface cellular automaton , 2000, Proceedings Seventh International Conference on Parallel and Distributed Systems: Workshops.

[43]  Marc Dellian,et al.  Acid production in glycolysis-impaired tumors provides new insights into tumor metabolism. , 2002, Clinical cancer research : an official journal of the American Association for Cancer Research.

[44]  John von Neumann,et al.  Theory Of Self Reproducing Automata , 1967 .

[45]  M. Chaplain,et al.  Free boundary value problems associated with the growth and development of multicellular spheroids , 1997, European Journal of Applied Mathematics.

[46]  M. L. Martins,et al.  Reaction-diffusion model for the growth of avascular tumor. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  R Hofmann-Wellenhof,et al.  Pattern interpretation by cellular automata (PICA)--evaluation of tumour cell adhesion in human melanomas. , 1994, Analytical cellular pathology : the journal of the European Society for Analytical Cellular Pathology.

[48]  E. Rofstad,et al.  Oxygen consumption rate and mitochondrial density in human melanoma monolayer cultures and multicellular spheroids , 1994, International journal of cancer.

[49]  L. D. de Pillis,et al.  A cellular automata model of tumor-immune system interactions. , 2006, Journal of theoretical biology.

[50]  H Lyng,et al.  Vascular density in human melanoma xenografts: relationship to angiogenesis, perfusion and necrosis. , 1998, Cancer letters.

[51]  Ami Radunskaya,et al.  The dynamics of an optimally controlled tumor model: A case study , 2003 .

[52]  Stephen Turner Using cell potential energy to model the dynamics of adhesive biological cells. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[53]  P. Leedman,et al.  Contribution by different fuels and metabolic pathways to the total ATP turnover of proliferating MCF-7 breast cancer cells. , 2002, The Biochemical journal.

[54]  G. Nicolson,et al.  The role of tumor cell adhesion as an important factor in formation of distant colorectal metastasis , 2001, Diseases of the colon and rectum.

[55]  P. Maini,et al.  The role of acidity in solid tumour growth and invasion. , 2005, Journal of theoretical biology.

[56]  Franco Celada,et al.  IMMSIM, a flexible model for in machina experiments on immune system responses , 2002, Future Gener. Comput. Syst..

[57]  Adrian L. Harris,et al.  Hypoxia — a key regulatory factor in tumour growth , 2002, Nature Reviews Cancer.

[58]  A. Radunskaya,et al.  Mixed Immunotherapy and Chemotherapy of Tumors: Modeling, Applications and Biological Interpretations , 2022 .

[59]  P. Vaupel,et al.  Current status of knowledge and critical issues in tumor oxygenation. Results from 25 years research in tumor pathophysiology. , 1998, Advances in experimental medicine and biology.

[60]  Richard P. Hill,et al.  The hypoxic tumour microenvironment and metastatic progression , 2004, Clinical & Experimental Metastasis.

[61]  Robert J Gillies,et al.  Contributions of cell metabolism and H+ diffusion to the acidic pH of tumors. , 2003, Neoplasia.

[62]  S. McDougall,et al.  Multiscale modelling and nonlinear simulation of vascular tumour growth , 2009, Journal of mathematical biology.

[63]  Zhongliang Ru,et al.  Cellular Automata to Simulate Rock Failure , 2006, ICAT Workshops.

[64]  Trachette L. Jackson,et al.  Vascular tumor growth and treatment: Consequences of polyclonality, competition and dynamic vascular support , 2002, Journal of mathematical biology.

[65]  E. Rofstad,et al.  Radiobiological and immunohistochemical assessment of hypoxia in human melanoma xenografts: acute and chronic hypoxia in individual tumours. , 1999, International journal of radiation biology.

[66]  Sebastiaan A.L.M. Kooijman,et al.  Dynamic Energy and Mass Budgets in Biological Systems , 2000 .

[67]  Michael A Henson,et al.  Incorporating energy metabolism into a growth model of multicellular tumor spheroids. , 2006, Journal of theoretical biology.

[68]  H. Frieboes,et al.  Predictive oncology: A review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth , 2007, NeuroImage.

[69]  Dirk Koczan,et al.  Mechanisms of Hypoxic Gene Regulation of Angiogenesis Factor Cyr61 in Melanoma Cells* , 2003, Journal of Biological Chemistry.

[70]  K. Rejniak,et al.  Current trends in mathematical modeling of tumor–microenvironment interactions: a survey of tools and applications , 2010, Experimental biology and medicine.

[71]  S. Kooijman,et al.  The embedded tumour: host physiology is important for the evaluation of tumour growth , 2003, British Journal of Cancer.

[72]  P. Vaupel,et al.  Tumor hypoxia: definitions and current clinical, biologic, and molecular aspects. , 2001, Journal of the National Cancer Institute.