A four-compartment multiscale model of fluid and drug distribution in vascular tumours

A three-dimensional, discrete-continuum model of blood pressure in microvascular networks

We present a 3D discrete-continuum model to simulate blood pressure in large microvascular tissues in the absence of known capillary network architecture. Our hybrid approach combines a 1D Poiseuille flow description for large, discrete arteriolar and venular networks coupled to a continuum-based Darcy model, point sources of flux, for transport in the capillary bed. We evaluate our hybrid approach using a vascular network imaged from the mouse brain medulla/pons using multi-fluorescence high-resolution episcopic microscopy (MF-HREM). We use the fully-resolved vascular network to predict the hydraulic conductivity of the capillary network and generate a fully-discrete pressure solution to benchmark against. Our results demonstrate that the discrete-continuum methodology is a computationally feasible and effective tool for predicting blood pressure in real-world microvascular tissues when capillary microvessels are poorly defined.

A multiscale orchestrated computational framework to reveal emergent phenomena in neuroblastoma

Neuroblastoma is a complex and aggressive type of cancer that affects children. Current treatments involve a combination of surgery, chemotherapy, radiotherapy, and stem cell transplantation. However, treatment outcomes vary due to the heterogeneous nature of the disease. Computational models have been used to analyse data, simulate biological processes, and predict disease progression and treatment outcomes. While continuum cancer models capture the overall behaviour of tumours, and agent-based models represent the complex behaviour of individual cells, multiscale models represent interactions at different organisational levels, providing a more comprehensive understanding of the system. In 2018, the PRIMAGE consortium was formed to build a cloud-based decision support system for neuroblastoma, including a multi-scale model for patient-specific simulations of disease progression. In this work we have developed this multi-scale model that includes data such as patient's tumour geometry, cellularity, vascularization, genetics and type of chemotherapy treatment, and integrated it into an online platform that runs the simulations on a high-performance computation cluster using Onedata and Kubernetes technologies. This infrastructure will allow clinicians to optimise treatment regimens and reduce the number of costly and time-consuming clinical trials. This manuscript outlines the challenging framework's model architecture, data workflow, hypothesis, and resources employed in its development.

Human brain solute transport quantified by glymphatic MRI-informed biophysics during sleep and sleep deprivation

Whether you are reading, running or sleeping, your brain and its fluid environment continuously interacts to distribute nutrients and clear metabolic waste. Yet, the precise mechanisms for solute transport within the human brain have remained hard to quantify using imaging techniques alone. From multi-modal human brain MRI data sets in sleeping and sleep-deprived subjects, we identify and quantify CSF tracer transport parameters using forward and inverse subject-specific computational modelling. Our findings support the notion that extracellular diffusion alone is not sufficient as a brain-wide tracer transport mechanism. Instead, we show that human MRI observations align well with transport by either by an effective diffusion coefficent 3.5[Formula: see text] that of extracellular diffusion in combination with local clearance rates corresponding to a tracer half-life of up to 5 h, or by extracellular diffusion augmented by advection with brain-wide average flow speeds on the order of 1-9 [Formula: see text]m/min. Reduced advection fully explains reduced tracer clearance after sleep-deprivation, supporting the role of sleep and sleep deprivation on human brain clearance.

Multi-compartmental model of glymphatic clearance of solutes in brain tissue

The glymphatic system is the subject of numerous pieces of research in biology. Mathematical modelling plays a considerable role in this field since it can indicate the possible physical effects of this system and validate the biologists' hypotheses. The available mathematical models that describe the system at the scale of the brain (i.e. the macroscopic scale) are often solely based on the diffusion equation and do not consider the fine structures formed by the perivascular spaces. We therefore propose a mathematical model representing the time and space evolution of a mixture flowing through multiple compartments of the brain. We adopt a macroscopic point of view in which the compartments are all present at any point in space. The equations system is composed of two coupled equations for each compartment: One equation for the pressure of a fluid and one for the mass concentration of a solute. The fluid and solute can move from one compartment to another according to certain membrane conditions modelled by transfer functions. We propose to apply this new modelling framework to the clearance of 14C-inulin from the rat brain.

Computational model of silica nanoparticle penetration into tumor spheroids: Effects of methoxy and carboxy PEG surface functionalization and hyperthermia

Drug delivery to tumors suffers from poor solubility, specificity, diffusion through the tumor micro-environment and nonoptimal interactions with components of the extracellular matrix and cell surface receptors. Nanoparticles and drug-polymer complexes address many of these problems. However, large size exasperates the problem of slow diffusion through the tumor. Three-dimensional tumor spheroids are good models to evaluate approaches to mitigate these difficulties and aid in design strategies to improve the delivery of drugs to treat cancer effectively. Diffusion of drug carriers is highly dependent on cell uptake rate parameters (association/dissociation) and temperature. Hyperthermia increases molecular transport and is known to act synergistically with chemotherapy to improve treatment. This study presents a new inverse estimation approach based on Bayesian probability for estimating nanoparticle cell uptake rates from experiments. The parameters were combined with a finite element computational model of nanoparticle transport under hyperthermia conditions to explore its effect on tumor porosity, diffusion and particle binding (association and dissociation) at cell surfaces. Carboxy-PEG-silane (cPEGSi) nanoparticles showed higher cell uptake compared to methoxy-PEG-silane (mPEGSi) nanoparticles. Simulations were consistent with experimental results from Skov-3 ovarian cancer spheroids. Amorphous silica (cPEGSi) nanoparticles (58 nm) concentrated at the periphery of the tumor spheroids at 37°C but mild hyperthermia (43°C) increased nanoparticle penetration. Thus, hyperthermia may enhance cancer treatment by improving blood delivery to tumors, enhancing extravasation and penetration into tumors, trigger release of drug from the carrier at the tumor site and possibly lead to synergistic anti-cancer activity with the drug.

Biologically-Based Mathematical Modeling of Tumor Vasculature and Angiogenesis via Time-Resolved Imaging Data

Tumor-associated vasculature is responsible for the delivery of nutrients, removal of waste, and allowing growth beyond 2-3 mm3. Additionally, the vascular network, which is changing in both space and time, fundamentally influences tumor response to both systemic and radiation therapy. Thus, a robust understanding of vascular dynamics is necessary to accurately predict tumor growth, as well as establish optimal treatment protocols to achieve optimal tumor control. Such a goal requires the intimate integration of both theory and experiment. Quantitative and time-resolved imaging methods have emerged as technologies able to visualize and characterize tumor vascular properties before and during therapy at the tissue and cell scale. Parallel to, but separate from those developments, mathematical modeling techniques have been developed to enable in silico investigations into theoretical tumor and vascular dynamics. In particular, recent efforts have sought to integrate both theory and experiment to enable data-driven mathematical modeling. Such mathematical models are calibrated by data obtained from individual tumor-vascular systems to predict future vascular growth, delivery of systemic agents, and response to radiotherapy. In this review, we discuss experimental techniques for visualizing and quantifying vascular dynamics including magnetic resonance imaging, microfluidic devices, and confocal microscopy. We then focus on the integration of these experimental measures with biologically based mathematical models to generate testable predictions.