Transport of Liposome Encapsulated Drugs in Voxelized Computational Model of Human Brain Tumors

Image-based predictive modelling frameworks for personalised drug delivery in cancer therapy

Personalised drug delivery enables a tailored treatment plan for each patient compared to conventional drug delivery, where a generic strategy is commonly employed. It can not only achieve precise treatment to improve effectiveness but also reduce the risk of adverse effects to improve patients' quality of life. Drug delivery involves multiple interconnected physiological and physicochemical processes, which span a wide range of time and length scales. How to consider the impact of individual differences on these processes becomes critical. Multiphysics models are an open system that allows well-controlled studies on the individual and combined effects of influencing factors on drug delivery outcomes while accommodating the patient-specific in vivo environment, which is not economically feasible through experimental means. Extensive modelling frameworks have been developed to reveal the underlying mechanisms of drug delivery and optimise effective delivery plans. This review provides an overview of currently available models, their integration with advanced medical imaging modalities, and code packages for personalised drug delivery. The potential to incorporate new technologies (i.e., machine learning) in this field is also addressed for development.

Predicting intratumoral fluid pressure and liposome accumulation using physics informed deep learning

Liposome-based anticancer agents take advantage of the increased vascular permeability and transvascular pressure gradients for selective accumulation in tumors, a phenomenon known as the enhanced permeability and retention(EPR) effect. The EPR effect has motivated the clinical use of nano-therapeutics, with mixed results on treatment outcome. High interstitial fluid pressure (IFP) has been shown to limit liposome drug delivery to central tumour regions. Furthermore, high IFP is an independent prognostic biomarker for treatment efficacy in radiation therapy and chemotherapy for some solid cancers. Therefore, accurately measuring spatial liposome accumulation and IFP distribution within a solid tumour is crucial for optimal treatment planning. In this paper, we develop a model capable of predicting voxel-by-voxel intratumoral liposome accumulation and IFP using pre and post administration imaging. Our approach is based on physics informed machine learning, a novel technique combining machine learning and partial differential equations. through application to a set of mouse data and a set of synthetically-generated tumours, we show that our approach accurately predicts the spatial liposome accumulation and IFP for an individual tumour while relying on minimal information. This is an important result with applications for forecasting tumour progression and designing treatment.

Evaluating the reliability of tumour spheroid-on-chip models for replicating intratumoural drug delivery: considering the role of microfluidic parameters

Several tumour spheroid-on-chip models have already been proposed in the literature to conduct high throughput drug screening assays. The microfluidic configurations in these models generally depend on the strategies adopted for spheroid formation and entrapment. However, it is not clear how successful they are to mimic in vivo transport mechanisms. In this study, drug transport in different tumour spheroid-on-chip models is numerically investigated under static and dynamic conditions using porous media theory. Moreover, the treatment of a solid tumour at the initial stage of development is modelled using bolus injection and continuous infusion methods. Then, the results of tumour spheroid-on-chip, including drug concentration, cell viability, as well as pressure and fluid shear stress distributions, are compared with those of the solid tumour, assuming identical transport properties in all models. Finally, a new configuration of the microfluidic device along with the optimal drug concentrations is proposed, which can well imitate a given in vivo situation.

Drug delivery to solid tumors with heterogeneous microvascular networks: Novel insights from image-based numerical modeling

The present study examines chemotherapy by incorporating multi-scale mathematical modeling to predict drug delivery and its effects. This approach leads to a more-realistic physiological tumor model than is possible with previous approaches, as it obtains the capillary network geometry from an image, and also considers the tumor's necrotic core, drug binding, and cellular uptake. Modeling of the fluid flow and drug transport is then performed in the extracellular matrix. The results demonstrate a 10% drop in the fraction of killed cancer cells 69% rather than the 79% reported earlier for a tumor of similar geometry a more-accurate value. This study examines how tumor-related parameters including the necrotic core size and tumor size, and also drug-related parameters drug dosage, binding affinity of drug, and drug degradation can affect the delivery of the drug to solid tumors. Results indicate that concentration of drug are high in the tumor, low in normal tissue, and remarkably low in the necrotic core. Results also offer a treatment of tumors with smaller necrotic core. Tumor size, which implies the tumor progression, has a considerable impact on treatment outcomes, so to be more effective, treatment should be applied at a specific size of tumor. It is demonstrated that binding affinity of drugs to cell-surface receptors and drug dosage have significant impact on treatment efficacy, so they should be regulated based on a balanced quantification between maximum treatment efficacy and minimum side effects. On the other hand, considering the effects of drug degradation in the model has not significant effect on treatment efficacy. The findings of the present study provide insight into the mechanism of drug delivery to solid tumors based on analyzing the effective parameters and modeling how their behavior in the tumor microenvironment affects treatment efficacy.

