Computational Modeling of Interstitial Fluid Pressure and Velocity in Head and Neck Cancer Based on Dynamic Contrast-Enhanced Magnetic Resonance Imaging: Feasibility Analysis

Image-guided patient-specific prediction of interstitial fluid flow and drug transport in solid tumors

Tumor fluid dynamics and drug delivery simulations in solid tumors are highly relevant topics in clinical oncology. The current study introduces a novel method combining computational fluid dynamics (CFD) modeling, quantitative magnetic resonance imaging (MRI; including dynamic contrast-enhanced (DCE) MRI and diffusion-weighted (DW) MRI), and a novel ex-vivo protocol to generate patient-specific models of solid tumors in four patients with peritoneal metastases. DCE-MRI data were analyzed using the extended Tofts model to estimate the spatial distribution of tumor capillary permeability using the Ktrans parameter. DW-MRI data analysis provided a 3D representation of drug diffusivity, and DW-MRI coupled to an ex-vivo measurement protocol informed the spatial heterogeneity of the hydraulic conductivity of tumor tissue. The patient-specific data were subsequently incorporated into a computational fluid dynamics (CFD) model to simulate individualized tumor perfusion and drug transport maps. The results on interstitial fluid flow demonstrated noticeable heterogeneity of interstitial fluid pressure and velocity within the tumor, along with heterogeneous drug penetration profiles among different tumors, even with a similar drug administration regimen.

A Comparative Analysis of Alpha and Beta Therapy in Prostate Cancer Using a 3D Image-Based Spatiotemporal Model

PURPOSE

In treating prostate cancer, distinguishing alpha and beta therapies is vital for efficient radiopharmaceutical delivery. Our study introduces a 3D image-based spatiotemporal computational model that utilizes MRI-derived images to evaluate the efficacy of 225Ac-PSMA and 177Lu-PSMA therapies. We examine the impact of tumor size, diffusion, interstitial fluid pressure (IFP), and interstitial fluid velocity (IFV) on the absorbed doses.

METHODS

An MRI-based geometric model of the tumor and its surrounding environment is initially developed. Subsequently, COMSOL Multiphysics software is utilized to solve convection-diffusion-reaction equations and conduct numerical analyses of blood pressure distribution. This computational methodology provides valuable insights into interstitial fluid patterns and the spatiotemporal distribution of extracellular and intracellular concentrations of 225Ac-PSMA and 177Lu-PSMA. In addition, our study investigates the impacts of increasing tumor size on absorbed doses and mechanisms involved in radiopharmaceutical transport and delivery.

RESULTS

Larger tumors have diminished absorbed doses, highlighting the need for customized treatments according to tumor size. Diffusion significantly influences the transport and delivery of radiopharmaceuticals. Additionally, alpha therapy was observed to consistently yield higher absorbed doses within the tumor than beta therapy.

CONCLUSIONS

This study reveals the complex interplay between radiopharmaceutical properties, the tumor microenvironment, and treatment outcomes. It highlights the potential of 225Ac-PSMA in prostate cancer treatment, advocating for personalized treatment strategies tailored to the specific characteristics of each patient and their tumor.

A Futile Cycle?: Tissue Homeostatic Trans-Membrane Water Co-Transport: Kinetics, Thermodynamics, Metabolic Consequences

The phenomenon of active trans-membrane water cycling (AWC) has emerged in little over a decade. Here, we consider H2O transport across cell membranes from the origins of its study. Historically, trans-membrane water transport processes were classified into: A) compensating bidirectional fluxes ("exchange"), and B) unidirectional flux ("net flow") categories. Recent literature molecular structure determinations and molecular dynamic (MD) simulations indicate probably all the many different hydrophilic substrate membrane co-transporters have membrane-spanning hydrophilic pathways and co-transport water along with their substrates, and that they individually catalyze category A and/or B water flux processes, although usually not simultaneously. The AWC name signifies that, integrated over the all the cell's co-transporters, the rate of homeostatic, bidirectional trans-cytolemmal water exchange (category A) is synchronized with the metabolic rate of the crucial Na+,K+-ATPase (NKA) enzyme. A literature survey indicates the stoichiometric (category B) water/substrate ratios of individual co-transporters are often very large. The MD simulations also suggest how different co-transporter reactions can be kinetically coupled molecularly. Is this (Na+,K+-ATPase rate-synchronized) cycling futile, or is it consequential? Conservatively representative literature metabolomic and proteinomic results enable comprehensive free energy analyses of the many transport reactions with known water stoichiometries. Free energy calculations, using literature intracellular pressure (Pi) values reveals there is an outward trans-membrane H2O barochemical gradient of magnitude comparable to that of the well-known inward Na+ electrochemical gradient. For most co-influxers, these gradients are finely balanced to maintain intracellular metabolite concentration values near their consuming enzyme Michaelis constants. The thermodynamic analyses include glucose, glutamate-, gamma-aminobutyric acid (GABA), and lactate- transporters. 2%-4% Pi alterations can lead to disastrous concentration levels. For the neurotransmitters glutamate- and GABA, very small astrocytic Pi changes can allow/disallow synaptic transmission. Unlike the Na+ and K+ electrochemical steady-states, the H2O barochemical steady-state is in (or near) chemical equilibrium. The analyses show why the presence of aquaporins (AQPs) does not dissipate the trans-membrane pressure gradient. A feedback loop inherent in the opposing Na+ electrochemical and H2O barochemical gradients regulates AQP-catalyzed water flux as an integral AWC aspect. These results also require a re-consideration of the underlying nature of Pi. Active trans-membrane water cycling is not futile, but is inherent to the cell's "NKA system" - a new, fundamental aspect of biology.

