Cell adhesion mechanisms and stress relaxation in the mechanics of tumours

A morpho-viscoelasticity theory for growth in proliferating aggregates

Despite significant research efforts in the continuum modeling of biological growth, certain aspects have been overlooked. For instance, numerous investigations have examined the influence of morphogenetic cell behaviors, like division and intercalation, on the mechanical response of passive (non-growing) tissues. Yet, their impact on active growth dynamics remains inadequately explored. A key reason for this inadequacy stems from challenges in the continuum treatment of cell-level processes. While some coarse-grained models have been proposed to address these shortcomings, a focus on cell division and cell expansion has been missing, rendering them unusable when it comes to modeling growth. Moreover, existing studies are limited to two-dimensional tissues and are yet to be formally extended to three-dimensional multicellular systems. To address these limitations, we here present a generalized multiscale model for three-dimensional aggregates that accounts for complex morphogenetic movements that include division, expansion, and intercalation. The proposed continuum theory thus allows for a comprehensive exploration into the growth and dissipation mechanics of proliferating aggregates, such as spheroids and organoids.

Selected aspects of avascular tumor growth reproduced by a hybrid model of cell dynamics and chemical kinetics

We here propose a hybrid computational framework to reproduce and analyze aspects of the avascular progression of a generic solid tumor. Our method first employs an individual-based approach to represent the population of tumor cells, which are distinguished in viable and necrotic agents. The active part of the disease is in turn differentiated according to a set of metabolic states. We then describe the spatio-temporal evolution of the concentration of oxygen and of tumor-secreted proteolytic enzymes using partial differential equations (PDEs). A differential equation finally governs the local degradation of the extracellular matrix (ECM) by the malignant mass. Numerical realizations of the model are run to reproduce tumor growth and invasion in a number scenarios that differ for cell properties (adhesiveness, duplication potential, proteolytic activity) and/or environmental conditions (level of tissue oxygenation and matrix density pattern). In particular, our simulations suggest that tumor aggressiveness, in terms of invasive depth and extension of necrotic tissue, can be reduced by (i) stable cell-cell contact interactions, (ii) poor tendency of malignant agents to chemotactically move upon oxygen gradients, and (iii) presence of an overdense matrix, if coupled by a disrupted proteolytic activity of the disease.

Minimal Morphoelastic Models of Solid Tumour Spheroids: A Tutorial

Tumour spheroids have been the focus of a variety of mathematical models, ranging from Greenspan's classical study of the 1970 s through to contemporary agent-based models. Of the many factors that regulate spheroid growth, mechanical effects are perhaps some of the least studied, both theoretically and experimentally, though experimental enquiry has established their significance to tumour growth dynamics. In this tutorial, we formulate a hierarchy of mathematical models of increasing complexity to explore the role of mechanics in spheroid growth, all the while seeking to retain desirable simplicity and analytical tractability. Beginning with the theory of morphoelasticity, which combines solid mechanics and growth, we successively refine our assumptions to develop a somewhat minimal model of mechanically regulated spheroid growth that is free from many unphysical and undesirable behaviours. In doing so, we will see how iterating upon simple models can provide rigorous guarantees of emergent behaviour, which are often precluded by existing, more complex modelling approaches. Perhaps surprisingly, we also demonstrate that the final model considered in this tutorial agrees favourably with classical experimental results, highlighting the potential for simple models to provide mechanistic insight whilst also serving as mathematical examples.

Mechanotransduction in tumor dynamics modeling

Mechanotherapy is a groundbreaking approach to impact carcinogenesis. Cells sense and respond to mechanical stimuli, translating them into biochemical signals in a process known as mechanotransduction. The impact of stress on tumor growth has been studied in the last three decades, and many papers highlight the role of mechanics as a critical self-inducer of tumor fate at the in vitro and in vivo biological levels. Meanwhile, mathematical models attempt to determine laws to reproduce tumor dynamics. This review discusses biological mechanotransduction mechanisms and mathematical-biomechanical models together. The aim is to provide a common framework for the different approaches that have emerged in the literature from the perspective of tumor avascularity and to provide insight into emerging mechanotherapies that have attracted interest in recent years.

Evaluation of innate and adaptive immune system interactions in the tumor microenvironment via a 3D continuum model

Immune cells in the tumor microenvironment (TME) are known to affect tumor growth, vascularization, and extracellular matrix (ECM) deposition. Marked interest in system-scale analysis of immune species interactions within the TME has encouraged progress in modeling tumor-immune interactions in silico. Due to the computational cost of simulating these intricate interactions, models have typically been constrained to representing a limited number of immune species. To expand the capability for system-scale analysis, this study develops a three-dimensional continuum mixture model of tumor-immune interactions to simulate multiple immune species in the TME. Building upon a recent distributed computing implementation that enables efficient solution of such mixture models, major immune species including monocytes, macrophages, natural killer cells, dendritic cells, neutrophils, myeloid-derived suppressor cells (MDSC), cytotoxic, helper, regulatory T-cells, and effector and regulatory B-cells and their interactions are represented in this novel implementation. Immune species extravasate from blood vasculature, undergo chemotaxis toward regions of high chemokine concentration, and influence the TME in proportion to locally defined levels of stimulation. The immune species contribute to the production of angiogenic and tumor growth factors, promotion of myofibroblast deposition of ECM, upregulation of angiogenesis, and elimination of living and dead tumor species. The results show that this modeling approach offers the capability for quantitative insight into the modulation of tumor growth by diverse immune-tumor interactions and immune-driven TME effects. In particular, MDSC-mediated effects on tumor-associated immune species' activation levels, volume fraction, and influence on the TME are explored. Longer term, linking of the model parameters to particular patient tumor information could simulate cancer-specific immune responses and move toward a more comprehensive evaluation of immunotherapeutic strategies.

Digital twin demonstrates significance of biomechanical growth control in liver regeneration after partial hepatectomy

