Mathematical modeling links Wnt signaling to emergent patterns of metabolism in colon cancer

Cell-intrinsic metabolic reprogramming is a hallmark of cancer that provides anabolic support to cell proliferation. How reprogramming influences tumor heterogeneity or drug sensitivities is not well understood. Here, we report a self-organizing spatial pattern of glycolysis in xenograft colon tumors where pyruvate dehydrogenase kinase (PDK1), a negative regulator of oxidative phosphorylation, is highly active in clusters of cells arranged in a spotted array. To understand this pattern, we developed a reaction- diffusion model that incorporates Wnt signaling, a pathway known to upregulate PDK1 and Warburg metabolism. Partial interference with Wnt alters the size and intensity of the spotted pattern in tumors and in the model. The model predicts that Wnt inhibition should trigger an increase in proteins that enhance the range of Wnt ligand diffusion. Not only was this prediction validated in xenograft tumors but similar patterns also emerge in radiochemotherapy-treated colorectal cancer. The model also predicts that inhibitors that target glycolysis or Wnt signaling in combination should synergize and be more effective than each treatment individually. We validated this prediction in 3D colon tumor spheroids.