Evaluating the influence of mechanical stress on anticancer treatments through a multiphase porous media model

Toward innovative approaches for exploring the mechanically regulated tumor-immune microenvironment

Within the complex tumor microenvironment, cells experience mechanical cues-such as extracellular matrix stiffening and elevation of solid stress, interstitial fluid pressure, and fluid shear stress-that significantly impact cancer cell behavior and immune responses. Recognizing the significance of these mechanical cues not only sheds light on cancer progression but also holds promise for identifying potential biomarkers that would predict therapeutic outcomes. However, standardizing methods for studying how mechanical cues affect tumor progression is challenging. This challenge stems from the limitations of traditional in vitro cell culture systems, which fail to encompass the critical contextual cues present in vivo. To address this, 3D tumor spheroids have been established as a preferred model, more closely mimicking cancer progression, but they usually lack reproduction of the mechanical microenvironment encountered in actual solid tumors. Here, we review the role of mechanical forces in modulating tumor- and immune-cell responses and discuss how grasping the importance of these mechanical cues could revolutionize in vitro tumor tissue engineering. The creation of more physiologically relevant environments that better replicate in vivo conditions will eventually increase the efficacy of currently available treatments, including immunotherapies.

Bio-Mechanical Model of Osteosarcoma Tumor Microenvironment: A Porous Media Approach

Osteosarcoma is the most common malignant bone tumor in children and adolescents with a poor prognosis. To describe the progression of osteosarcoma, we expanded a system of data-driven ODE from a previous study into a system of Reaction-Diffusion-Advection (RDA) equations and coupled it with Biot equations of poroelasticity to form a bio-mechanical model. The RDA system includes the spatio-temporal information of the key components of the tumor microenvironment. The Biot equations are comprised of an equation for the solid phase, which governs the movement of the solid tumor, and an equation for the fluid phase, which relates to the motion of cells. The model predicts the total number of cells and cytokines of the tumor microenvironment and simulates the tumor's size growth. We simulated different scenarios using this model to investigate the impact of several biomedical settings on tumors' growth. The results indicate the importance of macrophages in tumors' growth. Particularly, we have observed a high co-localization of macrophages and cancer cells, and the concentration of tumor cells increases as the number of macrophages increases.

Mathematical characterization of population dynamics in breast cancer cells treated with doxorubicin

The development of chemoresistance remains a significant cause of treatment failure in breast cancer. We posit that a mathematical understanding of chemoresistance could assist in developing successful treatment strategies. Towards that end, we have developed a model that describes the cytotoxic effects of the standard chemotherapeutic drug doxorubicin on the MCF-7 breast cancer cell line. We assume that treatment with doxorubicin induces a compartmentalization of the breast cancer cell population into surviving cells, which continue proliferating after treatment, and irreversibly damaged cells, which gradually transition from proliferating to treatment-induced death. The model is fit to experimental data including variations in drug concentration, inter-treatment interval, and number of doses. Our model recapitulates tumor cell dynamics in all these scenarios (as quantified by the concordance correlation coefficient, CCC > 0.95). In particular, superior tumor control is observed with higher doxorubicin concentrations, shorter inter-treatment intervals, and a higher number of doses (p < 0.05). Longer inter-treatment intervals require adapting the model parameterization after each doxorubicin dose, suggesting the promotion of chemoresistance. Additionally, we propose promising empirical formulas to describe the variation of model parameters as functions of doxorubicin concentration (CCC > 0.78). Thus, we conclude that our mathematical model could deepen our understanding of the cytotoxic effects of doxorubicin and could be used to explore practical drug regimens achieving optimal tumor control.

A continuum mechanical framework for modeling tumor growth and treatment in two- and three-phase systems

