Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method

Patient-Specific, Mechanistic Models of Tumor Growth Incorporating Artificial Intelligence and Big Data

Despite the remarkable advances in cancer diagnosis, treatment, and management over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires the integration of patient-specific information with an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous yet practical mathematical theory of tumor initiation, development, invasion, and response to therapy. We begin this review with an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on big data and artificial intelligence. We then present illustrative examples of mathematical models manifesting their utility and discuss the limitations of stand-alone mechanistic and data-driven models. We then discuss the potential of mechanistic models for not only predicting but also optimizing response to therapy on a patient-specific basis. We describe current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.

Numerical investigation of the pharmacokinetics and pharmacodynamics of the chemotherapeutic drug in avascular and vascular stages of a brain tumor

One of the most commonly used approaches for treating solid tumors is the systemic delivery of chemotherapeutic drugs. However, our understanding of the factors influencing treatment efficacy through this method is still limited. This study presents a comprehensive and realistic mathematical model that incorporates the dynamics of tumor growth, capillary network extension, and drug delivery in a coupled and simultaneous manner. The model covers two stages of tumor growth: avascular and vascular. For tumor growth, a continuum model is employed using the phase field interface-capturing method. The neo-vascularization process is modeled using a hybrid discrete-continuum approach. Additionally, a multi-scale model is used to describe the pharmacokinetics of doxorubicin, considering various agents. The study investigates the effect of haptotaxis and reveals that a higher haptotaxis coefficient leads to faster tumor growth (up to 2.6 times) and a quicker progression to angiogenesis. The impact of tumor-related and drug-related parameters is also examined, including tumor size, tumor sensitivity to the drug, chemotherapy initialization, treatment cycle duration, drug affinity to cells, and drug dose. The findings indicate that chemotherapy is more effective during the angiogenesis stage when active loops have formed. Other clinical methods such as radiotherapy and surgery may be more appropriate during the avascular stage or the transition period between angiogenesis initialization and loop formation. The penetration depth of the drug decreases by approximately 50% with an increase in the drug binding rate to surface-cell receptors. As a result, high-associate-rate drugs are preferred for chemotherapy after active loops have formed, while low-associate-rate drugs are suitable for earlier stages.

A phase-field model for non-small cell lung cancer under the effects of immunotherapy

Formulating mathematical models that estimate tumor growth under therapy is vital for improving patient-specific treatment plans. In this context, we present our recent work on simulating non-small-scale cell lung cancer (NSCLC) in a simple, deterministic setting for two different patients receiving an immunotherapeutic treatment. At its core, our model consists of a Cahn-Hilliard-based phase-field model describing the evolution of proliferative and necrotic tumor cells. These are coupled to a simplified nutrient model that drives the growth of the proliferative cells and their decay into necrotic cells. The applied immunotherapy decreases the proliferative cell concentration. Here, we model the immunotherapeutic agent concentration in the entire lung over time by an ordinary differential equation (ODE). Finally, reaction terms provide a coupling between all these equations. By assuming spherical, symmetric tumor growth and constant nutrient inflow, we simplify this full 3D cancer simulation model to a reduced 1D model. We can then resort to patient data gathered from computed tomography (CT) scans over several years to calibrate our model. Our model covers the case in which the immunotherapy is successful and limits the tumor size, as well as the case predicting a sudden relapse, leading to exponential tumor growth. Finally, we move from the reduced model back to the full 3D cancer simulation in the lung tissue. Thereby, we demonstrate the predictive benefits that a more detailed patient-specific simulation including spatial information as a possible generalization within our framework could yield in the future.

A morphological study of bicontinuous concentric lamellar silica synthesized at atmospheric pressure and its application as an internal micro-reflector in dye-sensitized solar cells

KCC-1, a nanostructured silica material with a bicontinuous concentric lamellar (bcl) morphology, provides plenty of functional characteristics, such as an open channel structure, excellent accessibility, and a large surface area. Although bcl silica exhibits various superior properties, studies on its morphology and its application in dye-sensitized solar cells (DSSCs) are still limited. Therefore, this work aims to study the influence of the synthesis time on the morphology of bcl silica. Moreover, we used the synthesized bcl silica as internal micro-reflectors in DSSCs. The bcl silica was synthesized using the reflux method by varying synthesis times. The morphology of bcl silica was observed using FESEM and HRTEM. FESEM images show that bcl silica has bicontinuous lamellar walls arranged concentrically to form spherical particles. As the synthesis time increases, the average particle size of bcl silica increases. The quantization of bcl silica binary images shows that the average lamellar cross-sectional area ratio decreases with increasing synthesis time. The simulation of the Cahn-Hilliard's spinodal decomposition model using MATLAB also describes the lamellar cross-sectional area ratio of bcl silica. In addition, to characterize the FESEM image's texture, a Shannon entropy calculation was performed. The line and circular gray value intensity profiles of the HRTEM image show that bcl silica has a denser core than the outer part. The denser core proves that the lamellae in bcl silica are concentrically arranged towards the particle core. Furthermore, we added bcl silica to a photoanode to see the effect of bcl characteristics on the DSSC performance. The results show that the bcl silica significantly improves the light-harvesting efficiency in DSSCs due to its low refractive index and open channel structure.

A multiscale chemical-mechanical model predicts impact of morphogen spreading on tissue growth

The exact mechanism controlling cell growth remains a grand challenge in developmental biology and regenerative medicine. The Drosophila wing disc tissue serves as an ideal biological model to study mechanisms involved in growth regulation. Most existing computational models for studying tissue growth focus specifically on either chemical signals or mechanical forces. Here we developed a multiscale chemical-mechanical model to investigate the growth regulation mechanism based on the dynamics of a morphogen gradient. By comparing the spatial distribution of dividing cells and the overall tissue shape obtained in model simulations with experimental data of the wing disc, it is shown that the size of the domain of the Dpp morphogen is critical in determining tissue size and shape. A larger tissue size with a faster growth rate and more symmetric shape can be achieved if the Dpp gradient spreads in a larger domain. Together with Dpp absorbance at the peripheral zone, the feedback regulation that downregulates Dpp receptors on the cell membrane allows for further spreading of the morphogen away from its source region, resulting in prolonged tissue growth at a more spatially homogeneous growth rate.

