Impact of metabolic heterogeneity on tumor growth, invasion, and treatment outcomes.

Histopathologic knowledge that extensive heterogeneity exists between and within tumors has been confirmed and deepened recently by molecular studies. However, the impact of tumor heterogeneity on prognosis and treatment remains as poorly understood as ever. Using a hybrid multiscale mathematical model of tumor growth in vascularized tissue, we investigated the selection pressures exerted by spatial and temporal variations in tumor microenvironment and the resulting phenotypic adaptations. A key component of this model is normal and tumor metabolism and its interaction with microenvironmental factors. The metabolic phenotype of tumor cells is plastic, and microenvironmental selection leads to increased tumor glycolysis and decreased pH. Once this phenotype emerges, the tumor dramatically changes its behavior due to acid-mediated invasion, an effect that depends on both variations in the tumor cell phenotypes and their spatial distribution within the tumor. In early stages of growth, tumors are stratified, with the most aggressive cells developing within the interior of the tumor. These cells then grow to the edge of the tumor and invade into the normal tissue using acidosis. Simulations suggest that diffusible cytotoxic treatments, such as chemotherapy, may increase the metabolic aggressiveness of a tumor due to drug-mediated selection. Chemotherapy removes the metabolic stratification of the tumor and allows more aggressive cells to grow toward blood vessels and normal tissue. Antiangiogenic therapy also selects for aggressive phenotypes due to degradation of the tumor microenvironment, ultimately resulting in a more invasive tumor. In contrast, pH buffer therapy slows down the development of aggressive tumors, but only if administered when the tumor is still stratified. Overall, findings from this model highlight the risks of cytotoxic and antiangiogenic treatments in the context of tumor heterogeneity resulting from a selection for more aggressive behaviors.

[1]  Jacob G Scott,et al.  Microenvironmental variables must influence intrinsic phenotypic parameters of cancer stem cells to affect tumourigenicity , 2014, PLoS Comput. Biol..

[2]  K. Schaller,et al.  'Go or grow': the key to the emergence of invasion in tumour progression? , 2012, Mathematical medicine and biology : a journal of the IMA.

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

[4]  C. Maley,et al.  Dispersal Evolution in Neoplasms: The Role of Disregulated Metabolism in the Evolution of Cell Motility , 2011, Cancer Prevention Research.

[5]  John S. Lowengrub,et al.  A New Ghost Cell/Level Set Method for Moving Boundary Problems: Application to Tumor Growth , 2008, J. Sci. Comput..

[6]  A. Anderson,et al.  A hybrid cellular automaton model of clonal evolution in cancer: the emergence of the glycolytic phenotype. , 2008, Journal of theoretical biology.

[7]  M G Nichols,et al.  Oxygen diffusion and reaction kinetics in the photodynamic therapy of multicell tumour spheroids. , 1994, Physics in medicine and biology.

[8]  A. Anderson,et al.  Evolution, regulation and disruption of homeostasis and its role in carcinogenesis , 2010 .

[9]  T. Ohtsubo,et al.  Acidic environment causes apoptosis by increasing caspase activity , 1999, British Journal of Cancer.

[10]  Jacob G. Scott,et al.  Mathematical Modeling of the Metastatic Process , 2013, 1305.4622.

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

[12]  P. Gerlee,et al.  Diffusion-limited tumour growth: simulations and analysis. , 2010, Mathematical biosciences and engineering : MBE.

[13]  S. Torquato,et al.  A Cellular Automaton Model for Tumor Dormancy: Emergence of a Proliferative Switch , 2014, PloS one.

[14]  Alexander R. A. Anderson,et al.  Modelling evolutionary cell behaviour using neural networks: Application to tumour growth , 2009, Biosyst..

[15]  K. Groebe,et al.  Glucose diffusion coefficients determined from concentration profiles in EMT6 tumor spheroids incubated in radioactively labeled L-glucose. , 1994, Advances in experimental medicine and biology.