Hydraulic conductivity of human cancer tissue: A hybrid study

Background

Elevated tumor tissue interstitial fluid pressure (IFP) is an adverse biomechanical biomarker that predicts poor therapy response and an aggressive phenotype. Advances in functional imaging have opened the prospect of measuring IFP non-invasively. Image-based estimation of the IFP requires knowledge of the tissue hydraulic conductivity (K), a measure for the ease of bulk flow through the interstitium. However, data on the magnitude of K in human cancer tissue are not available.

Methods

We measured the hydraulic conductivity of tumor tissue using modified Ussing chambers in surgical resection specimens. The effect of the tumor microenvironment (TME) on K was investigated by quantifying the collagen content, cell density, and fibroblast density of the tested samples using quantitative immune histochemistry. Also, we developed a computational fluid dynamics (CFD) model to evaluate the role of K on interstitial fluid flow and drug transport in solid tumors.

Results

The results show that the hydraulic conductivity of human tumor tissues is very limited, ranging from approximately 10-15 to 10-14 m2/Pa∙s. Moreover, K values varied significantly between tumor types and between different samples from the same tumor. A significant inverse correlation was found between collagen fiber density and hydraulic conductivity values. However, no correlation was detected between K and cancer cell or fibroblast densities. The computational model demonstrated the impact of K on the interstitial fluid flow and the drug concentration profile: higher K values led to a lower IFP and deeper drug penetration.

Conclusions

Human tumor tissue is characterized by a very limited hydraulic conductivity, representing a barrier to effective drug transport. The results of this study can inform the development of realistic computational models, facilitate non-invasive IFP estimation, and contribute to stromal targeting anticancer therapies.

Characterization of Microvascular Invasion in Hepatocellular Carcinoma Using Computational Modeling of Interstitial Fluid Pressure and Velocity

BACKGROUND

Most solid tumors show increased interstitial fluid pressure (IFP), and this increased IFP is an obstacle to treatment. A noninvasive model for measuring IFP in hepatocellular carcinoma (HCC) is an unresolved issue.

PURPOSE

To develop a noninvasive model to measure IFP and interstitial fluid velocity (IFV) in HCC and to characterize the microvascular invasion (MVI) status by using this model.

STUDY TYPE

Retrospective.

POPULATION

A total of 97 HCC patients (mean age 57.6 ± 10.9 years, 77.3% males), 53 of them with MVI and 44 of them without MVI.

FIELD STRENGTH/SEQUENCE

A 3-T, three-dimensional spoiled gradient-recalled echo.

ASSESSMENT

MVI was defined as microscopic vascular invasion of small vessels within the peritumoral liver tissue. The volumes of interest (VOIs) were manually delineated and enclosed the tumor lesion and healthy liver parenchyma, respectively. The extended Tofts model (ETM) was used to estimate permeability parameters from all the VOIs. Subsequently, the continuity partial differential equation (PDE) was implemented and IFP and IFV were acquired.

STATISTICAL TESTS

Wilcoxon signed-ranks tests, histogram analysis, Mann-Whitney U test, Fisher's exact test, least absolute shrinkage and selection operator (LASSO) logistic regression, receiver operating characteristic (ROC) curve analysis with the area under the curve (AUC), Youden index, DeLong test, and Benjamini-Hochberg correction. A P value <0.05 was considered statistically significant.

RESULTS