Morphological instability of solid tumors in a nutrient-deficient environment

A phenomenological reaction-diffusion model that includes a nutrient-regulated growth rate of tumor cells is proposed to investigate the morphological instability of solid tumors during the avascular growth. We find that the surface instability could be induced more easily when tumor cells are placed in a harsher nutrient-deficient environment, while the instability is suppressed for tumor cells in a nutrient-rich environment due to the nutrient-regulated proliferation. In addition, the surface instability is shown to be influenced by the growth moving speed of tumor rims. Our analysis reveals that a larger growth movement of the tumor front results in a closer proximity of tumor cells to a nutrient-rich region, which tends to inhibit the surface instability. A nourished length that represents the proximity is defined to illustrate its close relation to the surface instability.

Tumor Evolution Models of Phase-Field Type with Nonlocal Effects and Angiogenesis

In this survey article, a variety of systems modeling tumor growth are discussed. In accordance with the hallmarks of cancer, the described models incorporate the primary characteristics of cancer evolution. Specifically, we focus on diffusive interface models and follow the phase-field approach that describes the tumor as a collection of cells. Such systems are based on a multiphase approach that employs constitutive laws and balance laws for individual constituents. In mathematical oncology, numerous biological phenomena are involved, including temporal and spatial nonlocal effects, complex nonlinearities, stochasticity, and mixed-dimensional couplings. Using the models, for instance, we can express angiogenesis and cell-to-matrix adhesion effects. Finally, we offer some methods for numerically approximating the models and show simulations of the tumor's evolution in response to various biological effects.

Predicting the effect of a combination drug therapy on the prostate tumor growth via an improvement of a direct radial basis function partition of unity technique for a diffuse-interface model

Chemotherapy is usually applied to treat advanced prostate cancer that cancer cells spread outside the prostate gland. The treatment uses cytotoxic drugs to target cells that grow and divide quickly. On the other hand, the growth of such cancerous tumors depends on angiogenesis. In this paper, we numerically study a diffuse-interface model in a two-dimensional space related to the therapies of prostate cancer. The proposed model describes the tumor growth driven by a generic nutrient and producing the prostate-specific antigen. More precisely, the effect of cytotoxic chemotherapy in the model is evaluated by considering a time-dependent function in the tumor dynamics. Also, another function related to the antiangiogenic therapy is considered to show the reducing intratumoral nutrient supply in the nutrient dynamics. Here, a meshless approximation, i.e., a generalized form of the direct radial basis function partition of unity (D-RBF-PU) method is presented for finding the numerical simulations of this model utilizing in medical oncology. The method uses the lower number of trial points in each patch than the original D-RBF-PU scheme for approximating the trial function per test point. Hence, the time complexity of the method is less than the D-RBF-PU technique. Besides, a semi-implicit time discretization of order 1 has been used to deal with the time variable. Consequently, a linear system of algebraic equations could be solved iteratively per time step by the use of the biconjugate gradient stabilized method with zero-fill incomplete lower-upper preconditioner. Finally, the obtained results without using any adaptive algorithm demonstrate the response of the prostate tumor growth to the chemotherapy, antiangiogenic therapy and a combined therapy.

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

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.

The Dynamics of a Fractional-Order Mathematical Model of Cancer Tumor Disease

This article explores the application of the reduced differential transform method (RDTM) for the computational solutions of two fractional-order cancer tumor models in the Caputo sense: the model based on cancer chemotherapeutic effects which explain the relation between chemotherapeutic drugs, tumor cells, normal cells, and immune cells using a fractional partial differential equations, and the model that describes the different cases of killing rate K of cancer cells (the killing percentage of cancer cells K (I) is dependent on the number of cells, (II) is a function of time only, and (III) is a function of space only). The solutions are presented using Mathematica software as a convergent power series with elegantly computed terms using the suggested technique. The proposed method gives new series form results for various values of gamma. To clarify the complexity of the models, we plot the two- and three-dimensional and contour graphics of the obtained solutions at varied values of fractional-order gamma and the selected system parameters. The solutions are analyzed with fractional and reduced differential transform methods to obtain an idea of invariance regarding the computed solution of the designed mathematical model. The obtained results demonstrate the efficiency and preciseness of the proposed method to achieve a better understanding of chemotherapy effects. It is observed that chemotherapy drugs boost immunity against the specific cancer by decreasing the number of tumor cells, and the killing rate K of cancerous cells depend on the cells concentration.

Tumor growth towards lower extracellular matrix conductivity regions under Darcy's Law and steady morphology

We study a classic Darcy's law model for tumor cell motion with inhomogeneous and isotropic conductivity. The tumor cells are assumed to be a constant density fluid flowing through porous extracellular matrix (ECM). The ECM is assumed to be rigid and motionless with constant porosity. One and two dimensional simulations show that the tumor mass grows from high to low conductivity regions when the tumor morphology is steady. In the one-dimensional case, we proved that when the tumor size is steady, the tumor grows towards lower conductivity regions. We conclude that this phenomenon is produced by the coupling of a special inward flow pattern in the steady tumor and Darcy's law which gives faster flow speed in higher conductivity regions.

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.

Mechanics of 3D Cell-Hydrogel Interactions: Experiments, Models, and Mechanisms