Partial liver removal is an important therapy option for liver cancer. In most patients within a few weeks, the liver is able to fully regenerate. In some patients, however, regeneration fails with often severe consequences. To better understand the control mechanisms of liver regeneration, experiments in mice were performed, guiding the creation of a spatiotemporal 3D model of the regenerating liver. The model represents cells and blood vessels within an entire liver lobe, a macroscopic liver subunit. The model could reproduce the experimental data only if a biomechanical growth control (BGC)-mechanism, inhibiting cell cycle entrance at high compression, was taken into account and predicted that BGC may act as a short-range growth inhibitor minimizing the number of proliferating neighbor cells of a proliferating cell, generating a checkerboard-like proliferation pattern. Model-predicted cell proliferation patterns in pigs and mice were found experimentally. The results underpin the importance of biomechanical aspects in liver growth control.

Simulation of 3D centimeter-scale continuum tumor growth at sub-millimeter resolution via distributed computing

Simulation of cm-scale tumor growth has generally been constrained by the computational cost to numerically solve the associated equations, with models limited to representing mm-scale or smaller tumors. While the work has proven useful to the study of small tumors and micro-metastases, a biologically-relevant simulation of cm-scale masses as would be typically detected and treated in patients has remained an elusive goal. This study presents a distributed computing (parallelized) implementation of a mixture model of tumor growth to simulate 3D cm-scale vascularized tissue at sub-mm resolution. The numerical solving scheme utilizes a two-stage parallelization framework. The solution is written for GPU computation using the CUDA framework, which handles all Multigrid-related computations. Message Passing Interface (MPI) handles distribution of information across multiple processes, freeing the program from RAM and the processing limitations found on single systems. On each system, Nvidia's CUDA library allows for fast processing of model data using GPU-bound computing on fewer systems. The results show that a combined MPI-CUDA implementation enables the continuum modeling of cm-scale tumors at reasonable computational cost. Further work to calibrate model parameters to particular tumor conditions could enable simulation of patient-specific tumors for clinical application.

Impact of tumor-parenchyma biomechanics on liver metastatic progression: a multi-model approach

Colorectal cancer and other cancers often metastasize to the liver in later stages of the disease, contributing significantly to patient death. While the biomechanical properties of the liver parenchyma (normal liver tissue) are known to affect tumor cell behavior in primary and metastatic tumors, the role of these properties in driving or inhibiting metastatic inception remains poorly understood, as are the longer-term multicellular dynamics. This study adopts a multi-model approach to study the dynamics of tumor-parenchyma biomechanical interactions during metastatic seeding and growth. We employ a detailed poroviscoelastic model of a liver lobule to study how micrometastases disrupt flow and pressure on short time scales. Results from short-time simulations in detailed single hepatic lobules motivate constitutive relations and biological hypotheses for a minimal agent-based model of metastatic growth in centimeter-scale tissue over months-long time scales. After a parameter space investigation, we find that the balance of basic tumor-parenchyma biomechanical interactions on shorter time scales (adhesion, repulsion, and elastic tissue deformation over minutes) and longer time scales (plastic tissue relaxation over hours) can explain a broad range of behaviors of micrometastases, without the need for complex molecular-scale signaling. These interactions may arrest the growth of micrometastases in a dormant state and prevent newly arriving cancer cells from establishing successful metastatic foci. Moreover, the simulations indicate ways in which dormant tumors could "reawaken" after changes in parenchymal tissue mechanical properties, as may arise during aging or following acute liver illness or injury. We conclude that the proposed modeling approach yields insight into the role of tumor-parenchyma biomechanics in promoting liver metastatic growth, and advances the longer term goal of identifying conditions to clinically arrest and reverse the course of late-stage cancer.

A three dimensional model of multicellular aggregate compression

Multicellular aggregates are an excellent model system to explore the role of tissue biomechanics, which has been demonstrated to play a crucial role in many physiological and pathological processes. In this paper, we propose a three-dimensional mechanical model and apply it to the uniaxial compression of a multicellular aggregate in a realistic biological setting. In particular, we consider an aggregate of initially spherical shape and describe both its elastic deformations and the reorganisation of the cells forming the spheroid. The latter phenomenon, understood as remodelling, is accounted for by assuming that the aggregate undergoes plastic-like distortions. The study of the compression of the spheroid, achieved by means of two parallel, compressive plates, needs the formulation of a contact problem between the living spheroid itself and the plates, and is solved with the aid of the augmented Lagrangian method. The results of the performed numerical simulations are in qualitative agreement with the biological observations reported in the literature and can also be used to estimate quantitatively some fundamental aggregate mechanical parameters.

Multiscale modeling of solid stress and tumor cell invasion in response to dynamic mechanical microenvironment

Mathematical models can provide a quantitatively sophisticated description of tumor cell (TC) behaviors under mechanical microenvironment and help us better understand the role of specific biophysical factors based on their influences on the TC behaviors. To this end, we propose an off-lattice cell-based multiscale mathematical model to describe the dynamic growth-induced solid stress during tumor progression and investigate the influence of the mechanical microenvironment on TC invasion. At the cellular level, cell-cell and cell-matrix interactive forces depend on the mechanical properties of the cells and the cancer-associated fibroblasts in the stroma, respectively. The constitutive relationship between the interactive forces and cell migrations obeys the Hooke's law and damping effects. At the tissue level, the integrated growth-induced forces caused by proliferating cells within the simulation region are balanced by the external forces applied by the surrounding host tissues. Then, the cell movements are calculated according to the Newton's second law of motion, and the morphology of TC invasion is updated. The simulation results reveal the continuous changes of the macroscopic mechanical forces due to the interactions among the structural components and the microscopic environmental factors. Moreover, the simulation results demonstrate the adverse effect of the stiffness of tumor tissue on tumor growth and invasion. A decrease in the stiffness of tumor and matrix can promote TCs to proliferate at a much faster rate and invade into the surrounding healthy tissue more easily, whereas an increase in the stiffness can lead to an aggressive morphology of tumor invasion. We envision that the proposed model can be served as a quantitative theoretical platform to study the underlying biophysical role of the mechanical microenvironmental factors during tumor invasion and metastasis.