[16]  V. P. Collins,et al.  Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics , 2013, Proceedings of the National Academy of Sciences.

[17]  H. Frieboes,et al.  An integrated computational/experimental model of tumor invasion. , 2006, Cancer research.

[18]  Robert J Gillies,et al.  Acidity generated by the tumor microenvironment drives local invasion. , 2013, Cancer research.

[19]  Peter Kuhn,et al.  Spreaders and sponges define metastasis in lung cancer: a Markov chain Monte Carlo mathematical model. , 2013, Cancer research.

[20]  P Hahnfeldt,et al.  Migration rules: tumours are conglomerates of self-metastases , 2009, British Journal of Cancer.

[21]  A. Anderson,et al.  A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion. , 2005, Mathematical medicine and biology : a journal of the IMA.

[22]  K. Zierler,et al.  Whole body glucose metabolism. , 1999, The American journal of physiology.

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

[24]  J F Gross,et al.  Morphologic and hemodynamic comparison of tumor and healing normal tissue microvasculature. , 1989, International journal of radiation oncology, biology, physics.

[25]  K. Macleod,et al.  Abstract 3938: BNip3 limits HIF-1α stabilization and metastasis in a mouse model of breast cancer through effects on mitochondrial integrity and ROS generation , 2012 .

[26]  Sunwoo Park,et al.  A computational approach to resolve cell level contributions to early glandular epithelial cancer progression , 2009, BMC Systems Biology.

[27]  Robert A. Weinberg,et al.  Metastasis genes: A progression puzzle , 2002, Nature.

[28]  Christopher DuBois,et al.  A multiple time-scale computational model of a tumor and its micro environment. , 2012, Mathematical biosciences and engineering : MBE.

[29]  Sarah DiGiulio ASCO13: What ‘Building Bridges’ Means for Oncology , 2013 .

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

[31]  Salvatore Torquato,et al.  Modeling the effects of vasculature evolution on early brain tumor growth. , 2006, Journal of theoretical biology.

[32]  Robert Axelrod,et al.  Ecological therapy for cancer: defining tumors using an ecosystem paradigm suggests new opportunities for novel cancer treatments. , 2008, Translational oncology.

[33]  C. Vogel,et al.  Bevacizumab in the Treatment of Metastatic Breast Cancer: Friend or Foe? , 2012, Current Oncology Reports.

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

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

[36]  A. Anderson,et al.  Evolution of cell motility in an individual-based model of tumour growth. , 2009, Journal of theoretical biology.

[37]  R. Jain,et al.  Intussusceptive microvascular growth in a human colon adenocarcinoma xenograft: a novel mechanism of tumor angiogenesis. , 1996, Microvascular research.

[38]  A. Anderson,et al.  Senescent fibroblasts in melanoma initiation and progression: an integrated theoretical, experimental, and clinical approach. , 2013, Cancer research.

[39]  Gibin G Powathil,et al.  Modelling the effects of cell-cycle heterogeneity on the response of a solid tumour to chemotherapy: biological insights from a hybrid multiscale cellular automaton model. , 2012, Journal of theoretical biology.

[40]  Robert B. Heckendorn,et al.  Cellular automaton simulation examining progenitor hierarchy structure effects on mammary ductal carcinoma in situ. , 2007, Journal of theoretical biology.

[41]  Jie H. S. Zhou,et al.  Stretched cell cycle model for proliferating lymphocytes , 2014, Proceedings of the National Academy of Sciences.

[42]  David Basanta,et al.  The role of transforming growth factor-beta-mediated tumor-stroma interactions in prostate cancer progression: an integrative approach. , 2009, Cancer research.

[43]  P. A. Futreal,et al.  Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. , 2012, The New England journal of medicine.

[44]  Yi Jiang,et al.  A cell-based model exhibiting branching and anastomosis during tumor-induced angiogenesis. , 2007, Biophysical journal.