Hydrogels are highly water-swollen molecular networks that are ideal platforms to create tissue mimetics owing to their vast and tunable properties. As such, hydrogels are promising cell-delivery vehicles for applications in tissue engineering and have also emerged as an important base for ex vivo models to study healthy and pathophysiological events in a carefully controlled three-dimensional environment. Cells are readily encapsulated in hydrogels resulting in a plethora of biochemical and mechanical communication mechanisms, which recapitulates the natural cell and extracellular matrix interaction in tissues. These interactions are complex, with multiple events that are invariably coupled and spanning multiple length and time scales. To study and identify the underlying mechanisms involved, an integrated experimental and computational approach is ideally needed. This review discusses the state of our knowledge on cell-hydrogel interactions, with a focus on mechanics and transport, and in this context, highlights recent advancements in experiments, mathematical and computational modeling. The review begins with a background on the thermodynamics and physics fundamentals that govern hydrogel mechanics and transport. The review focuses on two main classes of hydrogels, described as semiflexible polymer networks that represent physically cross-linked fibrous hydrogels and flexible polymer networks representing the chemically cross-linked synthetic and natural hydrogels. In this review, we highlight five main cell-hydrogel interactions that involve key cellular functions related to communication, mechanosensing, migration, growth, and tissue deposition and elaboration. For each of these cellular functions, recent experiments and the most up to date modeling strategies are discussed and then followed by a summary of how to tune hydrogel properties to achieve a desired functional cellular outcome. We conclude with a summary linking these advancements and make the case for the need to integrate experiments and modeling to advance our fundamental understanding of cell-matrix interactions that will ultimately help identify new therapeutic approaches and enable successful tissue engineering.

Mathematical Modeling and Numerical Simulation of Atherosclerosis Based on a Novel Surgeon's View

This paper deals with the mathematical modeling of atherosclerosis based on a novel hypothesis proposed by a surgeon, Prof. Dr. Axel Haverich (Circulation 135(3):205-207, 2017). Atherosclerosis is referred as the thickening of the artery walls. Currently, there are two schools of thoughts for explaining the root of such phenomenon: thickening due to substance deposition and thickening as a result of inflammatory overgrowth. The hypothesis favored here is the second paradigm stating that the atherosclerosis is nothing else than the inflammatory response of of the wall tissues as a result of disruption in wall nourishment. It is known that a network of capillaries called vasa vasorum (VV) accounts for the nourishment of the wall in addition to the natural diffusion of nutrient from the blood passing through the lumen. Disruption of nutrient flow to the wall tissues may take place due to the occlusion of vasa vasorums with viruses, bacteria and very fine dust particles such as air pollutants referred to as PM 2.5. They can enter the body through the respiratory system at the first place and then reach the circulatory system. Hence in the new hypothesis, the root of atherosclerotic vessel is perceived as the malfunction of microvessels that nourish the vessel. A large number of clinical observation support this hypothesis. Recently and highly related to this work, and after the COVID-19 pandemic, one of the most prevalent disease in the lungs are attributed to the atherosclerotic pulmonary arteries, see Boyle and Haverich (Eur J Cardio Thorac Surg 58(6):1109-1110, 2020). In this work, a general framework is developed based on a multiphysics mathematical model to capture the wall deformation, nutrient availability and the inflammatory response. For the mechanical response an anisotropic constitutive relation is invoked in order to account for the presence of collagen fibers in the artery wall. A diffusion-reaction equation governs the transport of the nutrient within the wall. The inflammation (overgrowth) is described using a phase-field type equation with a double well potential which captures a sharp interface between two regions of the tissues, namely the healthy and the overgrowing part. The kinematics of the growth is treated by classical multiplicative decomposition of the gradient deformation. The inflammation is represented by means of a phase-field variable. A novel driving mechanism for the phase field is proposed for modeling the progression of the pathology. The model is 3D and fully based on the continuum description of the problem. The numerical implementation is carried out using FEM. Predictions of the model are compared with the clinical observations. The versatility and applicability of the model and the numerical tool allow.

A multi-scale agent-based model for avascular tumour growth

In this paper, we have developed a multi-scale, lattice-free, agent based model of avascular tumour growth in epithelial tissue. The model integrates different events to identify the underlying diversity within intracellular, cellular, and extracellular layer dynamics. The model considers every cell as an agent. A cellular agent may proliferate, spawns two identical daughter agents, or it may be transformed into other phenotypes during its life time depending on its internal proteins' activity as well as its external microenvironment. In this context, a simplified age-structured cell cycle model is adopted from the existing literature. The model considers that the intracellular events are regulated by p27 gene expression. In this model, p27 protein controls the overall tumour growth dynamics. Moreover, p27 is controlled by the external oxygen and nutrients that are modelled with the reaction-diffusion equations. The model also considers several biophysical forces which directly effect on the tumour growth dynamics. This modelling framework offers biologically realistic outcomes and also covers important criteria of the hallmarks of cancer which include oxygen and nutrient consumptions, micro-environmental heterogeneity, tumour cell proliferation by avoiding growth suppressor signals, replication of tumour cells at an abnormally faster rate, and resistance of apoptosis. The avascular tumour growth model is validated with immunohistochemistry and histopathology data. The outcome of the proposed model is very close to the range of the patient data, which concludes that the model is capable enough to mimic these complex biophysical phenomena.

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.

Mathematical simulation of tumour angiogenesis: angiopoietin balance is a key factor in vessel growth and regression

Excessive tumour growth results in a hypoxic environment around cancer cells, thus inducing tumour angiogenesis, which refers to the generation of new blood vessels from pre-existing vessels. This mechanism is biologically and physically complex, with various mathematical simulation models proposing to reproduce its formation. However, although temporary vessel regression is clinically known, few models succeed in reproducing this phenomenon. Here, we developed a three-dimensional simulation model encompassing both angiogenesis and tumour growth, specifically including angiopoietin. Angiopoietin regulates both adhesion and migration between vascular endothelial cells and wall cells, thus inhibiting the cell-to-cell adhesion required for angiogenesis initiation. Simulation results showed a regression, i.e. transient decrease, in the overall length of new vessels during vascular network formation. Using our model, we also evaluated the efficacy of administering the drug bevacizumab. The results highlighted differences in treatment efficacy: (1) earlier administration showed higher efficacy in inhibiting tumour growth, and (2) efficacy depended on the treatment interval even with the administration of the same dose. After thorough validation in the future, these results will contribute to the design of angiogenesis treatment protocols.

Generalized moving least squares approximation for the solution of local and non-local models of cancer cell invasion of tissue under the effect of adhesion in one- and two-dimensional spaces

