In silicomodeling of radiobiological effects using tumor spheroids

Objective.Radiotherapy remains one of the most widely used cancer treatments. However, optimizing its effectiveness requires a deeper understanding of the tumor architecture and its influence on radiation response. Three-dimensional spheroid cultures offer a more physiologically relevant model than traditional 2D assays, and their integration within silicotechniques enable the simulation of both growth dynamics and radiobiological response across multiple scales.Approach.In this study, we developed a multi-scalein silicoframework capable of simulating tumor spheroid growth, using a cellular potts model, and transferring that information to OpenTOPAS to simulate irradiation with any given particle type, where in this work we apply it specifically to a 6 MV photon beam. The dose deposited at each cell was used to determine cell fate via a two-lesion kinetic model (TLKM), which accounts for cellular oxygen concentration, the number of radiation-induced double-strand breaks, their type (simple or complex), and their origin (direct or indirect). The model was applied to two rectal cancer cell lines (SW837 and SW1463) subjected to a short-course radiotherapy regimen of 5 Gy/fraction over 5 consecutive days, replicating the fractionation schedule used in short-course radiotherapy for rectal cancer patients.Main results.The calibrated model reproduces experimental spheroid growth and viability data across fractionated treatment, with oxygen dynamics and hypoxic and necrotic zonation behaving as expected under the model's consumption-diffusion coupling, consistent with, though not independently validated, against spatially-resolved measurements. To our knowledge, this is also the first practical application of the oxygen fixation parameter within the TLKM.Significance.This framework provides a flexible platform for investigating alternative treatment regimens, combination therapies, and particle-based modalities. Expanding the experimental dataset, particularly with spatially resolved measurements of proliferative, hypoxic, and necrotic zones, will be essential to further constrain and validate the model. Incorporating subcellular damage mechanisms and explicit dose-rate effects into the TLKM formulation represents key directions for future development.