Kinetics of surface growth with coupled diffusion and the emergence of a universal growth path

Surface growth by association or dissociation of material on the boundary of a body is ubiquitous in both natural and engineering systems. It is the fundamental mechanism by which biological materials grow, starting from the level of a single cell, and is increasingly applied in engineering processes for fabrication and self-assembly. A significant challenge in modelling such processes arises due to the inherent coupled interaction between the growth kinetics, the local stresses and the diffusing constituents needed to sustain the growth. Moreover, the volume of the body changes not only due to surface growth but also by variation in solvent concentration within the bulk. In this paper, we present a general theoretical framework that captures these phenomena and describes the kinetics of surface growth while accounting for coupled diffusion. Then, by the combination of analytical and numerical tools, applied to a simple growth geometry, we show that the evolution of such growth processes tends towards a universal path that is independent of initial conditions. This path, on which surface growth and diffusion act harmoniously, can be extended to analytically portray the evolution of a body from inception up to a treadmilling state, in which addition and removal of material are balanced.

Simulation of Multispecies Desmoplastic Cancer Growth via a Fully Adaptive Non-linear Full Multigrid Algorithm

A fully adaptive non-linear full multigrid (FMG) algorithm is implemented to computationally simulate a model of multispecies desmoplastic tumor growth in three spatial dimensions. The algorithm solves a thermodynamic mixture model employing a diffuse interface approach with Cahn-Hilliard-type fourth-order equations that are coupled, non-linear, and numerically stiff. The tumor model includes extracellular matrix (ECM) as a major component with elastic energy contribution in its chemical potential term. Blood and lymphatic vasculatures are simulated via continuum representations. The model employs advection-reaction-diffusion partial differential equations (PDEs) for the cell, ECM, and vascular components, and reaction-diffusion PDEs for the elements diffusing from the vessels. This study provides the details of the numerical solution obtained by applying the fully adaptive non-linear FMG algorithm with finite difference method to solve this complex system of PDEs. The results indicate that this type of computational model can simulate the extracellular matrix-rich desmoplastic tumor microenvironment typical of fibrotic tumors, such as pancreatic adenocarcinoma.

Accumulation of mechanical forces in tumors is related to hyaluronan content and tissue stiffness

Hyaluronan is abundant in the extracellular matrix of many desmoplastic tumors and determines in large part the tumor biochemical and mechanical microenvironment. Additionally, it has been identified as one of the major physiological barriers to the effective delivery of drugs to solid tumors and its targeting with the use of pharmaceutical agents has shown to decompress tumor blood vessels, and thus improve tumor perfusion and efficacy of cytotoxic drugs. In this study, we investigated the contribution of hyaluronan to the accumulation of mechanical forces in tumors. Using experimental data from two orthotopic breast tumor models and treating tumors with two clinically approved anti-fibrotic drugs (tranilast and pirfenidone), we found that accumulation of growth-induced, residual forces in tumors are associated with hyaluronan content. Furthermore, mechanical characterization of the tumors revealed a good correlation of the accumulated forces with the elastic modulus of the tissue. Our results provide important insights on the mechano-pathology of solid tumors and can be used for the design of therapeutic strategies that target hyaluronan.

Predicting the role of microstructural and biomechanical cues in tumor growth and spreading

A continuum porous media model is developed for elucidating the role of the mechanical cues in regulating tumor growth and spreading. It is shown that stiffer matrices and higher cell-matrix adhesion limit tumor growth and spreading toward the surrounding tissue. Higher matrix porosities, conversely, favor the growth and the dissemination of tumor cells. This model could be used for predicting the response of malignant masses to novel therapeutic agents affecting directly the tumor microenvironment and its micromechanical cues.

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.

Solid tumors are poroelastic solids with a chemo-mechanical feedback on growth

The experimental evidence that a feedback exists between growth and stress in tumors poses challenging questions. First, the rheological properties (the "constitutive equations") of aggregates of malignant cells are still a matter of debate. Secondly, the feedback law (the "growth law") that relates stress and mitotic-apoptotic rate is far to be identified. We address these questions on the basis of a theoretical analysis of in vitro and in vivo experiments that involve the growth of tumor spheroids. We show that solid tumors exhibit several mechanical features of a poroelastic material, where the cellular component behaves like an elastic solid. When the solid component of the spheroid is loaded at the boundary, the cellular aggregate grows up to an asymptotic volume that depends on the exerted compression. Residual stress shows up when solid tumors are radially cut, highlighting a peculiar tensional pattern. By a novel numerical approach we correlate the measured opening angle and the underlying residual stress in a sphere. The features of the mechanobiological system can be explained in terms of a feedback of mechanics on the cell proliferation rate as modulated by the availability of nutrient, that is radially damped by the balance between diffusion and consumption. The volumetric growth profiles and the pattern of residual stress can be theoretically reproduced assuming a dependence of the target stress on the concentration of nutrient which is specific of the malignant tissue.

Model of vascular desmoplastic multispecies tumor growth

We present a three-dimensional nonlinear tumor growth model composed of heterogeneous cell types in a multicomponent-multispecies system, including viable, dead, healthy host, and extra-cellular matrix (ECM) tissue species. The model includes the capability for abnormal ECM dynamics noted in tumor development, as exemplified by pancreatic ductal adenocarcinoma, including dense desmoplasia typically characterized by a significant increase of interstitial connective tissue. An elastic energy is implemented to provide elasticity to the connective tissue. Cancer-associated fibroblasts (myofibroblasts) are modeled as key contributors to this ECM remodeling. The tumor growth is driven by growth factors released by these stromal cells as well as by oxygen and glucose provided by blood vasculature which along with lymphatics are stimulated to proliferate in and around the tumor based on pro-angiogenic factors released by hypoxic tissue regions. Cellular metabolic processes are simulated, including respiration and glycolysis with lactate fermentation. The bicarbonate buffering system is included for cellular pH regulation. This model system may be of use to simulate the complex interactions between tumor and stromal cells as well as the associated ECM and vascular remodeling that typically characterize malignant cancers notorious for poor therapeutic response.

