The Microvascular Network Density and Morphology Affect the Radiotherapy Outcome at the Microscale Level: A Computational Analysis.

Abstract

Half of the cancer patients are treated with radiotherapy to obtain tumor control while minimizing radiation toxicity. In this context, hypoxia is known to determine treatment resistance, driving tumor relapse after treatment. Therefore, we investigated the role of microvasculature density and morphology in shaping the tissue oxygen distribution and consequently affecting the treatment outcome. We developed an advanced computational model to describe oxygen delivery in the vascular network and the surrounding tissue. A peculiar aspect of the model is its mesoscale mixed-dimensional approach, which allows the explicit inclusion of vascular network geometry and the description of the red blood cells' effect. The oxygen delivery is modeled considering both diffusive and convective phenomena. A 30 × 2GyRBE treatment is delivered in silico simulating photons, protons, and carbon ions. First, we estimated the surviving fraction by the classical linear quadratic model modified to account for the oxygen effect. Then, leveraging the 3D description of the surviving fraction at the microenvironment scale, we compute the local tumor control probability (LTCP) in different oxygenation states (reference, acute hypoxia, high oxygen consumption). We report correlations between the LTCP and the hypoxic volume fractions (with pO2 lower than 1 mmHg) starting with photons. These hypoxic regions are present locally, even in highly vascularized tissue, if the network is not uniformly distributed, as it might be in cancer. They are also present in tissue with low microvascular density, even with regular morphology. Interestingly, the domain considered is comparable to or smaller than the clinical imaging standard voxel dimension, and the average oxygen partial pressure in the tissue region fails to spot treatments with low LTCP, questioning whether these hypoxia areas are visible clinically via imaging. Protons have a similar effect, highlighting a similar behavior across the oxygenation levels at the microscale level. Finally, carbon ions seem more effective than photons and protons in the presence of hypoxia due to the lower oxygen effect at high LET. For this reason, the treatment with carbon ions results in high LTCPs, whatever vascular network is considered (density and morphology differences). These results show the effect of microvascular density and regularity on the radiotherapy outcome and they help us understand how the microvascular network morphology affects tumor oxygenation and the radiotherapy outcome.