Combining cellular automata and Lattice Boltzmann method to model multiscale avascular tumor growth coupled with nutrient diffusion and immune competition.

In the last decades the Lattice Boltzmann method (LB) has been successfully used to simulate a variety of processes. The LB model describes the microscopic processes occurring at the cellular level and the macroscopic processes occurring at the continuum level with a unique function, the probability distribution function. Recently, it has been tried to couple deterministic approaches with probabilistic cellular automata (probabilistic CA) methods with the aim to model temporal evolution of tumor growths and three dimensional spatial evolution, obtaining hybrid methodologies. Despite the good results attained by CA-PDE methods, there is one important issue which has not been completely solved: the intrinsic stochastic nature of the interactions at the interface between cellular (microscopic) and continuum (macroscopic) level. CA methods are able to cope with the stochastic phenomena because of their probabilistic nature, while PDE methods are fully deterministic. Even if the coupling is mathematically correct, there could be important statistical effects that could be missed by the PDE approach. For such a reason, to be able to develop and manage a model that takes into account all these three level of complexity (cellular, molecular and continuum), we believe that PDE should be replaced with a statistic and stochastic model based on the numerical discretization of the Boltzmann equation: The Lattice Boltzmann (LB) method. In this work we introduce a new hybrid method to simulate tumor growth and immune system, by applying Cellular Automata Lattice Boltzmann (CA-LB) approach.

[1]  Pier Paolo Delsanto,et al.  Nutrient competition as a determinant for cancer growth , 1999 .

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

[3]  F. Pappalardo,et al.  Analysis of vaccine's schedules using models. , 2006, Cellular immunology.

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

[5]  Francesco Pappalardo,et al.  Genetic Algorithm Against Cancer , 2005, WILF.

[6]  Filippo Castiglione,et al.  Modelling vaccination schedules for a cancer immunoprevention vaccine , 2005, Immunome research.

[7]  Pierre Leone,et al.  A Hybrid Lattice Boltzmann Finite Difference Scheme for the Diffusion Equation , 2006 .

[8]  Philip K Maini,et al.  Angiogenesis and vascular remodelling in normal and cancerous tissues , 2009, Journal of mathematical biology.

[9]  Helen M. Byrne,et al.  A theoretical study of the response of vascular tumours to different types of chemotherapy , 2008, Math. Comput. Model..

[10]  E. T. Gawlinski,et al.  A reaction-diffusion model of cancer invasion. , 1996, Cancer research.

[11]  G. J. Pettet,et al.  The migration of cells in multicell tumor spheroids , 2001, Bulletin of mathematical biology.

[12]  P. J. Ponzo,et al.  A model for the growth of a solid tumor with non-uniform oxygen consumption , 1977 .

[13]  Francesco Pappalardo,et al.  Modeling the competition between lung metastases and the immune system using agents , 2010, BMC Bioinformatics.

[14]  Salvatore Musumeci,et al.  Modeling immune system control of atherogenesis , 2008, Bioinform..

[15]  S. Chen,et al.  Lattice methods and their applications to reacting systems , 1993, comp-gas/9312001.

[16]  Massimo Bernaschi,et al.  ImmunoGrid: towards agent-based simulations of the human immune system at a natural scale† , 2010, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[17]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[18]  S. Jonathan Chapman,et al.  Mathematical Models of Avascular Tumor Growth , 2007, SIAM Rev..

[19]  F. Pappalardo,et al.  In silico modeling and in vivo efficacy of cancer-preventive vaccinations. , 2010, Cancer research.

[20]  Serdal Pamuk,et al.  The method of lines for the numerical solution of a mathematical model for capillary formation: The role of tumor angiogenic factor in the extra-cellular matrix , 2007, Appl. Math. Comput..

[21]  Filippo Castiglione,et al.  Vaccine protocols optimization: in silico experiences. , 2010, Biotechnology advances.

[22]  Andreas Deutsch,et al.  Cellular Automaton Modeling of Biological Pattern Formation - Characterization, Applications, and Analysis , 2005, Modeling and simulation in science, engineering and technology.

[23]  Francesco Pappalardo,et al.  Optimal vaccination schedules using simulated annealing , 2008, Bioinform..

[24]  John Carl Panetta,et al.  Molecular biology of breast cancer metastasis: The use of mathematical models to determine relapse and to predict response to chemotherapy in breast cancer , 2000, Breast Cancer Research.