Finite Element Modeling of Avascular Tumor Growth Using a Stress-Driven Model

Tumor growth being a multistage process has been investigated from different aspects. In the present study, an attempt is made to represent a constitutive-structure-based model of avascular tumor growth in which the effects of tensile stresses caused by collagen fibers are considered. Collagen fibers as a source of anisotropy in the structure of tissue are taken into account using a continuous fiber distribution formulation. To this end, a finite element modeling is implemented in which a neo-Hookean hyperelastic material is assigned to the tumor and its surrounding host. The tumor is supplied with a growth term. The growth term includes the effect of parameters such as nutrient concentration on the tumor growth and the tumor's solid phase content in the formulation. Results of the study revealed that decrease of solid phase is indicative of decrease in growth rate and the final steady-state value of tumor's radius. Moreover, fiber distribution affects the final shape of the tumor, and it could be used to control the shape and geometry of the tumor in complex morphologies. Finally, the findings demonstrated that the exerted stresses on the tumor increase as time passes. Compression of tumor cells leads to the reduction of tumor growth rate until it gradually reaches an equilibrium radius. This finding is in accordance with experimental data. Hence, this formulation can be deployed to evaluate both the residual stresses induced by growth and the mechanical interactions with the host tissue.

Quantifying the relationship between cell division angle and morphogenesis through computational modeling

When biological cells divide, they divide on a given angle. It has been shown experimentally that the orientation of cell division angle for a single cell can be described by a probability density function. However, the way in which the probability density function underlying cell division orientation influences population or tissue scale morphogenesis is unknown. Here we show that a computational approach, with thousands of stochastic simulations modeling growth and division of a population of cells, can be used to investigate this unknown. In this paper we examine two potential forms of the probability density function: a wrapped normal distribution and a binomial distribution. Our results demonstrate that for the wrapped normal distribution the standard deviation of the division angle, potentially interpreted as biological noise, controls the degree of tissue scale anisotropy. For the binomial distribution, we demonstrate a mechanism by which direction and degree of tissue scale anisotropy can be tuned via the probability of each division angle. We anticipate that the method presented in this paper and the results of these simulations will be a starting point for further investigation of this topic.

The Solid Mechanics of Cancer and Strategies for Improved Therapy

Tumor progression and response to treatment is determined in large part by the generation of mechanical stresses that stem from both the solid and the fluid phase of the tumor. Furthermore, elevated solid stress levels can regulate fluid stresses by compressing intratumoral blood and lymphatic vessels. Blood vessel compression reduces tumor perfusion, while compression of lymphatic vessels hinders the ability of the tumor to drain excessive fluid from its interstitial space contributing to the uniform elevation of the interstitial fluid pressure. Hypoperfusion and interstitial hypertension pose major barriers to the systemic administration of chemotherapeutic agents and nanomedicines to tumors, reducing treatment efficacies. Hypoperfusion can also create a hypoxic and acidic tumor microenvironment that promotes tumor progression and metastasis. Hence, alleviation of intratumoral solid stress levels can decompress tumor vessels and restore perfusion and interstitial fluid pressure. In this review, three major types of tissue level solid stresses involved in tumor growth, namely stress exerted externally on the tumor by the host tissue, swelling stress, and residual stress, are discussed separately and details are provided regarding their causes, magnitudes, and remedies. Subsequently, evidence of how stress-alleviating drugs could be used in combination with chemotherapy to improve treatment efficacy is presented, highlighting the potential of stress-alleviation strategies to enhance cancer therapy. Finally, a continuum-level, mathematical framework to incorporate these types of solid stress is outlined.

On the morphological stability of multicellular tumour spheroids growing in porous media

Multicellular tumour spheroids (MCTSs) are extensively used as in vitro system models for investigating the avascular growth phase of solid tumours. In this work, we propose a continuous growth model of heterogeneous MCTSs within a porous material, taking into account a diffusing nutrient from the surrounding material directing both the proliferation rate and the mobility of tumour cells. At the time scale of interest, the MCTS behaves as an incompressible viscous fluid expanding inside a porous medium. The cell motion and proliferation rate are modelled using a non-convective chemotactic mass flux, driving the cell expansion in the direction of the external nutrients' source. At the early stages, the growth dynamics is derived by solving the quasi-stationary problem, obtaining an initial exponential growth followed by an almost linear regime, in accordance with experimental observations. We also perform a linear-stability analysis of the quasi-static solution in order to investigate the morphological stability of the radially symmetric growth pattern. We show that mechano-biological cues, as well as geometric effects related to the size of the MCTS subdomains with respect to the diffusion length of the nutrient, can drive a morphological transition to fingered structures, thus triggering the formation of complex shapes that might promote tumour invasiveness. The results also point out the formation of a retrograde flow in the MCTS close to the regions where protrusions form, that could describe the initial dynamics of metastasis detachment from the in vivo tumour mass. In conclusion, the results of the proposed model demonstrate that the integration of mathematical tools in biological research could be crucial for better understanding the tumour's ability to invade its host environment.

Hydrodynamics and convection enhanced macromolecular fluid transport in soft biological tissues: Application to solid tumor