The HCC lesions exhibited elevated IFP and reduced IFV. There were no significant differences in any measured demographic and clinical features between the MVI-positive and MVI-negative groups, except for tumor size. Nine IFP histogram analysis-derived parameters and seven IFV histogram analysis-derived parameters could be used to characterize the MVI status. LASSO regression selected five features: IFP maximum, IFP 10th percentile, IFP 90th percentile, IFV SD, and IFV 10th percentile. The combination of these features showed the highest AUC (0.781) and specificity (77.3%).

DATA CONCLUSION

A noninvasive IFP and IFV measurement model for HCC was developed. Specific IFP- and IFV-derived parameters exhibited significant association with the MVI status.

EVIDENCE LEVEL

3.

TECHNICAL EFFICACY

Stage 2.

Drug transport modeling in solid tumors: A computational exploration of spatial heterogeneity of biophysical properties

Inadequate uptake of therapeutic agents by tumor cells is still a major barrier in clinical cancer therapy. Mathematical modeling is a powerful tool to describe and investigate the transport phenomena involved. However, current models for interstitial flow and drug delivery in solid tumors have not yet embedded the existing heterogeneity of tumor biomechanical properties. The purpose of this study is to introduce a novel and more realistic methodology for computational models of solid tumor perfusion and drug delivery accounting for these regional heterogeneities as well as lymphatic drainage effects. Several tumor geometries were studied using an advanced computational fluid dynamics (CFD) modeling approach of intratumor interstitial fluid flow and drug transport. Hereby, the following novelties were implemented: (i) the heterogeneity of tumor-specific hydraulic conductivity and capillary permeability; (ii) the effect of lymphatic drainage on interstitial fluid flow and drug penetration. Tumor size and shape both have a crucial role on the interstitial fluid flow regime as well as drug transport illustrating a direct correlation with interstitial fluid pressure (IFP) and an inverse correlation with drug penetration, except for large tumors having a diameter larger than 50 mm. The results also suggest that the interstitial fluid flow and drug penetration in small tumors depend on tumor shape. A parameter study on the necrotic core size illustrated that the core effect (i.e. fluid flow and drug penetration alteration) was only profound in small tumors. Interestingly, the impact of a necrotic core on drug penetration differs depending on the tumor shape from having no effect in ideally spherical tumors to a clear effect in elliptical tumors with a necrotic core. A realistic presence of lymphatic vessels only slightly affected tumor perfusion, having no substantial effect on drug delivery. In conclusion, our findings illustrated that our novel parametric CFD modeling strategy in combination with accurate profiling of heterogeneous tumor biophysical properties can provide a powerful tool for better insights into tumor perfusion and drug transport, enabling effective therapy planning.

Interstitial fluid pressure as an emerging biomarker in solid tumors

The physical microenvironment of cancer is characterized by elevated stiffness and tissue pressure, the main component of which is the interstitial fluid pressure (IFP). Elevated IFP is an established negative predictive and prognostic parameter, directly affecting malignant behavior and therapy response. As such, measurement of the IFP would allow to develop strategies aimed at engineering the physical microenvironment of cancer. Traditionally, IFP measurement required the use of invasive methods. Recent progress in dynamic and functional imaging methods such as dynamic contrast enhanced (DCE) magnetic resonance imaging and elastography, combined with numerical models and simulation, allows to comprehensively assess the biomechanical landscape of cancer, and may help to overcome physical barriers to drug delivery and immune cell infiltration. Here, we provide a comprehensive overview of the origin of elevated IFP, and its role in the malignant phenotype. Also, we review the methods used to measure IFP using invasive and imaging based methods, and highlight remaining obstacles and potential areas of progress in order to implement IFP measurement in clinical practice.

Interaction between tumor microenvironment, autophagy, and epithelial-mesenchymal transition in tumor progression

Tumor microenvironment (TME) is the ecosystem surrounding a tumor to influence tumor cells' growth, metastasis and immunological battlefield, in which the tumor systems fight against the body system. TME has been considered as the essential link between the tumorigenesis and development of neoplasm. Both nutrients intake and tumor progression to malignancy require the participation of components in TME. Epithelial-mesenchymal transition (EMT) is a key step in the metastasis of tumor cells. Cells that lost polarity and acquired migration ability are prone to metastasize. Autophagy is an important self-protective mechanism in tumor cells and a necessity for the tumor cells to respond to harmful stress. Protective autophagy benefits tumor cells while abnormal autophagy leads to cell injury or death. EMT and autophagy are directly regulated by TME. To date, there are numerous studies on TME, autophagy and EMT separately, but few on their complex interrelationships. This review aims to comprehensively analyze the existing mechanisms and convincing evidence so far to seek novel therapeutic strategies and research directions.

The Role of in silico Research in Developing Nanoparticle-Based Therapeutics