Computational Modeling of Interstitial Fluid Pressure and Velocity in Non-small Cell Lung Cancer Brain Metastases Treated With Stereotactic Radiosurgery

Background: Early imaging-based treatment response assessment of brain metastases following stereotactic radiosurgery (SRS) remains challenging. The aim of this study is to determine whether early (within 12 weeks) intratumoral changes in interstitial fluid pressure (IFP) and velocity (IFV) estimated from computational fluid modeling (CFM) using dynamic contrast-enhanced (DCE) MRI can predict long-term outcomes of lung cancer brain metastases (LCBMs) treated with SRS.

Methods: Pre- and post-treatment T₁-weighted DCE-MRI data were obtained in 41 patients treated with SRS for intact LCBMs. The imaging response was assessed using RANO-BM criteria. For each lesion, extravasation of contrast agent measured from Extended Tofts pharmacokinetic Model (volume transfer constant, K^trans) was incorporated into a computational fluid model to estimate tumor IFP and IFV. Estimates of mean IFP and IFV and heterogeneity (skewness and kurtosis) were calculated for each lesion from pre- and post-SRS imaging. The Wilcoxon rank-sum test was utilized to assess for significant differences in IFP, IFV, and IFP/IFV change (Δ) between response groups.

Results: Fifty-three lesions from 41 patients were included. Median follow-up time after SRS was 11 months. The objective response (OR) rate (partial or complete response) was 79%, with 21% demonstrating stable disease (SD) or progressive disease (PD). There were significant response group differences for multiple posttreatment and Δ CFM parameters: post-SRS IFP skewness (mean −0.405 vs. −0.691, p = 0.022), IFP kurtosis (mean 2.88 vs. 3.51, p = 0.024), and IFV mean (5.75e-09 vs. 4.19e-09 m/s, p = 0.027); and Δ IFP kurtosis (mean −2.26 vs. −0.0156, p = 0.017) and IFV mean (1.91e-09 vs. 2.38e-10 m/s, p = 0.013). Posttreatment and Δ thresholds predicted non-OR with high sensitivity (sens): post-SRS IFP skewness (−0.432, sens 84%), kurtosis (2.89, sens 84%), and IFV mean (4.93e-09 m/s, sens 79%); and Δ IFP kurtosis (−0.469, sens 74%) and IFV mean (9.90e-10 m/s, sens 74%).

Conclusions: Objective response was associated with lower post-treatment tumor heterogeneity, as represented by reductions in IFP skewness and kurtosis. These results suggest that early post-treatment assessment of IFP and IFV can be used to predict long-term response of lung cancer brain metastases to SRS, allowing a timelier treatment modification.

Comparison of transport of chemotherapeutic drugs in voxelized heterogeneous model of human brain tumor

Systemic administration of chemotherapeutic drugs is widely used in the treatment of cancer. However, a good understanding of drug transport barriers that influence the treatment efficacy is still lacking. In this study, a voxelized numerical model based on dynamic contrast enhanced-magnetic resonance imaging (DCE-MRI) and computational fluid dynamics (CFD) is employed to study the transport and efficacy of three different chemotherapeutic drugs, namely methotrexate, doxorubicin and cisplatin in human brain tumors. DCE-MRI data provides realistic heterogeneous vasculature of the tumor, the permeability of tissue to contrast agent, interstitial volume fraction (porosity) of the tissue and patient-specific arterial input function (AIF). The permeability of tissue to aforementioned drugs is determined by correlating it with the permeability of tissue to the contrast agent. The model is employed to simulate drug concentration in the tissue and compare the effect of heterogeneous vasculature on the distribution of the drugs in the tumor. The drug accumulation is observed to be higher in high permeability areas initially, and in higher porosity areas at later times. Furthermore, it is observed that methotrexate remains in the interstitial space of the tumor in higher concentration for a longer duration as compared to other two drugs, facilitating more tumor cell killing.