This work addresses a theoretical framework for transvascular exchange and extravascular transport of solute macromolecules through soft interstitial space inside a solid tumor. Most of the soft biological tissues show materialistic properties similar to deformable porous material. They exhibit mechanical behavior towards the fluid motion since the solid phase of the tumor tissue gets compressed by the drag force that is associated with the extracellular fluid flow. This paper presents a general view about the transvascular and interstitial transport of solute nutrients inside a tumor in the macroscopic level. Modified Starling׳s equation is used to describe transvascular nutrient transport. On the macroscopic level, motion of extracellular fluid within the tumor interstitium is modeled with the help of biphasic mixture theory and a spherical symmetry solution is given as a simpler case. This present model describes the average interstitial fluid pressure (IFP), extracellular fluid velocity (EFV) and flow rate of extracellular fluid, as well as the deformation of the solid phase of the tumor tissue as an immediate cause of extracellular fluid flow. When the interstitial transport is diffusion dominated, an analytical treatment of advection-diffusion-reaction equation finds the overall nutrient distribution. We propose suitable criteria for the formation of necrosis within the tumor interstitium. This study introduces some parameters that represent the nutrient supply from tumor blood vessels into the tumor extracellular space. These transport parameters compete with the reversible nutrient metabolism of the tumor cells present in the interstitium. The present study also shows that the effectiveness factor corresponding to a first order nutrient metabolism may reach beyond unity if the strength of the distributive solute source assumes positive non-zero values.

Predicting the growth of glioblastoma multiforme spheroids using a multiphase porous media model

Tumor spheroids constitute an effective in vitro tool to investigate the avascular stage of tumor growth. These three-dimensional cell aggregates reproduce the nutrient and proliferation gradients found in the early stages of cancer and can be grown with a strict control of their environmental conditions. In the last years, new experimental techniques have been developed to determine the effect of mechanical stress on the growth of tumor spheroids. These studies report a reduction in cell proliferation as a function of increasingly applied stress on the surface of the spheroids. This work presents a specialization for tumor spheroid growth of a previous more general multiphase model. The equations of the model are derived in the framework of porous media theory, and constitutive relations for the mass transfer terms and the stress are formulated on the basis of experimental observations. A set of experiments is performed, investigating the growth of U-87MG spheroids both freely growing in the culture medium and subjected to an external mechanical pressure induced by a Dextran solution. The growth curves of the model are compared to the experimental data, with good agreement for both the experimental settings. A new mathematical law regulating the inhibitory effect of mechanical compression on cancer cell proliferation is presented at the end of the paper. This new law is validated against experimental data and provides better results compared to other expressions in the literature.

Computer Simulations of the Tumor Vasculature: Applications to Interstitial Fluid Flow, Drug Delivery, and Oxygen Supply

Tumor vasculature, the blood vessel network supplying a growing tumor with nutrients such as oxygen or glucose, is in many respects different from the hierarchically organized arterio-venous blood vessel network in normal tissues. Angiogenesis (the formation of new blood vessels), vessel cooption (the integration of existing blood vessels into the tumor vasculature), and vessel regression remodel the healthy vascular network into a tumor-specific vasculature. Integrative models, based on detailed experimental data and physical laws, implement, in silico, the complex interplay of molecular pathways, cell proliferation, migration, and death, tissue microenvironment, mechanical and hydrodynamic forces, and the fine structure of the host tissue vasculature. With the help of computer simulations high-precision information about blood flow patterns, interstitial fluid flow, drug distribution, oxygen and nutrient distribution can be obtained and a plethora of therapeutic protocols can be tested before clinical trials. This chapter provides an overview over the current status of computer simulations of vascular remodeling during tumor growth including interstitial fluid flow, drug delivery, and oxygen supply within the tumor. The model predictions are compared with experimental and clinical data and a number of longstanding physiological paradigms about tumor vasculature and intratumoral solute transport are critically scrutinized.

Multiscale modelling of solid tumour growth: the effect of collagen micromechanics

Here we introduce a model of solid tumour growth coupled with a multiscale biomechanical description of the tumour microenvironment, which facilitates the explicit simulation of fibre-fibre and tumour-fibre interactions. We hypothesise that such a model, which provides a purely mechanical description of tumour-host interactions, can be used to explain experimental observations of the effect of collagen micromechanics on solid tumour growth. The model was specified to mouse tumour data, and numerical simulations were performed. The multiscale model produced lower stresses than an equivalent continuum-like approach, due to a more realistic remodelling of the collagen microstructure. Furthermore, solid tumour growth was found to cause a passive mechanical realignment of fibres at the tumour boundary from a random to a circumferential orientation. This is in accordance with experimental observations, thus demonstrating that such a response can be explained as purely mechanical. Finally, peritumoural fibre network anisotropy was found to produce anisotropic tumour morphology. The dependency of tumour morphology on the peritumoural microstructure was reduced by adding a load-bearing non-collagenous component to the fibre network constitutive equation.

Stress-mediated progression of solid tumors: effect of mechanical stress on tissue oxygenation, cancer cell proliferation, and drug delivery

Oxygen supply plays a central role in cancer cell proliferation. While vascular density increases at the early stages of carcinogenesis, mechanical solid stresses developed during growth compress tumor blood vessels and, thus, drastically reduce not only the supply of oxygen, but also the delivery of drugs at inner tumor regions. Among other effects, hypoxia and reduced drug delivery compromise the efficacy of radiation and chemo/nanotherapy, respectively. In the present study, we developed a mathematical model of tumor growth to investigate the interconnections among tumor oxygenation that supports cancer cell proliferation, the heterogeneous accumulation of mechanical stresses owing to tumor growth, the non-uniform compression of intratumoral blood vessels due to the mechanical stresses, and the insufficient delivery of oxygen and therapeutic agents because of vessel compression. We found that the high vascular density and increased cancer cell proliferation often observed in the periphery compared to the interior of a tumor can be attributed to heterogeneous solid stress accumulation. Highly vascularized peripheral regions are also associated with greater oxygenation compared with the compressed, less vascularized inner regions. We also modeled the delivery of drugs of two distinct sizes, namely chemotherapy and nanomedicine. Model predictions suggest that drug delivery is affected negatively by vessel compression independently of the size of the therapeutic agent. Finally, we demonstrated the applicability of our model to actual geometries, employing a breast tumor model derived from MR images.

