A mesoscopic simulator to uncover heterogeneity and evolutionary dynamics in tumors

Increasingly complex in-silico modeling approaches offer a way to simultaneously access cancerous processes at different spatio-temporal scales. High-level models, such as those based on partial differential equations, are computationally affordable and allow large tumor sizes and long temporal windows to be studied, but miss the discrete nature of many key underlying cellular processes. Individual-based approaches provide a much more detailed description of tumors, but have difficulties when trying to handle full-sized real cancers. Thus, there exists a trade-off between the integration of macroscopic and microscopic information, now widely available, and the ability to attain clinical tumor sizes. In this paper we put forward a stochastic mesoscopic simulation framework that incorporates key cellular processes during tumor progression while keeping computational costs to a minimum. Our framework captures a physical scale that allows both the incorporation of microscopic information, tracking the spatio-temporal emergence of tumor heterogeneity and the underlying evolutionary dynamics, and the reconstruction of clinically sized tumors from high-resolution medical imaging data, with the additional benefit of low computational cost. We illustrate the functionality of our modeling approach for the case of glioblastoma, a paradigm of tumor heterogeneity that remains extremely challenging in the clinical setting. Author summary Computer simulation based on mathematical models provides a way to improve the understanding of complex processes in oncology. In this paper we develop a stochastic mesoscopic simulation approach that incorporates key cellular processes while keeping computational costs to a minimum. Our methodology captures the development of tumor heterogeneity and the underlying evolutionary dynamics. The physical scale considered allows microscopic information to be included, tracking the spatio-temporal evolution of tumor heterogeneity and reconstructing clinically sized tumors from high-resolution medical imaging data, with a low computational cost. We illustrate the functionality of the modeling approach for the case of glioblastoma, an epitome of heterogeneity in tumors.