Cancer metabolism has received renewed interest as a potential target for cancer therapy. In this study, we use a multi-scale modeling approach to interrogate the implications of three metabolic scenarios of potential clinical relevance: the Warburg effect, the reverse Warburg effect and glutamine addiction. At the intracellular level, we construct a network of central metabolism and perform flux balance analysis (FBA) to estimate metabolic fluxes; at the cellular level, we exploit this metabolic network to calculate parameters for a coarse-grained description of cellular growth kinetics; and at the multicellular level, we incorporate these kinetic schemes into the cellular automata of an agent-based model (ABM), iDynoMiCS. This ABM evaluates the reaction-diffusion of the metabolites, cellular division and motion over a simulation domain. Our multi-scale simulations suggest that the Warburg effect provides a growth advantage to the tumor cells under resource limitation. However, we identify a non-monotonic dependence of growth rate on the strength of glycolytic pathway. On the other hand, the reverse Warburg scenario provides an initial growth advantage in tumors that originate deeper in the tissue. The metabolic profile of stromal cells considered in this scenario allows more oxygen to reach the tumor cells in the deeper tissue and thus promotes tumor growth at earlier stages. Lastly, we suggest that glutamine addiction does not confer a selective advantage to tumor growth with glutamine acting as a carbon source in the tricarboxylic acid (TCA) cycle, any advantage of glutamine uptake must come through other pathways not included in our model (e.g., as a nitrogen donor). Our analysis illustrates the importance of accounting explicitly for spatial and temporal evolution of tumor microenvironment in the interpretation of metabolic scenarios and hence provides a basis for further studies, including evaluation of specific therapeutic strategies that target metabolism.
Cancer metabolism is an emerging hallmark of cancer. In the past decade, a renewed focus on cancer metabolism has led to several distinct hypotheses describing the role of metabolism in cancer. To complement experimental efforts in this field, a scale-bridging computational framework is needed to allow rapid evaluation of emerging hypotheses in cancer metabolism. In this study, we present a multi-scale modeling platform and demonstrate the distinct outcomes in population-scale growth dynamics under different metabolic scenarios: the Warburg effect, the reverse Warburg effect and glutamine addiction. Within this modeling framework, we confirmed population-scale growth advantage enabled by the Warburg effect, provided insights into the symbiosis between stromal cells and tumor cells in the reverse Warburg effect and argued that the anaplerotic role of glutamine is not exploited by tumor cells to gain growth advantage under resource limitations. We point to the opportunity for this framework to help understand tissue-scale response to therapeutic strategies that target cancer metabolism while accounting for the tumor complexity at multiple scales.