A Cellular Automaton Model of Early Tumor Growth and Invasion: The Effects of Native Tissue Vascularity and Increased Anaerobic Tumor Metabolism

A hybrid cellular automaton model is described and used to simulate early tumor growth and examine the roles of host tissue vascular density and tumor metabolism in the ability of a small number of monoclonal transformed cells to develop into an invasive tumor. The model incorporates normal cells, tumor cells, necrotic or empty space, and a random network of native microvessels as components of a cellular automaton state vector. Diffusion of glucose and H(+)ions (the latter largely resulting from the tumor's excessive reliance on anaerobic metabolism) to and from the microvessels, and their utilization or production by cells, is modeled through the solution of differential equations. In this way, the cells and microvessels affect the extracellular concentrations of glucose and H(+)which, in turn, affect the evolution of the automaton. Simulations of the model demonstrate that: (i) high tumor H(+)ion production is favorable for tumor growth and invasion; however for every H(+)ion production rate, there exists a range of optimal microvessel densities (leading to a local pH favorable to tumor but not to normal cells) for which growth and invasion is most effective, (ii) at vascular densities below this range, both tumor and normal cells die due to excessively low pH, (iii) for vascular densities above the optimal range the microvessel network is highly efficient at removing acid and therefore the tumor cells lose their advantage over normal cells gained by high local H(+)concentration. While significant spatial gradients of glucose formed, no regions of detrimentally poor glucose perfusion (for either cell type) were observed, regardless of microvessel density. Depending on metabolic phenotype, a variety of tumor morphologies similar to those clinically observed were realized in the simulations. Lastly, a sharp transition (analogous to that of the adenoma-carcinoma sequence) between states of initial tumor confinement and efficient invasiveness was observed when H(+)production reached a critical value.

[1]  Otto Warburn,et al.  THE METABOLISM OF TUMORS , 1931 .

[2]  Norman Goldstein,et al.  Metastases in carcinoma. Analysis of 1000 autopsied cases , 1950, Cancer.

[3]  M E KLINGER,et al.  Secondary tumors of the genito-urinary tract. , 1951, The Journal of urology.

[4]  J E Newsam,et al.  Metastatic tumours in the kidney. , 1966, British journal of urology.

[5]  P. Gullino,et al.  Relationship between oxygen and glucose consumption by transplanted tumors in vivo. , 1967, Cancer research.

[6]  William Francis Ganong,et al.  Review of Medical Physiology , 1969 .

[7]  G. Murphy,et al.  Secondary carcinomas of the kidney. , 1975, The Journal of urology.

[8]  S. Harguindey,et al.  Effects of systemic acidification of mice with Sarcoma 180. , 1979, Cancer research.

[9]  R. Sutherland,et al.  A move for the better. , 1994, Environmental health perspectives.

[10]  J S Fowler,et al.  Increased accumulation of 2-deoxy-2-[18F]Fluoro-D-glucose in liver metastases from colon carcinoma. , 1982, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[11]  I. Fidler,et al.  Biological diversity in metastatic neoplasms: origins and implications. , 1982, Science.

[12]  Cellular oxygen utilization and radiation response of V-79 spheroids. , 1983, Advances in experimental medicine and biology.

[13]  R. Sutherland,et al.  Growth and cellular characteristics of multicell spheroids. , 1984, Recent results in cancer research. Fortschritte der Krebsforschung. Progres dans les recherches sur le cancer.

[14]  B. Vogelstein,et al.  Clonal analysis of human colorectal tumors. , 1987, Science.

[15]  L. Tertoolen,et al.  The Na+/H+ exchanger is constitutively activated in P19 embryonal carcinoma cells, but not in a differentiated derivative. Responsiveness to growth factors and other stimuli. , 1987, The Journal of biological chemistry.

[16]  P. Vaupel,et al.  Glucose uptake, lactate release, ketone body turnover, metabolic micromilieu, and pH distributions in human breast cancer xenografts in nude rats. , 1988, Cancer research.

[17]  [Spontaneous regression of lung metastasis of renal cell carcinoma: a case report]. , 1988, Hinyokika kiyo. Acta urologica Japonica.

[18]  W. Paschen,et al.  Metabolic imaging in microregions of tumors and normal tissues with bioluminescence and photon counting. , 1988, Journal of the National Cancer Institute.

[19]  R. Gatenby,et al.  Oxygen distribution in squamous cell carcinoma metastases and its relationship to outcome of radiation therapy. , 1988, International journal of radiation oncology, biology, physics.

[20]  P. Vaupel,et al.  pH distributions in spontaneous and isotransplanted rat tumours. , 1988, British Journal of Cancer.

[21]  P. Okunieff,et al.  Blood flow, oxygen and nutrient supply, and metabolic microenvironment of human tumors: a review. , 1989, Cancer research.

[22]  I. Tannock,et al.  Acid pH in tumors and its potential for therapeutic exploitation. , 1989, Cancer research.

[23]  D. Hanahan,et al.  Induction of angiogenesis during the transition from hyperplasia to neoplasia , 1989, Nature.

[24]  [A case of spontaneous disappearance of pulmonary metastasis of renal cell carcinoma following nephrectomy]. , 1990, Hinyokika kiyo. Acta urologica Japonica.

[25]  H. Yamasaki Aberrant expression and function of gap junctions during carcinogenesis. , 1991, Environmental Health Perspectives.

[26]  J. Griffiths Are cancer cells acidic? , 1991, British Journal of Cancer.

[27]  R. Gatenby,et al.  Population ecology issues in tumor growth. , 1991, Cancer research.

[28]  C. Jm,et al.  Spontaneous disappearance of pulmonary metastasis secondary to renal cell carcinoma after nephrectomy. Presentation of a case and review of the literature , 1992 .

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

[30]  Y. Kato,et al.  Induction of 103-kDa gelatinase/type IV collagenase by acidic culture conditions in mouse metastatic melanoma cell lines. , 1992, The Journal of biological chemistry.

[31]  J. Folkman,et al.  The role of angiogenesis in tumor growth. , 1992, Seminars in cancer biology.

[32]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[33]  E. Keshet,et al.  Vascular endothelial growth factor induced by hypoxia may mediate hypoxia-initiated angiogenesis , 1992, Nature.

[34]  E. Hoffman,et al.  The role of positron emission tomography in oncology and other whole-body applications. , 1992, Seminars in nuclear medicine.

[35]  X. Zheng,et al.  A cellular automaton model of cancerous growth. , 1993, Journal of theoretical biology.

[36]  M S Jacobson,et al.  Normal reference laboratory values? , 1978, The New England journal of medicine.

[37]  J. Smolle,et al.  Computer simulation of tumour cell invasion by a stochastic growth model. , 1993, Journal of theoretical biology.

[38]  W N Hanafee,et al.  Extracranial head and neck: PET imaging with 2-[F-18]fluoro-2-deoxy-D-glucose and MR imaging correlation. , 1993, Radiology.

[39]  G B Ermentrout,et al.  Cellular automata approaches to biological modeling. , 1993, Journal of theoretical biology.

[40]  Z Bajzer,et al.  Analysis of growth of multicellular tumour spheroids by mathematical models , 1994, Cell proliferation.

[41]  R K Jain,et al.  Noninvasive measurement of interstitial pH profiles in normal and neoplastic tissue using fluorescence ratio imaging microscopy. , 1994, Cancer research.

[42]  Bonnie F. Sloane,et al.  Pericellular pH affects distribution and secretion of cathepsin B in malignant cells. , 1994, Cancer research.

[43]  W. Boron,et al.  Long-term expression of c-H-ras stimulates Na-H and Na(+)-dependent Cl-HCO3 exchange in NIH-3T3 fibroblasts. , 1994, The Journal of biological chemistry.

[44]  R. Gillies,et al.  31P-MRS measurements of extracellular pH of tumors using 3-aminopropylphosphonate. , 1994, The American journal of physiology.

[45]  K. Jeong,et al.  Metabolic consequences of a reversed pH gradient in rat tumors. , 1994, Cancer research.

[46]  R. Gatenby,et al.  The potential role of transformation-induced metabolic changes in tumor-host interaction. , 1995, Cancer research.

[47]  M. Neeman,et al.  Induction of vascular endothelial growth factor expression by hypoxia and by glucose deficiency in multicell spheroids: implications for tumor angiogenesis. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[48]  H M Byrne,et al.  Growth of nonnecrotic tumors in the presence and absence of inhibitors. , 1995, Mathematical biosciences.

[49]  S. Hubchak,et al.  Multiple genetic alterations in hamster pancreatic ductal adenocarcinomas. , 1995, Cancer research.

[50]  Nicola Bellomo,et al.  A Survey of Models for Tumor-Immune System Dynamics , 1996 .

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

[52]  M. Chaplain,et al.  Mathematical modelling, simulation and prediction of tumour-induced angiogenesis. , 1996, Invasion & metastasis.

[53]  M. Kraus,et al.  Implications of acidic tumor microenvironment for neoplastic growth and cancer treatment: a computer analysis. , 1996, Tumour biology : the journal of the International Society for Oncodevelopmental Biology and Medicine.

[54]  P. Maini,et al.  A mathematical model for the capillary endothelial cell-extracellular matrix interactions in wound-healing angiogenesis. , 1997, IMA journal of mathematics applied in medicine and biology.

[55]  M Buyse,et al.  Non-small-cell lung carcinoma tumor growth without morphological evidence of neo-angiogenesis. , 1997, The American journal of pathology.

[56]  Rakesh K. Jain,et al.  Interstitial pH and pO2 gradients in solid tumors in vivo: High-resolution measurements reveal a lack of correlation , 1997, Nature Medicine.

[57]  J. King,et al.  Mathematical modelling of avascular-tumour growth. , 1997, IMA journal of mathematics applied in medicine and biology.

[58]  A. Kratz,et al.  Case records of the Massachusetts General Hospital. Weekly clinicopathological exercises. Normal reference laboratory values. , 1998, The New England journal of medicine.

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

[60]  Robert A. Gatenby,et al.  Cellular automaton model of tumor growth , 1998 .

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

[62]  J Smolle,et al.  Fractal tumor stromal border in a nonequilibrium growth model. , 1998, Analytical and quantitative cytology and histology.

[63]  S. Kumar,et al.  Breast carcinoma: vascular density determined using CD105 antibody correlates with tumor prognosis. , 1999, Cancer research.

[64]  Alexander R. A. Anderson,et al.  A mathematical analysis of a model for capillary network formation in the absence of endothelial cell proliferation , 1999 .

[65]  M. Markus,et al.  Simulation of vessel morphogenesis using cellular automata. , 1999, Mathematical biosciences.

[66]  G. Yancopoulos,et al.  Vessel cooption, regression, and growth in tumors mediated by angiopoietins and VEGF. , 1999, Science.

[67]  J. King,et al.  Mathematical modelling of avascular-tumour growth. II: Modelling growth saturation. , 1999, IMA journal of mathematics applied in medicine and biology.

[68]  A. Perumpanani,et al.  Numerical interactions of random and directed motility during cancer invasion , 1999 .

[69]  T. A. Smith,et al.  Mammalian hexokinases and their abnormal expression in cancer. , 2000, British journal of biomedical science.

[70]  H M Byrne,et al.  A mathematical model of the stress induced during avascular tumour growth , 2000, Journal of mathematical biology.

[71]  S Torquato,et al.  Cellular automaton of idealized brain tumor growth dynamics. , 2000, Bio Systems.

[72]  D A Rew,et al.  Cell production rates in human tissues and tumours and their significance. Part II: clinical data. , 2000, European journal of surgical oncology : the journal of the European Society of Surgical Oncology and the British Association of Surgical Oncology.

[73]  Ovide Arino,et al.  Mathematical modelling of the intravenous glucose tolerance test , 2000, Journal of mathematical biology.

[74]  S Torquato,et al.  Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. , 2000, Journal of theoretical biology.

[75]  H M Byrne,et al.  A mathematical model to study the effects of drug resistance and vasculature on the response of solid tumors to chemotherapy. , 2000, Mathematical biosciences.

[76]  E. Rofstad,et al.  High lactate levels predict likelihood of metastases, tumor recurrence, and restricted patient survival in human cervical cancers. , 2000, Cancer research.