The main aim of this study is to solve numerically the mathematical models showing cancer cell invasion of tissue with/without considering the effect of cell-cell and cell-matrix adhesion. The mathematical models studied here are the systems of time-dependent reaction-diffusion-taxis equations in one- and two-dimensional spaces, which are formulated in the local and non-local forms. There are some difficulties in finding their solutions via numerical methods. The main difficulty is to compute the non-local term appearing in one of the studied models, which causes more CPU time during simulations. The current paper aims to overcome this problem, where a new meshless method, namely generalized moving least squares (GMLS) approximation in space and a semi-implicit backward differential formula of first-order (SBDF1) in time have been applied. Based on GMLS theory, the non-local term is approximated without any difficulties. Moreover, a simple method based on the GMLS technique is presented to implement the boundary conditions. The obtained discrete scheme for both mathematical models is a linear system of algebraic equations per time step. The biconjugate gradient stabilized (BiCGSTAB) algorithm with zero-fill incomplete lower-upper (ILU) preconditioner is used to solve the obtained linear system at each time step. At the end of this paper, some simulation results are reported to show the behavior of cancer cell invasion in the local model, and the non-local model due to reduction of cell-cell adhesion and increasing cell-matrix adhesion in one- and two-dimensional spaces, where two different types of distribution points have been considered in the square domain. The computational algorithms of the GMLS approximation and the developed numerical method for solving the non-local (local) model are included in the Appendix.

Collective cell mechanics of epithelial shells with organoid-like morphologies

The study of organoids, artificially grown cell aggregates with the functionality and small-scale anatomy of real organs, is one of the most active areas of research in biology and biophysics, yet the basic physical origins of their different morphologies remain poorly understood. Here, we propose a mechanistic theory of epithelial shells which resemble small-organoid morphologies. Using a 3D surface tension-based vertex model, we reproduce the characteristic shapes from branched and budded to invaginated structures. We find that the formation of branched morphologies relies strongly on junctional activity, enabling temporary aggregations of topological defects in cell packing. To elucidate our numerical results, we develop an effective elasticity theory, which allows one to estimate the apico-basal polarity from the tissue-scale modulation of cell height. Our work provides a generic interpretation of the observed epithelial shell morphologies, highlighting the role of physical factors such as differential surface tension, cell rearrangements, and tissue growth.

The physical origins of man-made tissue morphologies with organ-like microscopic anatomy and functionality remain poorly understood. Here, authors propose a mechanistic theory of these structures, employing a 3D vertex model to reproduce the characteristic morphologies such as branched shapes.

A Coupled Mass Transport and Deformation Theory of Multi-constituent Tumor Growth

We develop a general class of thermodynamically consistent, continuum models based on mixture theory with phase effects that describe the behavior of a mass of multiple interacting constituents. The constituents consist of solid species undergoing large elastic deformations and compressible viscous fluids. The fundamental building blocks framing the mixture theories consist of the mass balance law of diffusing species and microscopic (cellular scale) and macroscopic (tissue scale) force balances, as well as energy balance and the entropy production inequality derived from the first and second laws of thermodynamics. A general phase-field framework is developed by closing the system through postulating constitutive equations (i.e., specific forms of free energy and rate of dissipation potentials) to depict the growth of tumors in a microenvironment. A notable feature of this theory is that it contains a unified continuum mechanics framework for addressing the interactions of multiple species evolving in both space and time and involved in biological growth of soft tissues (e.g., tumor cells and nutrients). The formulation also accounts for the regulating roles of the mechanical deformation on the growth of tumors, through a physically and mathematically consistent coupled diffusion and deformation framework. A new algorithm for numerical approximation of the proposed model using mixed finite elements is presented. The results of numerical experiments indicate that the proposed theory captures critical features of avascular tumor growth in the various microenvironment of living tissue, in agreement with the experimental studies in the literature.

Adhesion modulates cell morphology and migration within dense fibrous networks

One of the most fundamental abilities required for the sustainability of complex life forms is active cell migration, since it is essential in diverse processes from morphogenesis to leukocyte chemotaxis in immune response. The movement of a cell is the result of intricate mechanisms, that involve the coordination between mechanical forces, biochemical regulatory pathways and environmental cues. In particular, epithelial cancer cells have to employ mechanical strategies in order to migrate through the tissue's basement membrane and infiltrate the bloodstream during the invasion stage of metastasis. In this work we explore how mechanical interactions such as spatial restriction and adhesion affect migration of a self-propelled droplet in dense fibrous media. We have performed a systematic analysis using a phase-field model and we propose a novel approach to simulate cell migration with dissipative particle dynamics modelling. With this purpose we have measured in our simulation the cell's velocity and quantified its morphology as a function of the fibre density and of its adhesiveness to the matrix fibres. Furthermore, we have compared our results to a previousin vitromigration assay of fibrosarcoma cells in fibrous matrices. The results show good agreement between the two methodologies and experiments in the literature, which indicates that these minimalist descriptions are able to capture the main features of the system. Our results indicate that adhesiveness is critical for cell migration, by modulating cell morphology in crowded environments and by enhancing cell velocity. In addition, our analysis suggests that matrix metalloproteinases (MMPs) play an important role as adhesiveness modulators. We propose that new assays should be carried out to address the role of adhesion and the effect of different MMPs in cell migration under confined conditions.

Optimal Control Theory for Personalized Therapeutic Regimens in Oncology: Background, History, Challenges, and Opportunities

Optimal control theory is branch of mathematics that aims to optimize a solution to a dynamical system. While the concept of using optimal control theory to improve treatment regimens in oncology is not novel, many of the early applications of this mathematical technique were not designed to work with routinely available data or produce results that can eventually be translated to the clinical setting. The purpose of this review is to discuss clinically relevant considerations for formulating and solving optimal control problems for treating cancer patients. Our review focuses on two of the most widely used cancer treatments, radiation therapy and systemic therapy, as they naturally lend themselves to optimal control theory as a means to personalize therapeutic plans in a rigorous fashion. To provide context for optimal control theory to address either of these two modalities, we first discuss the major limitations and difficulties oncologists face when considering alternate regimens for their patients. We then provide a brief introduction to optimal control theory before formulating the optimal control problem in the context of radiation and systemic therapy. We also summarize examples from the literature that illustrate these concepts. Finally, we present both challenges and opportunities for dramatically improving patient outcomes via the integration of clinically relevant, patient-specific, mathematical models and optimal control theory.