The growth and treatment of tumors is an important problem to society that involves the manifestation of cellular phenomena at length scales on the order of centimeters. Continuum mechanical approaches are being increasingly used to model tumors at the largest length scales of concern. The issue of how to best connect such descriptions to smaller-scale descriptions remains open. We formulate a framework to derive macroscale models of tumor behavior using the thermodynamically constrained averaging theory (TCAT), which provides a firm connection with the microscale and constraints on permissible forms of closure relations. We build on developments in the porous medium mechanics literature to formulate fundamental entropy inequality expressions for a general class of three-phase, compositional models at the macroscale. We use the general framework derived to formulate two classes of models, a two-phase model and a three-phase model. The general TCAT framework derived forms the basis for a wide range of potential models of varying sophistication, which can be derived, approximated, and applied to understand not only tumor growth but also the effectiveness of various treatment modalities.

Combined Influence of Nutrient Supply Level and Tissue Mechanical Properties on Benign Tumor Growth as Revealed by Mathematical Modeling

A continuous mathematical model of non-invasive avascular tumor growth in tissue is presented. The model considers tissue as a biphasic material, comprised of a solid matrix and interstitial fluid. The convective motion of tissue elements happens due to the gradients of stress, which change as a result of tumor cells proliferation and death. The model accounts for glucose as the crucial nutrient, supplied from the normal tissue, and can reproduce both diffusion-limited and stress-limited tumor growth. Approximate tumor growth curves are obtained semi-analytically in the limit of infinite tissue hydraulic conductivity, which implies instantaneous equalization of arising stress gradients. These growth curves correspond well to the numerical solutions and represent classical sigmoidal curves with a short initial exponential phase, subsequent almost linear growth phase and a phase with growth deceleration, in which tumor tends to reach its maximum volume. The influence of two model parameters on tumor growth curves is investigated: tissue hydraulic conductivity, which links the values of stress gradient and convective velocity of tissue phases, and tumor nutrient supply level, which corresponds to different permeability and surface area density of capillaries in the normal tissue that surrounds the tumor. In particular, it is demonstrated, that sufficiently low tissue hydraulic conductivity (intrinsic, e.g., to tumors arising from connective tissue) and sufficiently high nutrient supply can lead to formation of giant benign tumors, reaching tens of centimeters in diameter, which are indeed observed clinically.

In silico investigations of intratumoral heterogeneous interstitial fluid pressure

Recent preclinical studies have shown that interstitial fluid pressure (IFP) within tumors can be heterogeneous Andersen et al. (2019). In that study tumors of two xenograft models, respectively, HL-16 cervical carcinoma and Panc-1 pancreatic carcinoma, were investigated. Significant heterogeneity in IFP was reported and it was proposed that this was associated with division of tissue into compartments separated by thick connective tissue bands for the HL-16 tumors and with dense collagen-rich extracellular matrix for the Panc-1 tumors. The purpose of the current work is to explore these experimental observations by using in silico generated tumor models. We consider a mathematical multiphase model which accounts for tumor cells, fibroblasts and interstitial fluid. The model has been trained to comply with experimental in vitro results reported in Shieh et al. (2011) which has identified autologous chemotaxis, ECM remodeling, and cell-fibroblast interaction as drivers for invasive tumor cell behavior. The in silico model is informed with parameters that characterize the leaky intratumoral vascular network, the peritumoral lymphatics which collect the fluid, and the density of ECM as represented through the hydraulic conductivity of the interstitial space. Heterogeneous distribution of solid stress may result in heterogeneous compression of blood vessels and, thus, heterogeneous vascular density inside the tumor. To mimic this we expose the in silico tumor to an intratumoral vasculature whose net effect of density of blood vesssels and vessel wall conductivity is varied through a 2D Gaussian variogram constrained such that the resulting IFPs lie within the range as reported from the preclinical study. The in silico cervical carcinoma model illustrates that sparse ECM was associated with uniform intratumoral IFP in spite of heterogeneous microvascular network, whereas compartment structures resulted in more heterogeneous IFP. Similarly, the in silico pancreatic model shows that heterogeneity in the microvascular network combined with dense ECM structure prevents IFP to even out and gives rise to heterogeneous IFP. The computer model illustrates how a heterogeneous invasive front might form where groups of tumor cells detach from the primary tumor and form isolated islands, a behavior which is natural to associate with metastatic propensity. However, unlike experimental studies, the current version of the in silico model does not show an association between metastatic propensity and elevated IFP.

