Bridging cell-scale simulations and radiologic images to explain short-time intratumoral oxygen fluctuations

Radiologic images provide a way to monitor tumor development and its response to therapies in a longitudinal and minimally invasive fashion. However, they operate on a macroscopic scale (average value per voxel) and are not able to capture microscopic scale (cell level) phenomena. Nevertheless, to examine the causes of frequent fast fluctuations in tissue oxygenation, the models simulating individual cells’ behavior are needed. Here, we provided a link between the average data value recorded for radiologic image voxels and the cellular and vascular architecture of the tissue that fills these voxels. Using hybrid agent-based modeling, we generated a set of tissue morphologies capable of reproducing tissue oxygenation levels observed in radiologic images. We applied this approach to investigate whether oxygen fluctuations can be explained by changes in vascular oxygen supply or by modulations in cellular oxygen absorption. Our studies showed that intravascular changes in oxygen supply can reproduce the observed fluctuations in tissue oxygenation in all considered regions of interest. However, large magnitude fluctuations cannot be recreated by modifications in cellular absorption of oxygen in biologically feasible manner. Additionally, we developed a procedure to identify plausible tissue morphologies for a given temporal series of average data from radiology images. In future applications this approach can be used to generate a set of tissues representative for radiology images and to simulate tumor response to various anti-cancer treatments on the tissue-scale level. Authors Summary Low levels of oxygen, called hypoxia, are observable in many solid tumors. They are associated with more aggressive malignant cells which are resistant to chemo-, radio- and immunotherapies. Recently developed imaging techniques provide a way to measure the magnitude of frequent short-term oxygen fluctuation, however they operate on a macro-scale voxel level. To examine the causes of rapid oxygen fluctuations on the cell level, we developed a hybrid agent-based mathematical model. We tested two different mechanisms that could be responsible for these cyclic effects in tissue oxygenation: variations in vascular influx of oxygen and modulations in cellular oxygen absorption. Additionally, we developed a procedure to identify plausible tissue morphologies from data collected from radiological images. This will also provide a bridge between the micro-scale simulations with individual cells and the longitudinal medical images containing average voxel values. In the future applications, this approach can be used to generate a set of tissues representative of radiology images and to simulate tumor response to various anticancer treatments on the cell-scale level.