Mathematical Modelling for the Role of CD4+T Cells in Tumor-Immune Interactions

Mathematical modelling has been used to study tumor-immune cell interaction. Some models were proposed to examine the effect of circulating lymphocytes, natural killer cells, and CD8+T cells, but they neglected the role of CD4+T cells. Other models were constructed to study the role of CD4+T cells but did not consider the role of other immune cells. In this study, we propose a mathematical model, in the form of a system of nonlinear ordinary differential equations, that predicts the interaction between tumor cells and natural killer cells, CD4+T cells, CD8+T cells, and circulating lymphocytes with or without immunotherapy and/or chemotherapy. This system is stiff, and the Runge-Kutta method failed to solve it. Consequently, the "Adams predictor-corrector" method is used. The results reveal that the patient's immune system can overcome small tumors; however, if the tumor is large, adoptive therapy with CD4+T cells can be an alternative to both CD8+T cell therapy and cytokines in some cases. Moreover, CD4+T cell therapy could replace chemotherapy depending upon tumor size. Even if a combination of chemotherapy and immunotherapy is necessary, using CD4+T cell therapy can better reduce the dose of the associated chemotherapy compared to using combined CD8+T cells and cytokine therapy. Stability analysis is performed for the studied patients. It has been found that all equilibrium points are unstable, and a condition for preventing tumor recurrence after treatment has been deduced. Finally, a bifurcation analysis is performed to study the effect of varying system parameters on the stability, and bifurcation points are specified. New equilibrium points are created or demolished at some bifurcation points, and stability is changed at some others. Hence, for systems turning to be stable, tumors can be eradicated without the possibility of recurrence. The proposed mathematical model provides a valuable tool for designing patients' treatment intervention strategies.

A level-set method for the evolution of cells and tissue during curvature-controlled growth

Most biological tissues grow by the synthesis of new material close to the tissue's interface, where spatial interactions can exert strong geometric influences on the local rate of growth. These geometric influences may be mechanistic or cell behavioural in nature. The control of geometry on tissue growth has been evidenced in many in vivo and in vitro experiments, including bone remodelling, wound healing, and tissue engineering scaffolds. In this paper, we propose a generalisation of a mathematical model that captures the mechanistic influence of curvature on the joint evolution of cell density and tissue shape during tissue growth. This generalisation allows us to simulate abrupt topological changes such as tissue fragmentation and tissue fusion, as well as three dimensional cases, through a level-set-based method. The level-set method developed introduces another Eulerian field than the level-set function. This additional field represents the surface density of tissue-synthesising cells, anticipated at future locations of the interface. Numerical tests performed with this level-set-based method show that numerical conservation of cells is a good indicator of simulation accuracy, particularly when cusps develop in the tissue's interface. We apply this new model to several situations of curvature-controlled tissue evolutions that include fragmentation and fusion.

Radiation protraction schedules for low-grade gliomas: a comparison between different mathematical models

We optimize radiotherapy (RT) administration strategies for treating low-grade gliomas. Specifically, we consider different tumour growth laws, both with and without spatial effects. In each scenario, we find the optimal treatment in the sense of maximizing the overall survival time of a virtual low-grade glioma patient, whose tumour progresses according to the examined growth laws. We discover that an extreme protraction therapeutic strategy, which amounts to substantially extending the time interval between RT sessions, may lead to better tumour control. The clinical implications of our results are also presented.

Continuum descriptions of spatial spreading for heterogeneous cell populations: Theory and experiment

Variability in cell populations is frequently observed in both in vitro and in vivo settings. Intrinsic differences within populations of cells, such as differences in cell sizes or differences in rates of cell motility, can be present even within a population of cells from the same cell line. We refer to this variability as cell heterogeneity. Mathematical models of cell migration, for example, in the context of tumour growth and metastatic invasion, often account for both undirected (random) migration and directed migration that is mediated by cell-to-cell contacts and cell-to-cell adhesion. A key feature of standard models is that they often assume that the population is composed of identical cells with constant properties. This leads to relatively simple single-species homogeneous models that neglect the role of heterogeneity. In this work, we use a continuum modelling approach to explore the role of heterogeneity in spatial spreading of cell populations. We employ a three-species heterogeneous model of cell motility that explicitly incorporates different types of experimentally-motivated heterogeneity in cell sizes: (i) monotonically decreasing; (ii) uniform; (iii) non-monotonic; and (iv) monotonically increasing distributions of cell size. Comparing the density profiles generated by the three-species heterogeneous model with density profiles predicted by a more standard single-species homogeneous model reveals that when we are dealing with monotonically decreasing and uniform distributions a simple and computationally efficient single-species homogeneous model can be remarkably accurate in describing the evolution of a heterogeneous cell population. In contrast, we find that the simpler single-species homogeneous model performs relatively poorly when applied to non-monotonic and monotonically increasing distributions of cell sizes. Additional results for heterogeneity in parameters describing both undirected and directed cell migration are also considered, and we find that similar results apply.

A Distributed Control Problem for a Fractional Tumor Growth Model

In this paper, we study the distributed optimal control of a system of three evolutionary equations involving fractional powers of three self-adjoint, monotone, unbounded linear operators having compact resolvents. The system is a generalization of a Cahn–Hilliard type phase field system modeling tumor growth that has been proposed by Hawkins–Daarud, van der Zee and Oden. The aim of the control process, which could be realized by either administering a drug or monitoring the nutrition, is to keep the tumor cell fraction under control while avoiding possible harm for the patient. In contrast to previous studies, in which the occurring unbounded operators governing the diffusional regimes were all given by the Laplacian with zero Neumann boundary conditions, the operators may in our case be different; more generally, we consider systems with fractional powers of the type that were studied in a recent work by the present authors. In our analysis, we show the Fréchet differentiability of the associated control-to-state operator, establish the existence of solutions to the associated adjoint system, and derive the first-order necessary conditions of optimality for a cost functional of tracking type.