Integrative models of vascular remodeling during tumor growth

Malignant solid tumors recruit the blood vessel network of the host tissue for nutrient supply, continuous growth, and gain of metastatic potential. Angiogenesis (the formation of new blood vessels), vessel cooption (the integration of existing blood vessels into the tumor vasculature), and vessel regression remodel the healthy vascular network into a tumor-specific vasculature that is in many respects different from the hierarchically organized arterio-venous blood vessel network of the host tissues. Integrative models based on detailed experimental data and physical laws implement in silico the complex interplay of molecular pathways, cell proliferation, migration, and death, tissue microenvironment, mechanical and hydrodynamic forces, and the fine structure of the host tissue vasculature. With the help of computer simulations high-precision information about blood flow patterns, interstitial fluid flow, drug distribution, oxygen and nutrient distribution can be obtained and a plethora of therapeutic protocols can be tested before clinical trials. In this review, we give an overview over the current status of integrative models describing tumor growth, vascular remodeling, blood and interstitial fluid flow, drug delivery, and concomitant transformations of the microenvironment. WIREs Syst Biol Med 2015, 7:113-129. doi: 10.1002/wsbm.1295 For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

The role of mechanical forces in tumor growth and therapy

Tumors generate physical forces during growth and progression. These physical forces are able to compress blood and lymphatic vessels, reducing perfusion rates and creating hypoxia. When exerted directly on cancer cells, they can increase cells' invasive and metastatic potential. Tumor vessels-while nourishing the tumor-are usually leaky and tortuous, which further decreases perfusion. Hypoperfusion and hypoxia contribute to immune evasion, promote malignant progression and metastasis, and reduce the efficacy of a number of therapies, including radiation. In parallel, vessel leakiness together with vessel compression causes a uniformly elevated interstitial fluid pressure that hinders delivery of blood-borne therapeutic agents, lowering the efficacy of chemo- and nanotherapies. In addition, shear stresses exerted by flowing blood and interstitial fluid modulate the behavior of cancer and a variety of host cells. Taming these physical forces can improve therapeutic outcomes in many cancers.

Magnetic flattening of stem-cell spheroids indicates a size-dependent elastocapillary transition

Cellular aggregates (spheroids) are widely used in biophysics and tissue engineering as model systems for biological tissues. In this Letter we propose novel methods for molding stem-cell spheroids, deforming them, and measuring their interfacial and elastic properties with a single method based on cell tagging with magnetic nanoparticles and application of a magnetic field gradient. Magnetic molding yields spheroids of unprecedented sizes (up to a few mm in diameter) and preserves tissue integrity. On subjecting these spheroids to magnetic flattening (over 150g), we observed a size-dependent elastocapillary transition with two modes of deformation: liquid-drop-like behavior for small spheroids, and elastic-sphere-like behavior for larger spheroids, followed by relaxation to a liquidlike drop.

A tumor growth model with deformable ECM

Existing tumor growth models based on fluid analogy for the cells do not generally include the extracellular matrix (ECM), or if present, take it as rigid. The three-fluid model originally proposed by the authors and comprising tumor cells (TC), host cells (HC), interstitial fluid (IF) and an ECM, considered up to now only a rigid ECM in the applications. This limitation is here relaxed and the deformability of the ECM is investigated in detail. The ECM is modeled as a porous solid matrix with Green-elastic and elasto-visco-plastic material behavior within a large strain approach. Jauman and Truesdell objective stress measures are adopted together with the deformation rate tensor. Numerical results are first compared with those of a reference experiment of a multicellular tumor spheroid (MTS) growing in vitro, then three different tumor cases are studied: growth of an MTS in a decellularized ECM, growth of a spheroid in the presence of host cells and growth of a melanoma. The influence of the stiffness of the ECM is evidenced and comparison with the case of a rigid ECM is made. The processes in a deformable ECM are more rapid than in a rigid ECM and the obtained growth pattern differs. The reasons for this are due to the changes in porosity induced by the tumor growth. These changes are inhibited in a rigid ECM. This enhanced computational model emphasizes the importance of properly characterizing the biomechanical behavior of the malignant mass in all its components to correctly predict its temporal and spatial pattern evolution.

Role of constitutive behavior and tumor-host mechanical interactions in the state of stress and growth of solid tumors

Mechanical forces play a crucial role in tumor patho-physiology. Compression of cancer cells inhibits their proliferation rate, induces apoptosis and enhances their invasive and metastatic potential. Additionally, compression of intratumor blood vessels reduces the supply of oxygen, nutrients and drugs, affecting tumor progression and treatment. Despite the great importance of the mechanical microenvironment to the pathology of cancer, there are limited studies for the constitutive modeling and the mechanical properties of tumors and on how these parameters affect tumor growth. Also, the contribution of the host tissue to the growth and state of stress of the tumor remains unclear. To this end, we performed unconfined compression experiments in two tumor types and found that the experimental stress-strain response is better fitted to an exponential constitutive equation compared to the widely used neo-Hookean and Blatz-Ko models. Subsequently, we incorporated the constitutive equations along with the corresponding values of the mechanical properties - calculated by the fit - to a biomechanical model of tumor growth. Interestingly, we found that the evolution of stress and the growth rate of the tumor are independent from the selection of the constitutive equation, but depend strongly on the mechanical interactions with the surrounding host tissue. Particularly, model predictions - in agreement with experimental studies - suggest that the stiffness of solid tumors should exceed a critical value compared with that of the surrounding tissue in order to be able to displace the tissue and grow in size. With the use of the model, we estimated this critical value to be on the order of 1.5. Our results suggest that the direct effect of solid stress on tumor growth involves not only the inhibitory effect of stress on cancer cell proliferation and the induction of apoptosis, but also the resistance of the surrounding tissue to tumor expansion.

Tumor growth in complex, evolving microenvironmental geometries: a diffuse domain approach