Mechanical Control of Cell Proliferation Increases Resistance to Chemotherapeutic Agents

While many cellular mechanisms leading to chemotherapeutic resistance have been identified, there is an increasing realization that tumor-stroma interactions also play an important role. In particular, mechanical alterations are inherent to solid cancer progression and profoundly impact cell physiology. Here, we explore the influence of compressive stress on the efficacy of chemotherapeutics in pancreatic cancer spheroids. We find that increased compressive stress leads to decreased drug efficacy. Theoretical modeling and experiments suggest that mechanical stress decreases cell proliferation which in turn reduces the efficacy of chemotherapeutics that target proliferating cells. Our work highlights a mechanical form of drug resistance and suggests new strategies for therapy.

Drug delivery: Experiments, mathematical modelling and machine learning

We address the problem of determining from laboratory experiments the data necessary for a proper modeling of drug delivery and efficacy in anticancer therapy. There is an inherent difficulty in extracting the necessary parameters, because the experiments often yield an insufficient quantity of information. To overcome this difficulty, we propose to combine real experiments, numerical simulation, and Machine Learning (ML) based on Artificial Neural Networks (ANN), aiming at a reliable identification of the physical model factors, e.g. the killing action of the drug. To this purpose, we exploit the employed mathematical-numerical model for tumor growth and drug delivery, together with the ANN - ML procedure, to integrate the results of the experimental tests and feed back the model itself, thus obtaining a reliable predictive tool. The procedure represents a hybrid data-driven, physics-informed approach to machine learning. The physical and mathematical model employed for the numerical simulations is without extracellular matrix (ECM) and healthy cells because of the experimental conditions we reproduce.

Investigating the Physical Effects in Bacterial Therapies for Avascular Tumors

Tumor-targeting bacteria elicit anticancer effects by infiltrating hypoxic regions, releasing toxic agents and inducing immune responses. Although current research has largely focused on the influence of chemical and immunological aspects on the mechanisms of bacterial therapy, the impact of physical effects is still elusive. Here, we propose a mathematical model for the anti-tumor activity of bacteria in avascular tumors that takes into account the relevant chemo-mechanical effects. We consider a time-dependent administration of bacteria and analyze the impact of bacterial chemotaxis and killing rate. We show that active bacterial migration toward tumor hypoxic regions provides optimal infiltration and that high killing rates combined with high chemotactic values provide the smallest tumor volumes at the end of the treatment. We highlight the emergence of steady states in which a small population of bacteria is able to constrain tumor growth. Finally, we show that bacteria treatment works best in the case of tumors with high cellular proliferation and low oxygen consumption.

Extension of a multiphase tumour growth model to study nanoparticle delivery to solid tumours

One of the main challenges in increasing the efficacy of conventional chemotherapeutics is the fact that they do not reach cancerous cells at a sufficiently high dosage. In order to remedy this deficiency, nanoparticle-based drugs have evolved as a promising novel approach to more specific tumour targeting. Nevertheless, several biophysical phenomena prevent the sufficient penetration of nanoparticles in order to target the entire tumour. We therefore extend our vascular multiphase tumour growth model, enabling it to investigate the influence of different biophysical factors on the distribution of nanoparticles in the tumour microenvironment. The novel model permits the examination of the interplay between the size of vessel-wall pores, the permeability of the blood-vessel endothelium and the lymphatic drainage on the delivery of particles of different sizes. Solid tumours develop a non-perfused core and increased interstitial pressure. Our model confirms that those two typical features of solid tumours limit nanoparticle delivery. Only in case of small nanoparticles is the transport dominated by diffusion, and particles can reach the entire tumour. The size of the vessel-wall pores and the permeability of the blood-vessel endothelium have a major impact on the amount of delivered nanoparticles. This extended in-silico tumour growth model permits the examination of the characteristics and of the limitations of nanoparticle delivery to solid tumours, which currently complicate the translation of nanoparticle therapy to a clinical stage.

