Modelling the transport of fluid through heterogeneous, whole tumours in silico

Cancers exhibit spatially heterogeneous, unique vascular architectures across individual samples, cell-lines and patients. This inherently disorganised collection of leaky blood vessels contribute significantly to suboptimal treatment efficacy. Preclinical tools are urgently required which incorporate the inherent variability and heterogeneity of tumours to optimise and engineer anti-cancer therapies. In this study, we present a novel computational framework which incorporates whole, realistic tumours extracted ex vivo to efficiently simulate vascular blood flow and interstitial fluid transport in silico for validation against in vivo biomedical imaging. Our model couples Poiseuille and Darcy descriptions of vascular and interstitial flow, respectively, and incorporates spatially heterogeneous blood vessel lumen and interstitial permeabilities to generate accurate predictions of tumour fluid dynamics. Our platform enables highly-controlled experiments to be performed which provide insight into how tumour vascular heterogeneity contributes to tumour fluid transport. We detail the application of our framework to an orthotopic murine glioma (GL261) and a human colorectal carcinoma (LS147T), and perform sensitivity analysis to gain an understanding of the key biological mechanisms which determine tumour fluid transport. Finally we mimic vascular normalization by modifying parameters, such as vascular and interstitial permeabilities, and show that incorporating realistic vasculatures is key to modelling the contrasting fluid dynamic response between tumour samples. Contrary to literature, we show that reducing tumour interstitial fluid pressure is not essential to increase interstitial perfusion and that therapies should seek to develop an interstitial fluid pressure gradient. We also hypothesise that stabilising vessel diameters and permeabilities are not key responses following vascular normalization and that therapy may alter interstitial hydraulic conductivity. Consequently, we suggest that normalizing the interstitial microenvironment may provide a more effective means to increase interstitial perfusion within tumours.
Author summary The structure of tumours varies widely, with dense and chaotically-formed networks of blood vessels that differ between each individual tumour and even between different regions of the same tumour. This atypical environment can inhibit the delivery of anti-cancer therapies. Computational tools are urgently required which facilitate a deeper understanding of the relationship between blood vessel architectures and therapeutic response. We have developed a computational framework which integrates the complex tumour vascular architecture to predict fluid transport across all lengths scales in whole tumours. We apply our model to two tumour cell-lines and show that differences in their inherent vascular structures influence flow through cancerous tissue. We also use our platform to predict the fluid dynamic response following vascular normalization therapy in realistic, static tumour networks and show that the response is dependent on tumour vascular architecture. We hypothesise that therapy may alter the permeability of interstitial tissue to fluid transport and show that lowering interstitial fluid pressure is not a necessary therapeutic outcome to increase tumour perfusion.