Nanoparticles (NPs) hold great potential as therapeutics, particularly in the realm of drug delivery. They are effective at functional cargo delivery and offer a great degree of amenability that can be used to offset toxic side effects or to target drugs to specific regions in the body. However, there are many challenges associated with the development of NP-based drug formulations that hamper their successful clinical translation. Arguably, the most significant barrier in the way of efficacious NP-based drug delivery systems is the tedious and time-consuming nature of NP formulation—a process that needs to account for downstream effects, such as the onset of potential toxicity or immunogenicity, in vivo biodistribution and overall pharmacokinetic profiles, all while maintaining desirable therapeutic outcomes. Computational and AI-based approaches have shown promise in alleviating some of these restrictions. Via predictive modeling and deep learning, in silico approaches have shown the ability to accurately model NP-membrane interactions and cellular uptake based on minimal data, such as the physicochemical characteristics of a given NP. More importantly, machine learning allows computational models to predict how specific changes could be made to the physicochemical characteristics of a NP to improve functional aspects, such as drug retention or endocytosis. On a larger scale, they are also able to predict the in vivo pharmacokinetics of NP-encapsulated drugs, predicting aspects such as circulatory half-life, toxicity, and biodistribution. However, the convergence of nanomedicine and computational approaches is still in its infancy and limited in its applicability. The interactions between NPs, the encapsulated drug and the body form an intricate network of interactions that cannot be modeled with absolute certainty. Despite this, rapid advancements in the area promise to deliver increasingly powerful tools capable of accelerating the development of advanced nanoscale therapeutics. Here, we describe computational approaches that have been utilized in the field of nanomedicine, focusing on approaches for NP design and engineering.

Longitudinal Monitoring of Simulated Interstitial Fluid Pressure for Pancreatic Ductal Adenocarcinoma Patients Treated with Stereotactic Body Radiotherapy

Simple Summary

High vessel permeability, poor perfusion, low lymphatic drainage, and dense abundant stroma elevate interstitial fluid pressures (IFP) in pancreatic ductal adenocarcinoma (PDAC). The present study aims to monitor longitudinal changes in simulated tumor IFP and velocity (IFV) values using a dynamic contrast-enhanced (DCE)-MRI-based computational fluid modeling (CFM) approach in PDAC. Nine PDAC patients underwent DCE-MRI acquisition on a 3-Tesla MRI scanner at pre-treatment (TX (0)), immediately after the first fraction of stereotactic body radiotherapy (SBRT, (D1-TX)), and six weeks post-TX (D2-TX). The partial differential equation of IFP formulated from the continuity equation using the Starling Principle of fluid exchange and Darcy velocity–pressure relationship was solved in COMSOL Multiphysics software to generate IFP and IFV parametric maps using relevant tumor tissue physiological parameters. Initial results suggest that after validation, IFP and IFV can be imaging biomarkers of early response to therapy that may guide precision medicine in PDAC.

Abstract

The present study aims to monitor longitudinal changes in simulated tumor interstitial fluid pressure (IFP) and velocity (IFV) values using dynamic contrast-enhanced (DCE)-MRI-based computational fluid modeling (CFM) in pancreatic ductal adenocarcinoma (PDAC) patients. Nine PDAC patients underwent MRI, including DCE-MRI, on a 3-Tesla MRI scanner at pre-treatment (TX (0)), after the first fraction of stereotactic body radiotherapy (SBRT, (D1-TX)), and six weeks post-TX (D2-TX). The partial differential equation of IFP formulated from the continuity equation, incorporating the Starling Principle of fluid exchange, Darcy velocity, and volume transfer constant (K^trans), was solved in COMSOL Multiphysics software to generate IFP and IFV maps. Tumor volume (V_t), K^trans, IFP, and IFV values were compared (Wilcoxon and Spearman) between the time- points. D2-TX K^trans values were significantly different from pre-TX and D1-TX (p < 0.05). The D1-TX and pre-TX mean IFV values exhibited a borderline significant difference (p = 0.08). The IFP values varying <3.0% between the three time-points were not significantly different (p > 0.05). V_t and IFP values were strongly positively correlated at pre-TX (ρ = 0.90, p = 0.005), while IFV exhibited a strong negative correlation at D1-TX (ρ = −0.74, p = 0.045). V_t, K^trans, IFP, and IFV hold promise as imaging biomarkers of early response to therapy in PDAC.

Application of Community Detection Algorithm to Investigate the Correlation between Imaging Biomarkers of Tumor Metabolism, Hypoxia, Cellularity, and Perfusion for Precision Radiotherapy in Head and Neck Squamous Cell Carcinomas