On the Impact of Chemo-Mechanically Induced Phenotypic Transitions in Gliomas

Tumor microenvironment is a critical player in glioma progression, and novel therapies for its targeting have been recently proposed. In particular, stress-alleviation strategies act on the tumor by reducing its stiffness, decreasing solid stresses and improving blood perfusion. However, these microenvironmental changes trigger chemo-mechanically induced cellular phenotypic transitions whose impact on therapy outcomes is not completely understood. In this work we analyze the effects of mechanical compression on migration and proliferation of glioma cells. We derive a mathematical model of glioma progression focusing on cellular phenotypic plasticity. Our results reveal a trade-off between tumor infiltration and cellular content as a consequence of stress-alleviation approaches. We discuss how these novel findings increase the current understanding of glioma/microenvironment interactions and can contribute to new strategies for improved therapeutic outcomes.

Feeling Stress: The Mechanics of Cancer Progression and Aggression

The tumor microenvironment is a dynamic landscape in which the physical and mechanical properties evolve dramatically throughout cancer progression. These changes are driven by enhanced tumor cell contractility and expansion of the growing tumor mass, as well as through alterations to the material properties of the surrounding extracellular matrix (ECM). Consequently, tumor cells are exposed to a number of different mechanical inputs including cell–cell and cell-ECM tension, compression stress, interstitial fluid pressure and shear stress. Oncogenes engage signaling pathways that are activated in response to mechanical stress, thereby reworking the cell's intrinsic response to exogenous mechanical stimuli, enhancing intracellular tension via elevated actomyosin contraction, and influencing ECM stiffness and tissue morphology. In addition to altering their intracellular tension and remodeling the microenvironment, cells actively respond to these mechanical perturbations phenotypically through modification of gene expression. Herein, we present a description of the physical changes that promote tumor progression and aggression, discuss their interrelationship and highlight emerging therapeutic strategies to alleviate the mechanical stresses driving cancer to malignancy.

Modeling mechanical inhomogeneities in small populations of proliferating monolayers and spheroids

Understanding the mechanical behavior of multicellular monolayers and spheroids is fundamental to tissue culture, organism development, and the early stages of tumor growth. Proliferating cells in monolayers and spheroids experience mechanical forces as they grow and divide and local inhomogeneities in the mechanical microenvironment can cause individual cells within the multicellular system to grow and divide at different rates. This differential growth, combined with cell division and reorganization, leads to residual stress. Multiple different modeling approaches have been taken to understand and predict the residual stresses that arise in growing multicellular systems, particularly tumor spheroids. Here, we show that by using a mechanically robust agent-based model constructed with the peridynamic framework, we gain a better understanding of residual stresses in multicellular systems as they grow from a single cell. In particular, we focus on small populations of cells (1-100 s) where population behavior is highly stochastic and prior investigation has been limited. We compare the average strain energy density of cells in monolayers and spheroids using different growth and division rules and find that, on average, cells in spheroids have a higher strain energy density than cells in monolayers. We also find that cells in the interior of a growing spheroid are, on average, in compression. Finally, we demonstrate the importance of accounting for stochastic fluctuations in the mechanical environment, particularly when the cellular response to mechanical cues is nonlinear. The results presented here serve as a starting point for both further investigation with agent-based models, and for the incorporation of major findings from agent-based models into continuum scale models when explicit representation of individual cells is not computationally feasible.