We develop a mathematical model of tumor growth in complex, dynamic microenvironments with active, deformable membranes. Using a diffuse domain approach, the complex domain is captured implicitly using an auxiliary function and the governing equations are appropriately modified, extended and solved in a larger, regular domain. The diffuse domain method enables us to develop an efficient numerical implementation that does not depend on the space dimension or the microenvironmental geometry. We model homotypic cell-cell adhesion and heterotypic cell-basement membrane (BM) adhesion with the latter being implemented via a membrane energy that models cell-BM interactions. We incorporate simple models of elastic forces and the degradation of the BM and ECM by tumor-secreted matrix degrading enzymes. We investigate tumor progression and BM response as a function of cell-BM adhesion and the stiffness of the BM. We find tumor sizes tend to be positively correlated with cell-BM adhesion since increasing cell-BM adhesion results in thinner, more elongated tumors. Prior to invasion of the tumor into the stroma, we find a negative correlation between tumor size and BM stiffness as the elastic restoring forces tend to inhibit tumor growth. In order to model tumor invasion of the stroma, we find it necessary to downregulate cell-BM adhesiveness, which is consistent with experimental observations. A stiff BM promotes invasiveness because at early stages the opening in the BM created by MDE degradation from tumor cells tends to be narrower when the BM is stiffer. This requires invading cells to squeeze through the narrow opening and thus promotes fragmentation that then leads to enhanced growth and invasion. In three dimensions, the opening in the BM was found to increase in size even when the BM is stiff because of pressure induced by growing tumor clusters. A larger opening in the BM can increase the potential for further invasiveness by increasing the possibility that additional tumor cells could invade the stroma.

A stable scheme for a nonlinear, multiphase tumor growth model with an elastic membrane

In this paper, we extend the 3D multispecies diffuse-interface model of the tumor growth, which was derived in Wise et al. (Three-dimensional multispecies nonlinear tumor growth-I: model and numerical method, J. Theor. Biol. 253 (2008) 524-543), and incorporate the effect of a stiff membrane to model tumor growth in a confined microenvironment. We then develop accurate and efficient numerical methods to solve the model. When the membrane is endowed with a surface energy, the model is variational, and the numerical scheme, which involves adaptive mesh refinement and a nonlinear multigrid finite difference method, is demonstrably shown to be energy stable. Namely, in the absence of cell proliferation and death, the discrete energy is a nonincreasing function of time for any time and space steps. When a simplified model of membrane elastic energy is used, the resulting model is derived analogously to the surface energy case. However, the elastic energy model is actually nonvariational because certain coupling terms are neglected. Nevertheless, a very stable numerical scheme is developed following the strategy used in the surface energy case. 2D and 3D simulations are performed that demonstrate the accuracy of the algorithm and illustrate the shape instabilities and nonlinear effects of membrane elastic forces that may resist or enhance growth of the tumor. Compared with the standard Crank-Nicholson method, the time step can be up to 25 times larger using the new approach.

Physiology-based model of cell viscoelasticity

The measured viscoelastic properties of biological tissues is the result of the passive and active response of the cells. We propose an evolution law of the remodeling process in the cytoskeleton which is able to mimic the viscous properties of biological cellular tissues. Our model is based on dynamical changes of the resting length. We show that under the small strain regime, the linear rheology models are recovered, with the relaxation time being replaced by the cell resistance to remodel. We implement the one-dimensional model into network systems of two and three dimensions, and show that the same conclusions may be drawn for those systems.

Coevolution of solid stress and interstitial fluid pressure in tumors during progression: implications for vascular collapse

The stress harbored by the solid phase of tumors is known as solid stress. Solid stress can be either applied externally by the surrounding normal tissue or induced by the tumor itself due to its growth. Fluid pressure is the isotropic stress exerted by the fluid phase. We recently showed that growth-induced solid stress is on the order of 1.3 to 13.0 kPa (10-100 mmHg)--high enough to cause compression of fragile blood vessels, resulting in poor perfusion and hypoxia. However, the evolution of growth-induced stress with tumor progression and its effect on cancer cell proliferation in vivo is not understood. To this end, we developed a mathematical model for tumor growth that takes into account all three types of stresses: growth-induced stress, externally applied stress, and fluid pressure. First, we conducted in vivo experiments and found that growth-induced stress is related to tumor volume through a biexponential relationship. Then, we incorporated this information into our mathematical model and showed that due to the evolution of growth-induced stress, total solid stress levels are higher in the tumor interior and lower in the periphery. Elevated compressive solid stress in the interior of the tumor is sufficient to cause the collapse of blood vessels and results in a lower growth rate of cancer cells compared with the periphery, independently from that caused by the lack of nutrients due to vessel collapse. Furthermore, solid stress in the periphery of the tumor causes blood vessels in the surrounding normal tissue to deform to elliptical shapes. We present histologic sections of human cancers that show such vessel deformations. Finally, we found that fluid pressure increases with tumor growth due to increased vascular permeability and lymphatic impairment, and is governed by the microvascular pressure. Crucially, fluid pressure does not cause vessel compression of tumor vessels.

Buckling instability in growing tumor spheroids

A growing tumor is subjected to intrinsic physical forces, arising from the cellular turnover in a spatially constrained environment. This work demonstrates that such residual solid stresses can provoke a buckling instability in heterogeneous tumor spheroids. The growth rate ratio between the outer shell of proliferative cells and the inner necrotic core is the control parameter of this instability. The buckled morphology is found to depend both on the elastic and the geometric properties of the tumor components, suggesting a key role of residual stresses for promoting tumor invasiveness.

In vivo mathematical modeling of tumor growth from imaging data: soon to come in the future?