[45]  Alissa M. Weaver,et al.  Tumor Morphology and Phenotypic Evolution Driven by Selective Pressure from the Microenvironment , 2006, Cell.

[46]  J. Salk Clonal evolution in cancer , 2010 .

[47]  Stephanie Forrest,et al.  Modeling Somatic Evolution in Tumorigenesis , 2006, PLoS Comput. Biol..

[48]  C. Swanton Intratumor heterogeneity: evolution through space and time. , 2012, Cancer research.

[49]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer-associated genes , 2013 .

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

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

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

[53]  R. Gillies,et al.  Impact of metabolic heterogeneity on tumor growth, invasion, and treatment outcomes , 2016 .

[54]  David Basanta,et al.  An integrated computational model of the bone microenvironment in bone-metastatic prostate cancer. , 2014, Cancer research.

[55]  S Torquato,et al.  Emergence of a Subpopulation in a Computational Model of Tumor Growth , 2022 .

[56]  R. Jain Normalization of Tumor Vasculature: An Emerging Concept in Antiangiogenic Therapy , 2005, Science.

[57]  Carlo C. Maley,et al.  Clonal evolution in cancer , 2012, Nature.

[58]  Laura Esserman,et al.  Development of a novel method for measuring in vivo breast epithelial cell proliferation in humans , 2005, Breast Cancer Research and Treatment.

[59]  Bonnie F. Sloane,et al.  Bicarbonate increases tumor pH and inhibits spontaneous metastases. , 2009, Cancer research.

[60]  Axel R. Pries,et al.  Angiogenesis: An Adaptive Dynamic Biological Patterning Problem , 2013, PLoS Comput. Biol..

[61]  D. McElwain,et al.  The role of mechanical host-tumour interactions in the collapse of tumour blood vessels and tumour growth dynamics. , 2006, Journal of theoretical biology.

[62]  G. Melillo,et al.  Overcoming disappointing results with antiangiogenic therapy by targeting hypoxia , 2012, Nature Reviews Clinical Oncology.

[63]  Vittorio Cristini,et al.  Patient-calibrated agent-based modelling of ductal carcinoma in situ (DCIS): from microscopic measurements to macroscopic predictions of clinical progression. , 2012, Journal of theoretical biology.

[64]  B. Grammaticos,et al.  A Model for Short- and Long-range Interactions of Migrating Tumour Cell , 2008, Acta biotheoretica.

[65]  A. Deutsch,et al.  Evolutionary game theory elucidates the role of glycolysis in glioma progression and invasion , 2008, Cell proliferation.

[66]  David Basanta,et al.  Exploiting ecological principles to better understand cancer progression and treatment , 2013, Interface Focus.

[67]  Robert J. Gillies,et al.  Multiscale Modelling of Vascular Tumour Growth in 3D: The Roles of Domain Size and Boundary Conditions , 2011, PloS one.

[68]  E. T. Gawlinski,et al.  Acid-mediated tumor invasion: a multidisciplinary study. , 2006, Cancer research.

[69]  P. Maini,et al.  Modelling aspects of cancer dynamics: a review , 2006, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[70]  Bonnie F. Sloane,et al.  Chronic autophagy is a cellular adaptation to tumor acidic pH microenvironments. , 2012, Cancer research.

[71]  R. Gillies,et al.  A quantitative theoretical model for the development of malignancy in ductal carcinoma in situ. , 2010, Journal of theoretical biology.

[72]  Robert J Gillies,et al.  Systemic buffers inhibit carcinogenesis in TRAMP mice. , 2012, The Journal of urology.

[73]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer genes , 2014 .

[74]  Andreas Deutsch,et al.  Computational analysis of the influence of the microenvironment on carcinogenesis. , 2011, Mathematical biosciences.

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

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

[77]  Glen E. P. Ropella,et al.  Essential operating principles for tumor spheroid growth , 2008, BMC Systems Biology.