The present study aimed to investigate the correlation at pre-treatment (TX) between quantitative metrics derived from multimodality imaging (MMI), including 18F-FDG-PET/CT, 18F-FMISO-PET/CT, DW- and DCE-MRI, using a community detection algorithm (CDA) in head and neck squamous cell carcinoma (HNSCC) patients. Twenty-three HNSCC patients with 27 metastatic lymph nodes underwent a total of 69 MMI exams at pre-TX. Correlations among quantitative metrics derived from FDG-PET/CT (SUL), FMSIO-PET/CT (K1, k3, TBR, and DV), DW-MRI (ADC, IVIM [D, D*, and f]), and FXR DCE-MRI [Ktrans, ve, and τi]) were investigated using the CDA based on a "spin-glass model" coupled with the Spearman's rank, ρ, analysis. Mean MRI T2 weighted tumor volumes and SULmean values were moderately positively correlated (ρ = 0.48, p = 0.01). ADC and D exhibited a moderate negative correlation with SULmean (ρ ≤ -0.42, p < 0.03 for both). K1 and Ktrans were positively correlated (ρ = 0.48, p = 0.01). In contrast, Ktrans and k3max were negatively correlated (ρ = -0.41, p = 0.03). CDA revealed four communities for 16 metrics interconnected with 33 edges in the network. DV, Ktrans, and K1 had 8, 7, and 6 edges in the network, respectively. After validation in a larger population, the CDA approach may aid in identifying useful biomarkers for developing individual patient care in HNSCC.

Nongaussian Intravoxel Incoherent Motion Diffusion Weighted and Fast Exchange Regime Dynamic Contrast-Enhanced-MRI of Nasopharyngeal Carcinoma: Preliminary Study for Predicting Locoregional Failure

The aim of the present study was to identify whether the quantitative metrics from pre-treatment (TX) non-Gaussian intravoxel incoherent motion (NGIVIM) diffusion weighted (DW-) and fast exchange regime (FXR) dynamic contrast enhanced (DCE)-MRI can predict patients with locoregional failure (LRF) in nasopharyngeal carcinoma (NPC). Twenty-nine NPC patients underwent pre-TX DW- and DCE-MRI on a 3T MR scanner. DW imaging data from primary tumors were fitted to monoexponential (ADC) and NGIVIM (D, D*, f, and K) models. The metrics Ktrans, ve, and τi were estimated using the FXR model. Cumulative incidence (CI) analysis and Fine-Gray (FG) modeling were performed considering death as a competing risk. Mean ve values were significantly different between patients with and without LRF (p = 0.03). Mean f values showed a trend towards the difference between the groups (p = 0.08). Histograms exhibited inter primary tumor heterogeneity. The CI curves showed significant differences for the dichotomized cutoff value of ADC ≤ 0.68 × 10-3 (mm2/s), D ≤ 0.74 × 10-3 (mm2/s), and f ≤ 0.18 (p < 0.05). τi ≤ 0.89 (s) cutoff value showed borderline significance (p = 0.098). FG's modeling showed a significant difference for the K cutoff value of ≤0.86 (p = 0.034). Results suggest that the role of pre-TX NGIVIM DW- and FXR DCE-MRI-derived metrics for predicting LRF in NPC than alone.

Integrating Quantitative Assays with Biologically Based Mathematical Modeling for Predictive Oncology

We provide an overview on the use of biological assays to calibrate and initialize mechanism-based models of cancer phenomena. Although artificial intelligence methods currently dominate the landscape in computational oncology, mathematical models that seek to explicitly incorporate biological mechanisms into their formalism are of increasing interest. These models can guide experimental design and provide insights into the underlying mechanisms of cancer progression. Historically, these models have included a myriad of parameters that have been difficult to quantify in biologically relevant systems, limiting their practical insights. Recently, however, there has been much interest calibrating biologically based models with the quantitative measurements available from (for example) RNA sequencing, time-resolved microscopy, and in vivo imaging. In this contribution, we summarize how a variety of experimental methods quantify tumor characteristics from the molecular to tissue scales and describe how such data can be directly integrated with mechanism-based models to improve predictions of tumor growth and treatment response.

Quantitative Imaging Enters the Clinical Arena: A Personal Viewpoint

The National Cancer Institute's Quantitative Imaging Network (QIN) has thrived over the past 12 years with an emphasis on the development of image-based decision support software tools for improving measurements of imaging metrics. An overarching goal has been to develop advanced tools that could be translated into clinical trials to provide for improved prediction of response to therapeutic interventions. This article provides an overview of the successes in development and translation of new algorithms into the clinical workflow by the many research teams of the Quantitative Imaging Network.