Three-dimensional tumor growth in time-varying chemical fields: a modeling framework and theoretical study

Background

Contemporary biological observations have revealed a large variety of mechanisms acting during the expansion of a tumor. However, there are still many qualitative and quantitative aspects of the phenomenon that remain largely unknown. In this context, mathematical and computational modeling appears as an invaluable tool providing the means for conducting in silico experiments, which are cheaper and less tedious than real laboratory experiments.

Results

This paper aims at developing an extensible and computationally efficient framework for in silico modeling of tumor growth in a 3-dimensional, inhomogeneous and time-varying chemical environment. The resulting model consists of a set of mathematically derived and algorithmically defined operators, each one addressing the effects of a particular biological mechanism on the state of the system. These operators may be extended or re-adjusted, in case a different set of starting assumptions or a different simulation scenario needs to be considered.

Conclusion

In silico modeling provides an alternative means for testing hypotheses and simulating scenarios for which exact biological knowledge remains elusive. However, finer tuning of pertinent methods presupposes qualitative and quantitative enrichment of available biological evidence. Validation in a strict sense would further require comprehensive, case-specific simulations and detailed comparisons with biomedical observations.

Towards a Physiological Scale of Vocal Fold Agent-Based Models of Surgical Injury and Repair: Sensitivity Analysis, Calibration and Verification

Agent based models (ABM) were developed to numerically simulate the biological response to surgical vocal fold injury and repair at the physiological level. This study aimed to improve the representation of existing ABM through a combination of empirical and computational experiments. Empirical data of vocal fold cell populations including neutrophils, macrophages and fibroblasts were obtained using flow cytometry up to four weeks following surgical injury. Random Forests were used as a sensitivity analysis method to identify model parameters that were most influential to ABM outputs. Statistical Parameter Optimization Tool for Python was used to calibrate those parameter values to match the ABM-simulation data with the corresponding empirical data from Day 1 to Day 5 following surgery. Model performance was evaluated by verifying if the empirical data fell within the 95% confidence intervals of ABM outputs of cell quantities at Day 7, Week 2 and Week 4. For Day 7, all empirical data were within the ABM output ranges. The trends of ABM-simulated cell populations were also qualitatively comparable to those of the empirical data beyond Day 7. Exact values, however, fell outside of the 95% statistical confidence intervals. Parameters related to fibroblast proliferation were indicative to the ABM-simulation of fibroblast dynamics in final stages of wound healing.

Pattern formation of skin cancers: Effects of cancer proliferation and hydrodynamic interactions

We study pattern formation of skin cancers by means of numerical simulation of a binary system consisting of cancer and healthy cells. We extend the conventional model H for macrophase separations by considering a logistic growth of cancer cells and also a mechanical friction between dermis and epidermis. Importantly, our model exhibits a microphase separation due to the proliferation of cancer cells. By numerically solving the time evolution equations of the cancer composition and its velocity, we show that the phase separation kinetics strongly depends on the cell proliferation rate as well as on the strength of hydrodynamic interactions. A steady-state diagram of cancer patterns is established in terms of these two dynamical parameters and some of the patterns correspond to clinically observed cancer patterns. Furthermore, we examine in detail the time evolution of the average composition of cancer cells and the characteristic length of the microstructures. Our results demonstrate that different sequence of cancer patterns can be obtained by changing the proliferation rate and/or hydrodynamic interactions.

A 3D Multiscale Model to Explore the Role of EGFR Overexpression in Tumourigenesis

The epidermal growth factor receptor (EGFR) signalling cascade is one of the main pathways that regulate the survival and division of mammalian cells. It is also one of the most altered transduction pathways in cancer. Acquired mutations in the EGFR/ERK pathway can cause the overexpression of EGFR on the surface of the cell, while others downregulate the inactivation of switched on intracellular proteins such as Ras and Raf. This upregulates the activity of ERK and promotes cell division. We develop a 3D multiscale model to explore the role of EGFR overexpression on tumour initiation. In this model, cells are described as individual objects that move, interact, divide, proliferate, and die by apoptosis. We use Brownian Dynamics to describe the extracellular and intracellular regulations of cells as well as the spatial and stochastic effects influencing them. The fate of each cell depends on the number of active transcription factors in the nucleus. We use numerical simulations to investigate the individual and combined effects of mutations on the intracellular regulation of individual cells. Next, we show that the distance between active receptors increase the level of EGFR/ERK signalling. We demonstrate the usefulness of the model by quantifying the impact of mutational alterations in the EGFR/ERK pathway on the growth rate of in silico tumours.

Tumor growth and calcification in evolving microenvironmental geometries

In this paper, we apply the diffuse domain framework developed in Chen and Lowengrub (Tumor growth in complex, evolving microenvironmental geometries: A diffuse domain approach, J. Theor. Biol. 361 (2014) 14-30) to study the effects of a deformable basement membrane (BM) on the growth of a tumor in a confined, ductal geometry, such as ductal carcinoma in situ (DCIS). We use a continuum model of tumor microcalcification and investigate the tumor extent beyond the microcalcification. In order to solve the governing equations efficiently, we develop a stable nonlinear multigrid finite difference method. Two dimensional simulations are performed where the adhesion between tumor cells and the basement membrane is varied. Additional simulations considering the variation of duct radius and membrane stiffness are also conducted. The results demonstrate that enhanced membrane deformability promotes tumor growth and tumor calcification. When the duct radius is small, the cell-BM adhesion is weak or when the membrane is slightly deformed, the mammographic and pathologic tumor extents are linearly correlated, as predicted by Macklin et al. (J. Theor. Biol. 301 (2012) 122-140) using an agent-based model that does not account for the deformability of the basement membrane and the active forces that the membrane imparts on the tumor cells. Interestingly, we predict that when the duct radius is large, there is strong cell-BM adhesion or the membrane is highly deformed, the extents of the mammographic and pathologic tumors are instead quadratically correlated. The simulations can help surgeons to measure DCIS surgical margins while removing less non-cancerous tissue, and can improve targeting of intra- and post-operative radiotherapy.