[25]  Taro Imamura,et al.  Acceleration of steady-state lattice Boltzmann simulations on non-uniform mesh using local time step method , 2005 .

[26]  M Scalerandi,et al.  Local interaction simulation approach for the response of the vascular system to metabolic changes of cell behavior. , 2000, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  I. Fatt,et al.  Oxygen Consumption Rate of Tissue Measured by a Micropolarographic Method , 1966, The Journal of general physiology.

[28]  Nicola Bellomo,et al.  On the foundations of cancer modelling: Selected topics, speculations, and perspectives , 2008 .

[29]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[30]  Davide Alemani,et al.  A lattice Boltzmann numerical approach for modelling reaction-diffusion processes in chemically and physically heterogeneous environments , 2007 .

[31]  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 .

[32]  C. Cercignani.  The Boltzmann equation and its applications , 1988 .

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

[34]  Dimitri J. Mavriplis,et al.  Multigrid solution of the steady-state lattice Boltzmann equation , 2006 .

[35]  M. Chaplain,et al.  Does breast cancer exist in a state of chaos? , 1999, European journal of cancer.

[36]  Olga Bandman,et al.  A Hybrid Approach to Reaction-Diffusion Processes Simulation , 2001, PaCT.

[37]  M. L. Martins,et al.  Morphology transitions induced by chemotherapy in carcinomas in situ. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  J. D. B. MACDOUGALL,et al.  Diffusion Coefficient of Oxygen through Tissues , 1967, Nature.

[39]  D L S McElwain,et al.  A history of the study of solid tumour growth: The contribution of mathematical modelling , 2004, Bulletin of mathematical biology.

[40]  Filippo Castiglione,et al.  Modeling and simulation of cancer immunoprevention vaccine , 2005, Bioinform..

[41]  J. Panetta,et al.  A mathematical model of periodically pulsed chemotherapy: tumor recurrence and metastasis in a competitive environment. , 1996, Bulletin of mathematical biology.

[42]  E. T. Gawlinski,et al.  The glycolytic phenotype in carcinogenesis and tumor invasion: insights through mathematical models. , 2003, Cancer research.

[43]  Andrea Cristofolini,et al.  A hybrid implicit numerical method for the analysis of the magneto-plasmadynamics in a gas discharge , 2001 .

[44]  Lance L. Munn,et al.  Lattice-Boltzmann simulation of blood flow in digitized vessel networks , 2008, Comput. Math. Appl..

[45]  Ernst Rank,et al.  A Multigrid-Solver for the Discrete Boltzmann Equation , 2002 .

[46]  Takashi Suzuki,et al.  Mathematical models of tumour angiogenesis , 2007 .

[47]  D. Wolf-Gladrow.  Lattice-Gas Cellular Automata and Lattice Boltzmann Models: An Introduction , 2000 .

[48]  Francesco Pappalardo,et al.  SimB16: Modeling Induced Immune System Response against B16-Melanoma , 2011, PloS one.

[49]  H. Greenspan.  On the growth and stability of cell cultures and solid tumors. , 1976, Journal of theoretical biology.

[50]  Massimo Bernaschi,et al.  ImmunoGrid, an integrative environment for large-scale simulation of the immune system for vaccine discovery, design and optimization , 2008, Briefings Bioinform..

[51]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

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

[53]  T. Foster,et al.  Dosimetry in photodynamic therapy: oxygen and the critical importance of capillary density. , 1992, Radiation research.

[54]  R Mastrangelo,et al.  Antitumour activity of oxaliplatin in neuroblastoma cell lines. , 1999, European journal of cancer.

[55]  Panayotis G. Kevrekidis,et al.  A hybrid model for tumor-induced angiogenesis in the cornea in the presence of inhibitors , 2007, Math. Comput. Model..

[56]  D. Hanahan,et al.  The Hallmarks of Cancer , 2000, Cell.

[57]  D. Kirschner,et al.  A mathematical model of tumor-immune evasion and siRNA treatment , 2003 .

[58]  A. Anderson,et al.  A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion , 2005 .

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

[60]  Marilisa Neri,et al.  Coarse-grained model of proteins incorporating atomistic detail of the active site. , 2005, Physical review letters.

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