The future challenges in oncology imaging are to assess the response to treatment even earlier. As an addition to functional imaging, mathematical modeling based on the imaging is an alternative, cross-disciplinary area of development. Modeling was developed in oncology not only in order to understand and predict tumor growth, but also to anticipate the effects of targeted and untargeted therapies. A very wide range of these models exist, involving many stages in the progression of tumors. Few models, however, have been proposed to reproduce in vivo tumor growth because of the complexity of the mechanisms involved. Morphological imaging combined with "spatial" models appears to perform well although functioning imaging could still provide further information on metabolism and the micro-architecture. The combination of imaging and modeling can resolve complex problems and describe many facets of tumor growth or response to treatment. It is now possible to consider its clinical use in the medium term. This review describes the basic principles of mathematical modeling and describes the advantages, limitations and future prospects for this in vivo approach based on imaging data.

Mechano-transduction in tumour growth modelling

The evolution of biological systems is strongly influenced by physical factors, such as applied forces, geometry or the stiffness of the micro-environment. Mechanical changes are particularly important in solid tumour development, as altered stromal-epithelial interactions can provoke a persistent increase in cytoskeletal tension, driving the gene expression of a malignant phenotype. In this work, we propose a novel multi-scale treatment of mechano-transduction in cancer growth. The avascular tumour is modelled as an expanding elastic spheroid, whilst growth may occur both as a volume increase and as a mass production within a cell rim. Considering the physical constraints of an outer healthy tissue, we derive the thermo-dynamical requirements for coupling growth rate, solid stress and diffusing biomolecules inside a heterogeneous tumour. The theoretical predictions successfully reproduce the stress-dependent growth curves observed by in vitro experiments on multicellular spheroids.

Soft matter models of developing tissues and tumors

Analogies with inert soft condensed matter--such as viscoelastic liquids, pastes, foams, emulsions, colloids, and polymers--can be used to investigate the mechanical response of soft biological tissues to forces. A variety of experimental techniques and biophysical models have exploited these analogies allowing the quantitative characterization of the mechanical properties of model tissues, such as surface tension, elasticity, and viscosity. The framework of soft matter has been successful in explaining a number of dynamical tissue behaviors observed in physiology and development, such as cell sorting, tissue spreading, or the escape of individual cells from a tumor. However, living tissues also exhibit active responses, such as rigidity sensing or cell pulsation, that are absent in inert soft materials. The soft matter models reviewed here have provided valuable insight in understanding morphogenesis and cancer invasion and have set bases for using tissue engineering within medicine.

Tissue dynamics with permeation

Animal tissues are complex assemblies of cells, extracellular matrix (ECM), and permeating interstitial fluid. Whereas key aspects of the multicellular dynamics can be captured by a one-component continuum description, cell division and apoptosis imply material turnover between different components that can lead to additional mechanical conditions on the tissue dynamics. We extend our previous description of tissues in order to account for a cell/ECM phase and the permeating interstitial fluid independently. In line with our earlier work, we consider the cell/ECM phase to behave as an elastic solid in the absence of cell division and apoptosis. In addition, we consider the interstitial fluid as ideal on the relevant length scales, i.e., we ignore viscous stresses in the interstitial fluid. Friction between the fluid and the cell/ECM phase leads to a Darcy-like relation for the interstitial fluid velocity and introduces a new characteristic length scale. We discuss the dynamics of a tissue confined in a chamber with a permeable piston close to the homeostatic state where cell division and apoptosis balance, and we calculate the rescaled effective diffusion coefficient for cells. For different mass densities of the cell/ECM component and the interstitial fluid, a treadmilling steady state due to gravitational forces can be found.

Dynamic density functional theory of solid tumor growth: Preliminary models

Cancer is a disease that can be seen as a complex system whose dynamics and growth result from nonlinear processes coupled across wide ranges of spatio-temporal scales. The current mathematical modeling literature addresses issues at various scales but the development of theoretical methodologies capable of bridging gaps across scales needs further study. We present a new theoretical framework based on Dynamic Density Functional Theory (DDFT) extended, for the first time, to the dynamics of living tissues by accounting for cell density correlations, different cell types, phenotypes and cell birth/death processes, in order to provide a biophysically consistent description of processes across the scales. We present an application of this approach to tumor growth.

Invasion from a cell aggregate--the roles of active cell motion and mechanical equilibrium

Cell invasion from an aggregate into a surrounding extracellular matrix (ECM) is an important process during development disease, e.g., vascular network assembly or tumor progression. To describe the behavior emerging from autonomous cell motility, cell-cell adhesion and contact guidance by ECM filaments, we propose a suitably modified cellular Potts model. We consider an active cell motility process in which internal polarity is governed by a positive feedback from cell displacements, a mechanism that can result in highly persistent motion when constrained by an oriented ECM structure. The model allows us to explore the interplay between haptotaxis, matrix degradation and active cell movement. We show that for certain conditions the cells are able to both invade the ECM and follow the ECM tracks. Furthermore, we argue that enforcing mechanical equilibrium within a bulk cell mass is of key importance in multicellular simulations.

Micrometer-scale oxygen delivery rearranges cells and prevents necrosis in tumor tissue in vitro

Oxygen availability plays a critical role in cancer progression and is correlated with poor prognosis. Despite this connection, the independent effects of oxygen gradients on tumor tissues have not been measured. To address this, we developed an oxygen delivery device that uses microelectrodes to generate oxygen directly underneath three-dimensional tumor cylindroids composed of colon carcinoma cells. The extent of cell death was measured using fluorescence staining. Supplying oxygen for 60 h eliminated the necrotic region typically found in the center of cylindroids despite the continued presence of other nutrient gradients. A mathematical model of cylindroid growth showed that the rate of cell death was more sensitive to oxygen than the growth rate. After oxygenation, a ring of dead cells was observed at the outside edge of cylindroids, and dead cells were observed moving outward from cylindroid centers. This movement suggests that dead cells were pushed by viable cells migrating in response to oxygen gradients, a mechanism that may connect transient oxygen gradients to metastasis formation. These measurements show that oxygen gradients are a primary factor governing cell viability and rearrange cells in tumors.