A Visually Apparent and Quantifiable CT Imaging Feature Identifies Biophysical Subtypes of Pancreatic Ductal Adenocarcinoma

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) is a heterogeneous disease with variable presentations and natural histories of disease. We hypothesized that different morphologic characteristics of PDAC tumors on diagnostic computed tomography (CT) scans would reflect their underlying biology.

EXPERIMENTAL DESIGN

We developed a quantitative method to categorize the PDAC morphology on pretherapy CT scans from multiple datasets of patients with resectable and metastatic disease and correlated these patterns with clinical/pathologic measurements. We modeled macroscopic lesion growth computationally to test the effects of stroma on morphologic patterns, hypothesizing that the balance of proliferation and local migration rates of the cancer cells would determine tumor morphology.

RESULTS

In localized and metastatic PDAC, quantifying the change in enhancement on CT scans at the interface between tumor and parenchyma (delta) demonstrated that patients with conspicuous (high-delta) tumors had significantly less stroma, higher likelihood of multiple common pathway mutations, more mesenchymal features, higher likelihood of early distant metastasis, and shorter survival times compared with those with inconspicuous (low-delta) tumors. Pathologic measurements of stromal and mesenchymal features of the tumors supported the mathematical model's underlying theory for PDAC growth.

CONCLUSIONS

At baseline diagnosis, a visually striking and quantifiable CT imaging feature reflects the molecular and pathological heterogeneity of PDAC, and may be used to stratify patients into distinct subtypes. Moreover, growth patterns of PDAC may be described using physical principles, enabling new insights into diagnosis and treatment of this deadly disease.

Calibration of Multi-Parameter Models of Avascular Tumor Growth Using Time Resolved Microscopy Data

Two of the central challenges of using mathematical models for predicting the spatiotemporal development of tumors is the lack of appropriate data to calibrate the parameters of the model, and quantitative characterization of the uncertainties in both the experimental data and the modeling process itself. We present a sequence of experiments, with increasing complexity, designed to systematically calibrate the rates of apoptosis, proliferation, and necrosis, as well as mobility, within a phase-field tumor growth model. The in vitro experiments characterize the proliferation and death of human liver carcinoma cells under different initial cell concentrations, nutrient availabilities, and treatment conditions. A Bayesian framework is employed to quantify the uncertainties in model parameters. The average difference between the calibration and the data, across all time points is between 11.54% and 14.04% for the apoptosis experiments, 7.33% and 23.30% for the proliferation experiments, and 8.12% and 31.55% for the necrosis experiments. The results indicate the proposed experiment-computational approach is generalizable and appropriate for step-by-step calibration of multi-parameter models, yielding accurate estimations of model parameters related to rates of proliferation, apoptosis, and necrosis.

The invasion speed of cell migration models with realistic cell cycle time distributions

Cell proliferation is typically incorporated into stochastic mathematical models of cell migration by assuming that cell divisions occur after an exponentially distributed waiting time. Experimental observations, however, show that this assumption is often far from the real cell cycle time distribution (CCTD). Recent studies have suggested an alternative approach to modelling cell proliferation based on a multi-stage representation of the CCTD. In this paper we investigate the connection between the CCTD and the speed of the collective invasion. We first state a result for a general CCTD, which allows the computation of the invasion speed using the Laplace transform of the CCTD. We use this to deduce the range of speeds for the general case. We then focus on the more realistic case of multi-stage models, using both a stochastic agent-based model and a set of reaction-diffusion equations for the cells' average density. By studying the corresponding travelling wave solutions, we obtain an analytical expression for the speed of invasion for a general N-stage model with identical transition rates, in which case the resulting cell cycle times are Erlang distributed. We show that, for a general N-stage model, the Erlang distribution and the exponential distribution lead to the minimum and maximum invasion speed, respectively. This result allows us to determine the range of possible invasion speeds in terms of the average proliferation time for any multi-stage model.

Nonlinear studies of tumor morphological stability using a two-fluid flow model

We consider the nonlinear dynamics of an avascular tumor at the tissue scale using a two-fluid flow Stokes model, where the viscosity of the tumor and host microenvironment may be different. The viscosities reflect the combined properties of cell and extracellular matrix mixtures. We perform a linear morphological stability analysis of the tumors, and we investigate the role of nonlinearity using boundary-integral simulations in two dimensions. The tumor is non-necrotic, although cell death may occur through apoptosis. We demonstrate that tumor evolution is regulated by a reduced set of nondimensional parameters that characterize apoptosis, cell-cell/cell-extracellular matrix adhesion, vascularization and the ratio of tumor and host viscosities. A novel reformulation of the equations enables the use of standard boundary integral techniques to solve the equations numerically. Nonlinear simulation results are consistent with linear predictions for nearly circular tumors. As perturbations develop and grow, the linear and nonlinear results deviate and linear theory tends to underpredict the growth of perturbations. Simulations reveal two basic types of tumor shapes, depending on the viscosities of the tumor and microenvironment. When the tumor is more viscous than its environment, the tumors tend to develop invasive fingers and a branched-like structure. As the relative ratio of the tumor and host viscosities decreases, the tumors tend to grow with a more compact shape and develop complex invaginations of healthy regions that may become encapsulated in the tumor interior. Although our model utilizes a simplified description of the tumor and host biomechanics, our results are consistent with experiments in a variety of tumor types that suggest that there is a positive correlation between tumor stiffness and tumor aggressiveness.

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.

Modeling persistence of motion in a crowded environment: The diffusive limit of excluding velocity-jump processes