Convection enhanced delivery of liposome encapsulated doxorubicin for brain tumour therapy

Convection enhanced delivery is promising to overcome the blood brain barrier. However, the treatment is less efficient in clinic due to the rapid elimination of small molecular drugs in brain tumours. In this study, numerical simulation is applied to investigate the convection enhanced delivery of liposome encapsulated doxorubicin under various conditions, based on a 3-D brain tumour model that is reconstructed from magnetic resonance images. Treatment efficacy is evaluated in terms of the tumour volume where the free doxorubicin concentration is above LD90. Simulation results denote that intracerebral infusion is effective in increasing the interstitial fluid velocity and inhibiting the fluid leakage from blood around the infusion site. Comparisons with direct doxorubicin infusion demonstrate the advantages of liposomes in enhancing the doxorubicin accumulation and penetration in the brain tumour. Delivery outcomes are determined by both the intratumoural environment and properties of therapeutic agents. The treatment efficacy can be improved by either increasing the liposome solution concentration and infusion rate, administrating liposomes in the tumour with normalised microvasculature density, or using liposomes with low vascular permeability. The delivery is less sensitive to liposome diffusivity in the examined range (E-11~E-7 cm²/s) as convective transport is dominative in determining the liposome migration. Drug release rate is able to be optimised by keeping a trade-off between enhancing the drug penetration and providing sufficient free doxorubicin for effective cell killing. Results from this study can be used to improve the regimen of CED treatments.

Computational modelling of drug delivery to solid tumour: Understanding the interplay between chemotherapeutics and biological system for optimised delivery systems

Drug delivery to solid tumour involves multiple physiological, biochemical and biophysical processes taking place across a wide range of length and time scales. The therapeutic efficacy of anticancer drugs is influenced by the complex interplays among the intrinsic properties of tumours, biophysical aspects of drug transport and cellular uptake. Mathematical and computational modelling allows for a well-controlled study on the individual and combined effects of a wide range of parameters on drug transport and therapeutic efficacy, which would not be possible or economically viable through experimental means. A wide spectrum of mathematical models has been developed for the simulation of drug transport and delivery in solid tumours, including PK/PD-based compartmental models, microscopic and macroscopic transport models, and molecular dynamics drug loading and release models. These models have been used as a tool to identify the limiting factors and for optimal design of efficient drug delivery systems. This article gives an overview of the currently available computational models for drug transport in solid tumours, together with their applications to novel drug delivery systems, such as nanoparticle-mediated drug delivery and convection-enhanced delivery.

Numerical Study of Transport of Anticancer Drugs in Heterogeneous Vasculature of Human Brain Tumors Using Dynamic Contrast Enhanced-Magnetic Resonance Imaging

Systemic administration of drugs in tumors is a challenging task due to unorganized microvasculature and nonuniform extravasation. There is an imperative need to understand the transport behavior of drugs when administered intravenously. In this study, a transport model is developed to understand the therapeutic efficacy of a free drug and liposome-encapsulated drugs (LED), in heterogeneous vasculature of human brain tumors. Dynamic contrast enhanced-magnetic resonance imaging (DCE-MRI) data is employed to model the heterogeneity in tumor vasculature that is directly mapped onto the computational fluid dynamics (CFD) model. Results indicate that heterogeneous vasculature leads to preferential accumulation of drugs at the tumor position. Higher drug accumulation was found at location of higher interstitial volume, thereby facilitating more tumor cell killing at those areas. Liposome-released drug (LRD) remains inside the tumor for longer time as compared to free drug, which together with higher concentration enhances therapeutic efficacy. The interstitial as well as intracellular concentration of LRD is found to be 2–20 fold higher as compared to free drug, which are in line with experimental data reported in literature. Close agreement between the predicted and experimental data demonstrates the potential of the developed model in modeling the transport of LED and free drugs in heterogeneous vasculature of human tumors.