Persistence of motion is the tendency of an object to maintain motion in a direction for short time scales without necessarily being biased in any direction in the long term. One of the most appropriate mathematical tools to study this behavior is an agent-based velocity-jump process. In the absence of agent-agent interaction, the mean-field continuum limit of the agent-based model (ABM) gives rise to the well known hyperbolic telegraph equation. When agent-agent interaction is included in the ABM, a strictly advective system of partial differential equations (PDEs) can be derived at the population level. However, no diffusive limit of the ABM has been obtained from such a model. Connecting the microscopic behavior of the ABM to a diffusive macroscopic description is desirable, since it allows the exploration of a wider range of scenarios and establishes a direct connection with commonly used statistical tools of movement analysis. In order to connect the ABM at the population level to a diffusive PDE at the population level, we consider a generalization of the agent-based velocity-jump process on a two-dimensional lattice with three forms of agent interaction. This generalization allows us to take a diffusive limit and obtain a faithful population-level description. We investigate the properties of the model at both the individual and population levels and we elucidate some of the models' key characteristic features. In particular, we show an intrinsic anisotropy inherent to the models and we find evidence of a spontaneous form of aggregation at both the micro- and macroscales.

Three-Dimensional Spatiotemporal Modeling of Colon Cancer Organoids Reveals that Multimodal Control of Stem Cell Self-Renewal is a Critical Determinant of Size and Shape in Early Stages of Tumor Growth

We develop a three-dimensional multispecies mathematical model to simulate the growth of colon cancer organoids containing stem, progenitor and terminally differentiated cells, as a model of early (prevascular) tumor growth. Stem cells (SCs) secrete short-range self-renewal promoters (e.g., Wnt) and their long-range inhibitors (e.g., Dkk) and proliferate slowly. Committed progenitor (CP) cells proliferate more rapidly and differentiate to produce post-mitotic terminally differentiated cells that release differentiation promoters, forming negative feedback loops on SC and CP self-renewal. We demonstrate that SCs play a central role in normal and cancer colon organoids. Spatial patterning of the SC self-renewal promoter gives rise to SC clusters, which mimic stem cell niches, around the organoid surface, and drive the development of invasive fingers. We also study the effects of externally applied signaling factors. Applying bone morphogenic proteins, which inhibit SC and CP self-renewal, reduces invasiveness and organoid size. Applying hepatocyte growth factor, which enhances SC self-renewal, produces larger sizes and enhances finger development at low concentrations but suppresses fingers at high concentrations. These results are consistent with recent experiments on colon organoids. Because many cancers are hierarchically organized and are subject to feedback regulation similar to that in normal tissues, our results suggest that in cancer, control of cancer stem cell self-renewal should influence the size and shape in similar ways, thereby opening the door to novel therapies.

Activation of the HGF/c-Met axis in the tumor microenvironment: A multispecies model

The tumor microenvironment is an integral component in promoting tumor development. Cancer-associated fibroblasts (CAFs), which reside in the tumor stroma, produce Hepatocyte Growth Factor (HGF), an important trigger for invasive and metastatic tumor behavior. HGF contributes to a pro-tumorigenic environment by activating its cognate receptor, c-Met, on tumor cells. Tumor cells, in turn, secrete growth factors that upregulate HGF production in CAFs, thereby establishing a dynamic tumor-host signaling program. Using a spatiotemporal multispecies model of tumor growth, we investigate how the development and spread of a tumor is impacted by the initiation of a dynamic interaction between tumor-derived growth factors and CAF-derived HGF. We show that establishment of such an interaction results in increased tumor growth and morphological instability, the latter due in part to increased cell species heterogeneity at the tumor-host boundary. Invasive behavior is further increased if the tumor lowers responsiveness to paracrine pro-differentiation signals, which is a hallmark of neoplastic development. By modeling anti-HGF and anti-c-Met therapy, we show how disruption of the HGF/c-Met axis can reduce tumor invasiveness and growth, thereby providing theoretical evidence that targeting tumor-microenvironment interactions is a promising avenue for therapeutic development.

Cell migration in microengineered tumor environments

Recent advances in microengineered cell migration platforms are discussed critically with a focus on how cell migration is influenced by engineered tumor microenvironments, the medical relevance being to understand how tumor microenvironments may promote or suppress the progression of cancer. We first introduce key findings in cancer cell migration under the influence of the physical environment, which is systematically controlled by microengineering technology, followed by multi-cues of physico-chemical factors, which represent the complexity of the tumor environment. Recognizing that cancer cells constantly communicate not only with each other but also with tumor-associated cells such as vascular, fibroblast, and immune cells, and also with non-cellular components, it follows that cell motility in tumor microenvironments, especially metastasis via the invasion of cancer cells into the extracellular matrix and other tissues, is closely related to the malignancy of cancer-related mortality. Medical relevance of forefront research realized in microfabricated devices, such as single cell sorting based on the analysis of cell migration behavior, may assist personalized theragnostics based on the cell migration phenotype. Furthermore, we urge development of theory and numerical understanding of single or collective cell migration in microengineered platforms to gain new insights in cancer metastasis and in therapeutic strategies.

In Vitro/In Silico Study on the Role of Doubling Time Heterogeneity among Primary Glioblastoma Cell Lines

The application of accurate cancer predictive algorithms validated with experimental data is a field concerning both basic researchers and clinicians, especially regarding a highly aggressive form of cancer, such as Glioblastoma. In an aim to enhance prediction accuracy in realistic patient-specific environments, accounting for both inter- and intratumoral heterogeneity, we use patient-derived Glioblastoma cells from different patients. We focus on cell proliferation using in vitro experiments to estimate cell doubling times and sizes for established primary Glioblastoma cell lines. A preclinically driven mathematical model parametrization is accomplished by taking into account the experimental measurements. As a control cell line we use the well-studied U87MG cells. Both in vitro and in silico results presented support that the variance between tumor staging can be attributed to the differential proliferative capacity of the different Glioblastoma cells. More specifically, the intratumoral heterogeneity together with the overall proliferation reflected in both the proliferation rate and the mechanical cell contact inhibition can predict the in vitro evolution of different Glioblastoma cell lines growing under the same conditions. Undoubtedly, additional imaging techniques capable of providing spatial information of tumor cell physiology and microenvironment will enhance our understanding regarding Glioblastoma nature and verify and further improve our predictability.